Skip to main content

Transcriptomic and physiological analysis of the response of Spirodela polyrrhiza to sodium nitroprusside

Abstract

Background

Spirodela polyrrhiza is a simple floating aquatic plant with great potential in synthetic biology. Sodium nitroprusside (SNP) stimulates plant development and increases the biomass and flavonoid content in some plants. However, the molecular mechanism of SNP action is still unclear.

Results

To determine the effect of SNP on growth and metabolic flux in S. polyrrhiza, the plants were treated with different concentrations of SNP. Our results showed an inhibition of growth, an increase in starch, soluble protein, and flavonoid contents, and enhanced antioxidant enzyme activity in plants after 0.025 mM SNP treatment. Differentially expressed transcripts were analysed in S. polyrrhiza after 0.025 mM SNP treatment. A total of 2776 differentially expressed genes (1425 upregulated and 1351 downregulated) were identified. The expression of some genes related to flavonoid biosynthesis and NO biosynthesis was upregulated, while the expression of some photosynthesis-related genes was downregulated. Moreover, SNP stress also significantly influenced the expression of transcription factors (TFs), such as ERF, BHLH, NAC, and WRKY TFs.

Conclusions

Taken together, these findings provide novel insights into the mechanisms of underlying the SNP stress response in S. polyrrhiza and show that the metabolic flux of fixed CO2 is redirected into the starch synthesis and flavonoid biosynthesis pathways after SNP treatment.

Peer Review reports

Background

Duckweed plants are the smallest and simplest flowering aquatic plants in the world, being only a few millimetres in size. These tiny plants belong to the Lemnaceae family and Lemnoideae subfamily, which includes five genera: Landoltia, Lemna, Spirodela, Wolffiella, and Wolffia [1]. The plant body consists of a leaf-like structure (the frond), with or without a root-like structure. The plants propagate mainly via asexual reproduction and have a long yearly production period with an almost exponential growth rate, producing biomass faster than most other plants [2] and with a high content of starch [3] and protein [4] and a low content of lignocellulose [5]. In addition, they have extremely wide adaptability to different environments [6].

In recent years, S. polyrrhiza, as a chassis plant, has provided breakthroughs in biosynthesis. The advantages are as follows: (1) In terms of biological structure, S. polyrrhiza, the largest duckweed in North America, is approximately 7 to 15 mm long and has enormous fronds, making it useful for laboratory molecular biology research and tissue culture operations [7]. (2) This plant has a clear genetic background, the duckweed family's first species with a fully sequences genome was S. polyrrhiza, and a high-confidence genetic map of this plant has been created based on the genomic data [8]. (3) This plant is a well established genetic modification platform: to date, the genetic transformation of S. polyrrhiza has been achieved successfully, and the transformation efficiency has reached 5%-13% [9,10,11]. (4) This plant contains abundant precursor substances, including a variety of beneficial secondary metabolites, such as terpenoids and phenolic compounds [12].

Sodium nitroprusside (SNP), an established NO donor used in plant science research, simultaneously releases NO, cyanide (CN) and Fe (II) in solution [13,14,15]. Research has demonstrated that SNP treatment improves the activity of antioxidants in blueberry fruit, tomato, soybean sprouts, okra, chard and alfalfa [16,17,18,19,20,21]. In recent years, research has shown that SNP participates in the biosynthesis of secondary metabolites such as alkaloids [22], flavonoids [23], terpenoids [24] and saponins [25]. The majority of these studies, however, concentrated on vegetables and herbaceous plants [26,27,28], and few studies have been conducted on aquatic plants [29].

Previous studies have shown that S. polyrrhiza copes with abiotic stress by fixing CO2 into the starch biosynthesis pathway; as a result, biomass rich in starch rapidly accumulates [30]. We are curious about how the energy distribution of S. polyrrhiza changes under these conditions; it remains unknown whether the metabolic flux is primarily directed into the starch or flavonoid biosynthesis branch after SNP treatment. In our work, S. polyrrhiza was treated with SNP, and the changes in biomass accumulation, starch content, flavonoid content, and antioxidant-related enzyme activity dynamics of S. polyrrhiza were studied. Moreover, we conducted research on the response of S. polyrrhiza to SNP stress using transcriptomic approaches. This provided molecular support for research on the metabolic flux in S. polyrrhiza under SNP treatment and the use of S. polyrrhiza as a bioreactor. This is crucial for the development of S. polyrrhiza as a model plant for synthetic biology.

Result

Photosynthetic pigment content and chlorophyll fluorescence

Quantification of chlorophyll is frequently used to estimate how plants react to various stresses. As the SNP concentration increased, the photosynthetic pigment content showed a rising trend after first decreasing. The chlorophyll a, chlorophyll b, and carotenoid contents were reduced by 48%, 39.7%, and 35.3% at 0.01 mM SNP in comparison to the control group, respectively (Fig. 1 A, B, C). As the SNP concentration increased (Fig. 1 D, E, F), the PSII maximum photochemical conversion efficiency (Fv/Fm) showed a downwards trend, being lowest at 0.5 mM SNP and decreasing by 76.2% compared with the control group. The absorption flux per reaction centre (ABS/RC) and dissipation energy per reaction centre (DIo/RC) change trends were contrary to that of Fv/Fm, peaking at 0.5 mM SNP and increasing by 2.94- and 16.9-fold compared with the control group, respectively.

Fig. 1
figure 1

The effect of SNP stress on the photosynthetic pigment content and chlorophyll fluorescence of S. polyrrhiza. A chlorophyll a content; B chlorophyll b content; C carotenoid content; D Fv/Fm; E ABS/RC; F DIo/RC. The values are the means ± standard errors of triplicate assays. Bars with different lowercase letters are significantly different by one-way ANOVA (P < 0.05)

Effects of SNP treatment on the physiological changes of S. polyrrhiza

As the SNP concentration increased, the fresh weight showed a significant downwards trend after rising first and peaked at 0.01 mM SNP, increasing by 12.58% compared with the control group (Fig. 2A). Similar to the fresh weight trend, the dry weight of S. polyrrhiza showed an increase of 27.87% in the 0.01 mM SNP treatment samples compared with the control group (Fig. 2B). The starch content change trend was consistent with those of the fresh and dry weight; the starch content increased and then decreased, peaking at 0.025 mM SNP, which was approximately 1.5 times more than that in the control (Fig. 2C). The soluble protein change trend was consistent with that of starch. The highest soluble protein content was demonstrated in plants grown on medium with 0.025 mM SNP, increasing by 42.06% over the control group (Fig. 2D). The malondialdehyde (MDA) content was lower in the SNP treatment group than in the control group, with the lowest significant MDA content observed at 0.3 mM SNP, decreasing by 55.63% (Fig. 2E), which implied that exogenous NO at a specific concentration effectively inhibited the accumulation of MDA. Compared with the control group, the hydrogen peroxide (H2O2) content of the 1.0 mM SNP-treated group increased from 20.25 mol/g FW to 31.41 mol/g FW. Excessive accumulation of H2O2 may further aggravate membrane lipid peroxidation (Fig. 2F). The superoxide dismutase (SOD) activity showed a significant downwards trend after first increasing. Compared to the control group, the SOD activity of the 0.01 mM SNP treatment group increased rapidly, from 149.67 U/g FW to 201.67 U/g (Fig. 2G). The catalase (CAT) activity peaked at 0.025 mM SNP and was 108% higher than that in the control group (Fig. 2H). The peroxidase (POD) activity showed a significant downwards trend after rising first with increasing SNP concentration and was significantly higher in the 0.1 mM SNP treatment group (3860.03 U/g) than in the control group (1015.43 U/g) (F ig. 2I).

Fig. 2
figure 2

Effects of different concentrations of SNP on changes in the physiological and phytochemical characteristics of S. polyrrhiza. A fresh weight; B dry weight; C starch content; D soluble protein content; E MDA content; F H2O2 content; G SOD activity; H CAT activity; I POD activity. The values are the means ± standard errors of triplicate assays. Bars with different lowercase letters are significantly different by one-way ANOVA (P < 0.05)

Effects of SNP treatment on the flavonoid content of S. polyrrhiza

