Genome-Wide Characterization of SPL Gene Family in Codonopsis pilosula Reveals the Functions of CpSPL2 and CpSPL10 in Promoting the Accumulation of Secondary Metabolites and Growth of C. pilosula Hairy Root

SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) transcription factors play critical roles in regulating diverse aspects of plant growth and development, including vegetative phase change, plant architecture, anthocyanin accumulation, lateral root growth, etc. In the present study, 15 SPL genes were identified based on the genome data of Codonopsis pilosula, a well-known medicinal plant. Phylogenetic analysis clustered CpSPLs into eight groups (G1-G8) along with SPLs from Arabidopsis thaliana, Solanum lycopersicum, Oryza sativa and Physcomitrella patens. CpSPLs in the same group share similar gene structure and conserved motif composition. Cis-acting elements responding to light, stress and phytohormone widely exist in their promoter regions. Our qRT-PCR results indicated that 15 CpSPLs were differentially expressed in different tissues (root, stem, leaf, flower and calyx), different developmental periods (1, 2 and 3 months after germination) and various conditions (NaCl, MeJA and ABA treatment). Compared with the control, overexpression of CpSPL2 or CpSPL10 significantly promoted not only the growth of hairy roots, but also the accumulation of total saponins and lobetyolin. Our results established a foundation for further investigation of CpSPLs and provided novel insights into their biological functions. As far as we know, this is the first experimental research on gene function in C. pilosula.


Introduction
Transcription factors (TFs) function in various physiological and developmental processes via activating and/or repressing transcription of multiple target genes [1]. They have been usually divided into different families according to the sequence of DNA-binding domains and other conserved motifs [2]. SQUAMOSA-promoter binding protein-like (SPL or SBP) TFs are exclusive to plant and characterized by a highly conserved SBP domain and a nuclear localization signal (NLS) at the C-terminus. The SBP domain is approximately 76 amino acids and includes two zinc-binding sites (one zinc finger is C3H or C4, and the other is C2H4) essential for DNA binding, and the NLS partially overlaps with the second zinc finger [3][4][5][6]. AmSBP1 and AmSBP2 from Antirrhinum majus were the first discovered SBP-domain proteins in plants and were found to bind to the floral meristem identity gene SQUAMOSA promoter, so named them [3]. Since then, SPL genes

Plant Materials and Treatments
Seeds of C. pilosula were collected from Gansu Province, China. The botanical origin of the materials was identified by Professor ZheZhi Wang in Shaanxi Normal University. The specimens of the seeds were deposited in the herbarium of National Engineering Laboratory for Resource Development of Endangered Crude Drugs in Northwest of China, Shaanxi Normal University, Xi'an, China. The seeds of C. pilosula were germinated and incubated according to the method that we described previously [39].
For gene spatiotemporal expression analysis, the seeds were germinated and incubated in mixed soil of nutrient soil, perlite and vermiculite in a glass greenhouse with a temperature regime of 24 ± 2 • C, 40-50% relative humidity. The leaves, stems and roots were collected separately from one-month-old seedlings, two-month-old seedlings, and three-month-old plantlets, and the flower and calyx were collected from the plants at the flowering stage. To test CpSPLs responses to hormonal and stress treatments, two-week-old seedings were treated with 200 mmol NaCl, 200 µmol MeJA, and 100 µmol ABA, respectively, and samples were gathered after 6 h as we have described previously [40]. The control group was treated with the same amount of ddH 2 O and all the samples were collected 6 h after treatment.

Gene Expression Analysis
For qRT-PCR analysis, total RNA was extracted and then reverse transcribed into cDNA as we described previously [40]. All the primer sequences used for qRT-PCR are listed in Supplementary Material Table S1 and CpGAPDH was used as the internal control [40]. The relative expression levels of 15 CpSPLs were calculated according to the method described by Livak and Schmittgen [41]. All the experiments included three biological and three technical replicates.

