Salt stress enhances bioactive compound accumulation in Glycyrrhiza inflata: integrated transcriptomics and physiological analysis reveals germplasm-specific adaptation mechanisms

Abstract

Introduction:

Glycyrrhiza inflata Batal., a halophytic plant predominantly found in the saline-alkali deserts of southern Xinjiang, China, is renowned for its abundance of bioactive compounds like flavonoids and triterpenoids. It demonstrates considerable potential for applications within the pharmaceutical, food, health product, and cosmetic industries. Additionally, its cultivation presents the dual advantage of generating economic returns and facilitating the remediation of saline-alkali soils.

Methods:

This study examined 29 distinct provenances of G. inflata collected from various locations across Xinjiang. Key agronomic traits and the content of bioactive compounds in the underground parts of one-year-old plants grown in severely saline-alkali soil were measured to assess inter-germplasm variation. Subsequently, four germplasms displaying contrasting quality and salt tolerance were selected for controlled salt stress treatment (150 mM NaCl) under laboratory conditions. The effects on seed germination, root bioactive compound content, endogenous hormone levels, and key physiological and biochemical indices were analyzed. An integrated analysis of salt stress transcriptomic data was conducted using Weighted Gene Co-expression Network Analysis (WGCNA). This involved expression clustering and enrichment analysis of differentially expressed genes (DEGs) to investigate the impact of salt stress on genes related to bioactive compound biosynthesis (particularly flavonoids), endogenous hormone pathways, and key flavonoid biosynthesis enzymes.

Results:

The findings indicate that germplasms with superior stress tolerance maintained higher and more stable levels of antioxidant enzymes. In response to stress, these resilient germplasms modulated hormone signaling, notably upregulating abscisic acid (ABA) and downregulating auxin (IAA), thereby reallocating resources towards defense mechanisms. Crucially, salt stress was identified as an effective means to enhance the accumulation of bioactive compounds in G. inflata. Transcriptomic analysis revealed substantial divergence in post-stress gene expression patterns among germplasms, implicating key pathways such as plant hormone signal transduction, flavonoid biosynthesis, and phenylpropanoid metabolism.

Discussion:

This research establishes a foundation for breeding high-quality G. inflata germplasms adapted to desert saline-alkali environments and provides valuable insights into the molecular mechanisms regulating the synthesis and accumulation of its valuable bioactive compounds.

1 Introduction

Soil salinization is a major contributor to soil degradation, reduced crop yields, and ecosystem imbalance. This process has significantly compromised the sustainability of agricultural production systems and jeopardizes global food security (Hassani et al., 2021) Globally, approximately one-third of arable land is located in arid and semi-arid regions. These areas are frequently affected by salinization. Soil salinization adversely impacts crops through multiple mechanisms, high salt concentrations impede seed water imbibition, suppressing germination rates; It also restricts root development, impairing water and nutrient uptake. This subsequently inhibits shoot growth, leading to stunted plants, chlorosis, wilting, and even mortality. Furthermore, salinization alters soil microbial community structure and compromises ecosystem functions, indirectly affecting crop growth and productivity (van Zelm et al., 2020)Confronted with this urgent challenge, researchers worldwide are dedicated to developing innovative remediation technologies for saline-alkali soils. Several approaches with demonstrated potential and theoretical foundations have been established. Hydraulic engineering, a conventional method (Du et al., 2023), regulates soil water-salt dynamics through optimized irrigation and drainage. Agricultural management focuses on agronomic interventions, such as deep plowing, application of soil amendments (Dai et al., 2025) and strategic crop rotation.

Biological remediation has attracted considerable research focus in recent years, with halophyte-based approaches demonstrating distinct advantages and significant potential (Duan et al., 2022), Halophytes, defined as plants capable of thriving naturally in high-salinity environments, adapt to these conditions through specialized physiological and biochemical mechanisms (Jiao et al., 2022), These include intracellular ion homeostasis and osmoregulation (Sun et al., 2021) Beyond extracting soil salts to reduce salinity levels, halophytes enhance soil physicochemical properties via root exudates. This process promotes soil aggregate formation, increases porosity, and thereby creates favorable microenvironments for establishing other plant species (Li et al., 2022)。Currently, halophytes demonstrate substantial remediation efficacy in saline-alkali soil restoration. projects (Rahman et al., 2021), positioning them as promising candidates in bioremediation. For example, Suaeda salsa ameliorates soil conditions by modulating microbial communities and upregulating carbon fixation and sulfite oxidation genes (Tang et al., 2023), Medicago sativa cultivation elevates soil nitrogen and organic matter content, modifies nutrient profiles, and restructures fungal/bacterial community composition (Qi et al., 2023).

Among halophytes, Glycyrrhiza inflata Batal. emerges as a high-value psammohalophyte predominantly inhabiting saline-alkaline deserts of southern Xinjiang, China. This species exhibits remarkable environmental resilience under extreme conditions. Beyond its exceptional salt tolerance enabling survival in highly saline-alkali soils, G. inflata accumulates medicinally significant triterpenoids and flavonoids—including glycyrrhizin, liquiritin, and licochalcone A (Jiang et al., 2020), demonstrating potent antioxidant (Liu et al., 2024), anti-inflammatory (Wu et al., 2024),and anti-cancer activities (Meng et al., 2024). Within the global pharmaceutical market, G. inflata-derived products hold significant market share across Asia, Europe, and North America (Ma et al., 2025a),attributed to their unique therapeutic properties. Consequently, research and utilization of this species offers a viable bioremediation strategy for decertified saline-alkali lands while simultaneously stimulating regional economic growth through industry development and job creation.

Recent trends reveal a critical conservation status of wild Glycyrrhiza inflata resources (Yan et al., 2023),with dwindling populations failing to meet current industrial and medicinal demands for both quantity and quality. While commercially cultivated alternatives provide partial substitution, their inconsistent quality remains a significant limitation (Wang et al., 2019),Consequently, elucidating the salt and drought tolerance mechanisms of G. inflata transcends academic interest—it represents an urgent research imperative. This work is fundamental to ensuring sustainable resource supply and advancing evidence-based cultivation practices for the licorice industry.

Given that glycyrrhizin and liquiritin—classified as triterpenoid and flavonoid markers respectively (Yang et al., 2017), serve as critical quality determinants for licorice, enhancing their biosynthesis presents an urgent research priority. Recent advances in plant stress physiology demonstrate that moderate salt stress effectively stimulates bioactive compound accumulation across diverse species (Wang et al., 2020; Gu et al., 2024; Li et al., 2024), This evidence-based approach offers a novel pathway for optimizing G. inflata phytochemical yield.

We initiated germplasm screening by conducting germination assays and quantifying root glycyrrhizin/liquiritin content in 29 G. inflata accessions from severely saline-alkali soils, selecting elite germplasms exhibiting superior salt tolerance, high germination rates, and enhanced phytochemical profiles. Two salt-tolerant and two salt-sensitive accessions with contrasting metabolite accumulation were then subjected to controlled salt stress (150 mM NaCl). Treatment effects on germination dynamics, root bioactive compounds, endogenous hormones, and physiological indices were analyzed. Subsequent transcriptomic investigation integrated differential gene expression profiling, co-expression clustering, and functional enrichment to elucidate salt stress impacts on bioactive compound biosynthesis, hormone signaling, and flavonoid pathway genes. Finally, Weighted Gene Co-expression Network Analysis (WGCNA) (Liu et al., 2019) identified hub genes co-regulated with stress tolerance and flavonoid accumulation, revealing preliminary mechanistic relationships between specialized metabolism and stress adaptation.