Previous studies have reported that the duckweed species S. polyrrhiza has been widely utilized as folk medicine in China, Korea and a few European nations [31]. Duckweeds are medicinal herbs that do not have severe side effects [32]. Flavonoids, such as apigenin, orientin, vitexin, and luteolin-7-O-glucoside, are the main pharmacological components of duckweed [33]. Orientin, vitexin and luteolin-7-O-glucoside have anti-inflammatory, antioxidant, antitumour and other effects [34,35,36]. We determined the levels of three bioactive compounds of the flavonoid biosynthesis pathway: orientin, vitexin, and luteolin-7-O-glucoside. The orientin content showed a significant downwards trend after rising first, peaking at 0.025 mM SNP and increasing by 18.9% compared with the control group (Fig. 3A). The vitexin content significantly increased at SNP treatment concentrations ranging from 0.01 mM to 0.1 mM compared to the control group (Fig. 3B). The luteolin-7-O-glucoside content significantly increased at SNP treatment concentrations ranging from 0.025 mM to 0.5 mM and peaked at 0.025 mM SNP, and the value was 15.8% higher than that in the control group (Fig. 3C).

Fig. 3
figure 3

Effects of different concentrations of SNP on the flavonoid content of S. polyrrhiza. A orientin content; B vitexin content; C luteolin-7-O-glucoside content; D proanthocyanidin content. The values are the means ± standard errors of triplicate assays. Bars with different lowercase letters are significantly different by one-way ANOVA (P < 0.05)

Proanthocyanidins are a kind of polyphenol. Their special structure determines their strong antioxidation activity, and they have anticancer, antimutagenic, cardiovascular protection and other effects [37]. The proanthocyanidin content showed a downwards trend after rising first with increasing SNP concentration. The maximum value was reached at 0.025 mM SNP (Fig. 3D). This result is consistent with the observation that at 0.01 and 0.025 mM SNP, leaves exhibited chlorosis and accumulation of purplish red pigment on the back (Additional file 1: Figure S1).

In this study, we wanted to explore the redirection of metabolic flux from fixed CO2 into starch or flavonoid biosynthesis branches after SNP treatment. Therefore, we focused on the changes in starch and flavonoid contents and found that under 0.025 mM SNP treatment, the starch and flavonoid contents increased significantly. Based on these results, S. polyrrhiza treated with 0.025 mM SNP was selected for use in subsequent RNA expression analyses.

Global gene analysis in S. polyrrhiza after 0.025 mM SNP treatment

In total, six cDNA libraries were constructed and sequenced. While robust data were collected, after data filtering, 43.50, 44.84, 44.32, 39.85, 48.41, and 39.08 million high-quality reads were obtained in the CK and SNP treatments, respectively. Subsequently, these clean reads were de novo assembled into 33,684 genes, with a mean length of 1,360 bp and an N50 length of 2,791 bp (Table 1).

Table 1 Characteristics of the de novo assembled transcriptome in S. polyrrhiza after 0.025 mM SNP treatment

The unigenes mapped to the database searches are shown in Table 2. Using the National Center for Biotechnology Information (NCBI) NR database, 15,127 unigenes were mapped to ten identified plant species and ‘other’ plant species (Additional file 2: Figure S2). Using the SwissProt database, 11,816 unigenes were annotated and reviewed with the UniProt Knowledgebase. Searches against the KOG database led to 9,341 unigenes being divided into 25 categories. Using the Kyoto Encyclopedia of Genes and Genomes (KEGG) database, 14,966 unigenes were enriched in 137 biochemical pathways, such as plant‒pathogen interaction, carbon metabolism, plant hormone signal transduction and phenylpropanoid biosynthesis. Transcriptome profiles were compared between the CK and SNP - treated groups. A total of 2,776 DE genes were identified, of which 1,425 were upregulated, while 1,351 were downregulated (Fig. 4).

Table 2 Functional annotation of unigenes in the NR, SwissProt, KEGG, and COG databases
Fig. 4
figure 4

Volcano plots of the differentially expressed genes (DEGs) in S. polyrrhiza after 0.025 mM SNP treatment. The red dots indicate upregulated genes with significant differences, and the blue dots indicate downregulated genes with significant differences

Investigation of differentially expressed transcription factors under SNP treatment

Transcription factors (TFs) are key regulators that temporarily and spatially turn on or off the transcription of their target genes by binding certain upstream elements [38]. The de novo assembly data identified 818 TFs. Among these factors, 208 TFs (125 upregulated and 83 downregulated) exhibited significant differences in expression levels in the 0.025 mM SNP treatments (Fig. 5). In this study, several TFs responded significantly to SNP stress (Fig. 6), including ERF (22), BHLH (20), NAC (17), and WRKY (15) (Additional file 4: Table S1). ERFs are specific to plants and play a role in a variety of developmental processes, such as flowering [39], seed development [40], and fruit ripening [41]. Among these TFs, the most abundant were those of the ERF family, which accounted for 11.7% of the total TFs. The RNA-Seq results showed that the expression of 15 ERF-related genes was upregulated and that of 8 genes was downregulated. BHLH TFs are involved in stress tolerance, pathogen defence, and nutrient uptake but also play an important role in the production of secondary metabolites such as flavonoids and anthocyanins [42, 43]. Twenty BHLH TFs were differentially expressed between the SNP concentrations; ten of these genes were upregulated, and ten were downregulated. A crucial part of plant biological and abiotic stress responses is played by NAC TFs, which are unique to plants [44]. Of the seventeen members of the NAC TF family, thirteen were upregulated, and four were downregulated. WRKY TFs participate in the response to injury, ageing, development, and disease [45]. Additionally, WRKY TFs control how plants grow and develop as well as how they react to stress. Fourteen of the fifteen members of the WRKY TF family were upregulated, while one gene was downregulated.

Fig. 5
figure 5

Classification of differentially expressed transcription factors in S. polyrrhiza after 0.025 mM SNP treatment

Fig. 6
figure 6

Expression pattern of genes associated with TFs under SNP treatment in S. polyrrhiza. A ERF family; B BHLH family; C WRKY family; D NAC family. The expression level was Z score normalized. The colour gradient (green to red) corresponds to the gene expression levels (low to high). Each row represents a gene, with its name given in parentheses

Investigation of differentially expressed photosynthesis-related genes under SNP treatment

Based on their functions, the genes involved in photosynthesis were selected for this study and classified into two groups (Fig. 7). Most of the photosynthesis-related genes belong to the first class, which is involved in the light reaction, such as the photosystem I and II reaction centre genes PsbP, PsbQ, PsbS, PsbW, PsbY, Psb27, PsaD, PsaE, PsaF, PsaG, PsaL, PsaO and cytochrome b/6 complex PetC, which were expressed at lower levels in the SNP treatment (Additional file 5: Table S2). The second category consists of antenna proteins, which are light-harvesting chlorophyll a/b binding (LHC) proteins. Most of the genes (including Lchas and Lchbs, which are light-harvesting protein genes) showed similar expression patterns to the light reaction genes (Additional file 5: Table S2).

Fig. 7
figure 7

Heatmap displaying the expression levels of differentially expressed genes (DEGs) related to photosynthesis in S. polyrrhiza. A Expression pattern of genes involved in photosynthesis-antenna proteins under SNP treatment in S. polyrrhiza; B Expression pattern of genes involved in photosynthesis under SNP treatment in S. polyrrhiza. The expression level was Z score normalized. The colour gradient (green to red) corresponds to the gene expression levels (low to high). Each row represents a gene, with its name given in parentheses

Investigation of differentially expressed starch metabolism-related genes under SNP treatment

