Comparison of S. chinensis regenerated plant rooting under dark and light conditions
Light is the primary energy source for plant growth and is a complex environmental factor. To explore the effects of light and dark on the rooting of S. chinensis regenerated plants, we first determined the growth status, rooting rate, and germination rate of somatic embryos. As shown in Table 2, the embryogenic calli of S. chinensis that were treated with light for 14 days underwent rapid expansion, and some of the somatic embryos exhibited green pigmentation. The somatic embryos that were exposed to light for 28 days displayed delicate, red hypocotyls. Conversely, after 14 days of dark treatment, most of the somatic embryos subjected to this darkness treatment expanded and elongated, while those treated with dark for 28 days exhibited thicker, well-grown, pale yellow hypocotyls. The rooting and germination rates were markedly higher under dark conditions than under light conditions (Fig. 1).
In plant tissue culture, light plays a crucial role in the proliferation and differentiation of culture materials, a development known as morphological plasticity, as a form of adaptation. Previous studies have shown that the natural development process of zygotic embryos occurs in a dark and relatively closed environment, suggesting that the induction and proliferation of somatic embryos may also require continuous dark or low light conditions (Wei 2017). Recent research demonstrated that short-term dark treatment significantly increased the germination rate of somatic embryos in Paeifloria lactiflora (Liu 2020). Similarly, dark treatment is also needed for embryogenic callus induction in the pedicel (Yu 2020). Zhou (2008) found that inducing embryogenic callus of peony in full darkness was more conducive to the induction of embryogenic cells and the formation of somatic embryos. In this study, somatic embryos of S. chinensis treated in the dark for 28 days showed a significantly higher rooting rate (95.56%) and germination rate (36.67%) than those treated with light, indicating that the dark environment is suitable for somatic embryo radicle development and indirectly promotes germination (Sun et al. 2013). These findings suggest that dark treatment can improve the rooting and germination rates of somatic embryos in various plant species.
Table 2
Effects of light and dark on the growth and development of somatic embryos in S. chinensis
Light condition
|
Rooting
(%)
|
Germination
(%)
|
Growth and development
|
14 d
|
28 d
|
Light
|
4.11b
|
15.56b
|
Some somatic embryos rapidly expanded and acquired a green color
|
The hypocotyl of germinating somatic embryos are predominantly slender and acquired a red color
|
Dark
|
95.56a
|
36.67a
|
Most somatic embryos undergo swelling and elongation
|
The hypocotyl of germinating somatic embryos appears pale yellow and exhibits increased thickness
|
Histological observation at different stages of somatic embryo development in S. chinensis
Histological analysis of somatic embryo development confirmed the period of occurrence and development of root primordia. The morphological characteristics of the S. chinensis somatic embryos cultured under dark or light conditions were observed at 7, 14, 21, and 28 days (Fig. 2). After the embryogenic calli (Fig. 2A) were inoculated in somatic embryo induction medium under dark conditions, radicle primordia developed during the globular embryo stage on day 7 (Fig. 2B), followed by development of the radicle during the heart-shaped embryo stage on day 14 (Fig. 2C). Subsequently, the radicle gradually elongated during the torpedo-shaped embryo stage on day 21 (Fig. 2D) and fully developed during the cotyledon embryo stage on day 28 (Fig. 2E). In contrast, under light treatment, the majority of the S. chinensis somatic embryos failed to develop normally, with a large proportion of the embryos failing to form radicle primordia at the globular embryo stage (Fig. 2F). A few radicle primordia emerged at the heart-shaped embryo stage on day 14 (Fig. 2G) but did not develop and elongate during the torpedo-shaped embryo stage on day 21 (Fig. 2H). Most of the embryos remained in the torpedo-shaped embryo stage and gradually turned brown, possibly due to radicle primordium defects. Only a few of the embryos developed into cotyledon embryos on day 28 and formed complete plants (Fig. 2I).
Quantification of Endogenous Hormones during Somatic Embryogenesis of S. chinensis
The endogenous hormone content in plants has an important influence on plant organ morphogenesis. Here, we measured the endogenous hormone content in somatic embryo development under light and dark conditions. As shown in Fig. 3, the concentration of indole-3-acetic acid (IAA) was significantly higher in the light treatment (597.565 ng•g− 1) than in the dark treatment (190.965 ng•g− 1), exhibiting a threefold increase. Under dark culture conditions, the IAA content increased significantly during the torpedo-shaped embryo stage (371.230 ng•g− 1), potentially promoting radicle elongation. The abscisic acid (ABA) contents from the globular embryo stage to the heart-shaped embryo stage were higher under the light treatment (749.188 ng•g− 1 and 1019.859 ng•g− 1) than under the dark treatment (525.152 ng•g− 1 and 824.755 ng•g− 1) but lower at the torpedo-shaped embryo stage (329.906 ng•g− 1 and 629.372 ng•g− 1). The gibberellic acid (GA3) content was consistently higher under the dark treatment (141.227 ng•g− 1, 81.744 ng•g− 1, 51.900 ng•g− 1 and 48.676 ng•g− 1) than under the light treatment (245.153 ng•g− 1, 176.713 ng•g− 1, 61.460 ng•g− 1 and 93.131 ng•g− 1) during somatic embryo development, suggesting that dark treatment may enhance GA3 content, thereby promoting further somatic embryo development and cotyledon embryo formation. The zeatin riboside (ZR) contents during the globular embryo stage, heart-shaped embryo stage, and cotyledon embryo stage were higher under the dark treatment (1203.354 ng•g− 1, 1173.912 ng•g− 1, and 2328.489 ng•g− 1) than under the light treatment (2216.410 ng•g− 1, 1595.285 ng•g− 1, and 3497.238 ng•g− 1), with both treatment groups exhibiting the lowest ZR value at the torpedo-shaped embryo stage and the highest value at the cotyledon embryo stage.
The synergistic or antagonistic effects of hormones regulate plant somatic embryo development, radicle development, and plant regeneration. In this study, the S. chinensis radicle primordium was initiated during the globular embryo stage and developed during the heart-shaped embryo stage. The IAA content during this period was lower under the dark treatment than under the light treatment. Under light treatment, the IAA content during the globular embryo period was up to 598.99 ng/g, which is similar to the IAA content in abnormal embryos of S. chinensis (Sun et al. 2013). We also observed that more deformed embryos were produced under light culture conditions, possibly due to the high concentration of auxin (597.565 ng•g− 1) caused by light in the early stage of somatic embryo development. The IAA content decreased to a normal level (202.628 ng•g− 1) during the heart-shaped embryo stage and induced root primordium formation. However, the IAA content increased sharply again during the torpedo-shaped embryo stage, which may have inhibited further radicle development. When there are defects in radicle development, somatic embryos cannot develop completely, resulting in gradual browning. Similar to our result, Yan et al. (2019) found that the highest IAA content was 112.44 ng/g during the rooting induction period of Pyrus betulifolia tissue culture seedlings. Therefore, we speculate that dark treatment may reduce the occurrence of abnormal embryos and that an appropriate IAA concentration has a positive regulatory effect on somatic radicle production. Studies have shown that ABA may act as an endogenous signal to negatively regulate the role of IAA in promoting root formation during root morphogenesis (Eisner et al. 2021). Li et al. (2019) also showed that ABA negatively regulates IAA-related gene expression and affects rice root growth. In this study, the ABA level was lower under the dark treatment than under the light treatment from the globular embryo formation stage to the heart-shaped embryo stage. Therefore, we hypothesized that dark treatment reduces ABA levels from the globular embryo formation stage (7 d) to the heart-shaped embryo stage (14 d), weakening the negative regulatory effect of ABA on IAA and promoting the initiation of root primordia and the induction of radicles in S. chinensis.
Sequencing data and quality control
A total of 228.69 GB of raw sequence data were generated by sequencing samples of S. chinensis at each stage of the embryonic development cycle. Clean reads (194.03 GB) were obtained by filtering the sequence data using FastQC software. The quality control results are presented in Table 3. The total transcript count was 164,109, with a GC content of 42.70% (Fig. 4).
Table 3
List of data output quality
Sample
|
Raw_Reads
|
Raw_Bases
|
Valid_Reads
|
Valid_Bases
|
Valid%
|
Q20%
|
Q30%
|
GC%
|
CK
|
144605222
|
21.70G
|
141552600
|
19.77G
|
97.88
|
98.20
|
94.31
|
46.82
|
GZ1
|
140059378
|
21.00G
|
136838952
|
19.10G
|
97.66
|
98.11
|
94.10
|
46.97
|
GZ2
|
124945708
|
18.74G
|
122886202
|
17.16G
|
98.36
|
98.08
|
94.02
|
46.72
|
GZ3
|
134389322
|
20.16G
|
123980604
|
17.33G
|
92.20
|
98.21
|
94.31
|
46.90
|
GZ4
|
127168264
|
19.08G
|
122010946
|
17.04G
|
95.93
|
98.13
|
94.14
|
46.52
|
GZ5
|
143330914
|
21.50G
|
114968252
|
16.03G
|
79.71
|
98.09
|
94.06
|
46.72
|
HA1
|
130046422
|
19.51G
|
127297976
|
17.77G
|
97.89
|
98.08
|
93.99
|
46.89
|
HA2
|
150043174
|
22.51G
|
138540892
|
19.35G
|
92.28
|
98.13
|
94.14
|
47.09
|
HA3
|
142873138
|
21.44G
|
127014986
|
17.73G
|
88.90
|
98.10
|
94.09
|
47.10
|
HA4
|
132080040
|
19.81G
|
127228978
|
17.77G
|
96.63
|
98.27
|
94.54
|
46.96
|
HA5
|
154934928
|
23.24G
|
152011756
|
21.22G
|
98.13
|
98.10
|
94.08
|
47.06
|
Functional annotation and differential expression of Unigenes
As a complete genome sequence of S. chinensis is not available, in this study, the gene sequences that were obtained were compared with those in five publicly available databases: GO, KEGG, Pfam, Swissprot, eggNOG, and NR (Table 4). The results revealed that 21994 (35.11%), 17538 (28.00%), 20717 (33.88%), 18173 (29.00%), 25385 (40.53%), and 26794 (42.78%) genes were identified in the GO, KEGG, Pfam, Swissprot, eggNOG, and NR databases, respectively.
Table 4
Statistics of unigene annotation results
Database
|
Num
|
Ratio(%)
|
NR
|
26794
|
42.78
|
GO
|
21994
|
35.11
|
KEGG
|
17538
|
28.00
|
Pfam
|
20717
|
33.08
|
SwissProt
|
18173
|
29.01
|
eggNOG
|
25385
|
40.53
|
All
|
62636
|
100.00
|
The transcriptome data for the somatic embryos of S. chinensis under the two treatments were analyzed, revealing 3538 differentially expressed genes under dark conditions compared to light. The GO enrichment analysis (Fig. 5) showed that dark treatment enriched 733 GO terms, with biological process (BP) accounting for 60.66%. These included redox reduction, response to auxin and cytokinin, and cell wall organization. Cellular component (CC) and molecular function (MF) terms accounted for 12.02% and 27.32%, respectively. Most of the differentially expressed genes under the dark treatment were related to biological processes, including hormone and cell wall composition. The KEGG enrichment analysis (Fig. 6) allocated the differentially expressed genes to 20 pathways, with plant hormone signal transduction enriched in 264 genes, phenylpropanoid biosynthesis enriched in 63 genes, starch and sucrose metabolism enriched in 155 genes, circadian rhythm-plant enriched in 66 genes, and flavonoid biosynthesis enriched in 87 genes.
Plant hormone signal transduction pathway
In this study, the plant hormone signal transduction pathway revealed that dark treatment induced the responses of various hormones, including auxin, cytokinin, and gibberellin (Fig. 7). As one of the earliest discovered and most extensively studied plant hormones, auxin plays a crucial role in multiple stages of growth and development, such as plant embryogenesis, root development and elongation, and vascular tissue differentiation. In this study, the genes involved in auxin signal transduction and root development were ARF, AUX/IAA, etc. Compared to light culture conditions, under dark conditions, ARF1 was significantly upregulated from the heart-shaped embryo stage to the cotyledon embryo stage, while ARF2 and AUX/IAA were significantly upregulated from the globular embryo stage to the cotyledon embryo stage. ARF18, an auxin response factor, is a key factor that affects the downstream auxin transcription process, including activation or inhibition, to control changes in auxin content. In this study, ARF18 was significantly upregulated from the heart-shaped embryo stage to the cotyledon-shaped embryo stage under light conditions. Phenotypic analysis of OsARF18 transgenic materials showed that OsARF18 could inhibit the elongation of rice roots (Xu 2021). Aux/IAA (auxin/indole-3-acetic acid) is an early auxin response gene, and its protein product can specifically bind to ARF. The C-terminal domain CTD (C-terminal domain) of ARF is highly homologous to Domains III and IV of the Aux/IAA protein, which form dimers through these two regions to regulate the transcription of auxin response genes (Tiwari et al. 2003). Therefore, it is speculated that ARF and Aux/IAA effectively regulate the activity of auxin by specific binding under dark conditions, thus promoting radicle development.
Several key genes enriched in the cytokinin signal transduction pathway were identified, including type B ARR, CRE1, and AHK5. CRE1, a membrane-bound protein, acts as a negative regulator of cytokinin signal transduction (Kim et al. 2006), which is highly upregulated during the globular embryo stage under dark conditions. B-ARR plays a positive regulatory role in mediating CTK signal transduction by connecting IAA and CTK (Muller and Sheen 2007) and is significantly upregulated during each period of dark culture. ARR1 can bind to the promoter of the negative regulator of the auxin signaling pathway SHY2/IAA3 to activate gene expression, thereby regulating root development (Li et al. 2015). In this study, ARR1 gene expression was downregulated under dark treatment and upregulated under light treatment. The downregulated expression of CRE1 and the significantly upregulated expression of B-ARR in the dark environment may alter cytokinin content distribution. The combination of ARR1 and auxin response factors also affects auxin level changes, regulating the formation and development of the radicle. AHK5, a cytokinin receptor (Arabidopsis AHK4 histidine kinase), transmits cytokinin signals across the cell membrane, affecting endogenous cytokinin content.
MMK1, a downstream transcription factor of gibberellin, directly targets and binds to promoter elements of light signals and modulates the gibberellin content in somatic embryos. In this study, MMK1 was upregulated under light treatment and negatively regulated radial growth. PIF3, a photosensitive interaction factor downstream of the photoreceptor, which is also a downstream transcription factor of the DELLA protein in the gibberellin pathway, negatively regulates auxin signaling. The downregulation of PIF3 expression from the globular to torpedo-shaped embryo stages induced the downregulation of AUX/IAA, which in turn affected the upregulation of the cytokinin negative regulator B-ARR and negative regulator CRE1. PIF3 is also a downstream transcription factor of the DELLA protein in the gibberellin pathway, indicating that the combined action of auxin, cytokinin signaling pathways, and gibberellin signal transduction promotes the development of somatic radicles in the later stage.
Circadian rhythm pathway
The KEGG enrichment analysis revealed that 66 genes were enriched in circadian rhythm pathways (Fig. 8), including photosensitive interaction transcription factor PIFs, cryptochrome CRY1, and downstream transcription factors HY5 and COP1. Most genes in this pathway exhibited lower expression in the dark environment than in the light treatment, suggesting that these photoreceptors regulate the growth and development pattern of somatic embryos by up- or downregulating gene expression. PIF3, a downstream factor of the photoreceptor, negatively regulates auxin signaling and is also a downstream transcription factor of the DELLA protein in the gibberellin pathway. PIF4 activates IAA19 and IAA29 by binding to the G-box (CACGTG) sequence in the promoter of auxin/indoleacetic acid genes, thereby regulating auxin and light signal transduction in plants (Kami et al. 2012). COP1 has been shown to bind to HY5 in the dark environment, leading to the inactivation or degradation of these transcription factors, and HY5 can inhibit lateral root development (Jing and Lin 2017). Therefore, it is speculated that genes related to circadian rhythm that are expressed in response to the dark environment may interact with hormone-related genes and participate in the formation and development of radicles in S. chinensis.
Phenylpropanoid biosynthesis pathway
During the initial phase of dark treatment, there was a significant increase in the expression of genes associated with phenylpropanoid synthesis (Fig. 9). A total of 63 differentially expressed genes were enriched in the phenylpropanoid biosynthesis pathway, including 17 cinnamyl alcohol dehydrogenase ELIs, 36 peroxidase PERs, and 10 shikimate hydroxycinnamoyl transferase SALATs. These findings suggest that S. chinensis somatic embryos produce numerous secondary metabolites in response to changes in environmental light conditions. Notably, peroxidases are the primary components that may influence or participate in the development of somatic embryos.
The phenylpropanoid biosynthesis pathway primarily generates lignin precursors, which are subsequently polymerized into lignin and incorporated into the cell wall (Scully et al. 2016). The transcriptome sequencing in this study revealed a significant enrichment of genes involved in the phenylpropanoid biosynthesis pathway, including peroxidase, cinnamate dehydrogenase (CAD), and cinnamyl coenzyme A reductase (CCR). Under chromium stress, upregulation of CCR and CAD genes has been found to result in increased wood deposition in the cell wall of two bean (Vicia sativa L.) cultivars (Rui et al. 2018). These findings suggest that genes encoding CCR, CAD, and peroxidase are involved in cell wall expansion by regulating lignin biosynthesis, particularly in vascular and fiber cells. As a result, they may play a role in regulating the development of S. chinensis radicles.
RT‒qPCR analysis
Nine genes (GH3, SAUR, ARF1, ARF18, AUX/IAA, MMK1, AHK4, AHK5 and PIF3) related to the plant hormone signal transduction pathway were chosen for RT‒qPCR analysis. The qPCR analysis resulted in an expression pattern of target genes that was consistent with the transcriptome data under both light and dark cultivation conditions (Fig. 10), thereby affirming the accuracy of the transcriptome. The formation and development of the radicle in S. chinensis is a complex physiological process, and changes in light conditions may lead to up- or downregulation of these genes. The results suggest that the expression regulation of these nine genes plays a crucial role in modulating endogenous hormone content.