Front. Plant Sci. Frontiers in Plant Science Front. Plant Sci. 1664-462X Frontiers Media S.A. 10.3389/fpls.2024.1383918 Plant Science Original Research Effects of elicitors from culture filtrate of Fusarium solani CL105 on flavonoid production of Scutellaria baicalensis calli Cui Xiaoxuan 1 2 3 4 † Zhang Xin 1 2 3 † Sun Huigai 1 2 3 Zheng Yuguang 2 3 5 * Su Chunyan 1 2 3 * 1 College of Pharmacy, Hebei University of Chinese Medicine, Shijiazhuang, China 2 Traditional Chinese Medicine Processing Technology Innovation Center of Hebei Province, School of Pharmacy, Hebei University of Chinese Medicine, Shijiazhuang, China 3 International Joint Research Center on Resource Utilization and Quality Evaluation of Traditional Chinese Medicine of Hebei Province, School of Pharmacy, Hebei University of Chinese Medicine, Shijiazhuang, China 4 State Key Laboratory for Quality Ensurance and Sustainable Use of Dao-di Herbs, National ResourceCenter for Chinese Materia Medica, China Academy of Chinese Medical Sciences, Beijing, China 5 Department of Pharmaceutical Engineering, Hebei Chemical and Pharmaceutical College, Shijiazhuang, China

Edited by: Ahmed Abdelkhalek, City of Scientific Research and Technological Applications, Egypt

Reviewed by: Eleonora Harries, National Scientific and Technical Research Council (CONICET), Argentina

Said Ibrahim Behiry, Alexandria University, Egypt

*Correspondence: Yuguang Zheng, zyg314@163.com; Chunyan Su, suchunyan_198587@163.com

†These authors have contributed equally to this work

 04 06 2024 2024 15 1383918 08 02 2024 13 05 2024 Copyright © 2024 Cui, Zhang, Sun, Zheng and Su 2024 Cui, Zhang, Sun, Zheng and Su

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.

 Introduction

Endophytic fungi can promote secondary metabolite accumulation in medicinal plants. Previously, we observed that the culture filtrate of Fusarium solani CL105 promoted flavonoid production in Scutellaria baicalensis calli. However, the active ingredients and mechanisms associated with this secondary metabolite accumulation remain unclear.

 Methods

This study evaluates the effects of different elicitors from the culture filtrate of F. solani CL105 namely, exopolysaccharide (EPS), exoprotein (EP), and other parts (OP), on the flavonoid production in S. baicalensis calli by HPLC. Subsequently, the underlying mechanism of EPS induced flavonoid production in S. baicalensis calli was revealed by transcriptomics and RT-PCR.

 Results and discussion

The results indicated a significant increase in flavonoid production in S. baicalensis calli following treatment with EPS. Baicalin (1.40 fold), wogonoside (1.91 fold), and wogonin (2.76 fold) were most significantly up-regulated compared with the control. Transcriptome analysis further revealed up-regulation of key enzyme genes (CHS, CHI, FNS, and F6H) involved in flavonoid synthesis after 5 days of EPS treatment. Moreover, the expression of GA2ox and CYP707A—genes involved in gibberellin acid (GA) and abscisic acid biosynthesis (ABA), respectively—were significantly up-regulated. The expression levels of certain transcription factors, including MYB3, MYB8, and MYB13, were also significantly higher than in controls. Our results indicated that EPS was a main active elicitor involved in promoting flavonoid production in S. baicalensis calli. We postulated that EPS might stimulate the expression of MYB3, MYB8, MYB13, GA2ox, and CYP707A, leading to markedly upregulated CHS, CHI, FNS, and F6H expression levels, ultimately promoting flavonoid synthesis. This study provides a novel avenue for large-scale in vitro production of flavonoids in S. baicalensis.

 elicitor endophyte gene expression flavonoids Scutellaria baicalensis transcription factors section-in-acceptance Plant Pathogen Interactions Introduction

Scutellaria baicalensis Georgi, a perennial herb belonging to the family Labiatae, is a well-known traditional Chinese medicinal herb grown in northern China. Flavonoids and their glycosides are the major compounds in S. baicalensis and are classified into two types, namely free flavonoids and flavonoid glycosides, of which baicalein and wogonin are the most abundant flavonoids, while baicalin and wogonoside are the most abundant flavonoid glycosides (Zhao et al., 2019a). These flavonoids have pharmacological effects, including antioxidative, anti-inflammatory, antitumor, antiviral, antimicrobial, hepatoprotective, and neuroprotective properties (Wang et al., 2018). Baicalin, in particular, has demonstrated potential in the treatment of COVID-19 (Miguel et al., 2023). Furthermore, baicalin is clinically used as a post-marketing drug to treat acute and chronic hepatitis. However, the wild resources of S. baicalensis are continuing to gradually decrease, and the plant quality is limited by growth years (Yuan et al., 2010; Bai et al., 2020). Thus, tissue culturing provides a potential strategy to overcome these issues regarding S. baicalensis availability.

In S. baicalensis, flavonoids are synthesized via two pathways. The initial compound, phenylalanine, is produced through the shikimate pathway, which generates cinnamic acid via phenylalanine ammonialyase (PAL). In the aerial parts of S. baicalensis, cinnamic acid is converted to scutellarein by a series of enzymes, including cinnamate 4-hydroxylase (C4H), 4-coumarate CoA ligase (4CL), chalcone synthase (CHS), chalcone isomerase (CHI), flavone synthase (FNS), and flavone 6-hydroxylase (F6H). Alternatively, in the roots, baicalein, wogonin, and their glycosides are synthesized. This pathway involves the conversion of cinnamic acid to cinnamoyl-CoA in the presence of cinnamate-CoA ligase (CLL-7). Cinnamoyl-CoA is then converted to pinocembrin by CHS and CHI. Next, pinocembrin is converted by FNS to chrysin, which is converted to baicalein via F6H or to wogonin by the combined action of flavone 8-hydroxylase (F8H) and O-methyltransferase (OMT). Baicalein and wogonin then combine with sugar molecules via flavonoid 7-O-glucuronosyltransferase (UBGAT) to generate baicalin and wogonoside (Zhao et al., 2016, 2018; Zhao et al., 2019b).

