Skip to main content

ORIGINAL RESEARCH article

Front. Plant Sci., 21 January 2025
Sec. Plant Bioinformatics
This article is part of the Research Topic Integrating Omics to Unveil Legume Stress Response Mechanisms View all 5 articles

Integrated metabolome and transcriptome analysis reveals potential mechanism during the bud dormancy transition of Astragalus membranaceus (Fisch.) Bge. var. mongholicus (Bge.) Hsiao

Huan Guan,Huan Guan1,2Yuhuan Zhao,Yuhuan Zhao1,2Qing ChenQing Chen3Qianqian Zhang,Qianqian Zhang1,2Pengpeng Yang,Pengpeng Yang1,2Shuying Sun,*Shuying Sun1,2*Guilin Chen,*Guilin Chen1,2*
  • 1Key Laboratory of Herbage and Endemic Crop Biology, Ministry of Education, School of Life Sciences, Inner Mongolia University, Hohhot, China
  • 2The Good Agriculture Practice Engineering Technology Research Center of Chinese and Mongolian Medicine in Inner Mongolia, Inner Mongolia University, Hohhot, China
  • 3Spallation Neutron Source Science Center, Institute of High Energy Physics, Chinese Academy of Science, Dongguan, China

Astragalus membranaceus (Fisch.) Bge. var. mongholicus (Bge.) Hsiao (AMM) is an important medicinal plant that is used for both medicine and food. It is widely used in Asia and South Asia. It is normally cultivated by transplanting the annual rhizomes. Understanding the dormancy of underground buds of AMM is essential for its harvest and transplantation. Despite thorough research on bud dormancy in perennial woody plants, perennial herbs, and especially medicinal plants, such as AMM, have rarely been studied. We analyzed the transcriptome and non-targeted metabolome of dormant buds stage-by-stage to investigate the regulatory mechanism of the transition from endo- to ecodormancy. A total of 1,069 differentially accumulated metabolites (DAMs) participated in amino acid and carbohydrate metabolism. Transcriptome analysis revealed 16,832 differentially expressed genes (DEGs). Functional enrichment analysis indicated that carbohydrate metabolism, hormone signaling pathways, and amino acid metabolism contributed to the transition from endo- to ecodormancy. Starch and sucrose metabolism and hormone signaling pathways were mainly analyzed in the transition between different dormancy states. During the transition from endo- to ecodormancy, the highest content of indole-3-acetic acid (IAA) and the highest number of DEGs enriched in the IAA signaling pathway demonstrated that IAA may play a key role in this process. We obtained candidate genes through co-expression network analysis, such as BGL, GN, glgC, and glgB, which are involved in starch and sucrose metabolism. The transcription factors MYB, ERF, bHLH, zinc finger, and MADS-box may regulate the genes involved in hormone signal transduction and starch and sucrose metabolism, which are critical for regulating the transition from endo- to ecodormancy in AMM buds. In summary, these results provide insights into the novel regulatory mechanism of the transition of endo- to ecodormancy in underground buds of AMM and offer new analytical strategies for breaking dormancy in advance and shortening breeding time.

1 Introduction

Forming buds and undergoing dormancy are crucial adaptations for perennials during their life cycle to survive harsh winter environments. Environmental factors, such as temperature and light, play a vital role in regulating the induction and release of dormancy (Horvath et al., 2010; Wang et al., 2022a; Rothkegel et al., 2020). Dormancy can be divided into paradormancy, endodormancy, and ecodormancy, depending on the physiological state of dormant tissues (Lang et al., 1987). In the fall, perennial plants enter a dormant period, protecting themselves from harsh winter conditions by temporarily ceasing growth and activity until more favorable conditions return. This process, known as winter dormancy, is essential to the seasonal cycle and supports normal plant growth the following year. Winter dormancy consists of two consecutive stages: endodormancy and ecodormancy (Li et al., 2022b). Endodormancy is mainly regulated by internal factors that inhibit growth, even under favorable environmental conditions (Yordanov et al., 2014; Zhu et al., 2020). Only after the chilling requirements for a specific period, as determined by the variety, are met, leading to the release of endodormancy (Fadón et al., 2020). Insufficient CRs lead to abnormal growth and flowering, resulting in economic loss (Zhang et al., 2017; Ladwig et al., 2019). Once endodormancy is released, bud growth potential is restored and the bud enters an ecodormancy state. However, during this phase, they are sensitive to external environmental factors (Charrier, 2022). When adverse environmental conditions subside, warmer temperatures can lead to visible bud cracking. Studies on woody plants and perennial herbs have shown that bud dormancy is a complex biological phenomenon that requires a series of physiological, biochemical, and molecular mechanisms involving changes in the transport of nutrients, biosynthesis of nucleic acids and proteins, metabolism of carbohydrates, generation of plant endogenous hormones, and other compound shift activities (Horikoshi et al., 2018; Shangguan et al., 2020; André et al., 2022; Li et al., 2022a; Sapkota et al., 2023). Hormone and carbohydrate metabolism may be the main factors that regulate dormancy release (Wang et al., 2018). The regulation of winter dormancy in plants is closely tied to ambient temperature, which means that global warming will significantly affect this process in various ways (Lempe et al., 2024). Therefore, it is essential to understand the stages of dormancy and the mechanisms that regulate them to ensure the sustainable development of production. However, there are few reports on herbaceous plants, particularly medicinal plants.

Plant hormones are considered to be responsible for the regulation of bud dormancy. Abscisic acid (ABA) is involved in the establishment and maintenance of bud dormancy (Chao et al., 2016; Tylewicz et al., 2018). A high ABA content can contribute to maintaining dormancy and can slow down the process of dormancy release (Yang et al., 2021; Li et al., 2018; Horvath et al., 2003). Key genes such as Nine-cis-polyphenolenoid dioxygenase (NCED) and ABA 8’-hydroxylase (CYP707A) are responsible for ABA biosynthesis and degradation and are associated with dormancy in perennial plants, such as peaches (Li et al., 2021), Iris japonica Thunb (Xu et al., 2022), and petunia (Luo et al., 2023). High ABA levels induce and maintain endodormancy, whereas endodormancy release depends on a sufficient concentration of gibberellin (GA) (Li et al., 2020; Zhang et al., 2018; Gao et al., 2023). GA20 oxidase (GA20ox) and GA3ox are key enzymes that regulate GA biosynthesis, and GA2ox can catalyze the conversion of biologically active GA to non-active forms (Liu and Hou, 2018; Liu et al., 2022). Studies have shown that highly expressed GA3ox and GA20ox play an important role in regulating bud dormancy release in plants such as tea (Yue et al., 2018) and tree peony (Yuan et al., 2024). IAA and cytokinin (CTK) are signaling molecules that participate in bud dormancy (Qiu et al., 2019; Liu and Sherif, 2019). IAA is one of the most abundant auxins, and the levels of IAA in Chinese fir (Cunninghamia lanceolata (Lamb.) Hook) significantly increase during dormancy release (Qiu et al., 2013). IAA levels in underground buds and potato tubers of leafy spurge (Euphorbia esula L.) also increased during endodormancy release but decreased during ecodormancy (Chao et al., 2017; Sorce et al., 2000). Studies have shown that the trend of IAA content is similar to that of GA content, indicating that IAA may cooperate with GA to regulate dormancy and promote budbreak (Zhang et al., 2020; Fan et al., 2021). CTK also plays a positive role in the release of bud dormancy (Ito et al., 2021). Hormones may regulate bud dormancy through complex, intertwined pathways (Zhang et al., 2018; Liu and Sherif, 2019); however, the details remain to be further explored.

Carbohydrates have been shown to be closely related to bud dormancy (Park et al., 2009; Khalil-Ur-Rehman et al., 2017; Signorelli et al., 2018). Tarancón et al. (2017) reported that buds were in a state of carbon starvation, induced by both external environmental and endogenous factors during dormancy. Therefore, it is necessary to establish a balance between the stored starch and soluble sugars to maintain their basic functions. Dormancy release is triggered by starch degradation of starch and an increase in soluble sugar content (Zhang et al., 2021a). When Hippeastrum vittatum ‘Red Lion’ bulbs and Oriental hybrid lily ‘Sorbonne’ undergo dormancy release, the soluble sugar content increases due to starch degradation (Wang et al., 2014; Liu et al., 2014). Soluble sugars can function as energy supplies or signaling molecules to promote dormancy release. Some studies have confirmed that the interaction between carbohydrates and hormones plays an important role in bud dormancy release (Jiang et al., 2020; Wang et al., 2022a). GA3 can significantly enhance the expression of Sucrose Phosphate Synthase (SPS) during bud dormancy release in Paeonia lactiflora (Jiang et al., 2020). The interaction between carbohydrates and hormones mainly focuses on sucrose and IAA during dormancy release in Polygonatum kingianum (Wang et al., 2022b). Energy metabolism is also closely related to dormancy in grapes (Rehman et al., 2019) and Phalaenopsis amabilis (Chen et al., 2020). Furthermore, glycolysis, tricarboxylic acid (TCA) cycle, and pentose phosphate pathway (PPP) are involved in dormancy release (Zhang et al., 2020). The production of metabolites induced by hydrogen cyanamide (HC) treatment is closely involved in the TCA cycle and PPP, thereby generating sufficient energy to promote endodormancy release (Wang et al., 2021). Similarly, during dormancy release of blueberry flower buds, the expression levels of almost all DEGs involved in the TCA cycle were upregulated, indicating that the energy metabolism process was activated (Li et al., 2022a).

Transcription factors (TFs), including WRKY, APETALA2/ethylene response factors (AP2/ERF), basic region-leucine zipper (bZIP), MADS, Teosinte branched1/Cycloidea/Proliferating cell factor (TCP), and other families are an important group of factors that regulate bud dormancy transition. Previous studies have shown that the WRKY family is an important and conserved TF that plays a crucial role in dormancy release and the cold stress response (Zhang et al., 2017; Shangguan et al., 2020; Zhou et al., 2020). WRKYs and MYBs are upregulated during dormancy release in lily bulbs, indicating that these genes induce bulb dormancy release (Fan et al., 2021). The AP2/ERF family also plays a crucial role in regulating plant growth, development, and freezing stress responses (Yu et al., 2022). The DREB members within the AP2/ERF family are most closely associated with low-temperature stress responses (Liu et al., 2019) and are expressed under the induction of low temperatures (Li et al., 2019b). MADS-box genes are closely related to bud dormancy in various species (Niu et al., 2016), and SOC1 has been shown to affect the dormancy of kiwifruits (Voogd et al., 2015) and poplars (Gómez-Soto et al., 2021).

Astragalus membranaceus (Fisch.) Bge. var. mongholicus (Bge.) Hsiao (abbreviated as AMM in the following of the manuscript) is a perennial herbaceous medicinal plant that grows in northwest China (Chen et al., 2015). The dried roots of AMM have demonstrated antioxidant, anti-inflammatory, and immunoregulatory functions (Sheik et al., 2021; Zheng et al., 2020) and have a long history of use in traditional Chinese medicine (Dong et al., 2023). It is also one of the most commonly used medicinal materials in South and East Asia (Li et al., 2019a). Additionally, AMM has extensive applications in medicine, food, health products, cosmetics, etc (Ji et al., 2023). Similar to many perennial herbaceous plants, the underground buds on the crowns (the underground part of the stem) of AMM remain dormant during winter to ensure low temperatures. AMM is primarily cultivated on a large scale by transplanting high-quality annual rhizomes in either autumn or spring. Therefore, understanding the dormant period and potential mechanisms for dormancy release of underground buds is crucial for effective transplantation, growth, and harvesting, as well as for reducing breeding cycles. Although research on bud dormancy in perennial woody plants has been relatively thorough, the regulatory mechanism of perennial herbaceous plants, particularly medicinal plants such as AMM, remains not fully understood.