2 Materials and methods

2.1 Plant materials and treatments

All 29 Glycyrrhiza inflata accessions used were pure Xinjiang germplasms previously collected and authenticated by our laboratory. Seeds underwent 50-min 98% sulfuric acid treatment to break dormancy and enhance germination, followed by thorough rinsing to remove residues. Field cultivation occurred at Shihezi University Experimental Station (44°16’N, 86°00’E) under temperate continental climate (≈170-day frost-free period; July mean: 25–26°C; annual precipitation: 125–207 mm). Seeds were planted in May using drip irrigation (2.05-m wide mulch, six rows per bed, three drip lines) with alternating wide-narrow row spacing (18,000 plants/ha; irrigation: 4,200 m³/ha). Roots (20 cm segments from 5 plants/accession) were harvested during late September–early October for yield assessment (fresh/dry weight) and phytochemical analysis.

Saline-alkali field trials were conducted in Minfeng County (82°22’E, 35°20’N) featuring typical temperate desert climate (mean annual temp: stable; precipitation: 30.5 mm; evaporation: 2,756 mm; frost-free period: 194 days; sunshine: 2,842.2 hr). Identical root sampling protocols were followed.

Germination assays used plump, uniform seeds (50 seeds/dish on filter paper; 3 treatments × 3 replicates) with 10 mL test solution replenished periodically using original concentrations to maintain osmolarity (distilled water control). Salt stress experiments proceeded in climate-controlled chambers (200 μmol·m⁻²·s⁻¹ PPFD; 16/8-h light/dark; 50–55% RH; 28°C/25°C day/night). Seeds were sown in vermiculite-filled pots (10×10×10 cm; 6 plants/pot), watered biweekly and fertilized weekly with half-strength MS medium. At 20-cm height, plants received 150 mM NaCl (pre-optimized concentration). Aboveground and root tissues (harvested at root-shoot junction) were sampled at 0, 2, 6, 12, and 24 h post-treatment (three biological replicates). Samples were rinsed, flash-frozen in liquid nitrogen, and stored at −80 °C for subsequent phytochemical quantification and RNA extraction (transcriptome sequencing).

2.2 Determination of soil salinity and major elements

The measurement of this part was carried out by a third-party institution, the Agricultural Science Research Institute of the 7th Division of Xinjiang Production and Construction Corps. They took five samples from the 4,500 mu of land in Area 1 of the Salayuketek Township, Minfen County, Xinjiang Province, at depths of 0–30 cm and 30–60 cm. The pH measurement results of the 0–30 cm soil samples were 8.27, with water-soluble salts at 5.33 g/kg, available nitrogen at 1.62 mg/kg, available phosphorus at mg/kg, and available potassium at 106 mg/kg. The pH measurement results of the 30 –60 cm soil samples were 8.15, with water-soluble salts at 6.00 g/kg, available nitrogen at 0.82 mg/kg, available phosphorus at 0.09 mg/kg, and available potassium at 66 mg/kg.

2.3 Extraction and determination of active pharmaceutical components

Powdered samples (1.000 g) were ultrasonically extracted with 5 mL 100% methanol at 25°C for 60 min (1000 W), followed by centrifugation at 12,000 rpm (4°C, 5 min). Supernatants were collected, and residues were re-extracted with 5 mL methanol. Combined supernatants were diluted to volume in a volumetric flask with methanol. Aliquots (2.5 mL) were syringe-filtered through 0.22 µm membranes for phytochemical analysis. Standard solutions (1.0 mg/mL) were serially diluted in methanol to construct calibration curves (1 – 100 ng/mL). Quantitative analysis was performed using an ACQUITY UPLC H-Class system (Waters) coupled to a Xevo TQS triple quadrupole mass spectrometer. Separation was achieved on an ACQUITY UPLC BEH C18 column (50 × 2.1 mm, 1.7 μm; Waters) maintained at 40°C with a 0.3 mL/min flow rate. The mobile phase consisted of (A) 0.1% formic acid in water and (B) acetonitrile, using the following gradient: 0–2 min 5% B, 2–8 min 5-95% B, 8–10 min 95% B, 10-10.1 min 95-5% B, 10.1–12 min 5% B.

2.4 Extraction and determination of total flavonoids from licorice

For total flavonoid quantification, 2.00 g of licorice root powder was accurately weighed, wrapped in filter paper, and subjected to Soxhlet extraction with ethanol under 80°C water bath reflux until colorless effluent was observed. The extract was concentrated under reduced pressure, transferred to a 100-mL volumetric flask, and diluted to volume with ethanol. A 0.5-mL aliquot of this solution was mixed with 2.5 mL distilled water in a 10-mL volumetric flask. Sequentially, 0.5 mL of 5% NaNO₂ and 0.5 mL of 10% Al(NO₃)₃ were added with vortex mixing after each addition. After 5-min incubation at room temperature, 2.5 mL of 5% NaOH was added and mixed thoroughly. Following additional 5-min incubation, the solution was diluted to volume with distilled water. Absorbance was measured at 520 nm against a reagent blank using UV-Vis spectrophotometer.

2.5 Antioxidant enzyme content determination

Fresh leaf samples (0.5 g) were homogenized in 4.5 mL ice-cold 50 mM phosphate buffer (pH 7.0 - 7.4) containing 0.1 mM EDTA and 2 mM DTT. The homogenate was centrifuged at 4,000 ×g for 50 min at 4°C, and the supernatant was retained for antioxidant enzyme assays. Content of superoxide dismutase (SOD), peroxidase (POD), catalase (CAT), free proline, and malondialdehyde (MDA) were quantified in G. inflata roots using commercial assay kits (Nanjing Jiancheng Bioengineering Institute, China) according to manufacturer protocols.

2.6 Chlorophyll content measurement

Fresh tissue samples (0.1 g) were finely chopped and transferred to 50-mL conical flasks containing 10 mL of 95% (v/v) ethanol. After 3 – 5 h of dark incubation at room temperature until complete decolorization occurred, extracts were filtered into 25-mL amber volumetric flasks. The mortar and pestle were rinsed three times with small volumes of 95% ethanol, with washings filtered into the same flask. The combined filtrate was diluted to volume with ethanol and homogenized. Extract aliquots were transferred to cuvettes against a 95% ethanol blank. Absorbance was measured at 665 nm and 649 nm using UV-Vis spectrophotometer.

2.7 Transcriptome analysis