The formation and accumulation of plant secondary metabolites are regulated by a complex interplay of plant hormones and transcription factors (TFs). Plant hormones, such as gibberellic acid (GA) and abscisic acid (ABA), reportedly promote flavonoid production in S. baicalensis (Yuan et al., 2013; He et al., 2023). Meanwhile, TFs, such as the NAC and MYB families, are key regulators of flavonoid synthesis (Yuan et al., 2013, 2015; Qian et al., 2020; Fang et al., 2023, 2023; He et al., 2023). The crosstalk between plant hormones and TFs adds another layer of regulation to secondary metabolite formation and accumulation (Bang et al., 2008; Yuan et al., 2013; Zhou et al., 2022).

RNA-seq is a robust quantitative tool for investigating gene regulation and function. In plants, this method can be applied to analyze cell and tissue transcriptional profiles and identify the key genes associated with secondary metabolite biosynthesis and their complex networks involving various hormones and TFs. In fact, Cao et al. (2022) adopted transcriptomics to demonstrate that the endophytic fungus KL27 enhances taxol accumulation in Taxus chinensis by regulating specific hormones and TFs to promote the expression of key genes in taxol biosynthesis. Similarly, through transcriptomic analysis, Guo et al. (2018) identified the key plant hormones and TFs involved in light-regulated anthocyanin biosynthesis in cherries.

Endophytic fungi are microorganisms that reside in the internal tissues of living plants and do not cause any immediate or overt adverse effects (Hardoim et al., 2015). Endophytic fungi have various functions in host plants, including promoting host growth and secondary metabolite accumulation (Jia et al., 2016). Previous studies have found that several material bases, such as polysaccharide, oligosaccharide, protein, and polypeptide (Li et al., 2013; Yu et al., 2016; Zhang et al., 2019; Wu et al., 2019a), from fungal mycelial extracts or culture filtrates, can promote host growth and secondary metabolite accumulation. Xu et al. (2023) used the endophytic fungus Colletotrichum sp. AP12 to increase the expression of key genes in the andrographolide biosynthesis pathway, promoting its synthesis and accumulation in Andrographis paniculata. Additionally, Ming et al. (2013) found that the endophytic fungus Trichoderma atroviride promotes the expression of genes related to the tanshinone biosynthesis pathway, leading to its synthesis in Salvia miltiorrhiza. However, the precise mechanisms through which endophytic fungi facilitate flavonoid synthesis in S. baicalensis remain unclear.

We previously found that treating S. baicalensis calli with the culture filtrate of Fusarium solani CL105 enhanced cell growth and flavonoid production (Zhang et al., 2023). Accordingly, the current study analyzes the effects of three elicitors from the F. solani CL105 culture filtrate, namely exopolysaccharide (EPS), exoprotein (EP), and other parts (OP), on the growth and flavonoid production of S. baicalensis calli. We then employ transcriptome analysis and quantitative real-time reverse transcription polymerase chain reaction (qRT-PCR) to explore how the active elicitor regulates flavonoid production in S. baicalensis calli. This study provides new insights into the mechanism by which endophytic fungi enhance flavonoid accumulation in S. baicalensis. Additionally, it presents new ideas for large-scale flavonoid production using biotechnological methods.

 Materials and methods Calli induction

S. baicalensis seeds, which were identified as S. baicalensis Georgi by Professor Yuguang Zheng of Hebei University of Chinese Medicine, were purchased from Anguo Herbal Market, Hebei Province, China. The seeds were surface sterilized according to the methods described in our previous study (Zhang et al., 2023), inoculated on a MS medium (Hopebio, Qingdao, China) (Murashige and Skoog, 1962), and cultured at 25 ± 1°C for 14 days to obtain sterile seedlings. Calli induction in S. baicalensis was performed as previously described, with some modifications (Wan et al., 2012). Briefly, stem explants (grown in vitro) were cut aseptically into approximately 10 mm segments and placed on MS medium supplemented with 2.0 mg/L 6-BA (6-Benzylaminopurine) (Xiya Reagent, Linyi, China) and 1.0 mg/L NAA (-naphthalene acetic acid) (Solarbio, Beijing, China), which could promote the induction and proliferation of calli (Xu et al., 2022). After 2 weeks, S. baicalensis calli were subcultured in the same culture medium at 25 ± 1°C.

 Preparation of elicitors

The endophytic fungus F. solani CL105 was grown in 100 mL Potato Dextrose Broth (PDB) medium in a 250 mL Erlenmeyer flask. The fungal culture was incubated at 25 °C with shaking at 180 rpm for 7 days. After fermentation, the fermented broth was filtered, the filtrate was collected and concentrated to a third of the original volume, then mixed with four volumes of 95% ethanol and incubated at 4°C for 48 h, then filtered to obtain the supernatants and precipitates. The precipitate from the ethanol dispersion was collected as EPS, and crude protein was collected from the supernatant using the Sevage method (supernatant:liquid:chloroform:n-Butyl alcohol = 25:4:1, v:v:v), as EP; the crude protein elicitor was removed from the supernatant, as was OP. All the elicitors were dried using a vacuum freeze dryer (TF-FD-1, Shanghai Tianfeng Industrial Co., Ltd., Shanghai, China) and stored at 4°C.

 Elicitor treatment

We have previously shown that 5% (v/v) culture filtrate elicitors of F. solani CL105 could induce flavonoid production in S. baicalensis calli (Zhang et al., 2023). The EPS, EP, and OP were respectively dissolved in 25 mL MS medium to prepare a series of concentrations based on their yields (760 mg/L, 192 mg/L and 1040 mg/L) in culture filtrate and the active concentrate (5%, v/v) of culture filtrate (EPS: 20, 40, 80 and 160 mg/L; EP: 5, 10 and 15 mg/L; OP: 20, 50 and 80 mg/L). MS medium without an elicitor was used as a control. The media was sterilized at 121°C for 30 min. Approximately 1.50 g of S. baicalensis calli were subcultured into media and incubated at 25 ± 1°C. The experiments were repeated thrice for each elicitor concentration. S. baicalensis calli were harvested on 0, 5, 10, and 15 days following EPS, EP, and OP treatment, and the fresh weight of S. baicalensis calli was measured. One part of S. baicalensis was frozen in liquid nitrogen and sent to Shanghai Majorbio Bio-pharm Biotechnology Co., Ltd. (Shanghai, China) for RNA sequencing. While the other part was dried at 50°C in a dry heat oven until a constant weight, and the dry weight (DW) was measured.

 Quantification of flavonoids with HPLC