Transcriptome data indicated that carbon fixation pathways were significantly affected by 0.025 mM SNP treatment. Following SNP treatment, the expression of a few genes related to the Calvin cycle was reduced. The SNP-downregulated genes included those encoding sedoheptulose-bisphosphatase (SEBP), fructose-bisphosphate aldolase (ALDO), fructose-1,6-bisphosphatase (FBP), ribose 5-phosphate isomerase A (RPI), ribulose-bisphosphatecarboxylase (RBCS), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), and phosphoribulokinase (RPK). Starch metabolism is important for the supply of energy during S. polyrrhiza growth and development. Twenty-three DEGs were enriched in starch and sucrose metabolism pathways (Additional file 6: Table S3). Numerous important enzymes are involved in starch biosynthesis, such as ADP-glucose pyrophosphorylase (AGP). Half the transcripts encoding AGPase were upregulated, while the other half were downregulated (Fig. 8). The results showed that under SNP treatment, the expression of a transcript encoding AGPase (Unigene0007838) was upregulated from 35 to 90 FPKM, while the expression of another transcript encoding AGPase (Unigene0019277) was downregulated from 99 to 38 FPKM. This may be due to differences in gene expression regulation, which lays the foundation for the selection of target enzymes for further research in the future. The degradation of starch is catalysed by α-amylase and β-amylase [46], which are responsible for the degradation of starch to smaller hydrocarbons. β-amylase was downregulated under SNP treatment, while the expression of α-amylase did not change significantly. A transcript encoding β-amylase exhibited an expression level of 69.8 FPKM in the control group, and the value decreased to 28.5 FPKM under SNP treatment.

Fig. 8
figure 8

Expression pattern of genes involved in starch metabolism under SNP treatment in S. polyrrhiza. The arrows indicate the directions of catalytic reactions or transport. Red indicates upregulated expression, green indicates downregulated expression, and blue indicates up/downregulated expression. Dotted arrows indicate omitted steps. G-3P, glyceraldehyde 3-phosphate; fructose-6P, fructose 6-phosphate; glucose-6P, glucose 6-phosphate; PGI, glucose-6-phosphate isomerase; glucose-1P, glucose 1-phosphate; AGPase, ADP-glucose pyrophosphorylase; ADP-glucose, adenosine- 5-diphosphoglucose; SBE, starch branching enzyme; BAM, beta-amylase

Investigation of differentially expressed reactive oxygen species (ROS) homeostasis -related genes under SNP treatment

ROS are chemically reactive molecules that contain oxygen, and carbohydrates, lipids, proteins and DNA and are all damaged by the overproduction of ROS [47, 48]. The data from our transcriptome analysis revealed that the transcript abundance of genes encoding different components of the ROS scavenging machinery differed dramatically between the control and SNP treatments (Fig. 9). A total of 34 putative genes were found encoding proteins with antioxidant properties (Additional file 7: Table S4). There were 15 upregulated unigenes that encoded proteins such as glutathione s-transferases (GST), glutathione peroxidase (GPX), ascorbate peroxidase (APX), lipoxygenase (LOX) , monodehydroascorbate reductase (MDAR), and alternative oxidase (AOX). There were 6 downregulated unigenes that encoded proteins such as SOD and GST. Nine of the seventeen POD-encoding transcripts were upregulated, while eight were downregulated. However, there were no appreciable variations in the expression of CAT-encoding transcripts under SNP treatment compared to the untreated samples.

Fig. 9
figure 9

Expression pattern of genes involved in antioxidant activity under SNP treatment in S. polyrrhiza. The expression level was Z score normalized. The colour gradient (green to red) corresponds to the gene expression levels (low to high). Each row represents a gene, with its name given in parentheses

Investigation of differentially expressed flavonoid biosynthesis-related genes under SNP treatment

The phenylpropanoid biosynthesis route in plants contributes to numerous different biosynthetic branches, such as flavonoid biosynthesis and anthocyanin biosynthesis (Fig. 10). Analysis of the S. polyrrhiza flavonoid synthesis pathway identified 13 differentially expressed structural genes (Additional file 8: Table S5). Phenylalanine ammonialyase (PAL), cinnamate 4-hydroxylase (C4H) and 4-hydroxycinnamoyl-CoA ligase (4CL) are universal factors involved in flavonoid biosynthesis [49]. The expression of transcripts encoding PAL, C4H, and 4CL was upregulated 3.5-, 1.8-, and 2.8-fold, respectively, after SNP treatment. Chalcone synthase (CHS), chalcone isomerase (CHI) and flavanone 3-hydroxylase (F3H) are the enzymes that catalyse the first three reactions of the flavonoid biosynthesis branch. The expression of transcripts encoding CHS, CHI and F3H was upregulated 1.5-, 1.4- and 1.8-fold, respectively, after SNP treatment. Meanwhile, the expression of flavonoid-3′-hydroxylase (F3′H), dihydroflavonol 4-reductase (DFR), anthocyanidin synthase (ANS) and anthocyanidin 3-O-glucosyltransferase (UFGT) was upregulated 1.1-, 3.2-, 2.6- and 2.8-fold, respectively, after SNP treatment.

Fig. 10
figure 10

Biosynthetic pathway of flavonoids in S. polyrrhiza. Heatmap showing the changes in transcripts involved in flavonoid biosynthesis. Upregulated genes are shown in red. PAL, phenylalanine ammonium lyase; C4H, cinnamate 4-hydroxylase; 4CL, 4-coumarate-CoA ligase; CHS, chalcone synthase; CHI, chalcone isomerase; F3H, flavanone-3-hydroxylase; F3′H, flavonoid-3′-hydroxylase; DFR, dihydroflavonol reductase; ANS, anthocyanin synthase; UFGT, anthocyanidin 3-O-glucosyltransferase

Quantitative real-time PCR (qRT-PCR) validation

To further confirm the veracity and reliability of the DEG data generated by RNA-Seq, 17 DEGs (Additional file 9: Table S6) were selected for qRT‒PCR verification, including genes related to the flavonoid biosynthesis pathway, NO biosynthesis pathway, photosynthesis and TFs. The results showed that gene expression (17) was consistent with the RNA-Seq data (Additional file 3: Figure S3).

Discussion

It has been reported that SNP stress probably triggers the NO biosynthesis pathway, resulting in a high endogenous NO level [19, 50]. To date, the main NO synthesis pathways identified are the NOS (nitric oxide synthase) and NR (nitrate reductase) pathways [51]. The NOS pathway is mainly dependent on NO synthase (NOS), which generates NO and L-citrulline by oxidizing L-arginine [52]. In the NR pathway, nitrate reductase (NR) can reduce NO3 to NO2, nitrite reductase (NiR) can reduce NO2 to NO, and NR can also oxidize NO to nitrate to scavenge NO, which shows that NR is an important regulator of NO homeostasis in plants [51]. In the transcriptomic analysis of the key enzymes for NO synthesis, the NRT and NR genes exhibited increased expression under SNP stress (Additional file 10: Table S7). There is evidence that NO also increases the cytosolic level of free Ca2+ through upregulation of Ca2+ channels, transporter proteins, cGMP and cADPR [53, 54]. In the transcriptomic analysis of calmodulin expression, the genes encoding calmodulin, calcium-binding protein and the cyclic-nucleotide gated channel were upregulated, indicating that Ca2+ was regulated in S. polyrrhiza under SNP stress. Notably, every signalling pathway engages protein kinases, which can activate a variety of downstream TF family proteins [55]. In this study, several TFs responded significantly to SNP stress, including ERF (22), BHLH (20), NAC (17), and WRKY (15). However, more research is needed to determine whether these TFs work independently or in concert. A series of downstream NO-responsive genes are expressed as a result of TF activation, and their transcription further controls the equilibrium of cellular metabolism.