Total RNA was isolated using the Polyphenolics-Polysaccharides RNA Kit (TIANGEN Biotech, China). cDNA library construction was performed by Majorbio Bio-Pharm Technology Co., Ltd. Poly(A)+ mRNA was enriched from qualified RNA samples using oligo(dT)-attached magnetic beads. Following fragmentation (~300 bp fragments compatible with Illumina HiSeq 6000 platform requirements), first-strand cDNA synthesis was conducted with reverse transcriptase. Second-strand cDNA synthesis generated stable double-stranded DNA, which underwent end repair, adapter ligation, and PCR amplification for library enrichment. Thirty-six cDNA libraries were constructed, representing four G. inflata accessions under salt stress across two timepoints with triplicate biological replicates. Libraries were sequenced on the Illumina HiSeq 6000 platform (150-bp paired-end). Raw reads underwent quality control using SeqPrep and Sickle (Chen et al., 2025) to obtain clean data. High-quality reads were aligned to the G. uralensis reference genome (v2.0) using TopHat2 (Kim et al., 2013)and HISAT2 (Kim et al., 2019).Gene expression quantification was performed with RSEM (Li and Dewey, 2011). Differential expression analysis employed DEGseq (Wang et al., 2010)and DESeq2 (Love et al., 2014) (|log₂FC| > 1, FDR < 0.05). Functional annotation utilized DIAMOND for GO terms and KEGG Orthology (KO) assignment. Enrichment analyses were conducted with GOATools (Klopfenstein et al., 2018)and KOBAS (Bu et al., 2021) (FDR < 0.05). Weighted Gene Co-expression Network Analysis (WGCNA) (Liu et al., 2019)was implemented in R v4.2.2 to identify salt-responsive genes associated with specialized metabolism. An adjacency matrix was constructed using soft thresholding (power = β = 14). Modules were detected via dynamic tree cutting (mergeCutHeight = 0.25, minModuleSize = 30). Module eigengenes (first principal components) represented expression profiles. Spearman correlations between module eigengenes and flavonoid abundances were visualized using ggplot2.

2.8 Real-time fluorescence quantitative PCR

To validate RNA-seq accuracy, six differentially expressed genes (DEGs) regulating flavonoid biosynthesis were selected for qRT-PCR verification(Supplementary Figure S3). Specific primers were designed to measure relative expression levels of these genes in root tissues of four G. inflata accessions after 6-h salt treatment. The Glycyrrhiza uralensis lectin gene (Glyur000100s00008376) served as the endogenous reference. Reactions utilized the FastStart Universal SYBR Green Master Mix (Roche) on an Applied Biosystems 7500 Real-Time PCR System (Monad Biotech, China). Relative transcript abundance was calculated via the 2^–ΔΔCt method (Livak and Schmittgen, 2001), Statistical significance was determined by one-way ANOVA with Tukey’s post-hoc test (p < 0.05).

2.9 Statistical analysis

All statistical analyses were performed using R packages and SPSS 27.0. Intergroup differences were assessed by one-way analysis of variance (ANOVA), with post-hoc pairwise comparisons conducted via the least significant difference (LSD) test at a significance threshold of P < 0.05.

3 Result

3.1 Quantification of bioactive compounds and saline-alkali germination rates in one-year-old cultivated G. inflata

Confronting the dual challenges of diminishing wild Glycyrrhiza inflata resources and declining cultivated quality, we evaluated 29 germplasms from the Tarim Basin periphery (coordinates in Supplementary Table S1), revealing significant inter-accession variation in bioactive compounds (Table 1) and salt tolerance, as shown in Figure 1,all accessions achieved >45% germination in saline-alkali soil (pH 8.15; mean 65.13% ± 8.24%), with 20 exceeding 60% (nine >70%, two >80%, #10 peaking at 87.1%), while plant heights ranged 1.74 – 3.24 m (mean 2.23 m). Critical metabolites varied substantially: total flavonoids (17.93 ± 2.81 to 28.97 ± 5.53 mg/g), glycyrrhizic acid (2.06 ± 0.23 to 7.51 ± 1.48 mg/g), liquiritin (0.72 ± 0.14 to 4.01 ± 1.23 mg/g), liquiritigenin (0.05 ± 0.01 to 0.22 ± 0.03 mg/g), and isoliquiritigenin (0.09 ± 0.01 to 0.29 ± 0.01 mg/g), with hyper-arid desert margin accessions (#1,9,15) showing superior accumulation (P<0.05), suggesting strategic salt stress enhances medicinal biosynthesis. Through triaxial selection ranking germplasms by glycyrrhizic acid content, liquiritin content, and germination rate (Supplementary Figure S1)—driven by dual-purpose application in remediating southern Xinjiang’s saline-alkali deserts while generating local economic returns via harvestable bioactive compounds—we prioritized accessions meeting critical germination thresholds due to severely compromised viability in extreme salinity-alkalinity soils, thus selecting: stress-tolerant accessions (#10=T1, #28=T2) with high germination and enriched bioactive compounds, and sensitive accessions (#22=S1, #13=S2) exhibiting moderate germination rates (relatively lower than T1/T2) yet divergent bioactive compound accumulation for comparative stress-response analysis. their controlled-environment metabolite profiles (Figure 2) indicated no significant differences in general flavonoids, or isoglycyrrhizin (P>0.05), but glabridin, licochalcone A,liquiritigenin, liquiritin, and glycyrrhizic acid, The contents of these five active ingredients have varied.

Figure 1

Figure 2

3.2 Physiological and biochemical responses of four G. inflata germplasms to salt stress

To elucidate salt adaptation mechanisms and select elite germplasms with enhanced stress resilience and bioactive compound production, we first assessed germination rates under varying NaCl concentrations (Supplementary Figure S2). After 15-day exposure to 2% NaCl, significant inter-germplasm divergence emerged: T1 (62%) > T2 (32%) > S2 (8%) > S1 (2%). When plants reached a height of 20 cm, salt stress was applied using a concentration of 150 mM NaCl, previously determined as optimal in preliminary laboratory screening. Root tissues (sampled below the root-shoot junction) were collected at 0, 2, and 6 hours after stress initiation. All samples were immediately washed with distilled water, surface-dried, flash-frozen in liquid nitrogen, and stored at -80°C, with three biological replicates per treatment. As shown in Figure 3, the control groups of the four Glycyrrhiza inflata varieties exhibited no significant differences in malondialdehyde (MDA), proline, or superoxide dismutase (SOD) levels. However, following salt stress, the tolerant varieties (T1, T2) accumulated significantly higher proline levels (Figure 3C) than the sensitive varieties (S1, S2). Conversely, MDA content (Figure 3E) was significantly lower in T1 and T2 compared to S1 and S2. Catalase (CAT) activity (Figure 3A) was consistently higher in T1 and T2 both before and after stress exposure. These differences likely contribute to the enhanced stress tolerance of T1 and T2.While chlorophyll content (Figure 3B) differed between tolerant and sensitive varieties, the absolute difference was relatively small compared to the other parameters (CAT, proline, SOD, MDA) and did not reach statistical significance. Interestingly, SOD activity (Figure 3D) in the sensitive varieties (S1, S2) was significantly higher than in the tolerant varieties (T1, T2) after stress and increased progressively with prolonged stress duration. This result contradicts our initial expectations, and the underlying mechanism remains unclear, representing a key question for future investigation.

Figure 3

3.3 Phytohormonal dynamics in four G. inflata germplasms under salt stress

Overall, the four Glycyrrhiza inflata varieties exhibited distinct baseline levels of the measured compounds even under non-stressed conditions, and employed divergent strategies in response to salt stress. As a key hormone in plant stress responses, abscisic acid (ABA) levels were highest in the tolerant variety T1 under control conditions and remained relatively stable after stress onset. In contrast, the other tolerant variety (T2) significantly upregulated ABA synthesis following stress exposure. Interestingly, the sensitive variety S2 also adopted this strategy of post-stress ABA induction. Conversely, ABA content in the sensitive variety S1 steadily declined after stress imposition (Figure 4A), an unexpected pattern that may be a key factor contributing to its stress sensitivity. Regarding indole-3-acetic acid (IAA), post-stress dynamics differed markedly between tolerant and sensitive varieties (Figure 4B). The tolerant varieties T1 and T2 exhibited a progressive downregulation of IAA content with prolonged stress duration. Conversely, the sensitive varieties S1 and S2 displayed an increase in IAA levels over the same period. In contrast to both ABA and IAA, the four varieties displayed a consistent pattern for zeatin (ZT) synthesis (Figure 4C). All varieties initially increased ZT levels shortly after stress initiation, followed by a subsequent decline as stress duration extended.

Figure 4

3.4 Salt stress modulates bioactive compound accumulation in four G. inflata germplasms

Overall, short-term salt stress significantly increased the content of major active compounds in all four Glycyrrhiza inflata varieties, although accumulation patterns differed substantially among varieties for the same compound (Figure 5).As shown in Figures 5A, D, the sensitive varieties (S1, S2) accumulated significantly higher levels of glabridin and licochalcone A at 0, 2, and 6 hours post-stress compared to the tolerant varieties (T1, T2). Conversely, the tolerant varieties (T1, T2) exhibited greater accumulation of glycyrrhizic acid and liquiritin following stress exposure Figures 5B, E). Notably, glycyrrhizic acid and liquiritin are two primary compounds used in quality assessment of G. inflata. For isoliquiritigenin Figure 5C), accumulation patterns were similar between T1 and S2, and between T2 and S1.Analysis of the biosynthetic pathways (Figure 6) indicates that, with the exception of glycyrrhizic acid, the other four active compounds (glabridin, licochalcone A, liquiritin, isoliquiritigenin) are synthesized through multi-step reactions originating from coumaroyl-CoA. Notably, isoliquiritigenin and liquiritin share a common pathway, while glabridin and licochalcone A utilize the same precursor. This shared biosynthetic origin likely explains the observed similarities in their accumulation trends. The distinct accumulation patterns of these active compounds between sensitive and tolerant varieties may contribute to their differential stress tolerance.