For high-performance liquid chromatography (HPLC) analysis, 0.01 g samples of S. baicalensis calli powder were extracted with methanol under ultrasound for 30 min. The extracts were centrifuged at 12,000 rpm for 10 min and the supernatants were filtered through a 0.22 μm organic membrane. Analyses were performed using a Shimadzu LC2030 instrument (SHIMADZU, Kyoto, Japan). A DiKMA Diamonsil C18 (4.6×250 mm, 5.0 µm) column was applied for all analyses. Baicalin, baicalein, wogonoside, and wogonin, were detected and quantified by comparison with authentic standards (Biopurify Company, Chengdu, China). The mobile phases comprised A (water containing 0.1% formic acid) and B (acetonitrile). The gradient elution program was as follows: 25% B for 0–10 min, 25–45% B for 10–30 min, and 45–55% B for 30–55 min; the detection wavelength was 254 nm. The sample injection volume was 10 μL and the column temperature was maintained at 25°C. The method was validated in terms of linearity, precision, repeatability, stability, and recovery ( Supplementary Tables 1 , 2 ).

 Transcriptome determination and analysis Sample preparation, RNA extraction, library preparation, and transcriptome sequencing

Total RNA was extracted from the S. baicalensis calli using TRIzol® Reagent following the manufacturer’s instructions. RNA quality was determined using a 5300 Bioanalyzer (Agilent, Santa Clara, CA, USA) and quantified using an ND-2000 (NanoDrop Technologies, Wilmington, USA). Only high-quality RNAsample (OD260/280 = 1.8–2.2, OD260/230≥2.0, RIN≥6.5, 28S:18S≥1.0, >1 μg) was used to construct sequencing library. mRNA was enriched using oligo (dT) magnetic beads and was randomly fragmented using a fragmentation buffer. Under the action of reverse transcriptase, random hexamers were added to synthesize first-stranded cDNA using mRNA as a template. Thereafter, under DNA polymerase I, second-strand cDNA was synthesized. A poly (A) tail was added and connected to the sequence adaptors after end repair. The 300-bp cDNA target fragments were selected from the libraries on 2% low range ultra agarose gel. cDNA libraries were enriched by PCR amplification and then quantified using Qubit (v.4.0). A NovaSeq 6000 sequencer (2 × 150 bp read length) was used to sequence the paired-end RNA-seq (Wu et al., 2019b).

 <italic>De novo</italic> transcriptome assembly and gene functional annotation

Raw reads were filtered for quality using FASTQ (v.0.19.5), and the resulting clean reads were de novo-assembled using Trinity (v.2.8.5). To increase assembly quality, all assembled sequences were filtered using CD-Hit and translated.