Due to the complexity of genomes and the lack of efficiency in genetic transformation techniques, research on dormancy induction and release in most perennial plants is limited and is a relatively new research field. AMM is an important plant used in traditional Chinese medicine. The underground bud dormancy of AMM hinders the growth and reproduction of these medicinal plants to a certain extent, but the molecular mechanism that regulates dormancy release remains unclear. Through transcriptomics and metabolomics analysis, along with co-expression network analysis, we clarified the changes in genes and metabolites associated with dormancy transition in AMM underground buds. Additionally, we identified molecular pathways and metabolic processes that may be involved in dormancy release. Our findings indicate that hormones and carbohydrate metabolism are crucial in regulating dormancy transition in AMM underground buds. This study confirms that the transition from endo- to ecodormancy in AMM involves the integration of various environmental and internal factors, resulting in a complex regulatory network. Additionally, this study provides new insights into the mechanisms underlying this transition in perennial herbaceous plants.

2 Materials and methods

2.1 Plant materials

The artificially cultivated annual Astragalus membranaceus (Fisch.) Bge. var. mongholicus (Bge.) Hsiao (AMM) was introduced from Harqin Banner of Chifeng City (E 118°07′-119°20′, N 41°33′-42°14′) on October 2021, and was identified by Professor Guilin Chen, School of Life Sciences, Inner Mongolia University. Roots with underground buds of AMM were transplanted into pots of 17 cm diameter and 25 cm height with a mixture of soil, peat, and perlite (2:1:1), and carried out near the greenhouse of Inner Mongolia University (E 111°38′, N 40°17′), China.

2.2 Experimental design and morphological observation

The underground buds of potted one-year-old AMM (three biological replicates and seven plants per replicate) were transferred to a glasshouse (26°C–15°C day/night, natural light and artificial supplementary lighting, regular watering) every week or so from 31 October 2021. The specific transfer dates were 31 October, 8 November, 15 November, 22 November, 1 December, and 6 December, respectively. The dormancy release phase was assessed based on the method proposed by Wang et al. (2022a) with some modifications. The five morphological indices were observed weekly to evaluate potted plants after being transferred to the glasshouse, including the days until the first plant sprouted in the glasshouse (DFS), the days of bud break percentage ≥60 (DBS), bud break percentage four weeks after being transferred into the glasshouse (BPF), the days until all plants sprouted in the glasshouse (DAS), and plant height four weeks after being transferred into the glasshouse (PHF). The date when most morphological indices showed no significant increase or decrease, even when longer low-temperature exposure occurred, was defined as endodormancy release (Wang et al., 2020). The standard of bud break was when bud scales opened and exposed the leaves. If a bud sprouted on a plant in the pot, it was considered a sprouted plant.

2.3 Sample collection and selection

Morphological analysis indicated that the morphological indicators of AMM moved into the greenhouse on 31 October, which showed significant increases or decreases compared with the other groups. This suggested that the plants on this date were in the endodormancy stage. For the plants transferred on 15 and 22 November, most morphological indicators gradually stabilized, leading to the inference that by 22 November, the plants had been released from endodormancy. Once endodormancy breaks, plants enter the ecodormancy stage. Thus, the plants transferred on 1 and 6 December were considered to be in the ecodormancy stage (Table 1). Consequently, the samples from 6 December were not included in subsequent research. The potted plants used for sample collection were from batches different from those used for morphological observations. Underground buds were collected on 31 October, 8 November, 15 November, 22 November, and 1 December, respectively (Supplementary Figure S1). Since there were at most two large buds on the surface of the crowns, only one plump and one strong bud were collected from each crown. Underground buds were cut, washed with distilled water, and immediately frozen in liquid nitrogen. These samples were stored at −80°C for subsequent transcriptomic, metabolomic, and physiological analyses. Each time point included three biological replicates corresponding to the buds from three distinct plants.

Table 1
www.frontiersin.org

Table 1. Morphological observations of AMM after natural chilling treatments in growth room.

2.4 Metabolite extraction and analysis

Metabolites in underground buds collected at five time points were extracted and analyzed by Shanghai Majorbio Bio-Pharm Co., Ltd. (Shanghai, China). The underground buds were crushed, and lyophilized powder (50 mg) was weighed and extracted using 400 µL of 80% aqueous methanol with 0.02 mg/mL L-2-chlorophenylalanin as the internal standard. The mixture was treated with a high-throughput tissue crusher Wonbio-96c (Shanghai Wanbo Biotechnology Co., Ltd.). The supernatant was filtered through a 0.22 μm filter membrane after centrifugation and carefully transferred to sample vials for LC–MS/MS analysis.