Figure 5

Figure 6

3.5 Differential gene expression profiling

To validate our hypothesis and identify key salt-responsive genes in flavonoid biosynthesis, we performed transcriptome sequencing on root tissues of two salt-tolerant (T1, T2) and two salt-sensitive (S1, S2) G. inflata germplasms treated with 150 mM NaCl for 0, 2, and 6 hours. Differential expression analysis compared post-stress (2h/6h) versus control (0h) samples across accessions (T1_2h/T1_0h, T1_6h/T1_0h, T2_2h/T2_0h, T2_6h/T2_0h, S1_2h/S1_0h, S1_6h/S1_0h, S2_2h/S2_0h, S2_6h/S2_0h). Bar plots and Venn diagrams visualized DEG counts and overlaps (Figures 5A, B). At 2h post-stress, all accessions showed upregulated DEGs, but by 6h, tolerant (T1/T2) and moderately sensitive (S2) lines exhibited ∼50% reduction in upregulated DEGs (Figure 5C). Strikingly, S1 maintained elevated DEG counts, suggesting impaired transcriptional regulation under prolonged stress. Although total DEGs ranged 3,161-5,681 per accession, only 26 genes (0.19% of total DEGs) were commonly deregulated across all four germplasms (Figure 5D), indicating germplasm-specific transcriptional reprogramming.

3.6 Gene ontology analysis of DEGs

GO enrichment analysis of 2h/0h comparisons revealed significant enrichment (FDR < 0.05) in:Oxidative stress response, Flavonoid biosynthesis (T1:16 genes; T2:19; S1:14; S2:16),Flavonoid metabolism, Naringenin 2-hydroxylase activity, Flavonoid 3’-monooxygenase activity (exclusive to T1/S1 with 4–5 genes),Oxidoreductase activity, Hormone-mediated signaling (T1:122; T2:131; S1:91; S2:108 genes).At 6h/0h, pathways shifted toward: Isoflavonoid metabolism (S1-exclusive: 3 genes), Isoflavonoid biosynthesis, Flavonoid biosynthesis (T1:10; T2:15; S1:11; S2:10 genes),Glucosidase activity, Flavonoid 3’-monooxygenase activity (S1-exclusive: 5 genes),Hormone signaling (T1:99; S1:108; S2:81 genes; absent in T2).This temporal progression indicates transcriptional activation peaked at 2h post-stress, with flavonoid pathway genes and Hormone-mediated signaling showing maximal induction (Figure 7E). By 6h, attenuated activation coincided with specialized isoflavonoid pathway engagement, particularly in S1.

Figure 7

3.7 KEGG pathway enrichment of DEGs

To functionally characterize differentially expressed genes (DEGs), we performed KEGG enrichment analysis on 2h/0h and 6h/0h comparisons. The top 20 significantly enriched pathways (FDR < 0.01) shared across treatment groups revealed. At 2h post-stress:Isoflavonoid biosynthesis (ko00943), Phenylalanine metabolism (ko00360),Plant MAPK signaling (ko04016), Nitrogen metabolism (ko00910), Ubiquinone/terpenoid-quinone biosynthesis (ko00130), Phenylpropanoid biosynthesis (ko00940). Notably enriched flavonoid pathways: Flavone and flavonol biosynthesis (ko00944), Flavonoid biosynthesis (ko00941),Isoflavonoid biosynthesis (ko00943), Phenylpropanoid biosynthesis (ko00940). Concurrent enrichment of Plant hormone signal transduction (ko04075) suggests coordinated regulation. The significant enrichment (P <0.04) of multiple flavonoid biosynthesis pathways correlates with observed phytochemical upregulation (Section 3.4), indicating transcriptional reprogramming drives stress-induced metabolite accumulation.

3.8 Weighted gene co-expression network analysis