Vector Construction and Hairy Root Transformation
The complete open reading frames of CpSPL2 and CpSPL10 were amplified through PCR using the specific primer pairs CpSPL2-F/R and CpSPL10-F/R (Supplementary  Material Table S1), respectively, with the following PCR conditions: 95 • C, 3 min; 30 cycles of 95 • C, 10 s, 58 • C, 30 s, 72 • C, 45 s; 72 • C, 5 min. Then the products were digested with Pac I and Asc I and ligated into pMDC85 to generate overexpression (OE) vectors pMDC85-CpSPL2 and pMDC85-CpSPL10.
Transgenic hairy roots overexpressing CpSPL2 or CpSPL10 were obtained by Agrobacterium-mediated method according to the protocol established in our lab [39]. Briefly, the roots of four-week-old C. pilosula aseptic seedlings were cut into an average length of 0.5 cm. The surface of the explants was gently scratched with a scalpel back, inoculated on MS medium, and pre-cultured at 24 ± 2 • C for 2 days in dark. The pre-cultured explants were immersed in the infection solution for 5 min and inoculated onto 1/2 MS solid medium for co-cultivation for 2 days. After co-cultivation for 2 days, the explants were transferred to 1/2 MS solid medium supplemented with 200 mg/L cefotaxime sodium (Cef) and 2 mg/L hygromycin. In parallel, pMDC85 was introduced into C. pilosula as the empty vector control (EV). Every transgenic line was excised and sub-cultured separately as we described previously [39]. Five independent CpSPL2-OE lines, seven CpSPL10-OE lines and four EV lines were obtained, and then confirmed by genomic DNA PCR using primers hptII-F/R (Supplementary Material Table S1) for hygromycin phosphotransferase II gene (hptII), followed by expression analysis of CpSPL2 or CpSPL10 by qRT-PCR.

Determination of Lobetyolin and Total Saponins
Transgenic hairy roots sub-cultured for one month were used for determination of lobetyolin and total saponins.
To determine the concentration of lobetyolin, we ground 50 mg dried hairy roots into powder, followed by extracted three times with 10 mL methanol in an ultrasonic bath (Kunshan Instrument Co., Ltd., Kunshan, China) for 30, 20 and 15 min, respectively. The extracts were put together and the methanol solution was evaporated, followed by dissolving with methanol to 5 mL volumetric flask. After filtration with 0.22 µm microporous membrane, the solution was used for HPLC analysis on a Shimadzu LC-20A instrument (Shimadzu, Kyoto, Japan) equipped with an Agilent 5 TC-C 18 column (250 × 4.6 mm, 5 µm). The mobile phase consisted of ultrapure water (A) and methanol (B) and the gradient condition was 0-5 min, 20-40% B; 5-10 min, 40-70% B; 10-12 min, 79-90% B; 12-25 min, 90% B. The separation was performed at 30 • C, with the flow rate of 1.0 mL/min and UV detector wavelength at 220 nm. The concentration of total saponins in transgenic hairy roots was determined as we described previously [39].

Statistical Analysis
All the experiments and data presented here involved at least three biological repeats. SPSS version 20.0 software (SPSS Inc., Chicago, IL, USA) was used for statistical evaluation. The error bars indicate standard deviation. Significant difference of the mean values was set at p < 0.05.

Genome-Wide Identification and Sequence Feature Analysis of CpSPLs
To identify possible SPL genes in C. pilosula genome sequences, we employed the SBP domain (PF03110) to search the databases by HMMER. A total of 15 SPLs containing complete SBP domain were identified based on the genome sequence of C. pilosula and their cDNA sequences were listed in Supplementary Material Table S2. Consulting the homologous AtSPLs in Arabidopsis, 15 CpSPLs were named from CpSPL1 to CpSPL15. The deduced CpSPLs exhibited great variations in terms of their molecular weight (MW), ranging from 17.98 (CpSPL5) to 119.80 KDa (CpSPL14). Similarly, the lengths of the CDS were found to be varied in the CpSPLs, from 480 (CpSPL5) to 3276 bp (CpSPL14). The detailed information, including the gene length, intron number, protein length, predicted MW and theoretical isoelectric point (pI), are listed in Table 1.

Gene Structure and Conserved Motif Analysis
To clarify the structural diversities of 15 CpSPLs, we performed gene exon/intron structure analysis. The result displayed that the number of introns had a high variation and ranged from 1 to 10 ( Figure 2). Interestingly, we found that most CpSPLs in the same group share a similar structure. For instance, CpSPL1 and CpSPL14, belonging to G8, have 10 introns, respectively ( Figure 2).