The most significant photosynthetic pigment in plant chloroplasts is chlorophyll (Chl), which is crucial for absorbing light and transmitting light energy in antenna systems [56]. The findings of this study showed that SNP treatment caused the chlorophyll content of S. polyrrhiza to decrease, and similar findings were seen in other investigations, such as in Torreya grandis [57], chard [20] and soybean seedlings [58]. This may be due to the damage to photosynthetic pigments caused by SNP stress or the degradation of chlorophyll induced by NO [59]. In addition, there are also reports indicating that chlorophyll content increases under SNP treatment, such as in wheat seedlings [60] and barley seedlings [61]. The measurement of chlorophyll fluorescence characteristics, which is practical, precise, and sensitive, is frequently employed to reflect the photosystem alterations that occur in plants under environmental stress [62]. Generally, plants grown under stress have a lower Fv/Fm than nonstressed plants [63], which is consistent with the performance of S. polyrrhiza (Fig. 1D). Compared with the control group, SNP treatment caused an increase in PSII activity in S. polyrrhiza, which is indicated by the lower Fv/Fm and higher ABS/RC and DIo/RC. The increase in ABS/RC in S. polyrrhiza under SNP treatment can be explained by a decrease in the number of active reactive centres (RCs) of PSII, which might serve as a defence mechanism to reduce the burden to its systems under stress. These findings are consistent with previous reports on T. grandis [57]. Antenna proteins are essential for absorbing solar energy and for photoprotection under stress [64]. These proteins assume a conformation that allows them to release extra energy from excitation as heat when plants are exposed to stress that could cause photooxidative damage [65]. We found that SNP treatment downregulated the expression levels of LHCA4, LHCA5, LHCA6, LHCB2.1, and LHCB5. Similar results were also reported by a previous study in alfalfa seedlings [21]. Therefore, we propose that SNP reduced photosynthesis by downregulating the expression levels of LHCA4, LHCA5, LHCA6, LHCB2.1, and LHCB5.

Starch is one of the major storage compounds of duckweed and usually accumulates under stress conditions, which is an adaptive strategy gained during the long evolutionary process to enable the plant to survive adverse environments [30]. The starch content reached 4.57 mg/g of fresh weight in the 0.025 mM SNP-treated samples and 1.83 mg/g of fresh weight in the control samples (Fig. 2C). As a result, the total starch accumulation in the treated samples was 1.5 times higher than that in the control samples. Starch accumulation is closely related to starch synthesis and metabolism. Two tactics are generally used to increase starch accumulation: increasing the amount of substrate from other metabolic pathways used for starch and sucrose metabolism or decreasing the pace at which starch is degraded [30]. AGPase is essential for controlling the accumulation of starch in plants [66]. In our study, the RNA-seq results showed that the expression of half of the key AGPase enzyme-encoding genes in the starch biosynthesis pathway was upregulated under SNP stress, while that of the other half was downregulated. Previous studies have shown that under conditions of nutritional starvation, the starch content and the expression of the AGPase gene in Landoltia punctata significantly increase [67]. We hypothesize that starch accumulation is also related to the decreased starch degradation in S. polyrrhiza. β-amylase is considered the main hydrolytic enzyme for breaking down starch granules. Previous studies have shown that the β-amylase gene is downregulated in Landoltia punctata, L. turionifera, kiwifruit and alfalfa [21, 68,69,70], which is consistent with our results, suggesting that it probably plays an important regulatory role in starch accumulation.

Overproduction of ROS during abiotic stress is typically observed, and this might result in some damage that eventually leads to oxidative stress [59]. In this study, we further explored the effects of SNP on the oxidative damage of S. polyrrhiza. As a result of the SNP treatments, the MDA levels in S. polyrrhiza increased, and ROS accumulation was observed, indicating oxidative damage. This is because elevated ROS levels may result in oxidative stress, which causes membrane damage and lipid peroxidation, and MDA can been used as a marker of this damage [71]. ROS-scavenging enzymes are crucial components of the plant defence mechanism in response to ROS production. Antioxidant enzymes and antioxidants synthesized in plants can remove excess ROS in response to adversity. In regard to enzymatic antioxidants, the enzyme SOD serves as the first line of defence by converting O2- into H2O2, which is then metabolized by the enzymes APX, POD, and CAT [72]. SOD activity significantly increased at SNP treatment concentrations ranging from 0 mM to 0.01 mM compared to the control group. From 0.025 mM to 1.0 mM SNP, the SOD activity showed a downwards trend. The RNA-seq results showed that the gene expression of S. polyrrhiza SOD was lower after 0.025 mM SNP treatment. Although SOD activity was unchanged after 3 days of treatment with 0.025 mM SNP, we speculate that SOD activity might also show a downwards trend after treatment with 0.025 mM SNP for more than 3 days. We speculate that the change in SOD gene expression precedes SOD activity. The POD activity significantly increased at SNP treatment concentrations ranging from 0.025 mM to 0.5 mM compared to the control group. The RNA-seq results showed that nine of the seventeen POD-encoding transcripts were upregulated, while eight genes were downregulated under 0.025 mM SNP treatment. Therefore, the increased in POD activity is very understandable. This effect may be attributed to differences in gene expression regulation. In addition, the enzymatic properties of POD isoenzymes may be differ. A study of S. polyrrhiza indicated that enhancing antioxidative enzyme activity, such as that of POD and CAT, confers resistance to salt stress [73, 74]. In addition, the expression of genes related to antioxidants, such as GPX, GST, MDAR, APX, LOX, and AOX, was upregulated by SNP treatment. Research has shown that the expression of genes such as GPX was upregulated in SNP-treated Ganoderma lucidum mycelia [47].

The content of flavonoids such as orientin, vitexin, luteolin-7-O-glucoside and proanthocyanidins showed a significant downwards trend after rising first under SNP treatment. Transcriptome analysis showed that after SNP treatment, DEGs were highly enriched in phenylpropanoid metabolism (Fig. 10). Naringenin is one of the crucial substrates for the synthesis of flavonoids. Moreover, DEGs that encode the key enzymes in naringenin synthesis, such as PAL, C4H, 4CL, CHS, and CHI had high expression levels under 0.025 mM SNP. Therefore, we speculate that these genes were upregulated, which may have led to an increase in naringenin biosynthesis. Naringenin can be catalysed by F3′H to produce flavonoids such as orientin, vitexin, and luteolin-7-O-glucoside. Meanwhile, the expression of F3′H was higher than 70 FPKM. This possibly explains why S. polyrrhiza under 0.025 mM SNP treatment accumulated flavonoids. In addition, naringenin can be catalysed by F3H, DFR and ANS to produce proanthocyanidins. Simultaneously, the expression of genes such as F3H, DFR, and ANS was also upregulated. This may be the main reason why S. polyrrhiza treated could accumulate proanthocyanidins under 0.025 mM SNP treatment. Overall, we speculate that flavonoids accumulated through the high expression of flavonoid-related genes, which was consistent with the observation in Landoltia punctata under nutrient starvation [75]. Furthermore, in the 0.01 and 0.025 mM SNP treatments, purple colouration accumulated on the underside of the frond, whereas the control group showed no discernible changes. It has been reported that the reddish-purple tint on the underside of its fronds is a result of anthocyanin production [75]. In our study, the genes encoding anthocyanin synthesis, such as F3H, DFR, ANS, and UFGT were upregulated. We speculate that SNP treatment may have contributed to anthocyanin accumulation in S. polyrrhiza, which was consistent with the observation in Landoltia punctata under nutrient starvation [75].

S. polyrrhiza is a promising chassis plant. It is important to monitor the metabolic flux changes in S. polyrrhiza. Previous studies have indicated that this plant responds to various abiotic stresses, such as nitrogen starvation [30], maleic hydrazide [41], and uniconazole [76], and the metabolic flux is mainly directed into the starch branch. Therefore, starch biosynthesis was enhanced. This may be explained as a stress escape or stress avoidance response to complete the life cycle in advance by storing most carbon nutrients and energy in starch [75]. Our study found that after treatment of S. polyrrhiza with 0.025 mM SNP, the metabolic flux was partly directed to the starch branch and partly to the flavonoid branch. The production of flavonoids, which are essential for plant defence against pathogens [77], can partially explain why S. polyrrhiza is rarely infected by pathogens.

Our study provides insights for further development of S. polyrrhiza as a chassis plant for the synthesis of secondary metabolites.

Conclusions

SNP (0.025 mM) can effectively induce starch and flavonoid accumulation in S. polyrrhiza. Concretely, SNP may induce starch accumulation in S. polyrrhiza by downregulating the expression of key genes involved in starch degradation, while flavonoid accumulation in S. polyrrhiza may be induced by upregulating the expression of key genes involved in flavonoid biosynthesis. Moreover, SNP could promote the expression of some genes participating in starch and flavonoid accumulation, probably through the regulation of some TFs. These results suggest that after treatment of S. polyrrhiza with 0.025 mM SNP, the metabolic flux may be partly directed to the starch biosynthesis branches and partly into the flavonoid biosynthesis branches.