Using correlation-based clustering, we grouped co-expressed genes into modules via hierarchical clustering dendrograms, where branch colors represent distinct expression patterns (Figure 8A). WGCNA partitioned G. inflata genes into 68 co-expression modules, with the turquoise module being largest (1,235 genes) and plum smallest (31 genes). Module-trait associations were visualized through:Gene-trait correlation heatmap (Figure 8B),Mean connectivity distribution (Figure 8C), Scale-free topology fit index (Figure 8D). Spearman correlation analysis revealed significant compound-module associations: The greenyellow and black modules have the strongest correlation with Glabridin, the darkgreen module has the strongest correlation with Licochalcone A, the cyan module has the strongest correlation with glycyrrhizic acid module, and the plum1 module has the strongest correlation with Liquiritin. To further understand the relationship between the relevant modules and active components, and to identify key genes, we conducted KEGG annotation analysis on all modules with a coefficient of > 0.7. Genes encoding enzymes for primary bioactive compound synthesis were identified across multiple co-expression modules (darkorange, darkgreen, cyan), including Phenylalanine ammonia-lyase (PAL, GinfChr2G00031740) and Chalcone isomerase (CHI, GinfChr1G00083500; GinfChr1G00083530). Within the black module (138 genes), we detected:Five PAL genes (GinfChr1G00091980, GinfChr2G00031750, GinfChr2G00031760, GinfChr5G00197500, GinfChr1G00120940)/One Chalcone reductase (CHR, GinfChr5G00197280),One 3-Hydroxy-3-methylglutaryl-CoA reductase* (HMGR, GinfChr5G00173040),One Cinnamate 4-hydroxylase* (C4H, GinfChr1G00085290), Two 4-Coumarate: CoA ligases* (4CL, GinfChr4G00082360, GinfChr5G00165360), Four Chalcone synthases (CHS, GinfChr3G00137770, GinfChr4G00071830,GinfChr4G00071850).KEGG annotation revealed significant enrichment (P<0.001) in five functional categories (Figure 9): Metabolism (104 genes), Genetic Information Processing, Environmental Information Processing, Cellular Processes, and Organismal Systems. Metabolic genes predominated, with 24 genes each in secondary metabolite biosynthesis and amino acid metabolism. Crucially, this module contained:Phloridzin synthase (flavonoid pathway),2-Hydroxyisoflavanone dehydratase* (HIDH),2,7,4’-Trihydroxyisoflavanone 4’-O-methyltransferase* (HI4OMT),cis-Zeatin O-glucosyltransferase (cis-ZOG, EC:2.4.1.-). These findings indicate coordinated regulation of flavonoid diversification and phytohormone signaling (Table 2). We propose that salt stress induces crosstalk between hormone pathways (zeatin, gibberellins, ethylene, jasmonate) and specialized metabolism, particularly glabridin biosynthesis. Competitive flux partitioning may occur, where synthesis of certain flavonoids/isoflavonoids potentially downregulates glabridin-associated genes.

Figure 8

Figure 9

3.9 Gene expression and metabolite accumulation in flavonoid biosynthesis

Based on WGCNA and transcriptome results, we screened genes associated with active ingredient synthesis from the target modules and enrichment analysis. The accumulation levels of major active ingredients across different time points are presented in Figures 10 and 6. As shown in Figure 10, in the T1 Glycyrrhiza inflata germplasm, 18 genes within the glycyrrhizin biosynthetic pathway (e.g., GiHMGS1, GiHMGS2, GiHMGR1, GiHMGR2) exhibited significant upregulation 6 hours post-stress. Combined with the accumulation results of active compounds across different stress durations, this pronounced upregulation likely contributed to T1’s substantially higher glycyrrhizin accumulation at 12 hours compared to other germplasms. While the T2 germplasm also showed upregulation of most genes in this pathway, the magnitude of expression change was less pronounced than in T1. In contrast, the stress-sensitive germplasms S1 and S2 reached peak expression levels for most glycyrrhizin pathway genes at 2 hours post-stress, followed by a decline at 6 hours. This expression pattern corresponds to their lower glycyrrhizin accumulation at 6 and 12 hours relative to T1 and T2.Conversely, within the Phenylpropanoid Metabolism Pathway (PMP), the T1 germplasm showed significant downregulation of gene expression at 2 and 6 hours post-stress (Figure 6), exhibiting an expression pattern opposite to that observed in the glycyrrhizin pathway. However, its constitutive (control) expression levels for PMP genes were notably higher than those in the other three germplasms. This elevated baseline likely explains T1’s sustained higher levels of flavonoid and isoflavonoid active ingredients during short-term stress, consistent with previous findings. Germplasms T2, S1, and S2 exhibited upregulation of PMP pathway genes at 2 and 6 hours post-stress, suggesting resource allocation towards synthesizing active compounds like Liquiritin, Licochalcone A, and Glabridin. Notably, S1 displayed higher expression of GiCHI1, GiCHI4, GiIFS1, and GiGT1 at 6 hours compared to the other germplasms, resulting in significantly greater Glabridin synthesis post-stress. T2, meanwhile, upregulated genes such as GiCHS3, GiCHS4, GiCHS5, GiCHR1, and GiCHR2, potentially facilitating increased Isoliquiritigenin and subsequent Liquiritin production.

Figure 10

4 Discussion

4.1 Analysis of quality traits and salt-alkali stress tolerance in Glycyrrhiza inflata germplasm resources from Tarim basin

This study revealed that all 29 tested accessions of Glycyrrhiza inflata exhibited germination rates exceeding 45% under saline-alkaline stress (pH 8.15), with a mean germination rate of 65.13%. Notably, several accessions achieved rates surpassing 80%. Plant height ranged from 1.74 m to 3.24 m, collectively demonstrating the inherent capability of these diverse G. inflata germplasms from the Tarim Basin periphery to thrive in desert saline-alkali soils. Furthermore, the accessions displayed substantial variation in key bioactive constituents: total flavonoid content ranged from 17.926 ± 2.810 mg/g to 28.966 ± 5.533 mg/g, glycyrrhizic acid content from 2.056 ± 0.231 mg/g to 7.506 ± 1.480 mg/g, and liquiritin content from 0.717 ± 0.139 mg/g to 4.008 ± 1.232 mg/g. These results, particularly the high levels of the primary quality indicators glycyrrhizic acid and liquiritin, highlight.

4.2 Changes in key bioactive compounds, endogenous hormones, and antioxidant enzyme activities in Glycyrrhiza inflata

The primary bioactive constituents of Glycyrrhiza inflata are triterpenoids and flavonoids (Yang et al., 2017), These compounds not only exhibit free radical scavenging capacity and antioxidant activity (Kumar and Pandey, 2013), but also function as plant growth regulators (Kumar and Pandey, 2013), playing crucial roles in plant adaptation to environmental conditions and defense against abiotic stress (Shen et al., 2022), Furthermore, studies in both model (Ma et al., 2019) and non-model organisms (Zhang et al., 2023) have demonstrated that appropriate salt stress can effectively enhance desirable traits in cultivated germplasm. Consequently, we subjected roots of four G. inflata germplasms to hydroponic treatment with 150 mM NaCl (the previously determined optimal concentration in our laboratory) for 0h,2h, 6h. Subsequently, we measured changes in physiological parameters, endogenous hormone levels, and key bioactive compound concentrations. Following salt stress exposure, all four germplasms exhibited a characteristic triphasic response in CAT, SOD, and POD enzyme activities: an initial increase, followed by a decline, and then a secondary rise, ultimately stabilizing at levels slightly above or comparable to the control group (Figure 3). Generally, the salt-tolerant germplasms (T1, T2) displayed higher POD, and CAT activities and lower MDA content across most time points compared to the salt-sensitive germplasms (S1, S2). Chlorophyll and free proline content showed no significant changes post-stress, consistent with experimental expectations. This dynamic pattern aligns with findings reported by Ma et al. (Ma et al., 2025b)and Zhang et al(Zhang et al., 2023).Concurrently, endogenous hormone levels in all four germplasms underwent varying degrees of fluctuation, representing a physiological effort to mitigate the detrimental effects of salt stress on plant growth. Changes in abscisic acid (ABA), indole-3-acetic acid (IAA), and zeatin (ZT) levels are shown in Figures 4A–C. Notably, S1 displayed a significant decrease in ABA post-stress contrary to its established role as a key stress-responsive hormone in plants (Baron et al., 2012) which may contribute to its lower stress resistance. Salt stress typically triggers resource reallocation from growth to defense mechanisms, reducing auxin (IAA) demand and synthesis (Aizaz et al., 2024). Consistent with this, tolerant germplasms T1 and T2 showed significant IAA downregulation at 2h and 6h post-stress. In contrast, sensitive types S1 and S2 exhibited initial IAA increases; this delayed or attenuated stress response likely contributes to their reduced tolerance. Similarly, ZT levels in T1 remained stable throughout the 6h stress period, mirroring its ABA and CAT response patterns. T2, S1, and S2 displayed significant ZT increases at 2h followed by declines at 6h, representing an adaptive effort to mitigate stress impact. Numerous studies confirm that plants enhance flavonoid synthesis under salt stress to counteract its effects (Xu et al., 2020; Feng et al., 2023), establishing controlled salt exposure as an effective strategy for boosting bioactive compound production. As shown in Figure 5, all four germplasms exhibited significant upregulation of key bioactive compounds post-stress. Intriguingly, under control conditions, one-year-old plants from all four germplasms—particularly T1, T2, and S2—showed no significant differences in major bioactive compound levels (Figure 2). However, root exposure to 150 mM NaCl induced distinct accumulation patterns: T1/T2 accumulated substantially higher Glycyrrhizic acid than S1/S2; S1/S2 accumulated substantially higher Glabridin and LicochalconeA than T1/T2; Integrating these findings with preceding results, we infer that differential stress-induced expression of genes within flavonoid and glycyrrhizin biosynthetic pathways across the four germplasms underlies these compound-specific accumulation patterns.