Gene Structure and Conserved Motif Analysis
To clarify the structural diversities of 15 CpSPLs, we performed gene exon/intron structure analysis. The result displayed that the number of introns had a high variation and ranged from 1 to 10 ( Figure 2). Interestingly, we found that most CpSPLs in the same group share a similar structure. For instance, CpSPL1 and CpSPL14, belonging to G8, have To explore the conserved motifs, 15 CpSPLs were subjected to analysis with MEME program. Among the 12 conserved motifs identified ( Figure 3 and Supplementary Material Table S3), motif 1, motif 2 and motif 3 existed in all the 15 CpSPLs and formed the conserved SBP domain. Similar motif composition existed in the same group. For example, CpSPL2 and CpSPL10 in G1 all consisted of five conserved motifs (motif 1/2/3/10/11). The motif composition in CpSPL9 was completely consistent with that in CpSPL15, suggesting that CpSPL9 and CpSPL15 probably have similar and redundant functions in plant development. To explore the conserved motifs, 15 CpSPLs were subjected to analysis with MEME program. Among the 12 conserved motifs identified ( Figure 3 and Supplementary Material Table S3), motif 1, motif 2 and motif 3 existed in all the 15 CpSPLs and formed the conserved SBP domain. Similar motif composition existed in the same group. For example, CpSPL2 and CpSPL10 in G1 all consisted of five conserved motifs (motif 1/2/3/10/11). The motif composition in CpSPL9 was completely consistent with that in CpSPL15, suggesting that CpSPL9 and CpSPL15 probably have similar and redundant functions in plant development. To explore the conserved motifs, 15 CpSPLs were subjected to analysis with MEME program. Among the 12 conserved motifs identified ( Figure 3 and Supplementary Material Table S3), motif 1, motif 2 and motif 3 existed in all the 15 CpSPLs and formed the conserved SBP domain. Similar motif composition existed in the same group. For example, CpSPL2 and CpSPL10 in G1 all consisted of five conserved motifs (motif 1/2/3/10/11). The motif composition in CpSPL9 was completely consistent with that in CpSPL15, suggesting that CpSPL9 and CpSPL15 probably have similar and redundant functions in plant development.

Spatiotemporal Expression Analysis of CpSPL Genes
We investigated the expression patterns of 15 CpSPLs in the leaves, stems and roots from one-month-old seedlings, two-month-old seedlings, and three-month-old plantlets, and the flower and calyx from the plants at the flowering stage by qRT-PCR assay. The results showed that most CpSPLs expressed in almost all the tissues (Figure 4). Compared with other genes, the expression level of CpSPL7 was more constant in all the tissues tested. CpSPL8 showed highest level in calyx. CpSPL5 was expressed at relatively higher levels in leaf and calyx. The expression levels of CpSPL3, CpSPL8, CpSPL10, CpSPL12 and CpSPL13 in the stems gradually decreased with the maturation of the seedlings. CpSPL1 and CpSPL14, two members in G1, showed similar expression patterns and their expression levels in the root increased gradually with the maturation of the seedlings. In addition, the expression patterns of CpSPL9 and CpSPL15 were highly similar, with higher levels in flowers and 3-month-old roots. In summary, spatiotemporal expression analysis results indicated that CpSPL genes exhibited various expression patterns, which provide preliminary information for understanding their potential functions in the development of C. pilosula.