This provides a method to increase S. polyrrhiza starch and flavonoid accumulation in a selective manner. The genes regulated by SNP treatment would provide good candidates for improving the flavonoid content by genetic engineering in S. polyrrhiza.

Methods

Plant materials and treatments

Spirodela polyrrhiza strain P143 plants were described in the previous paper [78] and preserved at the Prof. Yong Wang’ s Lab. Plants of S. polyrrhiza were grown aseptically on DATKO medium as described by Wang and Kandeler [78] under a 16 h light/8 h dark cycle with a light intensity of ~ 45 µmolm−2s−1 at the plant level. The temperature of the culture room was maintained at 22 ± 2 °C. S. polyrrhiza plants were transferred to fresh medium every 10 days to minimize the effect of nutrient shortage. Plants were precultured under identical growth conditions for 8 days and then divided into eight groups as follows: (1) control (CK), without any treatment; (2) 0.01 mM SNP treatment; (3) 0.025 mM SNP treatment; (4) 0.05 mM SNP treatment; (5) 0.1 mM SNP treatment; (6) 0.3 mM SNP treatment; (7) 0.5 mM SNP treatment; and (8) 1.0 mM SNP treatment. Three biological replicates were performed. SNP was dissolved in distilled water, filter-sterilized and then added to DATKO medium. Plant samples were collected 1, 2, and 3 days after SNP treatment for measurement of NO content. Only the samples collected on the 3rd day were used for the analysis of other physiological parameters. The starch content showed a significant upwards trend within the range of 0-0.025 mM SNP. However, starch content showed a significant downward trend within the range of 0.025-1.0 mM SNP. Therefore, we selected five representative concentrations, namely, 0, 0.01, 0.025, 0.1, and 0.5 mM SNP, for determination of the flavonoid content.

Photosynthetic pigments contents and chlorophyll fluorescence

The photosynthetic pigment content was determined according to the method described by Arnon with some modifications [79]. Approximately 0.05 g (fresh weight) of sample was macerated in 3 mL of absolute ethanol, placed in the dark for 24 h at 28 °C, and centrifuged at 5000 × g for 5 min. The supernatant was collected, and the absorbance of chlorophyll a, b and carotenoids was recorded at 663, 645, and 470 nm, respectively. The PSII parameters were analysed using a pulse-amplitude modulation chlorophyll fluorometer (Handy-PEA, Hansatech, England). Measurements were performed with well-developed leaves treated with different concentrations of SNP, with three measurements per replicate. Each leaf was measured after dark adaptation for 30 min.

Assays of physiological parameters

Fresh and dry weights were measured with a digital balance. The activities of SOD, CAT and POD were determined using enzyme activity kits (Beijing SolarBio Science & Technology Co., Ltd., Beijing, China). The contents of starch, soluble protein, MDA and H2O2 were determined using a content measurement kit (Beijing SolarBio Science & Technology Co., Ltd., Beijing, China); the proanthocyanidin content was determined by the Luthar and Kreft method [80]. Approximately 0.05 g of sample was macerated in 5 mL 60% ethanol, ultrasonicated at 30 °C for 30 min, and centrifuged at 5000 × g for 10 min at 25 °C. An aliquot of 1 mL of supernatant was mixed with 2 mL of 4% vanillin solution and 1 mL of concentrated hydrochloric acid for a 20 min reaction at room temperature, and the absorbance was measured at 530 nm.

Determination of orientin, vitexin and luteolin-7-O-glucoside content

The dry materials were ground to a fine powder. Approximately 50 mg of powder was extracted in 5 mL of methanol by ultrasonication at 50 °C for 50 min, followed by centrifugation. Samples were separated using ultrahigh-performance liquid chromatography (UHPLC) with an ACQUITYUPLCBEH-C18 (2.1 mm × 50 mm, 1.7 μm, Waters). In the solvent system, mobile phase A was 0.1% (v: v) formic acid, and mobile phase B was methyl-alcohol. The gradient conditions applied were as follows: 0~7 min, 25% ~ 40% B; 7~10 min, 40% ~ 50% B; 10~11 min, 50% ~ 40% B; 11~12 min, 40% ~ 25% B; 12~14 min, 25% B. The flow rate was set at 0.3 mL/min, and the injection volume was 1 μL.

RNA-seq and bioinformatics analysis

As described above under Plant materials and treatments, there were 3 flasks of repeat cultures for the control and 0.025 mM SNP treatment samples prepared for RNA-Seq. Three individual biological replicates (6 samples in total) were used for transcriptome analysis of S. polyrrhiza. Total RNA samples were extracted. The libraries were prepared using mRNA isolated from total RNA using Dynabeads oligo (dT). Then, the enriched mRNA fragments were randomly cleaved into short fragments of 200 nt-700 nt by adding fragmentation buffer, and the obtained short fragmented mRNA was used as the template. First-strand cDNA was generated using random hexamer primers. Second-strand cDNA was generated using buffer, dNTPs, RNaseH, and DNA polymerase I (Invitrogen). The cDNA fragments were then purified with a QiaQuick PCR extraction kit, the ends were repaired, poly(A) was added, and the fragments were ligated to Illumina sequencing adapters. Finally, the ligation products were size-selected using agarose gel electrophoresis, PCR amplified, and subjected to RNA-seq using an Illumina platform from Gene Denovo Biotechnology (Guangzhou, China).

Raw reads obtained from Illumina sequencing were further filtered to obtain high-quality clean reads by removing reads containing adapters, reads with an N ratio greater than 10%, reads containing a single base, and low-quality reads. De novo assembly of clean reads was performed using Trinity [81]. The fragments per kilobase of transcript per million mapped reads (FPKM) method was used to determine the expression levels of each gene. In this study, the level of differential expression for each transcript, with a criterion of |log2(fold-change)| >1 and FDR < 0.05 at different growth stages, was employed to identify DEGs. DEGs were subjected to enrichment analysis against the GO and KEGG databases to identify changes in biological functions [82,83,84].

Quantification of gene expression levels

Total RNA was isolated from the CK and 0.025 mM SNP-treated S. polyrrhiza samples using an Eastep Super Total RNA Extraction Kit (Promega, Shanghai) according to the manufacturer’s protocol. The integrity of the total RNA was examined using 1.0% agarose gel electrophoresis. Primer sequences for selected DEGs (Additional file 11: Table S8) were designed using Primer Premier 5 and synthesized by Sales Genomics CN (Tianjin, China). First-strand cDNA was synthesized using a FastKing RT Kit containing gDNase (TIANGEN Biotech, China) with one cycle at 42 °C for 15 min and then 95 °C for 3 min. qRT‒PCR for gene expression analysis was carried out using SuperRealPreMix Plus (SYBR Green) (TIANGEN Biotech, China) with one cycle at 95 °C for 15 min followed by 40 cycles of 95 °C for 10 s and 60 °C for 30 s. Actin was used as a reference control gene. The 2−ΔΔCt method [85] was used to calculate relative gene expression.

Statistical analysis

All measurements were performed using three biological replicates. Statistical analysis was performed with GraphPad Prism 8.0 software, and significance analysis was conducted with one-way ANOVA with multiple comparisons (P < 0.05 or P < 0.01).

Availability of data and materials

The transcriptome datasets supporting the conclusions of this study have been uploaded to the National Center for Biotechnology Information (accession number: PRJNA922734) and can be accessed using the following link https://www.ncbi.nlm.nih.gov/bioproject/?term=PRJNA922734. The other supporting data are included as supplemental files.

Abbreviations

SNP:

Sodium nitroprusside

NO:

Nitric oxide

Fv/Fm:

PSII maximum photochemical conversion efficiency

ABS/RC:

Absorption flux per reaction center

DIo/RC:

Dissipation energy per reaction center

MDA:

Malondialdehyde

H2O2 :

Hydrogen peroxide

SOD:

Superoxide dismutase

POD:

Peroxidase

CAT:

Catalase

NCBI:

National Center for Biotechnology Information

KEGG:

Kyoto Encyclopedia of Genes and Genomes

DEGs:

Differentially expressed genes

TF:

Transcription factor

FPKM:

Fragments per kilobase of transcript per million mapped reads

ROS:

Reactive oxygen species