4.3 Salt stress-responsive gene expression

Transcriptome analysis in this study revealed substantial variation in the number of differentially expressed genes (DEGs) across different treatment combinations (Figure 7). The T2_6h/T2_0h comparison exhibited the highest number of upregulated genes (3,093, representing 6.43% of expressed genes), alongside 2,567 downregulated genes (5.96%). Conversely, the S2_6h/S2_0h comparison showed relatively fewer DEGs. Overall, most treatment groups induced significant transcriptional reprogramming, characterized by substantial numbers of both upregulated and downregulated genes. This reflects the complexity and diversity of gene expression regulation in different Glycyrrhiza inflata germplasms as they adapt physiologically to varying environmental or experimental conditions. Notably, the number of DEGs at 2h post-stress was slightly higher than at 6h.Gene Ontology (GO) enrichment analysis (Figure 11) indicated that G. inflata orchestrates a multifaceted response to salt stress by differentially regulating genes involved in diverse functions. These include, but are not limited to, flavonoid biosynthesis, modulation of oxidoreductase activity, plant hormone signal transduction, and amino acid degradation. Consistent with the DEG count data, GO terms were enriched with a greater number of genes at 2h compared to 6h, suggesting a more extensive transcriptional response during the initial phase of stress acclimation. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment (Figure 12) further highlighted significant differential expression of genes within key metabolic pathways post-stress, notably flavonoid biosynthesis, flavone and flavonol biosynthesis, and phenylalanine metabolism. The number of differentially expressed genes within these pathways varied considerably among the germplasms, potentially contributing to the observed differences in bioactive compound accumulation.

Figure 11

Figure 12

4.4 Biosynthetic regulatory gene network for key bioactive compounds in Glycyrrhiza inflata

Recent years have witnessed growing research interest in the clinical efficacy (Liu et al., 2025) therapeutic applications (Liu et al., 2021), and cosmetic utilization (Jiang et al., 2015) of bioactive compounds derived from licorice (Glycyrrhiza spp.). However, the specific biosynthetic mechanisms governing key active constituents, such as flavonoids and glycyrrhizin, in Glycyrrhiza inflata remain incompletely elucidated (Zhao et al., 2023). Furthermore, cultivated licorice, the primary supply source, faces challenges of declining quality, which constrains the development of related industries. Advances in bioinformatics have positioned Weighted Gene Co-expression Network Analysis (WGCNA) as a powerful tool for identifying target genes and key transcription factors (TFs) from RNA-seq data (He et al., 2023).In this study, we employed transcriptome analysis coupled with WGCNA to identify genes associated with bioactive compound accumulation in G. inflata under salt stress and constructed a putative biosynthetic network. Intriguingly, we identified a module (Black) exhibiting a strong negative correlation with Glabridin accumulation. Surprisingly, this module encompassed not only TFs involved in plant hormone signaling pathways (e.g., JAZ, ERF1, PIF3) but also the glycyrrhizin pathway rate-limiting enzyme HMGR and numerous genes encoding enzymes in flavonoid biosynthesis (PAL, 4CL, CHS, CHR).Glycyrrhizin biosynthesis primarily occurs via the mevalonate (MVA) pathway. Within this pathway, 3-hydroxy-3-methylglutaryl-CoA reductase (HMGR), a key rate-limiting enzyme, catalyzes the conversion of HMG-CoA to mevalonate (MVA), leading to the formation of isopentenyl pyrophosphate (IPP) and dimethylallyl pyrophosphate (DMAPP) (Figure 10). Conversely, other flavonoid bioactive compounds originate from phenylalanine. Sequential catalysis by phenylalanine ammonia-lyase (PAL), cinnamate 4-hydroxylase (C4H), and 4-coumarate-CoA ligase (4CL) yields coumaroyl-CoA. This precursor is then channeled through various chalcone synthases (CHS) to synthesize distinct compounds like Glabridin, Liquiritin, and Licochalcone A (Figure 6).Current understanding suggests no direct regulatory interplay between these two pathways. Therefore, our WGCNA finding of a strong negative correlation between HMGR (a key enzyme in glycyrrhizin synthesis) and Glabridin accumulation was unexpected. However, integrating this result with phenotypic data (Figures 5A, B) revealed a compelling pattern: the T1 germplasm exhibited the highest glycyrrhizic acid accumulation post-stress but the lowest Glabridin levels among the four germplasms. Conversely, the S1 germplasm accumulated the highest Glabridin but showed no significant change in glycyrrhizic acid within the first 12h post-stress. This inverse relationship led us to hypothesize a potential regulatory mechanism linking these seemingly independent pathways, with HMGR potentially playing a pivotal role. Supporting this hypothesis, transcriptomic data revealed that in the T1 germplasm, genes within the glycyrrhizin pathway (including HMGR) were highly upregulated at 6h post-stress, significantly exceeding levels in other germplasms. Conversely, several genes in the Glabridin pathway (PAL, C4H, 4CL) were downregulated. This coordinated differential expression across multiple genes in both pathways provides corroborative evidence for potential cross-talk. Furthermore, the Glabridin-correlated Black module contained CHR genes, which catalyze the formation of isoliquiritigenin. Notably, isoliquiritigenin, Liquiritin, Glabridin, and Licochalcone A all share coumaroyl-CoA as a common precursor. Their divergence hinges on the specific CHS isoform involved: CHS enzymes producing naringenin chalcone lead to isoliquiritigenin and Liquiritin, while those producing chalcone lead to Glabridin and Licochalcone A. Consequently, we infer that specific CHS isoforms (e.g., GiCHS4, GiCHS5, GiCHS6, which were highly negatively correlated with Glabridin) may preferentially channel coumaroyl-CoA towards alternative flavonoid branches, thereby influencing Glabridin flux. Additionally, the module contained annotated TFs such as ERF1, JAZ, and MYC. ERF1, established in Arabidopsis thaliana, integrates environmental signals to promote local auxin (IAA) accumulation (Li et al., 2021). Elevated IAA typically diverts resources away from stress-related compound synthesis. JAZ proteins act as repressors in jasmonate (JA) signaling; their high expression dampens JA responses (He et al., 2023). Given that JA signaling enhances flavonoid accumulation in plants like tomato (Li et al., 2021), the presence of these TFs within the network further validates the biological relevance of our WGCNA findings.