Overexpression of CpSPL2 or CpSPL10 Promotes the Growth of C. pilosula Hairy Root
To investigate the function of CpSPL2 and CpSPL10 in root development, we generated CpSPL2-overexpressing or CpSPL10-overexpressing transgenic hairy roots. The expression level of CpSPL2 or CpSPL10 in the transgenics was examined by qRT-PCR ( Figure 6A,B). Two independent CpSPL2-overexpressing lines (CpSPL2-OE3 and CpSPL2-OE5) and CpSPL10overexpressing lines (CpSPL10-OE2 and CpSPL2-OE3) with dramatically elevated CpSPL2 or CpSPL10 expression were selected for further analysis. In comparison to the control, the hairy roots overexpressing CpSPL2 or CpSPL10 grew faster ( Figure 6C,D). When the transgenic hairy roots with the length about 1.0 cm were cultured for one month, the biomass of CpSPL2-OE3, CpSPL2-OE5, CpSPL10-OE2 and CpSPL10-OE3 was 2.19, 1.98, 3.15 and 2.83 times that of the control (EV2), respectively ( Figure 6E). Our results indicated that both CpSPL2 and CpSPL10 promote the growth of hairy roots.
week-old seedings were treated with 100 µmol abscisic acid (ABA). Change multiples are relative to the expression level of each gene treated with ddH2O (set to 1). One-way ANOVA (followed by Tukey's comparisons) tested for significant differences among means (* indicates the value is significantly different from the control at p < 0.05).

Overexpression of CpSPL2 or CpSPL10 Promotes the Growth of C. pilosula Hairy Root
To investigate the function of CpSPL2 and CpSPL10 in root development, we generated CpSPL2-overexpressing or CpSPL10-overexpressing transgenic hairy roots. The expression level of CpSPL2 or CpSPL10 in the transgenics was examined by qRT-PCR ( Figure  6A,B). Two independent CpSPL2-overexpressing lines (CpSPL2-OE3 and CpSPL2-OE5) and CpSPL10-overexpressing lines (CpSPL10-OE2 and CpSPL2-OE3) with dramatically elevated CpSPL2 or CpSPL10 expression were selected for further analysis. In comparison to the control, the hairy roots overexpressing CpSPL2 or CpSPL10 grew faster ( Figure  6C,D). When the transgenic hairy roots with the length about 1.0 cm were cultured for one month, the biomass of CpSPL2-OE3, CpSPL2-OE5, CpSPL10-OE2 and CpSPL10-OE3 was 2.19, 1.98, 3.15 and 2.83 times that of the control (EV2), respectively ( Figure 6E). Our results indicated that both CpSPL2 and CpSPL10 promote the growth of hairy roots. CpSPL2-OE lines when 1 cm hairy roots were cultured for one month. (D) Growth phenotype of EV and CpSPL10-OE lines when 1 cm hairy roots were cultured for one month. (E) The biomass of EV, CpSPL2-OE and CpSPL10-OE lines when 1 cm hairy roots were cultured for one month. One-way ANOVA (followed by Tukey's comparisons) tested for significant differences among means (indicated by different letters at p < 0.05). CpSPL10-OE lines when 1 cm hairy roots were cultured for one month. One-way ANOVA (followed by Tukey's comparisons) tested for significant differences among means (indicated by different letters at p < 0.05).

Overexpression of CpSPL2 or CpSPL10 Promotes Accumulation of Lobetyolin and Total Saponins in C. pilosula Hairy Root
To evaluate the impact of CpSPL2 or CpSPL10 on active ingredients, HPLC and UV spectrophotometer were used to determine the concentrations of lobetyolin and total saponins in those transgenic lines, respectively. It was surprising that the production of both lobetyolin and total saponins were greatly increased in CpSPL2-OE or CpSPL10-OE lines. The concentration of lobetyolin in CpSPL2-OE3, CpSPL2-OE5, CpSPL10-OE2 and CpSPL10-OE3 was 6.43, 6.25, 6.29 and 7.03 times that of the control (EV2), respectively ( Figure 7A). The concentration of total saponins in CpSPL2-OE3, CpSPL2-OE5, CpSPL10-OE2 and CpSPL10-OE3 was 3.18, 2.72, 1.81 and 1.94 times that of the control (EV2), respectively ( Figure 7B). In summary, CpSPL2 and CpSPL10 promote not only the growth of hairy roots but also accumulations of lobetyolin and total saponins.

Overexpression of CpSPL2 or CpSPL10 Promotes Accumulation of Lobetyolin and Total Saponins in C. pilosula Hairy Root
To evaluate the impact of CpSPL2 or CpSPL10 on active ingredients, HPLC and UV spectrophotometer were used to determine the concentrations of lobetyolin and total saponins in those transgenic lines, respectively. It was surprising that the production of both lobetyolin and total saponins were greatly increased in CpSPL2-OE or CpSPL10-OE lines. The concentration of lobetyolin in CpSPL2-OE3, CpSPL2-OE5, CpSPL10-OE2 and CpSPL10-OE3 was 6.43, 6.25, 6.29 and 7.03 times that of the control (EV2), respectively ( Figure 7A). The concentration of total saponins in CpSPL2-OE3, CpSPL2-OE5, CpSPL10-OE2 and CpSPL10-OE3 was 3.18, 2.72, 1.81 and 1.94 times that of the control (EV2), respectively ( Figure 7B). In summary, CpSPL2 and CpSPL10 promote not only the growth of hairy roots but also accumulations of lobetyolin and total saponins.

Identification of SPL Genes in C. pilosula
SPLs are plant-specific TFs and characterized by a highly conserved SBP domain [5,6]. They play critical roles in regulating diverse aspects of plant growth and development, including vegetative phase change, plant architecture, anthocyanin accumulation, lateral root growth, etc. [18][19][20][21][22][23][24]. Since its first discovery in A. majus [3], the SPL gene family from various plants has been isolated and identified. For instance, there are 16 SPL gene family members in Arabidopsis thaliana [8], 19 in Oryza sativa [9], 15 in Solanum lycopersicum [11], 15 in Salvia miltiorrhiza [13] and 15 in Ricinus communis [17]. However, little information is known about SPL gene family in C. pilosula, a famous species with important medical and edible values. Here, we identified 15 CpSPL genes in C. pilosula genome.
We constructed the phylogenetic tree of 15 from C. pilosula, 19 from O. sativa, 17 from S. lycopersicum, 14 from Physcomitrella patens and 16 from A. thaliana. In total, 81 SPL genes were divided into eight groups and each group had at least one CpSPL (Figure 1). CpSPL family members in the same group showed similar gene structure and motif composition (Figures 2 and 3), which was consistent with previous report [17]. In Arabidopsis, members in the same group often have the same or similar function. For instance, AtSPL3, AtSPL4 and AtSPL5, clustered in G7, synergistically induce flowering under long-day photoperiod [43]. Most recently, AtSPL2, AtSPL10 and AtSPL11, members in G1, have been reported to inhibit root regeneration by dampening auxin biosynthesis [42]. We speculate that CpSPLs in the same group maybe have the same function, such as CpSPL2 and CpSPL10 in G1, CpSPL5 in G7 and so on.

Identification of SPL Genes in C. pilosula
SPLs are plant-specific TFs and characterized by a highly conserved SBP domain [5,6]. They play critical roles in regulating diverse aspects of plant growth and development, including vegetative phase change, plant architecture, anthocyanin accumulation, lateral root growth, etc. [18][19][20][21][22][23][24]. Since its first discovery in A. majus [3], the SPL gene family from various plants has been isolated and identified. For instance, there are 16 SPL gene family members in Arabidopsis thaliana [8], 19 in Oryza sativa [9], 15 in Solanum lycopersicum [11], 15 in Salvia miltiorrhiza [13] and 15 in Ricinus communis [17]. However, little information is known about SPL gene family in C. pilosula, a famous species with important medical and edible values. Here, we identified 15 CpSPL genes in C. pilosula genome.
We constructed the phylogenetic tree of 15 from C. pilosula, 19 from O. sativa, 17 from S. lycopersicum, 14 from Physcomitrella patens and 16 from A. thaliana. In total, 81 SPL genes were divided into eight groups and each group had at least one CpSPL (Figure 1). CpSPL family members in the same group showed similar gene structure and motif composition (Figures 2 and 3), which was consistent with previous report [17]. In Arabidopsis, members in the same group often have the same or similar function. For instance, AtSPL3, AtSPL4 and AtSPL5, clustered in G7, synergistically induce flowering under long-day photoperiod [43]. Most recently, AtSPL2, AtSPL10 and AtSPL11, members in G1, have been reported to inhibit root regeneration by dampening auxin biosynthesis [42]. We speculate that CpSPLs in the same group maybe have the same function, such as CpSPL2 and CpSPL10 in G1, CpSPL5 in G7 and so on.

CpSPL Genes' Expression Patterns in C. pilosula
Gene expression patterns, to a large extent, will provide valuable information for its potential function [44]. In this study, the spatiotemporal expression patterns of 15 CpSPLs in the leaves, stems and roots from one-month-old seedlings, two-month-old seedlings, and three-month-old plantlets, and the flower and calyx from the plants at the flowering stage were detected by qRT-PCR (Figure 4). The results showed that CpSPL1 and CpSPL14 in G8 exhibited similar expression patterns, and the expression patterns of paralogous CpSPL9 and CpSPL15 in G6 showed high similarity. Our results were consistent with previous conclusion that paralogous SPL genes in the same group often showed similar expression profiles [45,46]. AtSPL9, AtSPL10 and AtSPL15 contribute to the vegetative to reproductive phase transition [8]. Here, CpSPL9, CpSPL10 and CpSPL15 expressed predominantly in the flower, suggesting they might function in the development of flower in C. pilosula.
Some SPL genes have been proved to be involved in abiotic stress. For example, in Arabidopsis, AtSPL1 and AtSPL12 function redundantly in thermotolerance and overexpression of AtSPL1 or AtSPL12 increased plant thermotolerance [47]. In alfalfa, silencing MsSPL13 enhanced tolerance to drought and heat stress (40 • C) [26,48], and downregulation of MsSPL8 led to enhanced salt and drought tolerance [49]. In the present study, we investigated the expression levels of 15 CpSPLs under various stress conditions, including NaCl, MeJA, or ABA treatment. We found that the expression levels of most CpSPL genes significantly changed under NaCl, MeJA and ABA treatment ( Figure 5). Among those genes with significant change, CpSPL4, CpSPL6, CpSPL14 and CpSPL15 positively responded to all the treatments, while CpSPL5, CpSPL8 and CpSPL10 negatively responded to all the treatments. Compared with other genes, CpSPL5 and CpSPL8 showed higher fold change under different treatments. We speculate that those two genes are potential candidates involved in abiotic stress.

Functional Study of the CpSPL2 and CpSPL10 Genes
Since the medicinal and edible part of C. pilosula is the root, increasing root yield is one of the main goals of breeding for this species. In Arabidopsis, AtSPL2, AtSPL10 and AtSPL11, inhibit root regeneration by dampening auxin biosynthesis [42]. The miR156targeted SPL10 is involved in regulating not only lateral root growth but also primary root growth [21,23]. Recently, it was reported that overexpression of FvSPL10, a SPL gene from Fragaria vesca, resulted in increased organs size, including longer root, larger floral organ and seeds [30]. We speculated that CpSPL2 and CpSPL10, two members clustered in the same group with AtSPL2/10/11 (Figure 1), were probably involved in the regulation of root development. To investigate the function of CpSPL2 and CpSPL10, we generated transgenic hairy roots overexpressing CpSPL2 or CpSPL10. Compared with the control, transgenic lines overexpressing CpSPL2 or CpSPL10 grew faster and the biomass of CpSPL2-OE3, CpSPL2-OE5, CpSPL10-OE2 and CpSPL10-OE3 was 2.19, 1.98, 3.15 and 2.83 times that of the control when the transgenic hairy roots with the length about 1.0 cm were cultured for one month ( Figure 6). Our results indicated that overexpression of CpSPL2 or CpSPL10 significantly promote the growth of hairy root.
Furthermore, we determined the concentration of lobetyolin and total saponins in those transgenic lines. Unexpectedly, we found that overexpressing CpSPL2 or CpSPL10 dramatically promoted the accumulation of lobetyolin and total saponins in the hairy roots ( Figure 7). Among 16 AtSPLs in Arabidopsis, AtSPL9 is the only one that has been reported to regulate biosynthesis of secondary metabolites [24,25]. AtSPL9 negatively regulates anthocyanin accumulation by preventing the formation of MBW complex [24], and it positively regulates the formation of (E)-β-caryophyllene by binding to the promoter of sesquiterpene synthase gene TPS21 and activates its expression [25]. Our results indicated that CpSPL2 and CpSPL10 are potential candidates for genetic improvement of C. pilosula because they can significantly promote not only the growth of hairy roots, but also accumulation of lobetyolin and total saponins. The molecular mechanism that CpSPL2 and CpSPL10 function in hairy roots needs to be addressed in the future.

Conclusions
In this study, we identified 15 CpSPL genes, which were supported by confirmation of the SBP domain, based on the genome data of C. pilosula. All CpSPLs were clustered into eight groups and members in the same group share similar gene structure and conserved motif composition. The spatiotemporal expression analysis of 15 CpSPLs showed that the CpSPL gene family had various expression patterns. The expression levels of most CpSPLs significantly changed under NaCl, MeJA, or ABA treatment, and CpSPL5 and CpSPL8 showed higher fold change under different treatments. Overexpression of CpSPL2 or CpSPL10 significantly promoted not only the growth of hairy roots, but also the accumulation of lobetyolin and total saponins. CpSPL2 and CpSPL10 are potential candidates for genetic improvement of C. pilosula. These results established a foundation for further investigation of CpSPLs and provided novel insights into their biological functions.
Author Contributions: J.Y. performed the experiments, analyzed the data and wrote the draft manuscript. Z.G. and W.W. prepared the materials and assisted in the experiments. X.Y. and X.C. conceived the study and revised the manuscript. All authors have read and agreed to the published version of the manuscript.