References

  1. Bog M, Baumbach H, Schween U, Hellwig F, Landolt E, Appenroth K-J. Genetic structure of the genus Lemna L. (Lemnaceae) as revealed by amplified fragment length polymorphism. Planta. 2010;232:609–19.

    Article  CAS  PubMed  Google Scholar 

  2. Li J, Jain M, Vunsh R, Vishnevet sky J, Hanania U, Flaishman M, et al. Callus induction and regeneration in Spirodela and Lemna. Plant Cell Rep. 2004;22:457–64.

    Article  CAS  PubMed  Google Scholar 

  3. Guo L, Jin Y, Xiao Y, Tan L, Tian X, Ding Y, et al. Energy-efficient and environmentally friendly production of starch-rich duckweed biomass using nitrogen-limited cultivation. J Clean Prod. 2020;251(1):119726.

    Article  CAS  Google Scholar 

  4. Zhao YG, Fang Y, Jin YL, Huang J, Bao S, He ZM, et al. Effects of operation parameters on nutrient removal from wastewater and high-protein biomass production in a duckweed-based (Lemma aequinoctialis) pilot-scale system. Water Sci Technol. 2014;70(7):1195–204.

    Article  CAS  PubMed  Google Scholar 

  5. Su HF, Zhao Y, Jiang J, Lu QL, Li Q, Luo Y, et al. Use of duckweed (Landoltia punctata) as a fermentation substrate for the production of higher alcohols as biofuels. Energy Fuels. 2014;28:3206–16.

    Article  CAS  Google Scholar 

  6. Landolt E, Kandeler R. Biosystematic investigations in the family of duckweeds (Lemnaceae) = Biosystematische Untersuchungen in der Familie der Wasserlinsen (Lemnaceae). Biochemie Physiologie der Pflanzen. 1988;183(6):476.

    Google Scholar 

  7. Wang W, Haberer G, Gundlach H, Gläßer C, Nussbaumer T, Luo MC, et al. The Spirodela polyrhiza genome reveals insights into its neotenous reduction fast growth and aquatic lifestyle. Nat Commun. 2014;5:3311.

    Article  ADS  CAS  PubMed  Google Scholar 

  8. Cao HX, Vu GTH, Wang WQ, Appenroth KJ, Messing J, Schubert I. The map-based genome sequence of Spirodela polyrhiza aligned with its chromosomes, a reference for karyotype evolution. New Phytol. 2016;209:354–63.

    Article  CAS  PubMed  Google Scholar 

  9. Vunsh R, Li J, Hanania U, Edelman M, Flaishman M, Perl A, et al. High expression of transgene protein in Spirodela. Plant Cell Rep. 2007;26(9):1511–9.

    Article  CAS  PubMed  Google Scholar 

  10. Rival S, Wisniewski JP, Langlais A, Kaplan H, Freyssinet G, Vancanneyt G, et al. Spirodela (duckweed) as an alternative production system for pharmaceuticals: a case study, aprotinin. Transg Res. 2008;17(4):503–13.

    Article  CAS  Google Scholar 

  11. Yang J, Li G, Hu S, Bishopp A, Heenatigala PPM, Kumar S, et al. A protocol for efficient callus induction and stable transformation of Spirodela polyrhiza (L.) Schleiden using Agrobacterium tumefaciens. Aquatic Bot. 2018;151:80–6.

    Article  Google Scholar 

  12. Baek GY, Saeed M, Choi HK. Duckweeds: their utilization, metabolites and cultivation. Appl Biol Chem. 2021;73:64–73.

    Google Scholar 

  13. Filippou P, Antoniou C, Fotopoulos V. The nitric oxide donor sodium nitroprusside regulates polyamine and proline metabolism in leaves of Medicago truncatula plants. Free Radic Biol Med. 2012;56:172–83.

    Article  PubMed  Google Scholar 

  14. Floryszak-Wieczorek J, Milczarek G, Arasimowicz M, Ciszewski A. Do nitric oxide donors mimic endogenous NO-related response in plants? Planta. 2006;224:1363–72.

    Article  CAS  PubMed  Google Scholar 

  15. Keisham M, Jain P, Singh N, Von TC, Bhatla SC, Lindermayr C. Deciphering the nitric oxide, cyanide and iron-mediated actions of sodium nitroprusside in cotyledons of salt stressed sunflower seedlings. Nitric Oxide. 2019;88:10–26.

    Article  CAS  PubMed  Google Scholar 

  16. Ge Y, Li X, Li C, Tang Q, Duan B, Cheng Y, et al. Effect of sodium nitroprusside on antioxidative enzymes and the phenylpropanoid pathway in blueberry fruit. Food Chem. 2019;295:607–12.

    Article  CAS  PubMed  Google Scholar 

  17. Hayat S, Yadav S, Wani AS, Irfan M, Ahmad A. Nitric Oxide Effects on Photo-synthetic Rate, Growth, and Antioxidant Activity in Tomato. Int J Veg Sci. 2011;17:333–48.

    Article  Google Scholar 

  18. Jiao C, Yang R, Zhou Y, Gu Z. Nitric oxide mediates isoflavone accumulation and the antioxidant system enhancement in soybean sprouts. Food Chem. 2016;204:373–80.

    Article  CAS  PubMed  Google Scholar 

  19. Sun M, Yang XL, Zhu ZP, Xu QY, Wu KX, Kang YJ, et al. Comparative tran-scriptome analysis provides insight into nitric oxide suppressing lignin accumulation of postharvest okra (Abelmoschus esculentus L.) during cold storage. Plant Physiol Biochem. 2021;167:49–67.

    Article  CAS  PubMed  Google Scholar 

  20. Ekinci M, Ors S, Yildirim E, Turan M, Sahin U, Dursun A, et al. Determination of physiological indices and some antioxidant enzymes of chard exposed to nitric oxide under drought stress. Russ J Plant Physiol. 2020;67:740–9.

    Article  CAS  Google Scholar 

  21. Zhao Y, Wei XH, Long Y, Ji XZ. Transcriptional analysis reveals sodium nitroprusside affects alfalfa in response to PEG-induced osmotic stress at germination stage. Protoplasma. 2020;257:1345–58.

    Article  CAS  PubMed  Google Scholar 

  22. Xu MJ, Dong JF. Elicitor-induced nitric oxide burst is essential for triggering catharanthine synthesis in Catharanthus roseus suspension cells. Appl Microbiol Biotechnol. 2005;67:40–4.

    Article  CAS  PubMed  Google Scholar 

  23. Xu MJ, Dong JF, Zhu MY. Nitric oxide mediates the fungal elicitor-induced puerarin biosynthesis in Pueraria thomsonii Benth suspension cells through a salicylic acid (SA)-dependent and a jasmonic acid (JA)-dependent signal pathway. Sci China Series C Life Sci. 2006;49:379–89.

    Article  CAS  Google Scholar 

  24. Zheng LP, Guo YT, Wang JW, Tan RX. Nitric oxide potentiates oligosaccha-ride-induced artemisinin production in Artemisia annua hairy roots. J Integr Plant Biol. 2008;50:49–55.

    Article  CAS  PubMed  Google Scholar 

  25. Hao GP, Du XH, Zhao FX, Renjiu S, Wang JM. Role of nitric oxide in UV-B-induced activation of PAL and stimulation of flavonoid biosynthesis in Ginkgo biloba callus. Plant Cell Tissue Org Cult. 2009;97:175–85.

    Article  CAS  Google Scholar 

  26. Ren YF, Wang W, He JY, Zhang LY, Wei YJ, Yang M. Nitric oxide alleviates salt stress in seed germination and early seedling growth of pakchoi (Brassica chinensis L.) by enhancing physiological and biochemical parameters. Ecotoxicol Environ Saf. 2020;187:109785.

    Article  CAS  PubMed  Google Scholar 

  27. Valderrama R, Chaki M, Begara-Morales JC, Petrivalský M, Barroso JB. Editorial: Nitric Oxide in Plants. Front Plant Sci. 2021;12:705157.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Wani KI, Naeem M, Castroverde CDM, Kalaji HM, Albaqami M, Aftab T. Molecular Mechanisms of Nitric Oxide (NO) Signaling and Reactive Oxygen Species (ROS) Homeostasis during Abiotic Stresses in Plants. Int J Mol Sci. 2021;22:9656.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Liu LP. Effects of exogenous NO donor sodium nitroprusside (SNP) on physiological parameters of Chlorella. J Green Sci Technol. 2016;14:151–3.

    Google Scholar 

  30. Yu C, Zhao X, Qi G, Bai Z, Wang Y, Wang S, et al. Integrated analysis of transcriptome and metabolites reveals an essential role of metabolic flux in starch accumulation under nitrogen starvation in duckweed. Biotechnol Biofuels. 2017;10:167.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Ahn Y, Seong N, Choi HY. Study on the effect of medicinal herbs used as Bu-pyung (S. polyrhiza and L. paucicostata) on immune and anti-cancer. Korea J Herbol. 2004;19:117.

    Google Scholar 

  32. Al-Snafi AE. Lemna minor: Traditional uses, chemical constituents and pharmacological effects-a review. J Pharm. 2019;9:06–11.

    CAS  Google Scholar 

  33. Zhu XJ, Zhang F, Wang XX, Jing HL. Research progress on the pharmacological effects of Phyllostachys edulis. Guid J Tradit Chin Med Pharm. 2020;26(14):29–33.

    Google Scholar 

  34. Xu Q, Yang J, Li Q, Xie M, Qin M. Progress in the pharmacological study of lignan and its glycosides in Mimosa dioica. Chin Wild Plant Resour. 2019;38(04):53–62.

    Google Scholar 

  35. Wan SQ, Liu LY, Liu MS, Huang XL. Study on the pharmacological mechanism of action of Orientin. Jo Med Res. 2018;47(06):183–6.

    Google Scholar 

  36. Cheng J, Huang J. Study on the analgesic and anti-inflammatory immune effects of Mucuna pruriens. Guide China Med. 2016;14(31):29–30.

    Google Scholar 

  37. Gao CX, Huang Y, Sun WJ. Advances in the study of proanthocyanidins in tea. J Tea Sci. 2020;40(04):441–53.

    Google Scholar 

  38. Pérez-Rodríguez P, Riaño-Pachón DM, Corrêa LGG, Rensing SA, Kersten B, Mueller-Roeber B. PlnTFDB: updated content and new features of the plant transcription factor database. Nucleic Acids Res. 2010;38:822–7.

    Article  Google Scholar 

  39. Riechmann JL, Meyerowitz EM. The AP2/EREBP family of plant transcription factors. Biol Chem. 1998;379(6):633–46.

    CAS  PubMed  Google Scholar 

  40. El-Sharkawy I, Sherif S, Mila I, Bouzayen M, Jayasankar S. Molecular characterization of seven genes encoding ethylene-responsive transcriptional factors during plum fruit development and ripening. J Exp Bot. 2009;60:907–22.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Zhu YR, Li XX, Gao X, Sun JQ, Ji XY, Feng GD, et al. Molecular mechanism underlying the effect of maleic hydrazide treatment on starch accumulation in S. polyrhiza 7498 fronds. Biotechnol Biofuels Bioprod. 2021;14:99.

    Article  CAS  Google Scholar 

  42. Liu CX, Sun ZY, Luo YB, Zhu HL, Qu GQ. Progress in the study of phosphorylation of bHLH transcription factors regulating plant physiological functions. Biotechnol Bull. 2023;39(03):26–34.

    Google Scholar 

  43. Buti S, Hayes S, Pierik R. The bHLH network underlying plant shade-avoidance. Physiologia Plantarum. 2020;169:312–24.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Xu NT, Liu SA, Lu ZG, Pang SY, Wang L, Wang L, et al. Gene expression profiles and flavonoid accumulation during salt stress in Ginkgo biloba seedlings. Plants. 2020;9:1162.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Shi WY, Du YT, Ma J, Min DH, Jin LG, Chen J, et al. The WRKY transcription factor GmWRKY12 Confers drought and salt tolerance in soybean. Int J Mol Sci. 2018;19(12):4087.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Hedhly A, Vogler H, Schmid MW, Pazmino D, Gagliardini V, Santelia D, et al. Starch turnover and metabolism during flower and early embryo development. Plant Physiol. 2016;172:2388–402.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Zeng FS, Sun FK, Li LL, Liu K, Zhan YG. Genome-scale transcriptome analysis in re-sponse to nitric oxide in birch cells: implications of the triterpene biosynthetic pathway. Plos One. 2014;9(16):116157.

    Article  ADS  Google Scholar 

  48. Gill SS, Tuteja N. Reactive oxygen species and antioxidant machinery in abiotic stress tolerance in crop plants. Plant Physiol Biochem. 2010;48:909–30.

    Article  CAS  PubMed  Google Scholar 

  49. Hahlbrock K, Grisebach H. Enzymic controls in the biosynthesis of lignin and flavonoids. Ann Rev Plant Physiol. 1979;30(1):105–30.

    Article  CAS  Google Scholar 

  50. Li X, Zhang L, Ahammed GJ, Li ZX, Wei JP, Shen C, et al. Nitric oxide mediates brassinosteroid-induced flavonoid biosynthesis in Camellia sinensis L. J Plant Physiol. 2017;214:145–51.

    Article  CAS  PubMed  Google Scholar 

  51. Ma SY, Bao JY, Lu YQ, Lu X, Tian P, Zhang XL, et al. Glucoraphanin and sulforaphane biosynthesis by melatonin mediating nitric oxide in hairy roots of broccoli (Brassica oleracea L. var. italica Planch): insights from transcriptome data. BMC Plant Biol. 2022;22:403.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Pokora W, Tułodziecki S, Dettlaff-Pokora A, Aksmann A. Cross talk between hydrogen peroxide and nitric oxide in the unicellular green algae cell cycle: how does it work? Cells. 2022;11:2425.

    Article  PubMed  PubMed Central  Google Scholar 

  53. García-Mata C, Gay R, Sokolovski S, Hills A, Lamattina L, Blatt MR. Nitric oxide regulates K+ and Cl- channels in guard cells through a subset of abscisic acid-evoked signaling pathways. Natl Acad Sci. 2003;100:11116–21.

    Article  ADS  Google Scholar 

  54. Sokolovski S, Hills A, Gay R, Garcia-Mata C, Lamattina L, Blatt MR. Protein phosphorylation is a prerequisite for intracellular Ca2+ release and ion channel con-trol by nitric oxide and abscisic acid in guard cells. Plant J. 2005;43:520–9.

    Article  CAS  PubMed  Google Scholar 

  55. Deng SX, Ma J, Zhang LL, Chen FJ, Sang ZY, Jia ZK, et al. De novo transcriptome sequencing and gene expression profiling of Magnolia wufengensis in response to cold stress. BMC Plant Biol. 2019;19:321.

    Article  PubMed  PubMed Central  Google Scholar 

  56. Fromme P, Melkozernov A, Jordan P, Krauss N. Structure and function of photosystem I: interaction with its soluble electron carriers and external antenna systems. FEBS Lett. 2003;555(1):40–4.

    Article  CAS  PubMed  Google Scholar 

  57. Liu Y, Jiang ZK, Ye YT, Wang DH, Jin SH. Enhanced salt tolerance of Torreya grandis genders is related to nitric oxide level and antioxidant capacity. Front Plant Sci. 2022;13:906071.

    Article  PubMed  PubMed Central  Google Scholar 

  58. Jabeen Z, Fayyaz HA, Irshad F, Hussain N, Hassan MN, Li J, et al. Sodium nitro-prusside application improves morphological and physiological attributes of soybean (Glycine max L.) under salinity stress. Plos One. 2021;16:0248207.

    Article  Google Scholar 

  59. León J, Costa Á, Castillo M. Nitric oxide triggers a transient metabolic reprogramming in Arabidopsis. Sci Rep. 2016;6:37945.

    Article  ADS  PubMed  PubMed Central  Google Scholar 

  60. Iqbal N, Sehar Z, Fatma M, Umar S, Sofo A, Khan NA. Nitric oxide and abscisic acid mediate heat stress tolerance through regulation of osmolytes and antioxidants to protect photosynthesis and growth in wheat plants. Antioxidants. 2022;11:372.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  61. Alp K, Terzi H, Yildiz M. Proteomic and physiological analyses to elucidate nitric oxide-mediated adaptive responses of barley under cadmium stress. Physiol Mol Biol Plants. 2022;28:1467–76.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Pashkovskiy PP, Soshinkova TN, Korolkova DV, Kartashov AV, Zlobin IE, Lyubimov VY, et al. The effect of light quality on the pro-/antioxidant balance, activity of photosystem II, and expression of light-dependent genes in Eutrema salsugineum callus cells. Photosynth Res. 2018;136:199–214.

    Article  CAS  PubMed  Google Scholar 

  63. Baker NR. Chlorophyll fluorescence: a probe of photosynthesis in Vivo. Ann Rev Plant Biol. 2008;59:89–113.

    Article  CAS  Google Scholar 

  64. Rochaix JD. Regulation and dynamics of the light-harvesting system. Ann Rev Plant Biol. 2014;65:287–309.

    Article  CAS  Google Scholar 

  65. Adam BS, Frolow F, Nelson N. Light-harvesting features revealed by the structure of plant photosystem I. Photosynth Res. 2004;81:239–50.

    Article  Google Scholar 

  66. Wang W, Messing J. Analysis of ADP-glucose pyrophosphorylase expression during turion formation induced by abscisic acid in Spirodela polyrhiza (greater duckweed). BMC Plant Biol. 2012;12:5.

    Article  PubMed  PubMed Central  Google Scholar 

  67. Tao X, Fang Y, Xiao Y, Jin YL, Ma XR, Zhao Y, et al. Comparative transcriptome analysis to investigate the high starch accumulation of duckweed (Landoltia punctata) under nutrient starvation. Biotechnol Biofuels Bioprod. 2013;6:72.

    Article  CAS  Google Scholar 

  68. Tian XP, Fang Y, Jin YL, Yi ZL, Li J, Du A, et al. Ammonium detoxification mechanism of ammonium-tolerant duckweed (Landoltia punctata) revealed by carbon and nitrogen metabolism under ammonium stress. Environ Pollut. 2021;277(1):116834.

    Article  CAS  PubMed  Google Scholar 

  69. Wang L, Li SL, Sun L, Tong TN, Yang L, Zhu YR, et al. Over-Expression of Phosphoserine Aminotransferase-Encoding Gene (AtPSAT1) Prompts Starch Accumulation in L. turionifera under Nitrogen Starvation. Int J Mol Sci. 2022;23(19):11563.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  70. Yang R, Lin XC, Dou Y, Zhang W, Du H, Wan C, et al. Transcriptome profiling of postharvest kiwifruit in response to exogenous nitric oxide. Scientia Horticulturae. 2021;277:109788.

    Article  CAS  Google Scholar 

  71. Farmer EE, Mueller MJ. ROS-mediated lipid peroxidation and RES-activated signaling. Ann Rev Plant Biol. 2013;64:429–550.

    Article  CAS  Google Scholar 

  72. Gupta DK, Palma JM, Corpas FJ. Antioxidants and antioxidant enzymes in higher plants. 2018.

    Book  Google Scholar 

  73. Chang IH, Cheng KT, Huang PC, Lin YY, Cheng LJ, Cheng TS. Oxidative stress in greater duckweed (Spirodela polyrhiza) caused by long-term NaCl exposure. Acta Physiologiae Plantarum. 2012;34(3):1165–76.

    Article  CAS  Google Scholar 

  74. Cheng T. NaCl-induced responses in giant duckweed (Spirodela polyrhiza). J Aquatic Plant Manag. 2011;2:62–71.

    ADS  Google Scholar 

  75. Tao X, Fang Y, Huang MJ, Xiao Y, Liu Y, Ma XR, et al. High flavonoid accompanied with high starch accumulation triggered by nutrient starvation in bioenergy crop duckweed (Landoltia punctata). BMC Genomics. 2017;18(1):166.

    Article  PubMed  PubMed Central  Google Scholar 

  76. Liu Y, Fang Y, Huang M, Jin Y, Sun J, Tao X, et al. Uniconazole-induced starch accumulation in the bioenergy crop duckweed (Landoltia punctata) I :transcriptome analysis of the effects of uniconazole on chlorophyll and endogenous hormone biosynthesis. Biotechnol Biofuels. 2015;8:57.

    Article  PubMed  PubMed Central  Google Scholar 

  77. Huang JH, Li JJ, Zheng SB, Wu JZ, Zhang W, Sun T, et al. Epimedium flavonoids counteract the side effects of glucocorticoids on hypothalamic-pituitary-adrenal axis. Evid-Based Complement Altern Med. 2013;938425. https://doi.org/10.1155/2013/938425.

  78. Wang Y, Kandeler R. Promotion of flowering by a tumor promoter. J Plant Physiol. 1994;144:710–3.

    Article  CAS  Google Scholar 

  79. Arnon DI. Copper enzymes in isolated chloroplasts polyphenoloxidase in beta vulgaris. Plant Physiol. 1949;24:1–15.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  80. Luthar Z, Kreft I. Influence of temperature on tannin content in different ripening phases of buckwheat (Fagopyrum esculentum Moench) seeds. Fagopyrum. 1999;16:61–5.

    Google Scholar 

  81. Grabherr M, Haas B, Yassour M, Levin J, Thompson D, Amit I, et al. Full-Length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29:644–52.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  82. Kanehisa M, Goto S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28:27–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  83. Kanehisa M. Toward understanding the origin and evolution of cellular organisms. Protein Sci. 2019;28:1947–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  84. Kanehisa M, Furumichi M, Sato Y, Kawashima M, Ishiguro-Watanabe M. KEGG for taxonomy-based analysis of pathways and genomes. Nucleic Acids Res. 2023;51:587–92.

    Article  Google Scholar 

  85. Fu L, Ding Z, Han B, Hu W, Li Y, Zhang J. Physiological Investigation and Transcriptome Analysis of Polyethylene Glycol (PEG)-Induced Dehydration Stress in Cassava. Int J Mol Sci. 2016;17(3):283.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We are grateful to Guangzhou Genedenovo Biotechnology Co., Ltd for assisting in sequencing.