5 Conclusion

Glycyrrhiza inflata is a halophytic plant thriving in desert environments, valued for its diverse medicinal compounds. It offers dual benefits: remediating saline-alkaline soils while generating local economic value. Preliminary screening of 29 germplasms from the Tarim Basin periphery revealed significant variation in bioactive compound profiles and germination rates under saline-alkaline conditions. Consequently, we selected two salt-sensitive germplasms (S1, S2) and two stress-tolerant germplasms (T1, T2) for controlled salt-stress experiments. Key findings indicate that one-year-old plants from all four germplasms showed no significant differences in total flavonoids or five target bioactive compounds under optimal growth conditions. However, salt stress significantly enhanced the accumulation of all five bioactive compounds, with their content exhibiting a time-dependent increasing trend over the stress duration. Concurrently, short-term salt stress boosted antioxidant enzyme activities as an immediate defense response. All four germplasms underwent significant metabolic reprogramming, characterized by substantial alterations in transcriptomic profiles, flavonoid biosynthetic pathways, and phytohormone signaling cascades, ultimately leading to divergent stress tolerance phenotypes. This study not only advances our understanding of the molecular regulatory mechanisms governing flavonoid biosynthesis under salt stress but also provides crucial genetic insights for the targeted improvement of flavonoid production in saline environments. Furthermore, it establishes a foundation for identifying elite G. inflata germplasms with superior bioactive compound yields and enhanced stress resilience. Such germplasms hold significant promise for combating soil salinization in the southern Xinjiang deserts and other arid regions globally.

References