Raw data from LC/MS pretreated using Progenesis QI (Waters Corporation, Milford, USA) software, and a three-dimensional data matrix in CSV format was exported. Metabolite annotation was performed according to the main databases: HMDB (http://www.hmdb.ca/) and Metlin (https://metlin.scripps.edu/). Orthogonal partial least squares discriminant analysis (OPLS-DA) and PCA were performed with the R package (v1.6.2). Differentially accumulated metabolites (DAMs) among groups were determined by variable importance in the projection (VIP) >1 (P <0.05). These DAMs were further characterized through KEGG functional analysis (http://www.Genome.jp/keg/). K-means clustering of DAMs was performed using SciPy v1.0.0.

2.5 Transcriptomic profiling

Total RNA was extracted from 15 underground buds using the EasyPure Plant RNA Kit (TransGen, Beijing, China), according to the manufacturer’s instructions. The RNA quality was assessed using a 2100 Bioanalyzer (Agilent, USA).

RNA purification, reverse transcription, library construction, and sequencing were conducted by Shanghai Majorbio Bio-Pharm Co., Ltd. (Shanghai, China) according to the manufacturer’s instructions (Illumina, San Diego, CA, USA). RNA-seq libraries were sequenced using an Illumina NovaSeq 6000 sequencer (Illumina, San Diego, CA, USA). Raw reads were trimmed and quality controlled using SeqPrep (https://github.com/jstjohn/SeqPrep) and Sickle (https://github.com/najoshi/sickle) with default parameters. Clean reads were used for do de novo assembly with Trinity (http://trinityrnaseq.sourceforge.net/). The assembly unigenes were evaluated using BUSCO v3.0.2. The assembled unigenes were annotated using NR (https://www.ncbi.nlm.nih.gov/public/), Swiss-Prot (http://www.uniprot.org/), Pfam (http://pfam.xfam.org/), GO (http://www.geneontology.org), and KEGG (http://www.genome.jp/kegg/) databases. By analyzing the domain information contained in the transcript, gene families of transcription factors were predicted and analyzed using the PlantTFDB5.0 database (http://planttfdb.gao-lab.org/).

Fragments per kilobase of exon per million mapped fragments (FPKM) were used to quantify transcript or gene expression levels. Differential expression analysis between the two groups was performed using the DESeq2 v1.24.0, and DEGs screening criteria were |log2FC|>1 and false discovery rate (FDR) <0.05. FDR was estimated by correcting the P-value, and the Benjamini–Hochberg method was used to correct the P-value. GO functional enrichment of DEGs was carried out using Goatools v0.6.5, using Fisher’s exact test. KOBAS v2.1.1 was used for enrichment analysis of the DEGs in KEGG pathways based on the Fisher test.

2.6 qRT−PCR validation

The relative expression levels of the 16 selected genes were analyzed by qRT−PCR. Primers for candidate genes were designed using Primer-BLAST (https://www.ncbi.nlm.nih.gov/tools/primer-blast/) and are shown in Supplementary Table S1. The reference gene was 18S, and all reactions were performed in triplicate using the TransStart Tip Green qPCR SuperMix (TransGen, Beijing, China) on an ABI Prism7500 (Thermo Fisher, Singapore). The thermal cycle of SYBR Green qRT-PCR was as follows: 94°C for 30 s, followed by 40 cycles at 94°C for 30 s, 60°C for 34 s, and 72°C for 30 s. The relative gene expression was calculated using the 2−ΔΔCt method.

2.7 Weighted gene co-expression network analysis

After filtering the low-expression genes (mean FPKM of all samples <1), the remaining genes were used for co-expression network analysis with the WGCNA R package (v4.1.3). A hierarchical clustering dendrogram of the genes was created using the dynamic tree-cutting algorithm with mergeCutHeight = 0.25 and minModuleSize = 30. The correlations between gene modules and six sugar contents were analyzed by estimating module–trait relationships. Co-expression networks were formed using Cytoscape software v.3.10.1, with default settings.

2.8 Analysis of phytohormones in underground buds

The levels of ABA, GA3, IAA, and CTK in underground buds were measured by enzyme-linked immunosorbent assay (ELISA) (Yang et al., 2001) using an ELISA kit (MLBIO, Shanghai, China), repeated three times. Approximately 50 mg of buds was ground and extracted in 1 mL of cold 80%(v/v) methanol overnight at 4°C. After centrifugation for 10 min at 10,000×g and 4°C, the supernatant was collected. An additional 500 μL of 80% methanol was added to the precipitate for further extraction, which was performed for 2 h at 4°C. Following another round of centrifugation for 10 min at 10,000×g and 4°C, the supernatants were combined and dried with nitrogen. The pellet was then dissolved in 300 μL of 30% methanol and filtered through a filter membrane (0.22 µm). Hormone levels were determined according to the manufacturer’s instructions using an ELISA kit.

2.9 Carbohydrate concentration measurements

The soluble sugar and sucrose contents in the underground buds were determined using the approach of Du et al. (2020). Starch was extracted and quantified using a modified method described by Zhang et al. (2022).

2.10 Statistical analysis

SPSS software v19.0 (SPSS, Inc., Chicago, IL, USA) was used to conduct one-way analysis of the data with Duncan’s multiple comparison test. The results are expressed as mean ± standard deviation (SD) of three biological replicates. Line charts were generated using GraphPad Prism v9.0.

3 Results

3.1 Determining the dormancy status and physiological dynamics of underground buds of AMM

We followed the method used by Wang et al. (2022a) to determine the dormant states of P. lactiflora. Using this method, we transferred the one-year-old Astragalus membranaceus (Fisch.) Bge. Var. mongholicus (Bge.) Hsiao (AMM) under natural chilling conditions into a greenhouse on 31 October, 8 November, 15 November, 22 November, 1 December, and 6 December 2021, respectively, to monitor the development of underground buds. The bud break percentage four weeks after being transferred into the glasshouse (BPF) on 31 October was only 24%, and its BPF never reached 60% (Table 1). The days of the first plant to sprout in the glasshouse (DFS), the days of bud break percentage≧60 (DBS), and the days required for all plants to sprout in the glasshouse (DAS) gradually decreased. The BPF increased significantly on 15 November and reached 100% by 22 November; most of the indices were also stable. Based on these observations, 22 November was determined as the endodormancy release date. 1 and 6 December were classified as the ecodormancy stages. Because the samples collected on 1 and 6 December were in the same dormant stage, only the samples from the first five time points were selected for further analysis, and they were named 31 October, 8 November, 15 November, 22 November, and 1 December, respectively.

Plant hormones play an important role in the dormancy transition process; therefore, we measured hormone levels in underground buds at five time points. IAA had the highest concentration, exceeding 20 μg·g−1 FW, whereas CTK had the lowest concentration, recorded at just over 0.06 μg·g−1 FW (Figure 1A). The trends in the levels of the four hormones were similar: they gradually increased from 31 October to 15 November, reaching a peak on 15 November. Following this peak, the levels decreased sharply on 22 November, before experiencing a significant increase again.

Figure 1
www.frontiersin.org

Figure 1. Changes of hormone and sugar content in AMM underground buds during the transition from endo- to ecodormancy. (A) The contents of GA3, ABA, IAA, and CTK in the underground buds of AMM during different dormant stages. (B) Changes in the contents of soluble sugar, sucrose, and starch in AMM underground buds during the transition from endo- to ecodormancy.

Changes in carbohydrate levels are associated with dormancy and are often used as markers to divide bud dormancy. In underground buds, the contents of soluble sugar and sucrose consistently increased from 31 October to 1 December. However, starch levels sharply declined from 31 October to 8 November. Similarly, significant changes were observed in morphological indices, such as DFS, DBS, and DAS, during this period. On 1 December, the starch content decreased (Figure 1B). In conclusion, based on the morphological indices and trend in starch content, 31 October was divided into the endodormancy stage. The period from 31 October to 22 November was identified as the endodormancy release stage, and the ecodormancy stage was subsequently established.

3.2 Identifying the transition from endo- to ecodormancy of AMM underground buds by metabolic differences

To understand the metabolic profiles of AMM underground buds during the transition from endo-to ecodormancy, 15 samples from five transfer time points were collected and subjected to non-targeted metabolomic analysis using LC–MS/MS. In total, 1,816 metabolites were identified. Principal component analysis (PCA) showed that 31 October and 1 December were separated from the other samples (Figure 2A), indicating that the accumulation of metabolites was associated with different dormancy stages. After pairwise comparison, 1,069 DAMs were obtained, consisting of 21 classes, including 127 carboxylic acids and derivatives,117 fatty acyls, 71 organooxygen compounds, 71 prenol lipids, 52 glycerophospholipids, 24 steroids and steroid derivatives, 18 benzene and substituted derivatives, 12 flavonoids, 12 purine nucleosides, and 171 compounds from other classes (Supplementary Table S2; Figure 2B). After pairwise comparison, it was found that the highest number of DAMs were observed in the comparison of 1 December_vs_31 October, with 454 downregulated and 100 upregulated metabolites. This was followed by the 22 November_vs_31 October, which showed 356 downregulated and 193 upregulated metabolites. Additionally, a large number of DAMs were identified on 1 December_vs_ 22 November. These findings suggest a significant change in the metabolites of underground buds on 31 October, 22 November, and 1 December, which may indicate different dormancy states at these three time points (Figure 2C). During the endodormancy release in the underground buds of AMM, we analyzed 1,069 annotated DAMs using K-means cluster analysis to track changes in their relative abundance. The DAMs were divided into three subclusters (Figure 2D). Based on the accumulation patterns of the subclusters, we observed that the metabolites in subcluster 1 showed a continuously decreasing trend during the transition from endo- to ecodormancy (from 31 October to 1 December), with relatively high numbers of carboxylic acids and derivatives and fatty acyls, accounting for 21.3% and 21.04% of the total, respectively (Figure 2E; Supplementary Table S3). In contrast, the metabolites in subcluster 3 showed a gradually increasing trend during the transition from endo- to ecodormancy, with glycerophoripids, organooxygen compounds, carboxylic acids and derivatives, and fatty acyls dominating. The metabolites in subcluster 2 exhibited continuous accumulation during endodormancy release (from 31 October 31 to 22 November), with carboxylic acids and derivatives, fatty acyls, and organooxygen compounds being the most prevalent, accounting for 15.14%, 12.43%, and 11.89%, respectively.

Figure 2
www.frontiersin.org

Figure 2. Metabolomic analysis of the transition from endo- to ecodormancy in underground buds of AMM. (A) Principal component analysis (PCA) of metabolome data in different groups; (B) Classification and number of 1,069 DAMs; (C) Statistical analyses of the number of DAMs from different pairwise comparisons. (D) K-means cluster analysis of DAMs. (E) The number of different classes of DAMs in different subclusters. (F) KEGG enrichment analysis of DAMs in the three subclusters.

KEGG enrichment analysis showed that the clustered DAMs were involved in different metabolic and biosynthetic pathways (Figure 2F). The darker the color, the smaller is the P-value. During the transition from endo- to ecodormancy, DAMs are enriched in various amino acid metabolism pathways. In addition, the metabolic pathways of alpha-linolenic acid, glycerophospholipid, glycine, serine, threonine, phenylalanine, and galactose were weakened (subcluster 1). However, glycerophospholipid metabolism, zeatin biosynthesis, and starch and sucrose metabolism increased during this process (subcluster 3). The metabolites involved in flavonoid biosynthesis, tyrosine metabolism, and lysine biosynthesis were enriched during endodormancy release, indicating that these pathways were activated when endodormancy was released (subcluster 2).

3.3 Identifying the differently expressed genes in underground buds of AMM

As shown in Supplementary Table S4, the assembled unigenes were evaluated, revealing an N50 length of 1,598 bp and a BUSCO score of 65.7%. The transcriptome sequencing results showed a total of 43,271,254–72,963,680 clean reads with Q30 ranging from 91.36% to 92.82%, and a GC content of 43% in the AMM underground buds. The results showed that 72.31% to 75.14% of the clean reads could be mapped to the assembled reference sequence (Supplementary Table S5), which provided available transcript information for analysis. PCA of the transcriptome data revealed that the three biological repeats in each group were highly similar. Additionally, 31 October, 15 November, and 1 December were significantly different from those of the other groups. The first principal component (PC1) accounted for 15.04% of the total variance, separating 31 October, 8 November, and 15 November from 22 November and 1 December (Figure 3A). 1 December_vs_31 October had the largest number of DEGs, followed by 1 December_vs_15 November in which 3,026 and 2,874 were upregulated, and 5,065 and 3,173 downregulated DEGs were identified, respectively (Figure 3B). Unlike the transcriptome results, there were more DAMs on 1 December_vs_22 November, indicating significant changes in gene expression on 15 November, possibly in preparation for metabolite accumulation on 22 November. To investigate gene expression patterns during different dormancy stages, DEGs were further divided into five subclusters using K-means clustering analysis based on their relative expression, which was not the same as the clusters from the metabolic analysis. The gene expression levels in subclusters 1–5 peaked on 31 October, 8 November, 15 November, 22 November, and 1 December, respectively, which corresponded to the endodormancy stage, endodormancy release stage, early stage of complete release of endodormancy, endodormancy release date, and ecodormancy stage, respectively (Figure 3C).

Figure 3
www.frontiersin.org

Figure 3. Analysis of DEGs in underground buds of AMM during the transition from endo- to ecodormancy. (A) Principal component analysis PCA for transcriptomes; (B) Bar graph of up- and downregulated DEGs from different pairwise comparisons; (C) K-means cluster analysis and KEGG enrichment analysis of 16832 DEGs. The DEGs are divided into five subclusters. The top 6 KEGG pathways ranked by P-value in each subcluster are listed on the corresponding right panel. The font color represents the P-value, and the redder the color, the smaller the P-value.

To further determine the biological functions of the DEGs in the five subclusters, we performed GO and KEGG enrichment analyses. GO analysis of these subclusters showed that DEGs were enriched in molecular function, cellular component, and biological processes (Supplementary Table S6). Genes highly expressed during endodormancy were enriched in linoleic acid metabolism, flavonoid biosynthesis, and plant hormone signal transduction (subcluster 1). The genes in subcluster 2 were enriched in phenylpropanoid biosynthesis, flavonoid biosynthesis, and plant hormone signal transduction. The genes in subcluster 3 were enriched for tryptophan metabolism, DNA replication, and starch and sucrose metabolism. The KEGG pathways enriched by the DEGs in subcluster 4 included plant hormone signal transduction, amino sugar and nucleotide sugar metabolism, and phenylpropanoid biosynthesis, consistent with those observed in the metabolome analysis. Genes highly expressed during the ecodormancy phase were enriched in ribosome biogenesis in eukaryotes, circadian rhythm, autophagy, and glycerophospholipid metabolism (Figure 3C; Supplementary Table S6).

3.4 Validating DEGs by qRT-PCR

To validate the results of the transcriptomic analysis, we carried out qRT-PCR on 16 selected genes, which were mainly involved in hormone biosynthesis, signal transduction, starch and sucrose metabolism, and TFs (Supplementary Table S7). The results indicated that GA3 (TRINITY_DN34977_c0_g1), GID2 (TRINITY_DN18555_c0_g1), and ARF7 (TRINITY_DN6952_c1_g2) were upregulated and then downregulated during endodormancy release, whereas GA2ox2 (TRINITY_DN8660_c0_g1) and SAUR32 (TRINITY_DN10943_c0_g1) were downregulated (Supplementary Figure S2). Overall, the qRT-PCR results were consistent with the transcriptomic results, indicating the reliability of the RNA-sequencing data.

3.5 The responses of hormone signal transduction related genes during the dormancy transition

According to KEGG analysis, the genes involved in hormone signal transduction were significantly enriched during endodormancy release (Figure 3C). The profiles of the DEGs involved in hormone biosynthesis and metabolism were investigated (Supplementary Table S8). The key DEGs involved in ABA, GA, and IAA synthesis, such as NCEDs, GA20ox, GA3ox, YUCCA2, and YUCCA10, were highly upregulated from 31 October to 15 November and significantly decreased from 22 November to 1 December (Supplementary Figures S3A–C). Five ABA degradation genes, CYP707As were identified, in which the expression levels of three CYP707A2 genes were relatively low from 31 October to 15 November but increased significantly from 22 November to 1 December. Among the genes related to CTK metabolism, three CKXs (CKX1, CKX5, and CKX6) were downregulated from 31 October to 15 November (Supplementary Figure S3D), consistent with the changes observed in GA, ABA, IAA, and CTK contents.

To explore the role of hormone signal transduction in the transition from endo- to ecodormancy, we analyzed the genes in subclusters 1, 4, and 5. These subclusters exhibited high expression during the endodormancy stage, at the endodormancy release date, and during the ecodormancy stage. Additionally, we discovered that hormone levels peaked on 15 November, coinciding with the highest number of DEGs in the samples collected at this time. Therefore, we analyzed the genes in subcluster 3. A total of 78 DEGs related to ABA, GA, IAA, and CTK signaling were identified (Supplementary Table S9). There were more DEGs involved in IAA and ABA signaling pathways, 44 and 23, respectively (Figure 4). ABA-related genes PYLs, were downregulated from 31 October to 15 November. In contrast, SNRK2s and ABFs were significantly upregulated during the same period. DEGs associated with AUX/IAA and ARFs were upregulated gradually from 31 October to 15 November, peaked on 15 November, and then decreased significantly, whereas SAUR and GH3 displayed the opposite trend. Fewer genes were involved in GA and CTK signaling and GID2 was expressed significantly from 31 October to 15 November, followed by a significant decrease. Among the three annotated PIFs, the expression level of PIF3 was highest on 15 November. The expression patterns of genes in the CTK signaling pathway also varied, with two ARR-As gradually upregulated during the transition from endo- to ecodormancy.

Figure 4
www.frontiersin.org

Figure 4. Hormone signaling pathways during the transition from endo- to ecodormancy in the underground buds.

To identify the key genes involved in hormone signal transduction and the potential TFs regulating this pathway during dormancy release, we performed co-expression network analysis of DEGs involved in hormone signal transduction and differentially expressed TFs. In total, 187 differentially expressed TFs were identified within subclusters 1, 3, 4, and 5. Cytoscape was used to visualize the results (r >0.9 or <−0.9, P <0.05, Figure 5). The genes involved in the IAA and ABA signaling pathways were activated from 31 October to 22 November (Figure 5A). According to the gene connectivity, the top three genes in IAA signaling pathway were AUX22B (TRINITY_DN2156_c0_g1), IAA9 (TRINITY_DN12949_c0_g1), and IAA27 (TRINITY_DN1056_c0_g2). The top three genes in the ABA signaling pathway were PYL6 (TRINITY_DN51499_c0_g1), PYL4-like (TRINITY_DN38780_c0_g2), and PYL1 (TRINITY_DN20425_c0_g1). Additionally, DELLA (TRINITY_DN45150_c0_g1) in the GA signaling pathway also had a high level of connectivity. The AP2/ERF (18), WRKY (16), and zinc finger (8) families account for a large number of genes involved in hormone signal transduction. Among them, WRKY11 and WRKY28 within the WRKY family showed the same number of interaction pairs, and WRKY family genes were mainly co-expressed with IAA and ABA signaling pathway genes.

Figure 5
www.frontiersin.org

Figure 5. Co-expression network analysis of genes involved in hormone signal transduction and TFs. (A) Correlation network of structural genes involved in GA, ABA, IAA, and CTK signaling pathways and TFs in KEGG analysis during the process of endodormancy release. (B) Correlation network of structural genes involved in GA, ABA, IAA, and CTK signaling pathways and TFs in KEGG analysis during the transition from endodormancy to ecodormancy (r >0.9 or r <−0.9, P <0.05). The size of rectangles and circles indicates the degree of connectivity. The greater the degree of connectivity, the larger the shapes. The solid red line represents positive correlation and the dashed blue line represents negative correlation.

During the transition from endo- to ecodormancy, GID1C (TRINITY_DN37598_c0_g3), PIF4-like (TRINITY_DN2672_c0_g1), ARR-A (TRINITY_DN2166_c0_g1), SAPK3 (TRINITY_DN3171_c0_g2), AI5L2 (TRINITY_DN6632_c0_g1), and IAA8-like (TRINITY_DN12949_c0_g3) were highly correlated (Figure 5B). However, GID1C did not show high connectivity during endodormancy release (Figure 5A). This result suggests that genes from the same hormone signaling pathway may function at different stages of the transition from endo- to ecodormancy. Furthermore, WRKY, AP2/ERF, bHLH, and zinc finger were the key TF families regulating hormone signal transduction, with 12, 9, 6, and 6 annotated genes, respectively.

3.6 Analyzing the weighted gene correlation network

Metabolomic and transcriptomic analyses indicated that starch and sucrose metabolism was enriched during the transition from endo- to ecodormancy. We detected significant changes in the contents of soluble sugar, sucrose, and starch in underground buds. Additionally, three DAMs involved in starch and sucrose metabolism were identified through metabolomic analysis, including amylose, levan (a homopolysaccharide that contains fructose), and UDP-glucose. Thus, to obtain insights into the regulation of sugars, WGCNA was employed to identify DEGs with similar expression patterns related to six metabolites, including soluble sugar, sucrose, starch, UDP-glucose, amylose, and levan. DEGs with low expression (averaged FPKM <1) were filtered out, resulting in a total of 11,229 DEGs. A total of 21 distinct modules were identified with similar expression patterns (Figure 6A). Based on the Pearson correlation coefficients among the analyzed traits and gene modules, seven modules that were significantly associated with sugar content were identified. The purple module was positively correlated with the soluble sugar sucrose, and UDP-glucose contents (r = 0.786, 0.828, and 0.965, respectively) (Figure 6B). The red module was positively correlated with starch and amylose content (r = 0.79 and 0.616). The magenta module exhibited a positive relationship with the levan content (r = 0.826). These results suggest that the genes within these modules positively regulate the accumulation of these six sugars. Conversely, the blue module was negatively correlated with soluble sugar and levan contents (−0.721 and −0.78, respectively). The genes in the magenta, yellow, and green-yellow modules were highly negative correlated with the contents of starch, sucrose, and amylose (r = −0.868, −0.906, and −0.614), respectively. This suggests that these genes are involved in the negative regulation of sugar accumulation. In addition, we found no structural genes involved in starch and sucrose metabolism in the light yellow module. Consequently, genes from these six modules were selected for further analysis. KEGG enrichment analysis revealed that the DEGs were significantly enriched in starch and sucrose metabolism, and plant hormone signal transduction (Figure 6C).

Figure 6
www.frontiersin.org

Figure 6. WGCNA analysis of transcriptome data during dormancy transition in underground buds of AMM. (A) Dendrogram showing the co-expression modules of 16832 DEGs identified during the dormant phase transition through WGCNA, with the major tree branches consisting of 21 modules labeled with different colors. (B) Heat map showing the correlation between modules and metabolites in starch and sucrose metabolism pathways. (C) Top 20 of KEGG enrichment analysis of DEGs in the co-expression modules related to the sugar contents. Each column corresponds to a time point and each row corresponds to a module labeled with different colors. Pink and green represent positive and negative correlations, respectively.

3.7 Correlation analysis of DEGs and metabolites related to starch and sucrose metabolism

DEGs involved in starch and sucrose metabolism from the six modules related to sugar content were analyzed in association with the metabolites (Figure 7). A total of 36 DEGs involved in sucrose and starch metabolism were identified in the six modules (Supplementary Table S10). Metabolomic analysis revealed the presence of three DAMs involved in starch and sucrose metabolism: amylose, levan, and UDP-glucose. The glucan endo-1,3-beta-glucosidases (GNs) (TRINITY_DN3351_c0_g1, TRINITY_DN3351_c0_g3) were downregulated during the transition from endo- to ecodormancy, which was positively correlated with the content of amylose and starch, and negatively correlated with the levels of levan (Figures 7A, B). The expression levels of alpha-glucosidase (malZ, TRINITY_DN8503_c0_g1) and beta-glucosidase (BGL18, TRINITY_DN2095_c0_g1), which were negatively correlated with the contents of UDP-glucose and sucrose, were downregulated during endodormancy release and reached their lowest levels during ecodormancy. The level of BGL46-like (TRINITY_DN3816_c0_g1) decreased continuously during the transition from endo- to ecodormancy and was positively correlated with starch levels.

Figure 7
www.frontiersin.org

Figure 7. Starch and sucrose metabolism during the transition from endo- to ecodormancy in underground buds. (A) The heatmaps indicate the expression profiles of DEGs and metabolites involved in sucrose and starch metabolism during the transition from endo- to ecodormancy in the underground buds of AMM. (B) Correlation analysis between metabolites and structural genes involved in starch and sucrose metabolism. 31 October–22 November represents the endodormancy release stage, 31 October–1 December represents the transition from endo- to ecodormancy stage. The red, blue, and black fonts indicate accumulated differential metabolites, undifferentiated metabolites, and undetected metabolites, respectively. The ‘*’, ‘**’, and ‘***’ represent significant differences at the P <0.05, P <0.01, and P <0.001 levels according to Pearson’s correlation coefficient, respectively.

To better understand the key genes and TFs regulating starch and sucrose metabolism, differentially expressed TFs within the six modules were also screened, and a co-expression network of differential metabolites with DEGs and TFs was constructed (Figure 8). A total of 57 DEGs were highly correlated with four metabolites during endodormancy release (r >0.9 or <−0.9, P <0.05, (Figure 8A). The number of DEGs associated with amylose was the highest (28), and was positively correlated with amylose, including endoglucanase (EG11, TRINITY_DN51949_c0_g1) and two GNs. Additionally, 10 DEGs were related to UDP-glucose, and only BGL46-like was related to starch. The zinc finger (9), MYB (8), and TCP (5) families are likely to be TFs regulating starch and sucrose metabolism (Figure 8A). Furthermore, the two GNs were co-expressed with COL2 (TRINITY_DN8556_c0_g1) and DRE1F (TRINITY_DN3958_c0_g1), and these genes were positively related to amylose. These results demonstrated that COL2 and DRE1F may regulate the expression of GNs and, consequently, affect the amylose content which is consistent with the decrease in GN expression and amylose content during dormancy release (Figure 7B). During the transition from endo- to ecodormancy, 51 DEGs were identified to correlate with five metabolites (r >0.85 or <−0.85, P <0.05), and most of these DEGs were negatively correlated with metabolite contents (Figure 8B). In addition, five structural genes were identified: BGL18, BGL46-like, isoamylase 2 (ISA2, TRINITY_DN13858_c0_g1), malZ, and glgC. Compared to the endodormancy release state, BGL18, ISA2, and malZ may play an important role in endodormancy release. Furthermore, during the transition from endo- to ecodormancy, the zinc finger, ERF, TCP, and NAC families could be the key TFs involved in the regulation of starch and sucrose metabolism. Interestingly, we also identified that some MADS-box genes, AP1, SOC1, and AGL19, were all related to the concentration of sucrose, which suggested that MADS-box genes were mainly involved in regulating sucrose metabolism.

Figure 8
www.frontiersin.org

Figure 8. Correlation network of differential metabolites, structural genes, and TFs involved in starch and sucrose mechanism. (A) Co-expression networks of metabolites, structural genes, and TFs involved in starch and sucrose metabolic pathways during the process of endodormancy release. (B) Co-expression networks of metabolites, structural genes, and TFs involved in starch and sucrose metabolic pathways during the transition from endo- to ecodormancy. Hexagons and diamonds represent differential types of metabolites and the structural genes in starch and sucrose metabolism, respectively. TF families are represented by circles, and the color of the circles correspond to TFs of different families. The solid red line represents a positive correlation and the dashed blue line represents a negative correlation.

In addition, we identified genes involved in hormone signal transduction that regulates starch and sucrose metabolism. For example, a positive correlation between ARF5 (TRINITY_DN4693_c0_g1) and glgB as well as UDP-glucose was identified. Additionally, there was a correlation between glgB and UDP-glucose (Figure 8A), indicating that ARF5 may regulate the accumulation of UDP-glucose by positively regulating glgB expression. More hormone signal pathway-related genes were co-expressed with metabolites and structural genes related to starch and sucrose metabolism during the transition from endo- to ecodormancy (Figure 8B). PIF4-like (TRINITY_DN2672_c0_g1) negatively correlated with glgC and sucrose, whereas glgC positively correlated with sucrose. GAI (TRINITY_DN83684_c0_g1) also correlated with sucrose and malZ, and GAI and malZ were co-expressed. These findings indicate that genes related to the GA signaling pathway are involved in sucrose metabolism following endodormancy release. Overall, these results suggest that genes involved in hormone signal transduction also influence starch and sucrose metabolism by regulating the levels of UDP-glucose and sucrose.

3.8 Comprehensive analysis of energy metabolism and amino acid metabolism by transcriptome and metabolome

Omics analysis showed that energy and amino acid metabolism were activated during dormancy release in AMM. Therefore, we proposed a metabolic network including glycolysis, TCA cycle, shikimic acid pathway, and amino acid metabolism to visualize carbon flow (Figure 9; Supplementary Table S11). During endodormancy release, the expression of most genes involved in glycolysis was upregulated and peaked on 15 November. This indicated that glycolysis is highly active, preparing for the complete release of endodormancy. Genes associated with the TCA cycle exhibited similar expression patterns. Most genes were upregulated during endodormancy release and downregulated during ecodormancy. On 8 November, the ketoglutarate content dropped significantly but then gradually increased until endodormancy was fully lifted. In contrast, malate levels, which are the downstream products of ketoglutarate, showed opposite trends, aligning with the patterns of gene expression. The results showed that during endodormancy release, genes associated with energy metabolism reacted quickly, resulting in active metabolic activity and providing energy for transition into an ecodormancy state. Regarding amino acid metabolism, the contents of serine and isoleucine decreased with endodormancy release, while the contents of arginine, lysine, valine, tryptophan, and tyrosine increased. The expression patterns of DEGs involved in amino acid metabolism also showed a similar trend.

Figure 9
www.frontiersin.org

Figure 9. Combined analysis of the carbohydrate metabolism and amino acid metabolism pathways during the dormancy transition. (A) Metabolites involved in the carbohydrate metabolism and amino acid metabolism pathways. (B) Heatmap of the DEGs involved in the carbohydrate metabolism and amino acid metabolism pathways. The red, blue, and black fonts indicate accumulated differential metabolites, undifferentiated metabolites, and undetected metabolites, respectively.

4 Discussion

4.1 The combination of morphological analysis, transcriptome and metabolome provides the basis for determining the dormancy states

It is essential to have a sufficient CR for the transition from endo- to ecodormancy. We found that Astragalus membranaceus (Fisch.) Bge. var. mongholicus (Bge.) Hsiao (AMM) did not sprout completely on 31 October through the observation of several morphological indices (DFS, DBS, BPF, and DAS). The sprouting rates on 22 November reached 100% for the first time 4 weeks after transplanting, and several indices remained stable. Based on the research by Wang et al. (2022a) on the CR of P. lactiflora, we speculated that AMM emerged from endodormancy on 22 November. In addition, cluster analysis showed that differential metabolites accumulated significantly in underground buds collected on 31 October, 22 November, and 1 December (Figure 2D). These results indicate that the underground bud status of 31 October, 15 November, and 1 December may be different. Interestingly, unlike the accumulation pattern of DAMs, the number of DEGs was the highest on 15 November (Figure 3C). When comparing AMM that was transferred into the greenhouse on 8 November with AMM transferred on 15 November, we observed a significant increase in the BPF for the latter. However, there were no significant changes in DBS of the AMM transferred on 15 and 22 November. We speculate that, at this time, AMM may be in the early stages of fully breaking endodormancy, as it appears to have reached the minimum CRs for normal growth. This stage triggers the expression of a large number of genes, prepares the plant to completely break dormancy, and promotes various biosynthetic, degradative, and metabolic processes (Zhang et al., 2021b). In summary, it is determined that it is in endodormancy stage on 31 October, and 15–22 November is the stage when AMM underground buds completely break endodormancy. This study provides guidance for determining the period of dormancy release of AMM and clarifying the CR for dormancy release.

4.2 Endogenous hormones are actively involved in the dormancy transition of underground buds of AMM

The dormancy of perennial plants is primarily regulated by hormones that participate in the formation of a complex regulatory network during dormancy release (Liu and Sherif, 2019; Pan et al., 2021). IAA is an important hormone for plant growth and development and is involved in the regulation of bud dormancy (Hao et al., 2019). In this study, the IAA levels in underground buds were significantly higher than those of the other three hormones (Figure 1A), and a large number of DEGs related to IAA signaling pathways were identified (Figure 4). Similar findings have been reported in the Chinese fir (C. lanceolata) (Qiu et al., 2013) and plum (Prunus mume) (Zhang et al., 2018). Furthermore, co-expression network analysis showed that DEGs in the auxin signaling pathway exhibited greater connectivity than genes related to other hormone signaling pathways during endodormancy release (Figure 5A). These results highlight the key role of IAA in the endodormancy release of AMM underground buds. In most species, ABA and GA are considered key regulators of bud dormancy (Li et al., 2020; Zhang et al., 2018; Gao et al., 2023). The levels of ABA and GA also changed during underground bud dormancy release in AMM. Contrary to the common belief that a reduction in ABA content triggers dormancy release, our results showed that ABA levels remained high during endodormancy release, as observed in the upregulation of the biosynthesis gene NCEDs and the low level of the degradation gene CYP707A2 (Supplementary Figure S3), which was consistent with the results of Wang et al. (2019). This indicates that high levels of ABA may positively influence release from bud dormancy. However, further research is needed to verify the role of ABA in this process in AMM underground buds. Similar regulatory patterns were observed for GA synthesis genes GA20ox and GA3ox, and degradation gene GA2ox during the process, which allowed GA3 to maintain an upward trend throughout the release of endodormancy, which is consistent with the results observed in tree peony (Yuan et al., 2024).

In perennial plants, ABA, GA, IAA, and CK interactions regulate bud dormancy (Lv et al., 2018; Tosetti et al., 2021). The combined application of CTK and GA can effectively break the dormancy of tea buds within a short period (Thirugnanasambantham et al., 2020). This study found that IAA, ABA, GA, and CTK levels exhibited consistent trends during endodormancy release (Figure 1A). Additionally, co-expression network analysis revealed that genes associated with IAA, ABA, GA, and CTK signaling pathways were co-expressed during endodormancy release, with IAA signaling pathway genes showing stronger connectivity (Figure 5). In summary, we speculated that these hormones may form a synergistic network in regulating endodormancy release, with IAA potentially playing a leading role in this process. However, further research is needed to fully understand how these four hormones work together to regulate dormancy.

4.3 Carbohydrates participate in the regulation of the bud dormancy transition of AMM

It is well established that carbohydrates are closely related to dormancy (Gillespie and Volaire, 2017). Perennials’ dormant buds appear to be in an inactive state with limited metabolic activity, and reconfiguration of carbohydrate metabolism is associated with dormancy transition. The decrease in starch content and increase in soluble sugars trigger dormancy release (Marković et al., 2020). In our study, there was a significant decrease in the starch content at the onset of endodormancy release, followed by a gradual increase before dormancy was fully released. Subsequently, the starch content began to decline again, which was closely associated with the dormancy stage (Figure 2B). In addition, the level of levan increased during endodormancy release and reached a high level when endodormancy was released (Figure 1A), indicating that starch and levan could serve as potential markers for distinguishing the stages of dormancy (Ohanenye et al., 2019). As starch levels decreased, soluble sugar levels increased until ecodormancy was achieved. This observation aligns with metabolomic analyses showing that significantly abundant DAMs are enriched in starch and sucrose metabolism during ecodormancy. In addition, genes related to starch and sucrose metabolism were highly expressed during endodormancy release. We speculated that asynchrony in the expression of genes involved in starch and sucrose metabolism and metabolite accumulation may contribute to the promotion of endodormancy release. Genes involved in the synthesis and degradation of carbohydrates such as starch and sucrose were actively expressed. This activity was intended to provide adequate resources for respiration and energy metabolism. Therefore, there were few significant changes in metabolite content until endodormancy release. During the ecodormancy stage, metabolite levels peak in preparation for the subsequent regeneration stage (Zhang et al., 2020). However, more evidence is required to support this hypothesis.

Co-expression network analysis showed that GNs, EG, BGL, glgB, and glgC were associated with the differential metabolism of sugars during endodormancy release. We showed that GN plays an important role in starch degradation and the response to biotic and abiotic stresses (Nakatani et al., 2010). In our study, GN expression levels were highly correlated with starch and amylose during endodormancy release, and they had the highest levels of TF co-expression (Figure 7B). However, this correlation no longer exists following endodormancy release. We hypothesized that GNs may be important in regulating starch degradation during endodormancy release (Gao et al., 2021). Zeng et al. (2021) found that the decreased expression of BGL may be related to a decrease in sucrose content. In this study, most BGLs maintained a high level during endodormancy release, which was positively correlated with the sucrose content. Similar findings have also been reported for rhubarb (Wojtania et al., 2022) and P. lactiflora (Zhang et al., 2021a). In conclusion, GNs and BGL may play key roles in the release of endodormancy by regulating starch and sucrose metabolism.

Some studies have confirmed a relationship between carbohydrates and hormones during dormancy release (Wang et al., 2022a). This study found that the trends of the four hormone levels changed in the opposite direction to that of the starch content. Additionally, SNRK2 in the ABA signaling pathway was correlated with sucrose content during the transition from endo- to ecodormancy (Figure 8B). Additionally, the trend of the expression patterns of most SNRK2 genes was similar to that of the changes in ABA content. It has been speculated that the crosstalk between sucrose and ABA is mediated by SNRK2. AMY, which plays a key role in dormancy, also regulates the interaction between hormones and sugars (Xie et al., 2007). GA upregulates AMY expression through GAMYB, thereby promoting starch degradation (Li et al., 2020). The variation in GA3 levels in the underground buds of AMM during endodormancy release was opposite to the change in starch content. The expression pattern of GAMYB (TRINITY_DN15272_c0_g1) was consistent with GA3 levels, and GAMYB was positively and negatively correlated with GA3 and starch content, respectively (Supplementary Table S12). Therefore, we hypothesized that GA participates in starch degradation by GAMYB during the transition from endo- to ecodormancy (Chao et al., 2006).

The physiological activities involved in the transition to bud dormancy require significant energy and material resources, primarily obtained from the breakdown of stored substances. Metabolomic analysis revealed that, during the ecological dormancy stage, DAMs were notably enriched in starch and sucrose metabolism. This indicated that soluble sugars gradually accumulated, suggesting that the increase in soluble sugar content may fulfill the energy and material requirements for bud regrowth. In this study, most of the DEGs involved in glycolysis and the TCA cycle were upregulated during endodormancy release (Figure 9). This finding aligns with previous observations made during bud dormancy in peaches (Tang et al., 2019) and tree peonies (Zhang et al., 2020). This indicated that the activation of energy metabolism occurs through the upregulation of these genes during the endodormancy release stage, allowing them to serve as the primary energy source. Glycolysis and TCA cycle can also provide precursors for amino acid synthesis during budbreak. Metabolomic analysis showed that DAMs were significantly enriched in amino acid metabolism during the transition from endo- to ecodormancy, which is similar to the findings of Cheng et al. (2023). These results indicate that AMM may respond to endodormancy release through high levels of amino acids and their derivatives. Therefore, we speculated that amino acid metabolism might be necessary for endodormancy release in AMM.

4.4 The role of TFs in the release of dormancy

An increasing number of studies have demonstrated that the molecular mechanisms regulating dormancy are mediated by complex genetic networks composed of many TFs, such as WRKY, MYB, bZIP, NAC, and bHLH (Wang et al., 2018; Nakatani et al., 2010). Based on the co-expression network analysis, it was found that the dormancy of AMM underground buds was related to some key TF families, including WRKY, AP2/ERF, MYB, bHLH, zinc finger, HD-ZIP, and MADS-box. Interestingly, WRKY families were not screened for TFs regulating starch and sucrose metabolism, but Fan et al. (2021) found that WRKY32 and WRKY71 were closely related to sucrose metabolism during the dormancy transition of lily bulbs, which may vary from species to species. Moreover, our studies have shown that DREB1, DREB2, and DREB3 are co-expressed with genes associated with the IAA and ABA signaling pathways during endodormancy release. This result implies a significant role of DREB in the endodormancy release process, which aligns with the findings of Artlip et al. (2019). Genes encoding MADS-boxes are highly correlated with bud dormancy in many species (Wu et al., 2017), among which SOC1 has been proven to be crucial for breaking bud dormancy (Gómez-Soto et al., 2021). In this study, SOC1 was highly associated with PIF4-like proteins and sucrose during the transition from endo- to ecodormancy, indicating that SOC1 may be essential for regulating dormancy by simultaneously controlling the GA signaling pathway and sucrose accumulation (Yang et al., 2021). In addition, the flowering regulator AP1 is also involved in the expression of IAA signaling pathway genes and sucrose accumulation during dormancy release. These results suggest that the regulatory network mediated by TFs, including WRKY, AP2/ERF, bHLH, zinc finger, and MADS-box, has a potential function in regulating dormancy release in underground buds of AMM.

5 Conclusion

In this study, the regulation mechanism of underground bud dormancy of A. membranaceus (Fisch.) Bge. var. mongholicus (Bge.) Hsiao (AMM) was studied by a multi-omics approach for the first time. We comprehensively analyzed the metabolic and transcriptome data of the underground buds of one-year-old AMM at five successive dormancy stages and proposed a model for regulating the dormancy transition of underground buds of AMM induced by chilling (Figure 10). IAA may be a key hormone in the regulation of endodormancy release. Genes related to the IAA, ABA, and GA signaling pathways (IAA, ARF, SNRK2, and PIF4) were involved in regulating starch and sucrose metabolism during the transition from endo- to eco-dormancy. In addition to starch and sucrose, levan may also be used as a marker of dormancy. Moreover, the structural genes, including GN, EG, and glgB, for starch and sucrose metabolism may be associated with endodormancy release, whereas ISA and malZ are involved in regulating starch and sucrose metabolism following endodormancy release. Through co-expression network analysis, we found that TFs of the ERF, MYB, bHLH, zinc finger, and MADS-box families may control bud dormancy by regulating the expression of genes involved in hormone signal transduction and carbohydrate metabolism. Although the WRKY family only regulates hormone signal transduction, the interactions among molecules, cells, and metabolism need to be further explored. We hope that our study can contribute to a better understanding of the regulatory mechanism of the transition from endo- to ecodormancy in AMM and provide a new strategy for breaking dormancy and shortening the breeding time in advance.

Figure 10
www.frontiersin.org

Figure 10. A proposed model for regulating the transition from endo- to ecodormancy in underground buds of AMM induced by chilling. The blue circles represent genes that may play a role in the process of endodormancy release, the purple circles represent genes that may be active following endodormancy release, the yellow circles represent genes that are co-expressed with genes and metabolites related to starch and sucrose metabolism.

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

Author contributions

HG: Conceptualization, Data curation, Investigation, Methodology, Writing – original draft, Writing – review & editing. YZ: Data curation, Investigation, Methodology, Writing – review & editing. QC: Writing – review & editing. QZ: Writing – review & editing. PY: Investigation, Writing – review & editing. SS: Conceptualization, Methodology, Writing – review & editing. GC: Conceptualization, Funding acquisition, Methodology, Resources, 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 supported by grants from the National Natural Science Foundation of China (81660630) and Key R&D and Achievement Transformation Projects in the Inner Mongolia Autonomous Region (2023YFHH0105).

Acknowledgments

We would like to express our gratitude to Zhihe Liu, Yifan Dong, Ruijuan Wang, Zhi Wang, and Ding Zhang for their assistance with sample collection and transplantation.

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.

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.1483538/full#supplementary-material

References

André, D., Marcon, A., Lee, K. C., Goretti, D., Zhang, B., Delhomme, N., et al. (2022). Flowering Locus T paralogs control the annual growth cycle in Populus trees. Curr. Biol. 32, 2988–2996.e4. doi: 10.1016/j.cub.2022.05.023

PubMed Abstract | Crossref Full Text | Google Scholar

Artlip, T., McDermaid, A., Ma, Q., Wisniewski, M. (2019). Differential gene expression in non-transgenic and transgenic "M.26" apple overexpressing a peach CBF gene during the transition from eco-dormancy to bud break. Hortic. Res. 6, 86. doi: 10.1038/s41438-019-0168-9

PubMed Abstract | Crossref Full Text | Google Scholar

Chao, W. S., Doğramaci, M., Horvath, D. P., Anderson, J. V., Foley, M. E. (2016). Phytohormone balance and stress-related cellular responses are involved in the transition from bud to shoot growth in leafy spurge. BMC Plant Biol. 16, 47. doi: 10.1186/s12870-016-0735-2

PubMed Abstract | Crossref Full Text | Google Scholar

Chao, W. S., Doğramacı, M., Horvath, D. P., Anderson, J. V., Foley, M. E. (2017). Comparison of phytohormone levels and transcript profiles during seasonal dormancy transitions in underground adventitious buds of leafy spurge. Plant Mol. Biol. 94, 281–302. doi: 10.1007/s11103-017-0607-7

PubMed Abstract | Crossref Full Text | Google Scholar

Chao, W. S., Serpe, M. D., Anderson, J. V., Gesch, R. W., Horvath, D. P. (2006). Sugars, hormones, and environment affect the dormancy status in underground adventitious buds of leafy spurge (Euphorbia esula). Weed Sci. 54, 59–68. doi: 10.1614/WS-05-088R.1

Crossref Full Text | Google Scholar

Charrier, G. (2022). Is winter coming? Minor effect of the onset of chilling accumulation on the prediction of endodormancy release and budbreak. Physiologia plantarum 174, e13699. doi: 10.1111/ppl.13699

PubMed Abstract | Crossref Full Text | Google Scholar

Chen, C., Zeng, L., Zhao, H., Ye, Q. (2020). Proteomic analysis of the early development of the phalaenopsis amabilis flower bud under low temperature induction using the iTRAQ/MRM approach. Molecules. 25, 1244. doi: 10.3390/molecules25051244

PubMed Abstract | Crossref Full Text | Google Scholar

Chen, J., Wu, X. T., Xu, Y. Q., Zhong, Y., Li, Y. X., Chen, J. K., et al. (2015). Global transcriptome analysis profiles metabolic pathways in traditional herb Astragalus membranaceus Bge. var. mongolicus (Bge.) Hsiao. BMC Genomics 16, S15. doi: 10.1186/1471-2164-16-S7-S15

PubMed Abstract | Crossref Full Text | Google Scholar

Cheng, Y., Ban, Q., Mao, J., Lin, M., Zhu, X., Xia, Y., et al. (2023). Integrated metabolomic and transcriptomic analysis reveals that amino acid biosynthesis may determine differences in cold-tolerant and cold-sensitive tea cultivars. Int. J. Mol. Sci. 24, 1907. doi: 10.3390/ijms24031907

PubMed Abstract | Crossref Full Text | Google Scholar

Dong, M., Li, J., Yang, D., Li, M., Wei, J. (2023). Biosynthesis and pharmacological activities of flavonoids, triterpene saponins and polysaccharides derived from Astragalus membranaceus. Molecules. 28, 5018. doi: 10.3390/molecules28135018

PubMed Abstract | Crossref Full Text | Google Scholar

Du, Y., Zhao, Q., Chen, L., Yao, X., Zhang, W., Zhang, B., et al. (2020). Effect of drought stress on sugar metabolism in leaves and roots of soybean seedlings. Plant Physiol. Biochem. 146, 1–12. doi: 10.1016/j.plaphy.2019.11.003

PubMed Abstract | Crossref Full Text | Google Scholar

Fadón, E., Fernandez, E., Behn, H., Luedeling, E. (2020). A conceptual framework for winter dormancy in deciduous trees. Agronomy. 10, 241. doi: 10.3390/agronomy10020241

Crossref Full Text | Google Scholar

Fan, X., Yang, Y., Li, M., Fu, L., Zang, Y., Wang, C., et al. (2021). Transcriptomics and targeted metabolomics reveal the regulatory network of Lilium davidii var. unicolor during bulb dormancy release. Planta. 254, 59. doi: 10.1007/s00425-021-03672-7

PubMed Abstract | Crossref Full Text | Google Scholar

Gao, L., Niu, D., Chi, T., Yuan, Y., Liu, C., Gai, S., et al. (2023). PsRGL1 negatively regulates chilling- and gibberellin-induced dormancy release by PsF-box1-mediated targeting for proteolytic degradation in tree peony. Horticulture Res. 10, uhad044. doi: 10.1093/hr/uhad044

PubMed Abstract | Crossref Full Text | Google Scholar

Gao, Y., Yang, Q., Yan, X., Wu, X., Yang, F., Li, J., et al. (2021). High-quality genome assembly of 'Cuiguan' pear (Pyrus pyrifolia) as a reference genome for identifying regulatory genes and epigenetic modifications responsible for bud dormancy. Hortic. Res. 8, 197. doi: 10.1038/s41438-021-00632-w

PubMed Abstract | Crossref Full Text | Google Scholar

Gillespie, L. M., Volaire, F. A. (2017). Are winter and summer dormancy symmetrical seasonal adaptive strategies? The case of temperate herbaceous perennials. Ann. Bot. 119, 311–323. doi: 10.1093/aob/mcw264

PubMed Abstract | Crossref Full Text | Google Scholar

Gómez-Soto, D., Ramos-Sánchez, J. M., Alique, D., Conde, D., Triozzi, P. M., Perales, M., et al. (2021). Overexpression of a SOC1-related gene promotes bud break in ecodormant poplars. Front. Plant Sci. 12. doi: 10.3389/fpls.2021.670497

PubMed Abstract | Crossref Full Text | Google Scholar

Hao, X., Tang, H., Wang, B., Wang, L., Cao, H., Wang, Y., et al. (2019). Gene characterization and expression analysis reveal the importance of auxin signaling in bud dormancy regulation in tea plant. J. Plant Growth Regul. 38, 225–240. doi: 10.1007/s00344-018-9834-7

Crossref Full Text | Google Scholar

Horikoshi, H. M., Sekozawa, Y., Kobayashi, M., Saito, K., Kusano, M., Sugaya, S. (2018). Metabolomics analysis of 'Housui' Japanese pear flower buds during endodormancy reveals metabolic suppression by thermal fluctuation. Plant Physiol. Biochem. 126, 134–141. doi: 10.1016/j.plaphy.2018.02.028

PubMed Abstract | Crossref Full Text | Google Scholar

Horvath, D. P., Anderson, J. V., Chao, W. S., Foley, M. E. (2003). Knowing when to grow: Signals regulating bud dormancy. Trends Plant Sci. 8, 534–540. doi: 10.1016/j.tplants.2003.09.013

PubMed Abstract | Crossref Full Text | Google Scholar

Horvath, D. P., Sung, S., Kim, D., Chao, W., Anderson, J. (2010). Characterization, expression and function of DORMANCY ASSOCIATED MADS-BOX genes from leafy spurge. Plant Mol. Biol. 73, 169–179. doi: 10.1007/s11103-009-9596-5

PubMed Abstract | Crossref Full Text | Google Scholar

Ito, A., Tuan, P. A., Saito, T., Bai, S., Kita, M., Moriguchi, T. (2021). Changes in phytohormone content and associated gene expression throughout the stages of pear (Pyrus pyrifolia Nakai) dormancy. Tree Physiol. 41, 529–543. doi: 10.1093/treephys/tpz101

PubMed Abstract | Crossref Full Text | Google Scholar

Ji, B., Xuan, L., Zhang, Y., Zhang, G., Meng, J., Mu, W., et al. (2023). Advances in biotechnological production and metabolic regulation of Astragalus membranaceus. Plants (Basel). 12, 1858. doi: 10.3390/plants12091858

PubMed Abstract | Crossref Full Text | Google Scholar

Jiang, N. N., Xiang, Z. Q., Wang, Y., Sun, Y., Fang, Y. F., Xu, J. G. (2020). Effects of GA3 on dormancy release, endogenous hormones levels and sugar metabolism in Paeonia lactiflora ‘Da Fugui’. J. Nanjing Forestry Univ. 44, 26–32. doi: 10.3969/j.issn.1000-2006.201909019

Crossref Full Text | Google Scholar

Khalil-Ur-Rehman, M., Wang, W., Xu, Y. S., Haider, M. S., Li, C. X., Tao, J. M. (2017). Comparative study on reagents involved in grape bud break and their effects on different metabolites and related gene expression during winter. Front. Plant Sci. 8. doi: 10.3389/fpls.2017.01340

PubMed Abstract | Crossref Full Text | Google Scholar

Ladwig, L. M., Chandler, J. L., Guiden, P. W., Henn, J. J. (2019). Extreme winter warm event causes exceptionally early bud break for many woody species. Ecosphere. 10, e02542. doi: 10.1002/ecs2.2542

Crossref Full Text | Google Scholar

Lang, G. A., Early, J. D., Martin, G. C., Darnell, R. L. (1987). Endo-, para-, and ecodormancy: physiological terminology and classification for dormancy research. HortScience. 22, 371–377. doi: 10.21273/HORTSCI.22.3.371

Crossref Full Text | Google Scholar

Lempe, J., Moser, M., Asquini, E., Si-Ammour, A., Flachowsky, H. (2024). Functional evidence on the involvement of the MADS-box gene MdDAM4 in bud dormancy regulation in apple. Front. Plant Sci. 15. doi: 10.3389/fpls.2024.1433865

PubMed Abstract | Crossref Full Text | Google Scholar

Li, X., Fei, R., Chen, Z., Fan, C., Sun, X. (2020). Plant hormonal changes and differential expression profiling reveal seed dormancy removal process in double dormant plant-herbaceous peony. PloS One 15, e0231117. doi: 10.1371/journal.pone.0231117

PubMed Abstract | Crossref Full Text | Google Scholar

Li, Y., Guo, S., Zhu, Y., Yan, H., Qian, D. W., Wang, H. Q., et al. (2019a). Comparative analysis of twenty-five compounds in different parts of Astragalus membranaceus var. mongholicus and Astragalus membranaceus by UPLC-MS/MS. J. Pharm. Anal. 9, 392–399. doi: 10.1016/j.jpha.2019.06.002

PubMed Abstract | Crossref Full Text | Google Scholar

Li, Y., Ma, R., Li, R., Zhao, Q., Zhang, Z., Zong, Y., et al. (2022a). Comparative transcriptomic analysis provides insight into the key regulatory pathways and differentially expressed genes in blueberry flower bud endo- and ecodormancy release. Horticulturae. 8, 176. doi: 10.3390/horticulturae8020176

Crossref Full Text | Google Scholar

Li, D., Shao, L., Zhang, J., Wang, X., Zhang, D., Horvath, D. P., et al. (2022b). MADS-box transcription factors determine the duration of temporary winter dormancy in closely related evergreen and deciduous Iris spp. J. Exp. Bot. 73, 1429–1449. doi: 10.1093/jxb/erab484

PubMed Abstract | Crossref Full Text | Google Scholar

Li, S., Wang, Q., Wen, B., Zhang, R., Jing, X., Xiao, W., et al. (2021). Endodormancy release can be modulated by the GA4-GID1c-DELLA2 module in peach leaf buds. Front. Plant Sci. 12. doi: 10.3389/fpls.2021.713514

PubMed Abstract | Crossref Full Text | Google Scholar

Li, J., Xu, Y., Niu, Q., He, L., Teng, Y., Bai, S. (2018). Abscisic Acid (ABA) promotes the induction and maintenance of pear (Pyrus pyrifolia white pear group) flower bud endodormancy. Int. J. Mol. Sci. 19, 310. doi: 10.3390/ijms19010310

PubMed Abstract | Crossref Full Text | Google Scholar

Li, Y., Zhang, H., Zhang, Q., Liu, Q., Zhai, H., Zhao, N., et al. (2019b). An AP2/ERF gene, IbRAP2-12, from sweetpotato is involved in salt and drought tolerance in transgenic Arabidopsis. Plant Sci. 281, 19–30. doi: 10.1016/j.plantsci.2019.01.009

PubMed Abstract | Crossref Full Text | Google Scholar

Liu, Y., Dang, P., Liu, L., He, C. (2019). Cold acclimation by the CBF-COR pathway in a changing climate: Lessons from Arabidopsis thaliana. Plant Cell Rep. 38, 511–519. doi: 10.1007/s00299-019-02376-3

PubMed Abstract | Crossref Full Text | Google Scholar

Liu, X., Hou, X. (2018). Antagonistic Regulation of ABA and GA in metabolism and signaling pathways. Front. Plant Sci. 9. doi: 10.3389/fpls.2018.00251

PubMed Abstract | Crossref Full Text | Google Scholar

Liu, J., Sherif, S. M. (2019). Hormonal orchestration of bud dormancy cycle in deciduous woody perennials. Front. Plant Sci. 10. doi: 10.3389/fpls.2019.01136

PubMed Abstract | Crossref Full Text | Google Scholar

Liu, X., Wang, Q., Gu, J., Lü, Y. (2014). Vernalization of Oriental hybrid lily ‘Sorbonne’: changes in physiology metabolic activity and molecular mechanism. Mol. Biol. Rep. 41, 6619–6634. doi: 10.1007/s11033-014-3545-3

PubMed Abstract | Crossref Full Text | Google Scholar

Liu, X., Wang, J., Sabir, I. A., Sun, W., Wang, L., Xu, Y., et al. (2022). PavGA2ox-2L inhibits the plant growth and development interacting with PavDWARF in sweet cherry (Prunus avium L.). Plant Physiol. Biochem. 186, 299–309. doi: 10.1016/j.plaphy.2022.07.019

PubMed Abstract | Crossref Full Text | Google Scholar

Luo, Z., Jones, D., Philp-Wright, S., Putterill, J., Snowden, K. C. (2023). Transcriptomic analysis implicates ABA signaling and carbon supply in the differential outgrowth of petunia axillary buds. BMC Plant Biol. 23, 482. doi: 10.1186/s12870-023-04505-3

PubMed Abstract | Crossref Full Text | Google Scholar

Lv, L., Huo, X., Wen, L., Gao, Z., Khalil-Ur-Rehman, M. (2018). Isolation and role of PmRGL2 in GA-mediated floral bud dormancy release in Japanese Apricot (Prunus mume Siebold et Zucc.). Front. Plant Sci. 9. doi: 10.3389/fpls.2018.00027

PubMed Abstract | Crossref Full Text | Google Scholar

Marković, M., Trifunović Momčilov, M., Uzelac, B., Radulović, O., Milošević, S., Jevremović, S., et al. (2020). Breaking the dormancy of snake’s head fritillary (Fritillaria meleagris L.) in vitro bulbs-part 2: effect of GA3 soaking and chilling on sugar status in sprouted bulbs. Plants (Basel). 9, 1573. doi: 10.3390/plants9111573

PubMed Abstract | Crossref Full Text | Google Scholar

Nakatani, Y., Lamont, I. L., Cutfield, J. F. (2010). Discovery and characterization of a distinctive exo-1,3/1,4-{beta}-glucanase from the marine bacterium Pseudoalteromonas sp. strain BB1. Appl. Environ. Microbiol. 76, 6760–6768. doi: 10.1128/AEM.00758-10

PubMed Abstract | Crossref Full Text | Google Scholar

Niu, Q., Li, J., Cai, D., Qian, M., Jia, H., Bai, S., et al. (2016). Dormancy-associated MADS-box genes and microRNAs jointly control dormancy transition in pear (Pyrus pyrifolia white pear group) flower bud. J. Exp. Bot. 67, 239–257. doi: 10.1093/jxb/erv454

PubMed Abstract | Crossref Full Text | Google Scholar

Ohanenye, I. C., Alamar, M. C., Thompson, A. J., Terry, L. A. (2019). Fructans redistribution prior to sprouting in stored onion bulbs is a potential marker for dormancy break. Postharvest Biol. Technol. 149, 221–234. doi: 10.1016/j.postharvbio.2018.12.002

PubMed Abstract | Crossref Full Text | Google Scholar

Pan, W., Liang, J., Sui, J., Li, J., Liu, C., Xin, Y., et al. (2021). ABA and bud dormancy in perennials: current knowledge and future perspective. Genes (Basel). 12, 1635. doi: 10.3390/genes12101635

PubMed Abstract | Crossref Full Text | Google Scholar

Park, J. Y., Canam, T., Kang, K. Y., Unda, F., Mansfield, S. D. (2009). Sucrose phosphate synthase expression influences poplar phenology. Tree Physiol. 29, 937–946. doi: 10.1093/treephys/tpp028

PubMed Abstract | Crossref Full Text | Google Scholar

Qiu, Y., Guan, S. C., Wen, C., Li, P., Gao, Z., Chen, X. (2019). Auxin and cytokinin coordinate the dormancy and outgrowth of axillary bud in strawberry runner. BMC Plant Biol. 19, 528. doi: 10.1186/s12870-019-2151-x

PubMed Abstract | Crossref Full Text | Google Scholar

Qiu, Z., Wan, L., Chen, T., Wan, Y., He, X., Lu, S., et al. (2013). The regulation of cambial activity in Chinese fir (Cunninghamia lanceolata) involves extensive transcriptome remodeling. New Phytol. 199, 708–719. doi: 10.1111/nph.12301

PubMed Abstract | Crossref Full Text | Google Scholar

Rehman, M., Wang, W., Dong, Y., Faheem, M., Xu, Y., Gao, Z., et al. (2019). Comparative transcriptomic and proteomic analysis to deeply investigate the role of hydrogen cyanamide in grape bud dormancy. Int. J. Mol. Sci. 20, 3528. doi: 10.3390/ijms20143528

PubMed Abstract | Crossref Full Text | Google Scholar

Rothkegel, K., Sandoval, P., Soto, E., Ulloa, L., Riveros, A., Lillo-Carmona, V., et al. (2020). Dormant but Active: Chilling Accumulation Modulates the Epigenome and Transcriptome of Prunus avium During Bud Dormancy. Front. Plant Sci. 11. doi: 10.3389/fpls.2020.01115

PubMed Abstract | Crossref Full Text | Google Scholar

Sapkota, S., Salem, M., Jahed, K. R., Artlip, T. S., Sherif, S. M. (2023). From endodormancy to ecodormancy: the transcriptional landscape of apple floral buds. Front. Plant Sci. 14. doi: 10.3389/fpls.2023.1194244

PubMed Abstract | Crossref Full Text | Google Scholar

Shangguan, L., Chen, M., Fang, X., Xie, Z., Gong, P., Huang, Y., et al. (2020). Comparative transcriptome analysis provides insight into regulation pathways and temporal and spatial expression characteristics of grapevine (Vitis vinifera) dormant buds in different nodes. BMC Plant Biol. 20, 390. doi: 10.1186/s12870-020-02583-1

PubMed Abstract | Crossref Full Text | Google Scholar

Sheik, A., Kwanwoo, K., Ganji Lakshmi, V., Lee, H., Kim, S., Kim, E., et al. (2021). The anti-cancerous activity of adaptogenic herb Astragalus membranaceus. Phytomedicine. 91, 153698. doi: 10.1016/j.phymed.2021.153698

PubMed Abstract | Crossref Full Text | Google Scholar

Signorelli, S., Agudelo-Romero, P., Meitha, K., Foyer, C. H., Considine, M. J. (2018). Roles for light, energy, and oxygen in the fate of quiescent axillary buds. Plant Physiol. 176, 1171–1181. doi: 10.1104/pp.17.01479

PubMed Abstract | Crossref Full Text | Google Scholar

Sorce, C., Lorenzi, R., Ceccarelli, N., Ranalli, P. (2000). Changes in free and conjugated IAA during dormancy and sprouting of potato tubers. Aust. J. Plant Physiol. 27, 371–377. doi: 10.1016/j.jplph.2008.12.003

PubMed Abstract | Crossref Full Text | Google Scholar

Tang, L., Chhajed, S., Vashisth, T., Olmstead, M., Olmstead, J., Colquhoun, T. (2019). Transcriptomic study of early responses to the bud dormancy-breaking agent hydrogen cyanamide in 'TropicBeauty' peach. J. Am. Soc. Hortic. Sci. 144, 244–256. doi: 10.21273/JASHS04686-19

Crossref Full Text | Google Scholar

Tarancón, C., González-Grandío, E., Oliveros, J. C., Nicolas, M., Cubas, P. (2017). A conserved carbon starvation response underlies bud dormancy in woody and herbaceous species. Front. Plant Sci. 8. doi: 10.3389/fpls.2017.00788

PubMed Abstract | Crossref Full Text | Google Scholar

Thirugnanasambantham, K., Prabu, G., Mandal, A. K. A. (2020). Synergistic effect of cytokinin and gibberellins stimulates release of dormancy in tea (Camellia sinensis (L.) O. Kuntze) bud. Physiol. Mol. Biol. Plants. 26, 1035–1045. doi: 10.1007/s12298-020-00786-2

PubMed Abstract | Crossref Full Text | Google Scholar

Tosetti, R., Waters, A., Chope, G., Cools, K., Alamar, M., McWilliam, S., et al. (2021). New insights into the effects of ethylene on ABA catabolism, sweetening and dormancy in stored potato tubers. Postharvest Biol. Technol. 173, 111420. doi: 10.1016/j.postharvbio.2020.111420

PubMed Abstract | Crossref Full Text | Google Scholar

Tylewicz, S., Petterle, A., Marttila, S., Miskolczi, P., Azeez, A., Singh, R. K., et al. (2018). Photoperiodic control of seasonal growth is mediated by ABA acting on cell-cell communication. Science. 360, 212–215. doi: 10.1126/science.aan8576

PubMed Abstract | Crossref Full Text | Google Scholar

Voogd, C., Wang, T., Varkonyi-Gasic, E. (2015). Functional and expression analyses of kiwifruit SOC1-like genes suggest that they may not have a role in the transition to flowering but may affect the duration of dormancy. J. Exp. Bot. 66, 4699–4710. doi: 10.1093/jxb/erv234

PubMed Abstract | Crossref Full Text | Google Scholar

Wang, X., Li, D., Zhang, D., Shi, X., Wu, Y., Qi, Z., et al. (2020). Improving crucial details and selecting the optimal model for evaluating the chilling requirement of Paeonia lactiflora Pall. at low latitudes during four winters. Sci. Hortic. 265, 109175. doi: 10.1016/j.scienta.2020.109175

Crossref Full Text | Google Scholar

Wang, Y., Liu, T., Ma, C., Li, G., Wang, X., Wang, J., et al. (2022b). Carbohydrate regulation response to cold during rhizome bud dormancy release in Polygonatum kingianum. BMC Plant Biol. 22, 163. doi: 10.1186/s12870-022-03558-0

PubMed Abstract | Crossref Full Text | Google Scholar

Wang, Y., Liu, X., Su, H., Yin, S., Han, C., Hao, D., et al. (2019). The regulatory mechanism of chilling-induced dormancy transition from endo-dormancy to non-dormancy in Polygonatum kingianum Coll.et Hemsl rhizome bud. Plant Mol. Biol. 99, 205–217. doi: 10.1007/s11103-018-0812-z

PubMed Abstract | Crossref Full Text | Google Scholar

Wang, W., Su, X., Tian, Z., Liu, Y., Zhou, Y., He, M. (2018). Transcriptome profiling provides insights into dormancy release during hypothermia of Lilium pumilum. BMC Genomics 19, 196. doi: 10.1186/s12864-018-4536-x

PubMed Abstract | Crossref Full Text | Google Scholar

Wang, H., Xia, X., An, L. (2021). Metabolomics analysis reveals the mechanism of hydrogen cyanamide in promoting flower bud break in blueberry. Agronomy. 11, 102. doi: 10.3390/agronomy11010102

Crossref Full Text | Google Scholar

Wang, X., Zhang, R., Huang, Q., Shi, X., Li, D., Shao, L., et al. (2022a). Comparative study on physiological responses and gene expression of bud endodormancy release between two herbaceous peony cultivars (Paeonia lactiflora Pall.) with contrasting chilling requirements. Front. Plant Sci. 12. doi: 10.3389/fpls.2021.772285

PubMed Abstract | Crossref Full Text | Google Scholar

Wang, Q., Zhang, J. J., Teixeira da Silva, J. A., Yu, X. N. (2014). Temporal physiological and biochemical changes in Hippeastrum vittatum ‘Redlion’ bulbs stored at different temperatures. Maejo Int. J. Sci. 8, 114–121. doi: 10.14456/mijst.2014.10

Crossref Full Text | Google Scholar

Wojtania, A., Markiewicz, M., Waligórski, P. (2022). Regulation of the bud dormancy development and release in micropropagated rhubarb 'Malinowy'. Int. J. Mol. Sci. 23, 1480. doi: 10.3390/ijms23031480

PubMed Abstract | Crossref Full Text | Google Scholar

Wu, R., Tomes, S., Karunairetnam, S., Tustin, S. D., Hellens, R. P., Allan, A. C., et al. (2017). SVP-like MADS box genes control dormancy and budbreak in apple. Front. Plant Sci. 8. doi: 10.3389/fpls.2017.00477

PubMed Abstract | Crossref Full Text | Google Scholar

Xie, Z., Zhang, Z., Hanzlik, S., Cook, E., Shen, Q. (2007). Salicylic acid inhibits gibberellin-induced alpha-amylase expression and seed germination via a pathway involving an abscisic-acid-inducible WRKY gene. Plant Mol. Biol. 64, 293–303. doi: 10.1007/s11103-007-9152-0

PubMed Abstract | Crossref Full Text | Google Scholar

Xu, T., Zhang, J., Shao, L., Wang, X., Zhang, R., Ji, C., et al. (2022). Later growth cessation and increased freezing tolerance potentially result in later dormancy in evergreen iris compared with deciduous iris. Int. J. Mol. Sci. 23, 11123. doi: 10.3390/ijms231911123

PubMed Abstract | Crossref Full Text | Google Scholar

Yang, Q., Gao, Y., Wu, X., Moriguchi, T., Bai, S., Teng, Y. (2021). Bud endodormancy in deciduous fruit trees: Advances and prospects. Hortic. Res. 8, 139. doi: 10.1038/s41438-021-00575-2

PubMed Abstract | Crossref Full Text | Google Scholar

Yang, Y., Xu, C., Wang, B., Jia, J. (2001). Effects of plant growth regulators on secondary wall thickening of cotton fibres. Plant Growth Regul. 35, 233–237. doi: 10.1023/A:1014442015872

Crossref Full Text | Google Scholar

Yordanov, Y. S., Ma, C., Strauss, S. H., Busov, V. B. (2014). EARLY BUDBREAK 1 (EBB1) is a regulator of release from seasonal dormancy in poplar trees. Proc. Natl. Acad. Sci. U S A. 111, 10001–10006. doi: 10.1073/pnas.1405621111

PubMed Abstract | Crossref Full Text | Google Scholar

Yu, Z., Zhang, D., Hu, S., Liu, X., Zeng, B., Gao, W., et al. (2022). Genome-wide analysis of the almond AP2/ERF superfamily and its functional prediction during dormancy in response to freezing stress. Biol. (Basel). 11, 1520. doi: 10.3390/biology11101520

PubMed Abstract | Crossref Full Text | Google Scholar

Yuan, Y., Zeng, L., Kong, D., Mao, Y., Xu, Y., Wang, M., et al. (2024). Abscisic acid-induced transcription factor PsMYB306 negatively regulates tree peony bud dormancy release. Plant Physiol. 194, 2449–2471. doi: 10.1093/plphys/kiae014

PubMed Abstract | Crossref Full Text | Google Scholar

Yue, C., Cao, H., Hao, X., Zeng, J., Qian, W., Guo, Y., et al. (2018). Differential expression of gibberellin- and abscisic acid-related genes implies their roles in the bud activity-dormancy transition of tea plants. Plant Cell Rep. 37, 425–441. doi: 10.1007/s00299-017-2238-5

PubMed Abstract | Crossref Full Text | Google Scholar

Zeng, R., Chen, T., Wang, X., Cao, J., Li, X., Xu, X., et al. (2021). Physiological and expressional regulation on photosynthesis, starch and sucrose metabolism response to waterlogging stress in peanut. Front. Plant Sci. 12. doi: 10.3389/fpls.2021.601771

PubMed Abstract | Crossref Full Text | Google Scholar

Zhang, J., Li, D., Shi, X., Zhang, D., Qiu, S., Wei, J., et al. (2017). Mining and expression analysis of candidate genes involved in regulating the chilling requirement fulfillment of Paeonia lactiflora 'Hang Baishao'. BMC Plant Biol. 17, 262. doi: 10.1186/s12870-017-1205-1

PubMed Abstract | Crossref Full Text | Google Scholar

Zhang, J., Liu, H., Yang, B., Wang, H., Niu, S., El-Kassaby, Y. A., et al. (2022). Phytohormone profiles and related gene expressions after endodormancy release in developing Pinus tabuliformis male strobili. Plant Sci. 316, 111167. doi: 10.1016/j.plantsci.2021.111167

PubMed Abstract | Crossref Full Text | Google Scholar

Zhang, R., Wang, X., Shi, X., Shao, L., Xu, T., Xia, Y., et al. (2021a). Chilling requirement validation and physiological and molecular responses of the bud endodormancy release in Paeonia lactiflora 'Meiju'. Int. J. Mol. Sci. 22, 8382. doi: 10.3390/ijms22168382

PubMed Abstract | Crossref Full Text | Google Scholar

Zhang, Y., Yuan, Y., Liu, Z., Zhang, T., Li, F., Liu, C., et al. (2021b). GA3 is superior to GA4 in promoting bud endodormancy release in tree peony (Paeonia suffruticosa) and their potential working mechanism. BMC Plant Biol. 21, 323. doi: 10.1186/s12870-021-03106-2

PubMed Abstract | Crossref Full Text | Google Scholar

Zhang, T., Yuan, Y., Zhan, Y., Cao, X., Liu, C., Zhang, Y., et al. (2020). Metabolomics analysis reveals Embden Meyerhof Parnas pathway activation and flavonoids accumulation during dormancy transition in tree peony. BMC Plant Biol. 20, 484. doi: 10.1186/s12870-020-02692-x

PubMed Abstract | Crossref Full Text | Google Scholar

Zhang, Z., Zhuo, X., Zhao, K., Zheng, T., Han, Y., Yuan, C., et al. (2018). Transcriptome profiles reveal the crucial roles of hormone and sugar in the bud dormancy of Prunus mume. Sci. Rep. 8, 5090. doi: 10.1038/s41598-018-23108-9

PubMed Abstract | Crossref Full Text | Google Scholar

Zheng, Y., Ren, W., Zhang, L., Zhang, Y., Liu, D., Liu, Y. (2020). A review of the pharmacological action of Astragalus polysaccharide. Front. Pharmacol. 11. doi: 10.3389/fphar.2020.00349

PubMed Abstract | Crossref Full Text | Google Scholar

Zhou, H., He, Y., Zhu, Y., Li, M., Song, S., Bo, W., et al. (2020). Comparative transcriptome profiling reveals cold stress responsiveness in two contrasting Chinese jujube cultivars. BMC Plant Biol. 20, 240. doi: 10.1186/s12870-020-02450-z

PubMed Abstract | Crossref Full Text | Google Scholar

Zhu, H., Chen, P. Y., Zhong, S., Dardick, C., Callahan, A., An, Y. Q., et al. (2020). Thermal-responsive genetic and epigenetic regulation of DAM cluster controlling dormancy and chilling requirement in peach floral buds. Hortic. Res. 7, 114. doi: 10.1038/s41438-020-0336-y

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: Astragalus membranaceus (Fisch.) Bge. var. mongholicus (Bge.) Hsiao, bud dormancy transition, metabolome, transcriptome, co-expression network

Citation: Guan H, Zhao Y, Chen Q, Zhang Q, Yang P, Sun S and Chen G (2025) Integrated metabolome and transcriptome analysis reveals potential mechanism during the bud dormancy transition of Astragalus membranaceus (Fisch.) Bge. var. mongholicus (Bge.) Hsiao. Front. Plant Sci. 15:1483538. doi: 10.3389/fpls.2024.1483538

Received: 20 August 2024; Accepted: 18 December 2024;
Published: 21 January 2025.

Edited by:

Inês Maria Valente, LAQV Network of Chemistry and Technology, Portugal

Reviewed by:

Klára Kosová, Crop Research Institute (CRI), Czechia
Yanfei Cai, Flower Research Institute of Yunnan Academy of Agricultural Sciences, China

Copyright © 2025 Guan, Zhao, Chen, Zhang, Yang, Sun and Chen. 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: Shuying Sun, bm1nc3Vuc2h1eWluZ0AxNjMuY29t; Guilin Chen, Z3VpbGluY2hlbjYxQDE2My5jb20=

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.