Funding

This work was financially supported by National Natural Science Foundation of China (82074120) for the design of the study and the collection, analyses, and interpretation of data.

Author information

Authors and Affiliations

Authors

Contributions

BBX contributed to the experimental design and manuscript review. YH organized the experiment process. YMZ modified the picture and wrote the paper. RJ, TYH, YW and WJW participated in the experiments. YW, YRZ and LY also reviewed this manuscript All authors have read and approved the final manuscript.

Corresponding author

Correspondence to Beibei Xiang.

Ethics declarations

Ethics approval and consent to participate

Our research did not involve any human or animal subjects, material, or data. We declare that the plant material in the experiment was collected and studied in accordance with relevant institutional, national, and international guidelines and legislation.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1:

 Figure S1. Phenotypic changes of S. polyrrhiza after SNP treatment.

Additional file 2:

 Figure S2. Species distribution of Nr annotation.

Additional file 3:

 Figure S3. The relative expression levels of seventeen selected DEGs were compared by RNA-seq and qRT-PCR.

Additional file 4

Table S1. The DEGs were identified as transcription factor.

Additional file 5:

 Table S2. The DEGs were identified as photosynthesis.

Additional file 6:

 Table S3. The DEGs were identified as starch metabolism.

Additional file 7:

 Table S4. The DEGs were identified as antioxidant.

Additional file 8:

 Table S5. The DEGs were identified as flavonoids biosynthesis pathways.

Additional file 9:

 Table S6. DEGs involved in quantitative real-time PCR validation.

Additional file 10:

 Table S7. DEGs involved in NO biosynthetic and signal transduction.

Additional file 11:

 Table S8. Primers used for QRT-PCR analysis.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhang, Y., Jia, R., Hui, T. et al. Transcriptomic and physiological analysis of the response of Spirodela polyrrhiza to sodium nitroprusside. BMC Plant Biol 24, 95 (2024). https://doi.org/10.1186/s12870-024-04766-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12870-024-04766-6

Keywords