• 1

  AizazM.LubnaJanR.AsafS.BilalS.KimK. M.et al. (2024). Regulatory dynamics of plant hormones and transcription factors under salt stress. Biol. (Basel)13. doi: 10.3390/biology13090673

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 2

  BaronK. N.SchroederD. F.StasollaC. (2012). Transcriptional response of abscisic acid (ABA) metabolism and transport to cold and heat stress applied at the reproductive stage of development in Arabidopsis thaliana. Plant Sci.188-189, 48–59. doi: 10.1016/j.plantsci.2012.03.001

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 3

  BuD.LuoH.HuoP.WangZ.ZhangS.HeZ.et al. (2021). KOBAS-i: intelligent prioritization and exploratory visualization of biological functions for gene enrichment analysis. Nucleic Acids Res.49, W317–w325. doi: 10.1093/nar/gkab447

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 4

  ChenK. Y.Toro-MorenoM.SubramaniamA. R. (2025). GitHub enables collaborative and reproducible laboratory research. PLoS Biol.23, e3003029. doi: 10.1371/journal.pbio.3003029

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 5

  DaiS.FengW.SongF.LiT.TaoY.YangF.et al. (2025). Review of biological algal fertilizer technology: Alleviating salinization, sequestering carbon, and improving crop productivity. Bioresour Technol.429, 132507. doi: 10.1016/j.biortech.2025.132507

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 6

  DuY.LiuX.ZhangL.ZhouW. (2023). Drip irrigation in agricultural saline-alkali land controls soil salinity and improves crop yield: Evidence from a global meta-analysis. Sci. Total Environ.880, 163226. doi: 10.1016/j.scitotenv.2023.163226

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 7

  DuanQ.ZhuZ.WangB.ChenM. (2022). Recent progress on the salt tolerance mechanisms and application of tamarisk. Int. J. Mol. Sci.23. doi: 10.3390/ijms23063325

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 8

  FengS.YaoY. T.WangB. B.LiY. M.LiL.BaoA. K. (2023). Flavonoids are involved in salt tolerance through ROS scavenging in the halophyte Atriplex canescens. Plant Cell Rep.43, 5. doi: 10.1007/s00299-023-03087-6

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 9

  GuJ.JiaT.MaM. (2024). Influence of habitat and effects of salt stress on biochemical and physiological parameters of glycyrrhiza uralensis. Plants (Basel)13. doi: 10.3390/plants13152108

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 10

  HassaniA.AzapagicA.ShokriN. (2021). Global predictions of primary soil salinization under changing climate in the 21st century. Nat. Commun.12, 6663. doi: 10.1038/s41467-021-26907-3

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 11

  HeK.DuJ.HanX.LiH.KuiM.ZhangJ.et al. (2023). Phosphate starvation response1 (phr1) interacts with jasmonate zim-domain (JAZ) and MYC2 to modulate phosphate deficiency-induced jasmonate signaling in Arabidopsis. Plant Cell35, 2132–2156. doi: 10.1093/plcell/koad057

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 12

  JiangM.ZhaoS.YangS.LinX.HeX.WeiX.et al. (2020). An "essential herbal medicine"-licorice: A review of phytochemicals and its effects in combination preparations. J. Ethnopharmacol249, 112439. doi: 10.1016/j.jep.2019.112439

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 13

  JiangZ.ZhouX.LiR.MichalJ. J.ZhangS.DodsonM. V.et al. (2015). Whole transcriptome analysis with sequencing: methods, challenges and potential solutions. Cell Mol. Life Sci.72, 3425–3439. doi: 10.1007/s00018-015-1934-y

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 14

  JiaoY.ZhangJ.PanC. (2022). Integrated physiological, proteomic, and metabolomic analyses of pecan cultivar 'Pawnee' adaptation to salt stress. Sci. Rep.12, 1841. doi: 10.1038/s41598-022-05866-9

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 15

  KimD.PaggiJ. M.ParkC.BennettC.SalzbergS. L. (2019). Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat. Biotechnol.37, 907–915. doi: 10.1038/s41587-019-0201-4

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 16

  KimD.PerteaG.TrapnellC.PimentelH.KelleyR.SalzbergS. L. (2013). TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol.14, R36. doi: 10.1186/gb-2013-14-4-r36

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 17

  KlopfensteinD. V.ZhangL.PedersenB. S.RamírezF.Warwick VesztrocyA.NaldiA.et al. (2018). GOATOOLS: A Python library for Gene Ontology analyses. Sci. Rep.8, 10872. doi: 10.1038/s41598-018-28948-z

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 18

  KumarS.PandeyA. K. (2013). Chemistry and biological activities of flavonoids: an overview. ScientificWorldJournal2013, 162750. doi: 10.1155/2013/162750

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 19

  LiB.DeweyC. N. (2011). RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinf.12, 323. doi: 10.1186/1471-2105-12-323

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 20

  LiF.GuoY.WangZ.MuY. (2022). Influence of different phytoremediation on soil microbial diversity and community composition in saline-alkaline land. Int. J. Phytoremediation24, 507–517. doi: 10.1080/15226514.2021.1955240

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 21

  LiX.XuY.ZhangJ.XuK.ZhengX.LuoJ.et al. (2024). Integrative physiology and transcriptome reveal salt-tolerance differences between two licorice species: Ion transport, Casparian strip formation and flavonoids biosynthesis. BMC Plant Biol.24, 272. doi: 10.1186/s12870-024-04911-1

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 22

  LiZ.PengR.YaoQ. (2021). SlMYB14 promotes flavonoids accumulation and confers higher tolerance to 2,4,6-trichlorophenol in tomato. Plant Sci.303, 110796. doi: 10.1016/j.plantsci.2020.110796

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 23

  LiuW.FengY.YuS.FanZ.LiX.LiJ.et al. (2021). The flavonoid biosynthesis network in plants. Int. J. Mol. Sci.22. doi: 10.3390/ijms222312824

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 24

  LiuW.LinL.ZhangZ.LiuS.GaoK.LvY.et al. (2019). Gene co-expression network analysis identifies trait-related modules in Arabidopsis thaliana. Planta249, 1487–1501. doi: 10.1007/s00425-019-03102-9

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 25

  LiuX.OuX.ZhangT.LiX.QiaoQ.JiaL.et al. (2024). In situ neutrophil apoptosis and macrophage efferocytosis mediated by Glycyrrhiza protein nanoparticles for acute inflammation therapy. J. Control Release369, 215–230. doi: 10.1016/j.jconrel.2024.03.029

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 26

  LiuJ.XuX.JianM.GuoY.ZhaiL.SunG.et al. (2025). Glycyrrhiza glabra extract as a skin-whitening Agent: Identification of active components and CRTC1/MITF pathway-inhibition mechanism. J. Ethnopharmacol349, 119948. doi: 10.1016/j.jep.2025.119948

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 27

  LivakK. J.SchmittgenT. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta C(T)) Method. Methods25, 402–408. doi: 10.1006/meth.2001.1262

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 28

  LoveM. I.HuberW.AndersS. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol.15, 550. doi: 10.1186/s13059-014-0550-8

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 29

  MaF.HuangX.YangZ.NiuP.PangK.WangM.et al. (2025a). Suitable habitat prediction and decertified landscape remediation potential of three medicinal Glycyrrhiza species in China. Sci. Rep.15, 11198. doi: 10.1038/s41598-025-91874-4

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 30

  MaL.LianY.LiS.FahimA. M.HouX.LiuL.et al. (2025b). Integrated transcriptome and metabolome analysis revealed molecular regulatory mechanism of saline-alkali stress tolerance and identified bHLH142 in winter rapeseed (Brassica rapa). Int. J. Biol. Macromol295, 139542. doi: 10.1016/j.ijbiomac.2025.139542

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 31

  MaS.LvL.MengC.ZhouC.FuJ.ShenX.et al. (2019). Genome-wide analysis of abscisic acid biosynthesis, catabolism, and signaling in sorghum bicolor under saline-alkali stress. Biomolecules9. doi: 10.3390/biom9120823

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 32

  MengY.CaiX.CongS.SunJ.DuW.CuiH.et al. (2024). Diammonium glycyrrhizinate inhibited inflammatory response and modulated serum metabolism in poly(i:c)-induced pneumonia model mice. Shock61, 905–914. doi: 10.1097/shk.0000000000002353

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 33

  QiJ.FuD.WangX.ZhangF.MaC. (2023). The effect of alfalfa cultivation on improving physicochemical properties soil microorganisms community structure of grey desert soil. Sci. Rep.13, 13747. doi: 10.1038/s41598-023-41005-8

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 34

  RahmanM. M.MostofaM. G.KeyaS. S.SiddiquiM. N.AnsaryM. M. U.DasA. K.et al. (2021). Adaptive mechanisms of halophytes and their potential in improving salinity tolerance in plants. Int. J. Mol. Sci.22. doi: 10.3390/ijms221910733

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 35

  ShenN.WangT.GanQ.LiuS.WangL.JinB. (2022). Plant flavonoids: Classification, distribution, biosynthesis, and antioxidant activity. Food Chem.383, 132531. doi: 10.1016/j.foodchem.2022.132531

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 36

  SunT.MaN.WangC.FanH.WangM.ZhangJ.et al. (2021). A golgi-localized sodium/hydrogen exchanger positively regulates salt tolerance by maintaining higher K(+)/na(+) ratio in soybean. Front. Plant Sci.12. doi: 10.3389/fpls.2021.638340

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 37

  TangL.ZhanL.HanY.WangZ.DongL.ZhangZ. (2023). Microbial community assembly and functional profiles along the soil-root continuum of salt-tolerant Suaeda glauca and Suaeda salsa. Front. Plant Sci.14. doi: 10.3389/fpls.2023.1301117

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 38

  van ZelmE.ZhangY.TesterinkC. (2020). Salt tolerance mechanisms of plants. Annu. Rev. Plant Biol.71, 403–433. doi: 10.1146/annurev-arplant-050718-100005

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 39

  WangC.CaiZ.ShiJ.ChenS.TanM.ChenJ.et al. (2019). Comparative metabolite profiling of wild and cultivated licorice based on ultra-fast liquid chromatography coupled with triple quadrupole-time of flight tandem mass spectrometry. Chem. Pharm. Bull. (Tokyo)67, 1104–1115. doi: 10.1248/cpb.c19-00423

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 40

  WangC.ChenL.CaiZ. C.ChenC.LiuZ.LiuX.et al. (2020). Comparative Proteomic Analysis Reveals the Molecular Mechanisms Underlying the Accumulation Difference of Bioactive Constituents in Glycyrrhiza uralensis Fisch under Salt Stress. J. Agric. Food Chem.68, 1480–1493. doi: 10.1021/acs.jafc.9b04887

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 41

  WangL.FengZ.WangX.WangX.ZhangX. (2010). DEGseq: an R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics26, 136–138. doi: 10.1093/bioinformatics/btp612

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 42

  WuY.ZhangF.KuangD.LiD.YanJ.YangJ.et al. (2024). Efficacy of botanical lozenges in the treatment of chronic pharyngitis: a randomized controlled trial. Front. Pharmacol.15. doi: 10.3389/fphar.2024.1162883

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 43

  XuN.LiuS.LuZ.PangS.WangL.WangL.et al. (2020). Gene Expression Profiles and Flavonoid Accumulation during Salt Stress in Ginkgo biloba Seedlings. Plants (Basel)9. doi: 10.3390/plants9091162

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 44

  YanB.HouJ.LiW.LuoL.YeM.ZhaoZ.et al. (2023). A review on the plant resources of important medicinal licorice. J. Ethnopharmacol301, 115823. doi: 10.1016/j.jep.2022.115823

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 45

  YangR.YuanB. C.MaY. S.ZhouS.LiuY. (2017). The anti-inflammatory activity of licorice, a widely used Chinese herb. Pharm. Biol.55, 5–18. doi: 10.1080/13880209.2016.1225775

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 46

  ZhangZ.FangJ.ZhangL.JinH.FangS. (2023). Genome-wide identification of bHLH transcription factors and their response to salt stress in Cyclocarya paliurus. Front. Plant Sci.14. doi: 10.3389/fpls.2023.1117246

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 47

  ZhaoP.ZhangJ.ChenS.ZhangZ.WanG.MaoJ.et al. (2023). ERF1 inhibits lateral root emergence by promoting local auxin accumulation and repressing ARF7 expression. Cell Rep.42, 112565. doi: 10.1016/j.celrep.2023.112565

  • Pubmed Abstract
  • CrossRef
  • Google Scholar