All de novo-assembled unigenes were annotated using the following databases: NCBI non-redundant protein sequences (NR, http://www.ncbi.nlm.nih.gov); Protein family (Pfam, http://xfam.org/); Swiss-Prot (http://www.expasy.ch/sprot); Kyoto Encyclopedia of Genes and Genomes (KEGG, http://www.genome.jp) databases using Diamond to identify the proteins that had the highest sequence similarity with the given transcripts to retrieve their function annotations and a typical cut-off E-values less than 1.0×10–5 was set. The BLAST2GO program was used to obtain GO (http://www.geneontology.org) annotations of uniquely assembled transcripts to describe biological processes, molecular functions, and cellular components. Metabolic pathway analysis was performed using KEGG.

 Differentially expressed genes (DEGs) analysis

Gene expression levels were estimated using the RSEM (Li and Dewey, 2011) for each sample. Clean data were mapped onto assembled transcriptomes. The read count of each gene was obtained from the mapping results. Differential expression analysis of the two samples was performed using the DEseq2 package. Genes with P value < 0.05 and |log2 (fold change)| > 1 found by DEseq2 were assigned as differentially expressed.

 Quantitative real-time PCR (qRT-PCR) validation

To validate the RNAseq gene expression, total RNA was isolated from S. baicalensis calli and analyzed using quantitative real-time polymerase chain reaction (qRT-PCR). Four genes, namely CHS, CHI, FNS, and F6H, associated with flavonoid biosynthesis were selected from the RNA-Seq data for qRT-PCR analysis. The constitutively expressed 18S rRNA gene was used as the housekeeping gene (Zhou et al., 2014). All specific primer pairs used in this study are listed in Supplementary Table 3 . The 25 μL reaction mixture consisted of 1 μL Strand cDNA, 10µL qRT-PCR taq master mix (YEASEN, Shanghai, China), 0.7 µL forward primer (10 µM), 0.7 µL reverse primer (10 µM), and 7.6 µL highly pure water. The reaction was performed using an Thermofisher QuantStudio 1(Thermofisher, Waltham, MA, USA). The experiment was conducted in triplicates. The expression level of each target gene was determined using the 2-ΔΔCt method, normalized based on the CT value of the housekeeping gene.

 Statistical analysis

The fresh and dry weight, flavonoids contents, and relative-expression level of S. baicalensis calli were expressed as mean ± SD from three separate observations. Data were analyzed using one-way analysis of variance (ANOVA) and least significant difference (LSD) using SPSS19.0 software. p value < 0.05 was considered statistically significant.

 Results Effects of elicitors on <italic>S. baicalensis</italic> calli growth

As shown in Figure 1A , with 40 or 160 mg/L EPS treatments for 5 days, the fresh weight of S. baicalensis calli was significantly increased by ~1.18- and 1.17-fold compared with the control. Meanwhile, after treatment with 20, 40, 80, or 160 mg/L EPS for 5 days, the dry weight was significantly increased by ~1.20- to 1.39-fold compared with the control. After 10 days of culture, the fresh weight of S. baicalensis calli was significantly increased following treatment with 80 or 160 mg/L EPS by ~1.19- and 1.15-fold, respectively. However, reductions were observed in the fresh and dry weights at 15 days after culture. In particular, the dry weight of S. baicalensis calli was significantly reduced with 160 mg/L EPS treatment (~13% decrease) compared with the control.

Effect of EPS (A), EP (B) and OP (C) on growth in S. baicalensis calli. Data are presented as means ± SD, n = 3. *P<0.05, **P<0.01, ***P<0.001; #P<0.05, ##P<0.01, ###P<0.001.

The fresh weight of S. baicalensis calli continuously increased following EP treatment from 5 to 15 days, while the dry weight was significantly increased only after 5 days of treatment with 5 mg/L EP (1.12-fold; Figure 1B ). In contrast, the fresh and dry weights were significantly decreased following treatment with 10 mg/L (18% and 13%) or 15 mg/L (19% and 15%) EP for 5 days compared with the control. After 10 days of culture with 10 mg/L EP, the fresh weight of S. baicalensis calli significantly increased by ~1.30-fold compared with the control. After 15 days of treatment with 10 mg/L or 15 mg/L EP, the fresh weight significantly increased by ~1.15- and 1.22-fold, respectively, compared with the control.

Treatment with OP at 20, 50, or 80 mg/L did not significantly impact on the growth of S. baicalensis calli after 5, 10 or 15 days ( Figure 1C ).

 Effects of elicitors on <italic>S. baicalensis</italic> calli flavonoid production

The contents of baicalin, wogonoside, and wogonin in S. baicalensis calli were significantly increased following EPS treatment for 5–15 days ( Figure 2 ). After 5 days of culture with 20, 40, 80, or 160 mg/L EPS, baicalin and wogonin levels were significantly increased by ~1.17–2.76-fold. Meanwhile, the wogonoside content was significantly increased following treatment with 40, 80, or 160 mg/L EPS by ~1.51–1.91-fold compared with the control. After culturing for 10 days with any concentration of EPS, the contents of baicalin, wogonoside, and wogonin were significantly increased by ~1.10 to 2.03-fold compared with the control. Similarly, after 15 days, the baicalin content was significantly increased by ~1.19–1.22-fold following treatment with all EPS concentrations. However, after 15 days, the wogonoside content was only significantly increased following treatment with 20, 40, or 80 mg/L EPS (1.23–1.28-fold) and the wogonin content was significantly increased following 80 mg/L (2.11-fold) or 160 mg/L (1.64-fold) EPS treatment.

Effects of EPS on the accumulation of flavonoids in S. baicalensis calli on days 5, 10, and 15, respectively. The treatments were control, 20, 40, 80, and 160 mg/L. Data are presented as means ± SD, n = 3. *P<0.05, **P<0.01, ***P<0.001.

EP had a weak effect on flavonoids. After culturing for 5 days with 5 mg/L EP, the baicalin content was significantly increased by ~1.06-fold compared with the control ( Supplementary Figure 1 ). At the late stage of induction, OP significantly promoted flavonoid synthesis. After culturing for 15 days with 50 or 80 mg/L OP, the wogonoside contents were significantly increased by ~1.45- and 1.31-fold, respectively, while those of baicalein and wogonin were significantly increased following 80 mg/L OP treatment by ~2.48-fold and 1.61-fold, respectively, compared with the control ( Supplementary Figure 2 ).

 Transcriptomic analysis and differentially expressed genes

Since EPS is the main active elicitor, we performed transcriptome sequencing on EPS-treated S. baicalensis calli. In total, 212.11 Gb of clean data were generated from the samples. All clean reads were subjected to de novo assembly using Trinity, producing 94,597 transcripts and 94,597 unigenes. Sequencing data quality evaluations are listed in Supplementary Table 4 . The unigenes were functionally annotated based on the seven largest public databases. A total of 38,106 unigenes (40.28% of the total unigenes) were annotated; 19,209 unigenes (50.41% of the total annotated unigenes) were matched in the Gene Ontology (GO) database, and 9,942 (26.09%), 19,740 (51.80%), 22,154 (58.14%), 17,482 (45.88%), and 22,667 (45.27%) unigenes exhibited significant similarity to sequences in the GO, KEGG, eggNOG, NCBI non-redundant protein sequences (NR), Swiss-Prot, and Pfam databases, respectively.

All unigenes were searched against the GO database to classify their functions based on the NR annotation. The 19,209 unigenes assigned to one or more GO terms were classified into three main GO categories and 53 groups ( Supplementary Figure 3 ). Within the “biological process” domain, the most evident matches were “cellular process” (8,198), “metabolic process” (7,239), and “biological regulation” (3,074). In the “cellular component” domain, “cell part” (9,774), “membrane part” (7,431), and “organelle” (5,641) predominated. For the “molecular function” domain, the genes were primarily enriched in “binding” (10,479) and “catalytic activity” (10,123).

For further analysis, the unigenes were mapped onto the KEGG database to categorize gene function and identify biochemical pathways. A total of 9,942 unigenes were annotated and assigned to five main KEGG metabolic pathways, 18 sub-branches, and 5,814 KEGG pathways. The most common sub-branch was “translation” (954), followed by “carbohydrate metabolism” (862) and “folding, sorting, and degradation” (590) ( Supplementary Figure 4 ). An additional 147 unigenes matched “phenylpropanoid biosynthesis” (ko00940) and 54 unigenes matched “flavonoid biosynthesis” (ko00941).

A gene was designated as differentially expressed based on the following cutoff criteria: |Log2FC| ≥ 1 and P < 0.05. A total of 1,024, 1,715, 565, and 883 DEGs were identified in the EPS-treated for 0 d vs. control, EPS-treated for 5 days vs. control, EPS-treated for 10 days vs. control, and EPS-treated for 15 days vs. control comparisons, respectively. Of these DEGs, 150, 704, 181, and 310 were up-regulated, whereas 874, 1,011, 384, and 573 were down-regulated at 0, 5, 10, and 15 days, respectively ( Figure 3 ).

Number of DEGs at four comparisons. The abscissa represents EPS-treated versus untreated at different time points (0, 5, 10 and 15 days), and the ordinate represents the number of the up- and down-regulated DEGs.

A volcano plot was constructed to illustrate the distribution of the DEGs among the four comparison groups. The number of DEGs in the 5-day group was higher than in the 10- or 15-day groups. Moreover, the expression levels of a few genes in the 5-day group were significantly up-regulated with a higher degree of significance compared with the 10- or 15-day groups ( Supplementary Figure 5 ).

 DEGs involved in flavonoid biosynthesis and qRT-PCR validation

The flavonoid biosynthesis pathway in S. baicalensis has been previously described (Hu et al., 2022; Pei et al., 2022). To further analyze how the identified DEGs contribute to higher flavonoid accumulation following EPS treatment, their expression patterns in flavonoid biosynthesis pathways were analyzed ( Figure 4A ). As shown in Figure 4B , most genes involved in flavonoid biosynthesis were up-regulated in S. baicalensis calli following treatment with EPS after 5 days compared with the control. Particularly evident was the expression of genes encoding CHS, CHI, FNS, and F6H. However, ≤ 50% were up-regulated at 10 or 15 days. Hence, EPS promoted flavonoid accumulation by up-regulating genes encoding enzymes in the flavonoid biosynthesis pathway. Activation of the flavonoid biosynthesis pathway by EPS occurs in the early stages of stimulation and decreases over time.

DEGs involved flavonoid biosynthesis and qRT-PCR validation. (A) Biosynthetic pathways of the flavonoids in S. baicalensis. Enzymes abbreviations are: PAL: phenylalanine ammonialyase; C4H: cinnamate 4-hydroxylase; 4CL: 4-coumarate CoA ligase; CHS: chalcone synthase; CHI: chalcone isomerase; FNS: flavone synthase; F6H: flavone 6-hydroxylase; CLL-7: cinnamate-CoA ligase; F8H: flavone 8-hydroxylase; OMT: O-methyltransferase; UBGAT: flavonoid 7-O-glucuronosyltransferase. (B) Expression analysis of the flavonoid biosyntheisis-related unigenes at 5, 10 and 15days. The bar indicated the “log2 (fold change)”. (C) Effects of EPS on the expression of flavonoid biosynthesis key genes in S. baicalensis calli. Data are presented as means ± SD, n = 3. *P<0.05, ***P<0.001.

To assess the reliability of our transcriptomic data, four up-regulated genes associated with the flavonoid biosynthesis pathway in S. baicalensis calli during 5-day cultures were selected for validation via RNA-seq. The expression profiles for the four DEGs determined via qRT-PCR were consistent with those obtained via transcriptomics ( Figure 4C , Supplementary Table 5 ). Therefore, these findings suggest that the transcriptome information is reliable.

 DEGs involved in hormone biosynthesis

Phytohormones influence the expression of key enzymes in the flavonoid biosynthesis pathway (Yu et al., 2021). Our transcriptome data showed that the expression of genes encoding enzymes in GA, ABA, ethylene (ET), and brassinosteroid (BR) biosynthesis was significantly increased. The diterpenoid biosynthesis, carotenoid biosynthesis, cysteine and methionine metabolism, and brassinosteroid biosynthesis pathways occur in response to GA, ABA, ET, and BR biosynthesis. After EPS treatment for 5 days, the genes encoding GA2ox (K04125) in GA biosynthesis, CYP707A (K09843) in ABA biosynthesis, ERF1 (K14516) in ET biosynthesis, and CYP85A1 (K09590) and CYP92A6 (K20623) in BR biosynthesis were significantly differentially expressed compared to the control. GA2ox, CYP707A, and CYP85A1 were significantly up-regulated, while, CYP92A6 and ERF1 were down-regulated ( Supplementary Table 6 ).

 DEGs involved in transcription factors expression

A total of 857 unigenes were annotated to TFs, belonging to 34 families ( Supplementary Table 7 ). The top TF families included the MYB_superfamily (143 unigenes), AP2/ERF (93 unigenes), C2C2 (68 unigenes), bHLH (59 unigenes), WRKY (59 unigenes), NAC (56 unigenes), and GRAS (44 unigenes). Following EPS treatment for 5 days, 122 DEGs were significantly up-regulated and 106 down-regulated. After EPS treatment for 10 days, 24 DEGs were significantly up-regulated and 30 down-regulated. After EPS treatment for 15 days, 60 DEGs were significantly up-regulated and 55 down-regulated. Hence, the most TFs were significantly up-regulated after 5 days of EPS treatment; those related to flavonoid synthesis were primarily concentrated in the MYB superfamily (Yuan et al., 2015; Qian et al., 2020; Fang et al., 2023); however, they only included MYB3, MYB8, and MYB13.

 Discussion

The stimulatory effects of EPS and EP on the growth of S. baicalensis calli were observed. The fresh and dry weights of S. baicalensis calli were significantly increased following EPS treatment, consistent with previous findings (Li et al., 2013; Zhong et al., 2016). Moreover, the fresh and dry weights of S. baicalensis calli were promoted within the early and middle stages of treatment, i.e., 5–10 days, consistent with the results of Zhong et al. (2016). We observed similar results for EP treatment, whereas Xu et al. (2009) reported that the EP from Phytophthora boehmeriae culture filtrate had no effect on host growth. We speculate that the EP from endophytic fungi has different physiological effects due to differences in structures.

In our study, EPS and OP elicited stimulatory effects on S. baicalensis calli flavonoid production. Different concentrations promoted baicalin, baicalein, and wogonin biosynthesis across all time points. This strong effect elicited by EPS is consistent with previous findings (Li et al., 2013; Zhong et al., 2016; Chen et al., 2022). The maximum stimulatory effects on flavonoid production in S. baicalensis calli were observed after 5 days of EPS treatment. Zhong et al. (2016) showed that the highest rutin and quercetin contents in tartary buckwheat sprout cultures were obtained 9 days following EPS treatment, whereas the rutin and quercetin contents showed a decreasing trend at 10 days. We speculate that EPS elicits chemical defense responses in S. baicalensis calli and rapidly stimulates flavonoid production following exposure. With prolonged exposure, flavonoids may break down into compounds that are harmless to the host. In addition, the high concentration of OP showed mild promotional effects on baicalein, wogonoside, and wogonin production during the late stage of culture. However, the main components of OP remain unknown. Data on the effects of other elicitors (beyond exopolysaccharides and exoproteins) from the culture filtrates of endophytic fungi on the synthesis of plant secondary metabolites are lacking.

Interestingly, our results indicated that EPS exhibited a slight growth-promoting effect and, promoted flavonoid synthesis in S. baicalensis calli. Similar findings have been observed in previous studies. The polysaccharide derived from Pestalotiopsis sp. DO14 not only enhanced the accumulation of flavonoids in Dendrobium ferrugineum seedlings, but also showed a slight improvement in the growth (1.19–1.35-fold, fresh weight) (Zhu et al., 2018). Ming et al. (2013) found that polysaccharides from T. atroviride D16 not only promoted the growth of S. miltiorrhiza hairy roots, but also markedly promoted the synthesis of tanshinones (5–66-fold).

Transcriptome analysis was adopted to investigate the mechanisms underlying the cell growth and flavonoid production in S. baicalensis calli treated with 80 mg/L EPS and cultured for 0, 5, 10, and 15 days. The number of DEGs initially increases and then decreases over time. We speculate that EPS rapidly stimulates the expression of DEGs in the early stage of treatment. With prolonged exposure, the number of DEGs decreases due to attenuation of the stimulation. CHS, CHI, FNS, and F6H were identified as being involved in the early stages of EPS induction. qRT-PCR analysis confirmed that after culturing for 5 days, the transcription levels of CHS, CHI, FNS, and F6H were higher in EPS-treated S. baicalensis calli than in the control. This suggests that the promotion of flavonoid biosynthesis by EPS occurred primarily during the initial stages of induction. Moreover, the increased accumulation of baicalin, wogonoside, baicalein, and wogonin in EPS-treated S. baicalensis calli correlated with the increased expression of the above genes. Although no biotic elicitors of S. baicalensi have been reported, abiotic elicitors promote flavonoid accumulation in S. baicalensis by up-regulating key genes in the flavonoid biosynthesis pathway. Park et al. (2011) reported that following treatment of S. baicalensis cell suspensions with methyl jasmonate, the expression of CHI increased. Meanwhile, Xu et al. (2010) found that the expressions of PAL1, PAL2, PAL3, C4H, 4CL, and CHS in S. baicalensis cell suspensions were increased following methyl jasmonate induction. Hence, we speculate that different types of elicitors stimulate key gene expression in the flavonoid biosynthesis pathway in S. baicalensis.

ABA and GA can promote flavonoid synthesis in S. baicalensis (Yuan et al., 2013; He et al., 2023). Our RNA-seq results revealed significant up-regulation of GA2ox and CYP707A, which are involved in GA and ABA synthesis, respectively. Therefore, we speculate that EPS might promote the synthesis of GA and ABA by promoting the expression of GA2ox and CYP707A, ultimately promoting flavonoid synthesis. Although we also found that ET and BR were differentially expressed in this culture, it remains unclear whether they participate in the flavonoid synthesis process in S. baicalensis.

TFs play important roles in the growth and development of plants as they adapt to changes in the external environment (Strader et al., 2022). TFs with regulatory roles in the flavonoid synthesis pathway are concentrated in the NAC and MYB families (Yuan et al., 2013, 2015; Qian et al., 2020; Fang et al., 2023, 2023; He et al., 2023). After 5 days of EPS treatment, MYB3, MYB8, MYB13, FNS, CHS, and CHI expression was markedly up-regulated compared with the control. Hence, we posit that EPS might promote the expression of enzyme-encoding genes (CHS, CHI, and FNS) by up-regulating the expression of TFs (MYB3, MYB8, and MYB13), effectively promoting flavonoid synthesis ( Figure 5 ). Interestingly, studies have shown that phytohormones can also promote the expression of TFs. For example, Yuan et al. (2013) showed that GA metabolism contributes to flavonoid biosynthesis, with the expression of MYB8 significantly improving in S. baicalensis within 3 h of GA treatment. Furthermore, ABA can improve the expression of MYB3 and MYB8 in a relatively short period of time (Bang et al., 2008; Zhou et al., 2022). Accordingly, EPS not only stimulates the expression of ABA and GA key enzyme genes (GA2ox and CYP707A) in S. baicalensis calli but also promotes the expression of MYB3, MYB8, and MYB13. Based on our experimental results and relevant literature (Bang et al., 2008; Yuan et al., 2013, 2015; Qian et al., 2020; Zhou et al., 2022; Fang et al., 2023; He et al., 2023), EPS treatment induces flavonoid synthesis in S. baicalensis calli by regulating the expression of phytohormones, TFs, and enzymes. That is, EPS might stimulate the expression of genes encoding MYB3, MYB8, MYB13, GA2ox, and CYP707A, leading to significantly higher expression of CHS, CHI, FNS, and F6H and ultimately promoting the synthesis of flavonoids.

Schematic diagram of the EPS regulation of the flavonoid biosynthesis pathway in S. baicalensis calli. The solid arrow indicates a promoting effect. The dotted arrow indicates the biosynthetic pathways.

This study has certain limitations. First, although our transcriptomic analysis provided a large amount of data, we focused solely on genes related to flavonoid biosynthesis pathways. Thus, further analysis is required to identify other genes associated with the effect elicited by EPS on S. baicalensis calli. Second, we observed that EPS also promoted S. baicalensis calli growth. However, EPS is a crude polysaccharide and the active components are still unknown. It is, therefore, necessary to research the mechanisms promoting growth and the structural characterization of EPS in future studies.

 Conclusion

In this study, EPS from the culture filtrate of the endophytic fungus F. solani CL105 constituted the main active elicitor involved in the promotion of growth and flavonoid production in S. baicalensis calli. Transcriptome and qRT-PCR analyses further revealed that CHS, CHI, FNS, and F6H participate in regulating the flavonoid biosynthesis pathways following EPS treatment for 5 days. Moreover, the genes encoding GA2ox and CYP707A were significantly up-regulated and associated with GA and ABA biosynthesis, respectively. Meanwhile, the expression levels of TFs, including MYB3, MYB8, and MYB13, were significantly up-regulated. We speculate that EPS stimulates the expression of genes encoding MYB3, MYB8, MYB13, GA2ox, and CYP707A, which leads to significantly increased CHS, CHI, FNS, and F6H expression levels and promotion of baicalin, wogonoside, baicalein, and wogonin accumulation in S. baicalensis calli. Our findings revealed the mechanism by which F. solani CL105 promotes flavonoid accumulation in S. baicalensis calli. This study contributes to the growing body of evidence supporting the role of endophytic fungi in regulating the accumulation of active compounds in S. baicalensis and provides insights for the large-scale synthesis of flavonoids through biotechnological methods. However, the mechanisms promoting growth and the structural characterization of EPS require further investigation.

 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 below: BioProject, PRJNA999053.

 Author contributions

XC: Conceptualization, Data curation, Investigation, Methodology, Writing – original draft, Writing – review & editing. XZ: Conceptualization, Data curation, Investigation, Methodology, Writing – original draft, Writing – review & editing. HS: Investigation, Methodology, Writing – review & editing. YZ: Funding acquisition, Project administration, Supervision, Writing – review & editing. CS: Conceptualization, Funding acquisition, Project administration, Investigation, Supervision, Writing – review & editing.

 Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by the Natural Science Foundation of Hebei Province (No. H2020423254), Program of Hebei Administration of Traditional Chinese Medicine (No. 2019085 and No. 2024100), Fourth National Traditional Chinese Medicine Resources Survey of China (Z135080000022), Innovation Team of Hebei Province Modern Agricultural Industry Technology System (No. HBCT2018060205), and Basic Research Project for Outstanding Youth Teachers of Hebei University of Traditional Chinese Medicine (YQ2021008).

 Acknowledgments

We would like to thank all the members in Traditional Chinese Medicine Processing Technology Innovation Center of Hebei Province for fruitful discussions.

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

 References Bai C. Yang J. Cao B. Xue Y. Gao P. Liang H. . (2020). Growth years and post-harvest processing methods have critical roles on the contents of medicinal active ingredients of Scutellaria baicalensis . Ind. Crops Prod 158, 112985. doi: 10.1016/j.indcrop.2020.112985 Bang W. Y. Kim S. W. Jeong I. S. Koiwa H. Bahk J. D. (2008). The C-terminal region (640–967) of Arabidopsis CPL1 interacts with the abiotic stress-and ABA-responsive transcription factors. Biochem. Biophys. Res. Commun. 372, 907–912. doi: 10.1016/j.bbrc.2008.05.161 Cao X. Xu L. Wang J. Dong M. Xu C. Kai G. . (2022). Endophytic fungus Pseudodidymocyrtis lobariellae KL27 promotes taxol biosynthesis and accumulation in Taxus chinensis . BMC Plant Biol. 22, 1–18. doi: 10.1186/s12870–021-03396–6 Chen H. Chen J. Qi Y. Chu S. Ma Y. Xu L. . (2022). Endophytic fungus Cladosporium tenuissimum DF11, an efficient inducer of tanshinone biosynthesis in Salvia miltiorrhiza roots. Phytochemistry 194, 113021. doi: 10.1016/j.phytochem.2021.113021 Fang S. Qiu S. Chen K. Lv Z. Chen W. (2023). The transcription factors SbMYB45 and SbMYB86. 1 regulate flavone biosynthesis in Scutellaria baicalensis . Plant Physiol. Biochem. 200, 107794. doi: 10.1016/j.plaphy.2023.107794 Fang Y. Liu J. Zheng M. Zhu S. Pei T. Cui M. . (2023). SbMYB3 transcription factor promotes root-specific flavone biosynthesis in Scutellaria baicalensis . Hortic. Res. 10, uhac266. doi: 10.1093/hr/uhac266 Guo X. Wang Y. Zhai Z. Huang T. Zhao D. Peng X. . (2018). Transcriptomic analysis of light-dependent anthocyanin accumulation in bicolored cherry fruits. Plant Physiol. Biochem. 130, 663–677. doi: 10.1016/j.plaphy.2018.08.016 Hardoim P. R. Van Overbeek L. S. Berg G. Pirttilä A. M. Compant S. Campisano A. . (2015). The hidden world within plants: ecological and evolutionary considerations for defining functioning of microbial endophytes. Microbiol. Mol. Biol. Rev. 79, 293–320. doi: 10.1128/mmbr.00050–14 He H. Li Q. Fang L. Yang W. Xu F. Yan Y. . (2023). Comprehensive analysis of NAC transcription factors in Scutellaria baicalensis and their response to exogenous ABA and GA3. Int. J. Biol. Macromol. 244, 125290. doi: 10.1016/j.ijbiomac.2023.125290 Hu S. Wang D. Wang W. Zhang C. Li Y. Wang Y. . (2022). Whole genome and transcriptome reveal flavone accumulation in Scutellaria baicalensis roots. Front. Plant Sci. 13. doi: 10.3389/fpls.2022.1000469 Jia M. Chen L. Xin H. L. Zheng C. J. Rahman K. Han T. . (2016). A friendly relationship between endophytic fungi and medicinal plants: A systematic review. Front. Microbiol. 7. doi: 10.3389/fmicb.2016.00906 Li B. Dewey C. N. (2011). RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinform. 12, 323. doi: 10.1186/1471–2105-12–323 Li P. Luo C. Luo R. Mou Y. Sun W. Zhou L. (2013). Effects of polysaccharides and oligosaccharides from endophytic fungus Berkleasmium sp. dzf12 on diosgenin accumulation in Dioscorea zingiberensis cell and seedling cultures. Afr. J. Microbiol. Res. 7, 3049–3055. doi: 10.5897/AJMR12.599 Miguel V. Rey-Serra C. Tituaña J. Sirera B. Alcalde-Estévez E. Herrero J. I. . (2023). Enhanced fatty acid oxidation through metformin and baicalin as therapy for COVID-19 and associated inflammatory states in lung and kidney. Redox Biol. 68, 102957. doi: 10.1016/j.redox.2023.102957 Ming Q. Su C. Zheng C. Jia M. Zhang Q. Zhang H. . (2013). Elicitors from the endophytic fungus Trichoderma atroviride promote Salvia miltiorrhiza hairy root growth and tanshinone biosynthesis. J. Exp. Bot. 64, 5687–5694. doi: 10.1093/jxb/ert342 Murashige T. Skoog F. (1962). A Revised medium for rapid growth and bio assays with tobacco tissue cultures. Physiol. Plant 15, 473–497. doi: 10.1111/j.1399–3054.1962.tb08052.x Park N. I. Xu H. Li X. Kim S. J. Park S. U. (2011). Enhancement of flavone levels through overexpression of chalcone isomerase in hairy root cultures of Scutellaria baicalensis . Funct. Integr. Genomics 11, 491–496. doi: 10.1007/s10142–011-0229–0 Pei T. Yan M. Huang Y. Wei Y. Martin C. Zhao Q. (2022). Specific flavonoids and their biosynthetic pathway in scutellaria baicalensis . Front. Plant Sci. 13. doi: 10.3389/fpls.2022.866282 Qian C. Chen Z. Liu Q. Mao W. Chen Y. Tian W. . (2020). Coordinated transcriptional regulation by the UV-B photoreceptor and multiple transcription factors for plant UV-B responses. Mol. Plant 13, 777–792. doi: 10.1016/j.molp.2020.02.015 Strader L. Weijers D. Wagner D. (2022). Plant transcription factors—Being in the right place with the right company. Curr. Opin. Plant Biol. 65, 102136. doi: 10.1016/j.pbi.2021.102136 Wan G. X. Ma L. Zhang J. (2012). Effect on diferent doncentrations of exogenous hormones on baicalin in Scutellaria baicalensis callus. Zhongguo Zhong Yao Za Zhi 37, 3799–3802. doi: 10.4268/cjcmm20122424 Wang Z. L. Wang S. Kuang Y. Hu Z. M. Qiao X. Ye M. (2018). A comprehensive review on phytochemistry, pharmacology, and flavonoid biosynthesis of Scutellaria baicalensis . Pharm. Biol. 56, 465–484. doi: 10.1080/13880209.2018.1492620 Wu J. Ming Q. Zhai X. Wang S. Zhu B. Zhang Q. . (2019a). Structure of a polysaccharide from Trichoderma atroviride and its promotion on tanshinones production in Salvia miltiorrhiza hairy roots. Carbohydr. Polym. 223, 115125. doi: 10.1016/j.carbpol.2019.115125 Wu J. Peng L. Dong S. Xia X. Zhao L. (2019b). Transcriptome analysis of Chelidonium majus elaiosomes and seeds provide insights into fatty acid biosynthesis. PeerJ 7, e6871. doi: 10.7717/peerj.6871 Xu M. Dong J. Wang H. Huang L. (2009). Complementary action of jasmonic acid on salicylic acid in mediating fungal elicitor-induced flavonol glycoside accumulation of Ginkgo biloba cells. Plant Cell Environ. 32, 960–967. doi: 10.1111/j.1365–3040.2009.01976.x Xu D. Li N. Gu Y. Q. Huang J. Hu B. S. Zheng J. Y. . (2023). Endophytic fungus Colletotrichum sp. AP12 promotes growth physiology and andrographolide biosynthesis in Andrographis paniculata (Burm. f.) Nees. Front. Plant Sci. 14. doi: 10.3389/fpls.2023.1166803 Xu H. Park N. I. Li X. Kim Y. K. Lee S. Y. Park S. U. (2010). Molecular cloning and characterization of phenylalanine ammonia-lyase, cinnamate 4-hydroxylase and genes involved in flavone biosynthesis in Scutellaria baicalensis . Bioresour. Technol. 101, 9715–9722. doi: 10.1016/j.biortech.2010.07.083 Xu J. Shi R. Cheng Y. Xie G. Zhu Y. Qin M. (2022). Metabolomic association analysis reveals defense related secondary metabolic reprogramming in callus of Scutellaria baicalensis Georgi induced by endophytic Pseudomonas sp. 2B. Ind. Crops Prod 187, 115409. doi: 10.1016/j.indcrop.2022.115409 Yu X. Bi Y. Yan L. Liu X. Wang Y. Shen K. . (2016). Activation of phenylpropanoid pathway and PR of potato tuber against Fusarium sulphureum by fungal elicitor from Trichothecium roseum . World J. Microbiol. Biotechnol. 32, 142. doi: 10.1007/s11274–016-2108–2 Yu J. S. Zhang M. M. Shi J. Yang Y. Meng X. Xue J. P. . (2021). Research progress on mechanism of phytohormones in regulating flavonoid metabolism. Zhongguo Zhong Yao Za Zhi 46, 3806–3813. doi: 10.19540/j.cnki.cjcmm.20210522.103 Yuan Y. Qi L. Yang J. Wu C. Liu Y. Huang L. (2015). A Scutellaria baicalensis R2R3-MYB gene, SbMYB8, regulates flavonoid biosynthesis and improves drought stress tolerance in transgenic tobacco. Plant Cell Tissue Organ Cult. 120, 961–972. doi: 10.1007/s11240–014-0650-x Yuan Y. Wu C. Liu Y. Yang J. Huang L. (2013). The Scutellaria baicalensis R2R3-MYB transcription factors modulates flavonoid biosynthesis by regulating GA metabolism in transgenic tobacco plants. PloS One 8, e77275. doi: 10.1371/journal.pone.0077275 Yuan Q. J. Zhang Z. Y. Hu J. Guo L. P. Shao A. J. Huang L. Q. (2010). Impacts of recent cultivation on genetic diversity pattern of a medicinal plant, Scutellaria baicalensis (Lamiaceae). BMC Genom. 11, 1–13. doi: 10.1186/1471–2156-11–29 Zhang X. Cui X. X. Zheng Y. G. Su C. Y. (2023). Screening and molecular identification of endophytic fungi promoting the accumulation of flavonoids in callus of Scutellaria baicalensis . Zhongguo Zhong Yao Za Zhi 48, 4974–4980. doi: 10.19540/j.cnki.cjcmm.20230606.101 Zhang S. Guo S. L. Wang Q. B. Liu Y. Shen H. W. Wang Z. Y. (2019). Effects of fungi fraction on growth and anti-oxidative activity of Eleutherococcus senticosus . Zhongguo Zhong Yao Za Zhi 44, 1517–1523. doi: 10.19540/j.cnki.cjcmm.20190129.008 Zhao Q. Cui M. Y. Levsh O. Yang D. Liu J. Li J. . (2018). Two CYP82D enzymes function as flavone hydroxylases in the biosynthesis of root-specific 4′-deoxyflavones in Scutellaria baicalensis . Mol. Plant 11, 135–148. doi: 10.1016/j.molp.2017.08.009 Zhao T. Tang H. Xie L. Zheng Y. Ma Z. Sun Q. . (2019a). Scutellaria baicalensis Georgi.(Lamiaceae): a review of its traditional uses, botany, phytochemistry, pharmacology and toxicology. J. Pharm. Pharmacol. 71, 1353–1369. doi: 10.1111/jphp.13129 Zhao Q. Yang J. Cui M. Y. Liu J. Fang Y. Yan M. . (2019b). The reference genome sequence of Scutellaria baicalensis provides insights into the evolution of wogonin biosynthesis. Mol. Plant 12, 935–950. doi: 10.1016/j.molp.2019.04.002 Zhao Q. Zhang Y. Wang G. Hill L. Weng J. K. Chen X. Y. . (2016). A specialized flavone biosynthetic pathway has evolved in the medicinal plant, Scutellaria baicalensis . Sci. Adv. 2, e1501780. doi: 10.1126/sciadv.1501780 Zhong L. Niu B. Tang L. Chen F. Zhao G. Zhao J. (2016). Effects of Polysaccharide Elicitors from Endophytic Fusarium oxysporum Fat9 on the Growth, Flavonoid Accumulation and Antioxidant Property of Fagopyrum tataricum Sprout Cultures. Molecules 21, 1590. doi: 10.3390/molecules21121590 Zhou L. Li J. Zeng T. Xu Z. Luo J. Zheng R. . (2022). TcMYB8, a R3-MYB transcription factor, positively regulates pyrethrin biosynthesis in tanacetum cinerariifolium . Int. J. Mol. Sci. 23, 12186. doi: 10.3390/ijms232012186 Zhou X. Liu J. Zhuang Y. (2014). Selection of appropriate reference genes in eggplant for quantitative gene expression studies under different experimental conditions. Sci. Hortic. 176, 200–207. doi: 10.1016/j.scienta.2014.07.010 Zhu B. Wu L. Wan H. Yang K. Si J. Qin L. (2018). Fungal elicitors stimulate biomass and active ingredients accumulation in Dendrobium catenatum plantlets. Biologia 73, 917–926. doi: 10.2478/s11756–018-0091–9