Skip to main content
  • Research article
  • Open access
  • Published:

High flavonoid accompanied with high starch accumulation triggered by nutrient starvation in bioenergy crop duckweed (Landoltia punctata)

Abstract

Background

As the fastest growing plant, duckweed can thrive on anthropogenic wastewater. The purple-backed duckweed, Landoltia punctata, is rich in starch and flavonoids. However, the molecular biological basis of high flavonoid and low lignin content remains largely unknown, as does the best method to combine nutrients removed from sewage and the utilization value improvement of duckweed biomass.

Results

A combined omics study was performed to investigate the biosynthesis of flavonoid and the metabolic flux changes in L. punctata grown in different culture medium. Phenylalanine metabolism related transcripts were identified and carefully analyzed. Expression quantification results showed that most of the flavonoid biosynthetic transcripts were relatively highly expressed, while most lignin-related transcripts were poorly expressed or failed to be detected by iTRAQ based proteomic analyses. This explains why duckweed has a much lower lignin percentage and higher flavonoid content than most other plants. Growing in distilled water, expression of most flavonoid-related transcripts were increased, while most were decreased in uniconazole treated L. punctata (1/6 × Hoagland + 800 mg•L-1 uniconazole). When L. punctata was cultivated in full nutrient medium (1/6 × Hoagland), more than half of these transcripts were increased, however others were suppressed. Metabolome results showed that a total of 20 flavonoid compounds were separated by HPLC in L. punctata grown in uniconazole and full nutrient medium. The quantities of all 20 compounds were decreased by uniconazole, while 11 were increased and 6 decreased when grown in full nutrient medium. Nutrient starvation resulted in an obvious purple accumulation on the underside of each frond.

Conclusions

The high flavonoid and low lignin content of L. punctata appears to be predominantly caused by the flavonoid-directed metabolic flux. Nutrient starvation is the best option to obtain high starch and flavonoid accumulation simultaneously in a short time for biofuels fermentation and natural products isolation.

Background

Flavonoids, also known as vitamin P, constitute a vast class of secondary metabolites widely distributed in plants, which encompasses more than 10,000 structures [1]. They have a low molecular weight and a general structure of three rings, including two phenyl rings (A and B) and a heterocyclic ring (C). With different substituent groups, flavonoids can be divided into seven subgroups, including chalcones, flavones, flavonols, flavandiols, anthocyanins, condensed tannins and aurones [2]. Some specialized forms of flavonoids can be synthesized by some plant species, such as the isoflavonoids [3] and 3-deoxyanthocyanins [4, 5]. Different flavonoids usually play various roles in plants by regulating several developmental processes [610]. Furthermore, these secondary metabolites are well characterized as defense compounds and signaling molecules that can withstand a wide array of environmental stresses in plants and diseases in humans [11, 12], due to their capacity to absorb ultraviolet (UV) radiation, and inhibiting the generation of reactive oxygen species (ROS) [1315]. Through their ability to inhibit DNA gyrase, energy metabolism and cytoplasmic membrane function, flavonoids possess antifungal, antiviral and antibacterial activity [16, 17].

Duckweed (Lemnacecae family) is the smallest and simplest flowering aquatic plant in the world, and its growth highly adaptable across a broad range of climates [18]. It has a long yearly production period with an almost exponential growth rate, producing biomass faster than most other plants. It can thrive on eutrophic wastewater, through its ability to remove nutrients from sewage [19] and large amounts of CO2 from the atmosphere [2022]. In warm seasons, duckweed can remove up to 85% of total Kjehldahl nitrogen (TKN) and 78% of total phosphorous (TP) from sewage [22]. The value of duckweed as a test species for the registration of agrochemicals has been discussed worldwide [23]. A previous study indicated that this plant possesses negligible lignin content [24]. Depending on the duckweed species and the growing conditions, the starch content of duckweed ranges from 3% to 75% [2527]. Furthermore, it has been found that the purple-backed duckweed has high flavonoid content in crude plant form [28, 29], while only flavonoid-rich fractions of the most prevalent flavonoid sources, tartary buckwheat or ginkgo, can be used to extract these kinds of flavonoids [3032]. Together with a much higher biomass production, duckweed should be a more promising flavonoid resource plant than tartary buckwheat and gingko. These characteristics make purple-backed duckweed a potential sustainable source for bioenergy production [33, 34], animal feed [35] and even human food [36]. A company from Israel has found that duckweed can address the challenges of consumer health concerns, rising health care costs and food security issues by exploring the nutritional value, traditional consumption in Southeast Asia, and commercialization possibilities of duckweed [36]. However, a method to combine nutrients remove from sewage and the utilization value improvement of duckweed biomass remains as yet unknown.

In a previous study, it was found that the total flavonoid content of L. punctata increased from 4.51% to 5.56% following nutrient starvation for 168 h [28], accompanied by high starch accumulation for bioethanol fermentation [25]. Spraying with 800 mg · L-1 uniconazole is an alternative method to accumulate high levels of starch [26, 27], but whether it underwent the same physiological and molecular alteration remains unknown. In this study, the changes of flavonoids in full nutrient, starvation and uniconazole treated L. punctata groups were investigated and compared by a combined omics study. This provided molecular support for the simultaneous accumulation of high starch and high flavonoid levels in this bio-resource plant.

Results

Comprehensive transcriptome construction for L. punctata

In order to construct a comprehensive transcriptome for L. punctata, Illumina HiSeq 2000 paired-end (PE) reads of nutrient starvation (distilled water, NS) [25] and uniconazole (1/6 × Hoagland + 800 mg•L-1 uniconazole, UT) [26, 27] responsive transcriptomes, and also the full nutrient (1/6 × Hoagland, FN) transcriptome, were pooled together and de novo assembled using Trinity (v2012-06-08) [37]. All PE reads were deposited in Sequence Read Archive database (SRA) under accession number of PRJNA185389. A total of 543,912,936 PE 90 bp reads were obtained from the three RNA-Seq groups, corresponding to 48.95 Gbp in total (Table 1). Furthermore, 155,903 contigs with lengths ≥200 bp were assembled, corresponding to a transcriptome size of 170.34 Mb. The average length, N50 length and max length was 1093 bp, 2190 bp and 17,234 bp, respectively. Among these contigs, 51,873 were longer than 1000 bp and 26,931 were longer than 2000 bp. The results from scanning of the Open Reading Frames (ORFs) of all contigs showed that there were 67,061 ORFs with lengths ≥600 bp (from ATG to stop codon), and 37,797 ORFs with lengths ≥900 bp.

Table 1 RNA-Seq statistics for different duckweed samples

All PE reads were used separately for short-read alignment for each sample through the perl script provided with the Trinity package (v2012-06-08) [37]. The number of aligned reads for each contig was counted and used for expression profiling. To normalize the bias introduced by the sequencing library size and mRNA composition, edgeR (the Empirical analysis of Digital Gene Expression in R) [38] in the Trinity package (v2012-06-08) [37] was used to make an effective library size for each sample and normalize the number of aligned reads per transcript to generate a FPKM (Fragments Per Kilobase of transcripts per Million mapped fragments) value using the RESM-based algorithm. It was found that the number of expressed transcripts ranged from 36,950 to 60,854, with only 20,776 transcripts expressed in all samples (Table 1). Furthermore, the results showed that full nutrient conditions stimulated more transcript expression compared to the other two experimental groups.

A BlastX sequence similarity search against the non-redundant protein database (NR) of NCBI [http://www.ncbi.nlm.nih.gov/] was conducted by a locally-installed blast program to investigate functional annotation of each contig. BlastX results were then uploaded to the Blast2GO platform [39, 40] for annotation. A total of 98,106 (62.9%) contigs (transcripts) had significant BlastX hits. Of the 26,931 contigs ≥2000 bp in length, 26,273 were annotated, corresponding to an annotation rate of 97.6%. For the 51,873 contigs ≥1000 bp, the annotation rate was 93.0% (annotated 48,266). For contigs ≥900 bp and ≥600 bp, this rate was 91.7% (50,762 of 55,328) and 85.3% (59,595 of 69,878), respectively.

Biosynthetic network of phenylalanine metabolism

To construct the biosynthetic network related to phenylalanine metabolism, enzyme codes were extracted and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways retrieved from the KEGG web server [http://www.genome.jp/kegg/]. Transcripts were detected that corresponded to almost all of the enzymes involved in flavonoid and lignin biosynthesis, except for flavone synthase, aureusidin synthase, flavanone 7-O-beta-glucosyltransferase and flavanone 7-O-glucoside 2''-O-beta-L-rhamnosyltransferase (Fig. 1, Additional file 1: Table S1, Additional file 2: Table S2). Phenylalanine ammonialyase (EC: 4.3.1.24, PAL), cinnamate 4-hydroxylase (EC: 1.14.13.11, C4H) and 4-hydroxycinnamoyl-CoA ligase (EC: 6.2.1.12, 4CL) are the universal factors involved in flavonoid and lignin biosynthesis [41]. Transcripts encoding these enzymes were highly expressed with abundance higher than 100 FPKM. Cinnamoyl-CoA reductase (EC: 1.2.1.44, CCR) and hydroxycinnamoyl transferase (EC: 2.3.1.133, HCT) catalyze the initial reaction in the lignin biosynthesis branch. Almost all transcripts related to CCR and HCT had expression levels lower than 20 FPKM. Conversely, chalcone synthase (EC: 2.3.1.74, CHS) and chalcone isomerase (EC: 5.5.1.6, CHI), the enzymes that catalyze the first two reactions of flavonoid biosynthesis branch, were more highly expressed (Fig. 1). This possibly explains why duckweed has a much lower lignin percentage than most other plants [24, 34, 42, 43]. As p-coumaroyl CoA is the product of a 4CL or C4H catalyzed reaction, it can be converted into isoliquiritigenin or naringenin chalcone, then catalyzed by chalcone isomerase (EC: 5.5.1.6, CHI) to feed this product into the isoflavonoid biosynthesis pathway. Expression data showed that CHS and CHI were both highly expressed at over 100 FPKM, which may indicate that a large amount of p-coumaroyl CoA was directed into the isoflavonoid biosynthesis pathway. Meanwhile, the expression of flavanone 3-hydroxylase (EC: 1.14.11.9, F3H), dihydroflavonol 4-reductase (EC: 1.1.1.234, DFR) and anthocyanidin synthase (EC: 1.14.11.19, ANS), which direct the metabolic flux into the anthocyanin biosynthesis branch, were all higher than 100 FPKM. However, flavone synthase (EC: 1.14.11.22, FNS) was unable to be detected in this study. Of course, this may be because there was no expression of FNS.

Fig. 1
figure 1

Phenylalanine metabolism networks in L. punctata. The abbreviations correspond to enzymes involved in phenylalanine metabolic networks. Different colors represent different expression levels. PAL: phenylalanine ammonialyase, EC: 4.3.1.24. C4H: cinnamate 4-hydroxylase, EC: 1.14.13.11. 4CL: 4-coumarate-CoA ligase, EC: 6.2.1.12. HCT: hydroxycinnamoyl transferase, EC: 2.3.1.133. C3H: 4-coumarate 3-hydroxylase, EC: 1.14.14.9. CCoAOMT: caffeoyl-CoA O-methyl transferase, EC: 2.1.1.104. COMT: caffeic acid o-methyl transferase, EC: 2.1.1.68. F5H: ferulate 5-hydroxylase, EC:1.14.-.-. CCR: cinnamoyl-CoA reductase, EC: 1.2.1.44; CAD: cinnamyl-alcohol dehydrogenase, EC: 1.1.1.195; LACC: laccase, EC: 1.10. 3.2. CHS: chalcone synthase, EC: 2.3.1.74. CHI: chalcone isomerase, EC: 5.5.1.6. F3H: flavanone 3-hydroxylase, EC: 1.14.11.9. FLS: flavonol synthase EC:1.14.11.23. DFR: dihydroflavonol 4-reductase, EC: 1.1.1.234. F3’H: flavonoid 3'-hydroxylase, EC: 1.14.13.21. F3’5’H: EC: 1.14.13.88. FNS: flavone synthase, EC:1.14.11.22. ANS: anthocyanidin synthase, EC: 1.14.11.19. ANR: anthocyanidin reductase, EC:1.3.1.77. LAR: leucoanthocyanidin reductase, EC:1.17.1.3. AS1: aureusidin synthase, EC:1.21.3.6. The bold arrows show the main metabolic flux

RNA-Seq based flavonoid biosynthetic analyses of nutrient starvation or uniconazole treated L. punctata

Expression patterns of genes involved in specific pathways can affect the metabolic flux. All transcripts described above were quantified by RNA-Seq analyses (Fig. 2, Additional file 1: Table S1, Additional file 2: Table S2). Culturing duckweed in distilled water, the highest expressed PAL, comp39767_c0_seq1, showed no obvious change (132.07, 140.13 and 113.28 FPKM). The other highly expressed PAL, comp39767_c0_seq2, was increased in NS-2 (29.62 FPKM) and NS-24 (39.36 FPKM) when compared with that in NS-0 (17.70 FPKM). Expression level of comp46865_c0_seq1, a C4H encoding transcripts, was increased from 78.22 FPKM to 152.93 FPKM in 24 h. The expression level of comp46833_c0_seq1, a 4CL encoding transcript, was increased 2.22 times in NS-24 (294.61 FPKM) compared with that in NS-0 (132.43 FPKM). Spraying with 800 mg · L-1 uniconazole resulted in a slight initial increase in the expression of the two highly expressed PAL, comp39767_c0_seq1 and comp39767_c0_seq2, followed by a slight decrease after 72 h. The highest expressed C4H, comp46865_c0_seq1, was increased in UT-72 (197.19 FPKM) and UT-240 (189.65 FPKM) when compared with that in UT-0 (102.44 FPKM), UT-2 (105.64 FPKM) and UT-5 (108.16 FPKM). Moreover, two other poorly expressed C4H, comp45597_c0_seq1 and comp45597_c1_seq2, were increased during the first 240 h in the uniconazole treated group. However, the results also showed that the highest expressed 4CL encoding transcripts (comp46833_c0_seq1, comp46135_c0_seq2) were obviously suppressed by uniconazole. In addition, several other 4CL were also suppressed, including comp44243_c0_seq10, comp44243_c0_seq9, comp35216_c0_seq1 and comp35216_c1_seq1. 4-coumaroyl-CoA produced by 4CL can be catalyzed by HCT to produce 4-coumaroylshikimate or catalyzed by CHS to be converted into naringenin chalcone, which directs the substrate into lignin biosynthesis or flavonoid biosynthesis. Enzymes involved in the flavonoid metabolism branch were carefully analyzed. It was found that a CHS encoding transcript (comp46654_c0_seq1) was increased by starvation treatment from the original 401.10 FPKM to 500.20 FPKM at 2 h and 519.48 FPKM at 24 h. Similarly, comp46654_c0_seq1 showed a comparable upward tendency in the uniconazole treated group. However, its expression level was much lower than in NS-0, NS-2 and NS-24. Furthermore, expression level of several CHI, FLS, F3H, F3’H, ANS, ANR, LAR, IFR and GT4 encoding transcripts were increased by nutrient starvation. Conversely, the highest expressed DFR transcript (comp36170_c1_seq2) was suppressed (176.81, 118.47 and 36.75 FPKM in NS-0, NS-2 and NS-24), while several other lowly expressed DFR transcripts were increased (comp19960_c0_seq2, comp44884_c0_seq1, comp44884_c0_seq2 and comp44884_c0_seq3). In the uniconazole treated group, the highly expressed FLS, comp48210_c0_seq1, was increased in the first 5 h, but suppressed thereafter by uniconazole. Meanwhile, three F3’H encoding transcripts (comp43561_c0_seq1, comp43561_c0_seq3, comp43561_c0_seq4) were also suppressed by this plant growth regulator. Upon cultivation of L. punctata in full nutrient (1/6Hoagland) solution, the number of increased transcripts was larger than that of the suppressed transcripts. For example, comp39767_c0_seq4, comp39767_c0_seq6, comp46865_c0_seq1, comp44243_c0_seq10, comp44483_c0_seq2, comp35216_c0_seq1, comp35216_c1_seq1, comp35307_c1_seq1, comp46833_c0_seq1, comp18376_c0_seq1, comp41458_c0_seq1, comp46286_c0_seq1, comp46286_c0_seq2, comp19956_c0_seq1, comp43561_c0_seq6, comp46246_c0_seq1 were increased in full nutrient conditions. These results suggest that these three different treatments may trigger different molecular responses.

Fig. 2
figure 2

Expression changes of transcripts related to flavonoid biosynthesis based on RNA-Seq. A heatmap was drawn by HemI toolkit using log2FC values [75]. Most abbreviations correspond to the enzymes listed in Fig. 1. Transcripts with extremely low expression levels are not shown in this figure

iTRAQ based flavonoids biosynthetic analyses of L. punctata treated with nutrient starvation or uniconazole

RNA-Seq study provides a global expression pattern to reveal mRNA composition, but it cannot reveal information about the proteome. As the newest developed quantitative technology, iTRAQ is widely used for proteome characterization. In this study, iTRAQ data of previous studies [28] was re-analyzed using the transcriptome described above as a reference database. The abundance of the most detected flavonoid related proteins, including the PAL, C4H, 4CL, CHS, CHI, F3H and ANS (Fig. 3, Additional file 3: Table S3), was detected to be improved in L. punctata when grown in distilled water. CHS (comp46654_c0_seq1), the first pivotal catalyzer for the flavonoid biosynthetic branch, was increased 6.92 times at 2 h, and 4.37, 8.71 and 15.28 times at 5, 24 and 72 h, respectively, when compared with that of 0 h. Conversely, the first key enzyme of the other branch of the phenylalanine metabolic networks, the lignin biosynthetic branch, failed to be detected in starvation treated samples. Although 25 laccase transcripts were assembled by our de novo RNA-Seq study, none of the proteins corresponding to these transcripts were detected by iTRAQ. F3H, DFR and ANS are the enzymes directing metabolic flux into the anthocyanin biosynthesis branch. In starvation treated L. punctata, the expression of F3H (comp47195_c0_seq1) was increased by 1.14, 1.04, 1.42 and 3.73 times at 2, 5, 24 and 72 h, respectively, when compared with that of 0 h, whereas DFR (comp44884_c0_seq2) was increased by 0.83, 0.83, 0.97 and 1.32, and 1.33, 1.00, 1.39 and 3.13 for ANS (comp18732_c0_seq1), respectively. These results may indicate that the metabolic flux was regulated by nutrient starvation to direct more substrates toward the anthocyanin biosynthesis branch. Spraying with uniconazole resulted in the levels of most flavonoid related proteins being decreased. Three universal factors of phenylalanine metabolism, PAL (comp39767_c0_seq1), C4H (comp46865_c0_seq1) and 4CL (comp35307_c0_seq1), were suppressed immediately (at 2 h) by uniconazole. Although the expression of C4H was increased thereafter, the 4CL was suppressed in all of the four UT samples. The other 4CL, comp46135_c0_seq2, was also suppressed at 72 h and 240 h. Furthermore, CHI (comp41458_c0_seq1), IFR (comp46286_c0_seq1, comp46286_c0_seq2), COMT (comp44590_c0_seq3, comp44487_c0_seq1), 3GT (comp19783_c0_seq1, comp40167_c0_seq4), UF3GT (comp45471_c0_seq2), and 5AT (anthocyanin 5-aromatic acyltransferase, comp45766_c0_seq1) showed a downward trend by the application of uniconazole.

Fig. 3
figure 3

Expression changes of proteins involved in flavonoid biosynthesis based on iTRAQ. A heatmap was drawn by HemI toolkit according to log2FC values [75]. Most abbreviations correspond to the enzymes listed in Fig. 1

Expression of lignin biosynthesis related genes

Lignin provides mechanical support for plant growth, but is not necessary in duckweed, which floats on water surfaces. It has been found that most of the lignin biosynthesis related genes in L. punctata had lower expression than that seen in the flavonoid biosynthesis pathway, and the last key gene involved in lignin biosynthesis was only poorly expressed [25]. These results were supported by an iTRAQ study of the same test samples [28]. When treated with 800 mg · L-1 uniconazole, the first enzyme of the lignin biosynthetic branch, HCT, had limited expression (Additional file 1: Table S1, Additional file 2: Table S2). The highest expression level recorded, for HCT, was 14.64 FPKM in UT-5. For the other enzymes in this pathway, only three transcripts had expression levels ≥100 FPKM, including comp46970_c0_seq1 (CCoAOMT), comp44590_c0_seq3 (COMT) and comp36380_c1_seq1 (F5H). 25 laccase (LACC, EC: 1.10.3.2) encoding transcripts were assembled and all had lower expression levels than 10.00 FPKM (Additional file 1: Table S1, Additional file 2: Table S2). The highest expressed transcript exhibited an expression abundance of 2.97, 5.10, 6.76, 4.95 and 2.75 FPKM in five samples, while expression levels of the other transcripts were all lower than 1.5 FPKM. iTRAQ proteomics profiling results also strongly supported this view. When treated with nutrient starvation, almost all key enzymes involved in lignin biosynthesis were not detected [28]. When exposed to uniconazole and full nutrient medium, most of these assembled transcripts were not detected either.

Expression of flavone, flavonol, isoflavonoid and anthocyanin biosynthesis involved genes

The results described above suggest that the metabolic flux may be primarily directed to the isoflavonoid or anthocyanin biosynthesis branches in L. punctata (Fig. 1). To verify this, enzyme encoding genes involved in flavone, flavonol, isoflavonoid and anthocyanin biosynthesis were carefully analyzed. It was found that almost all genes involved in isoflavonoid biosynthesis, or the flavone and flavonol biosynthesis pathway were not detected in the transcriptome. Despite most anthocyanin biosynthesis related genes failing to be identified, most detected transcripts were increased by nutrient starvation and uniconazole (Additional file 1: Table S1, Additional file 2: Table S2). Five UDP-glucose flavonoid 3-O-glucosyltransferase (EC: 2.4.1.51, UF3GT) encoding sequences, including comp44420_c0_seq1, comp44420_c0_seq2, comp44420_c0_seq3, comp38450_c0_seq2 and comp47570_c0_seq1, were identified. Four of them (comp44420_c0_seq2, comp44420_c0_seq3, comp38450_c0_seq2 and comp47570_c0_seq1) were increased by nutrient starvation, while being slightly increased by uniconazole treatment and suppressed by full nutrient treatment. Moreover, expression levels of two anthocyanidin 3-o-glucosyltransferase (GT1, EC: 2.4.1.115, comp46472_c0_seq3, comp46472_c0_seq14) genes were also increased by the first two treatments described above, but decreased by full nutrient treatment. Whereas 5AT (comp19598_c1_seq1) was only increased by uniconazole, anthocyanin 3-o-beta-glucosyltransferase (3GT, 2.4.1.238, comp33618_c0_seq1) was increased by starvation. These observations support the hypothesis that the metabolic flux was mainly directed into the anthocyanin biosynthesis branch and not the others.

Flavonoid content of uniconazole and full nutrient treated L. punctata

In a previous study, it was found that the total flavonoid content of L. punctata increased from 4.51% to 5.56% during nutrient starvation for 168 h, of which seven of the 17 components showed an obvious increase [28]. Growing L. punctata under natural conditions, the same number of flavonoid compounds was separated by spectroscopic, chemical and biochemical methods, and four of these were identified as new flavonoids in duckweed [29]. However, whether these 17 flavonoid compounds are the same as those observed in the starvation or uniconazole treated L. punctata has not been verified. In this study, flavonoids were extracted and characterized from uniconazole treated L. punctata following the protocol described in the study of Wang, et al. [29]. The results showed that a total of 20 compounds were separated, including the additional compounds 1, 9, 14, 15, 16, 17, 18, 19 and 20 that did not separate in the previous study (Fig. 4). In contrast to the starvation treated L. punctata, all of the 20 compounds were decreased by uniconazole treatment. In addition, several compounds (compound 1, 3, 4, 5, 6, 7, 8, 9, 13, 14 and 17) were increased and several were decreased (compound 11, 12, 15, 18, 19 and 20) by full nutrient treatment. The total flavonoid content of entire plants changed from 2.83% to 0.94% and 3.37% in uniconazole and full nutrient treated L. punctata, respectively. Furthermore, it was found that purple coloration accumulated on the frond underside in starved L. punctata, whereas no obvious changes were observed in the full nutrient group, and only slight changes seen for uniconazole treatment group (Fig. 5). These results suggest that anthocyanin accumulation may be one of the main factors of flavonoid increase caused by nutrient starvation.

Fig. 4
figure 4

Flavonoid profiles of uniconazole or full nutrient treated L. punctata. a flavonoids of uniconazole treated L. punctata, 0342014-5-11, 0242014-5-11 and 0322014-5-11 corresponded to UT-0, UT-72 and UT-240, respectively. b flavonoids of full nutrient treated L. punctata, 0252014-5-10 and 0332014-5-11 corresponded to FN-72 and FN-240, respectively

Fig. 5
figure 5

Color change of frond underside under different growth conditions. L. punctata 0202 monoclonal was cultivated in 1/6 × Hoagland (FN), 1/6 × Hoagland and sprayed with 800 mg•L-1 uniconazole (Aoke Biotech Corp, Japan) solution on the surface (UT), or distilled water for 12 days (NS)

Discussion

Special phenylalanine metabolic flux resulted in high flavonoid and low lignin content

Duckweeds are the fastest growing and smallest flowering plants. A number of studies have shown the potential for duckweeds to be developed as feedstock, for biofuel production and as a natural purifier for swine wastewater treatment [21, 25, 4449], due to its high fermentable starch level (40-70% of dry weight), negligible lignin content and capacity to thrive on anthropogenic wastewater [24, 25, 33, 42, 5052]. More recently, its high flavonoid content (>4% of dry weight) in crude plant form [28, 29], has been found. As flavonoids play a crucial role in plant defense against pathogens [16, 53], they can be used to partially explain why duckweeds are rarely infected by pathogens. With near-exponential growth rates, duckweed can achieve a biomass of 13 to 38 metric tons/hectare/year dry weight [54], resulting in more than 520 kg/hectare/year flavonoid production. However, the molecular mechanism responsible for high flavonoid content remains largely uninvestigated. Newly-developed, high-throughput DNA sequencing technology provides an opportunity for genome-wide global transcriptome studies and metabolic pathway analyses. In this study, phenylalanine metabolism involved genes were carefully analyzed based on the RNA-Seq data of starvation, uniconazole and full nutrient treated L. punctata. Except flavone synthase, aureusidin synthase, flavanone 7-O-beta-glucosyltransferase, and flavanone 7-O-glucoside 2''-O-beta-L-rhamnosyltransferase, all of the other key enzymes involved in phenylalanine metabolism were successfully detected from the transcriptome (Fig. 1, Additional file 2: Table S2). p-coumaroyl CoA is the common substrate for the biosynthesis of flavonoid and lignin. The expression levels of CHS, HCT and CCR provided cues that p-coumaroyl CoA may be predominantly directed into the flavonoid branch and rarely into the lignin branch, resulting in the high flavonoid and low lignin content in L. punctata (Fig. 1). It is well known that lignin primarily provides mechanical support for plants to stand upright and enables xylems to withstand the negative pressure generated during water transport. Consequently, lignin is useless for L. punctata as these plants usually grow on the water’s surface with no need for mechanical support. To effectively avoid the accumulation of a helpless product, the metabolic flux is therefore mainly directed into the flavonoid branch. With this characteristic, L. punctata can be developed as a promising resource plant for biofuels fermentation and flavonoids extraction.

The following iTRAQ based proteomics analyses supported these results. The majority of lignin synthesis involved transcripts identified by RNA-Seq were not detected in the iTRAQ study. Although possibly due to technology bias, these lignin related enzymes were present in levels lower than the detection limit of this technology, as most enzymes involved in the other branch were successfully quantified using the same samples. Moreover, the global expression pattern of the phenylalanine metabolism pathway revealed that the metabolic flux was directed to the following anthocyanin biosynthesis branch with priority, but not the isoflavonoid biosynthesis or flavone and flavonol biosynthesis branches (Fig. 1). Since L. punctata also known as purple-backed duckweed due to the reddish-purple tint on the underside of its fronds as a result of anthocyanin production, the metabolic flux can be explained by its morphological characteristics. Furthermore, almost all of the enzymes involved in isoflavonoid biosynthesis or the flavone and flavonol biosynthesis pathway, failed to be de novo assembled using RNA-Seq reads, probably because of extremely low levels of expression. This study combined omics data to investigate flavonoid biosynthesis in L. punctata for the first time. The expression profiling not only gives a valuable insight into the molecular biological basis of the high flavonoid content in L. punctata, but also supports the morphological characteristics of this plant species by the analyses of metabolic flux.

Nutrient starvation is the optimized method to accumulate high starch and flavonoid content simultaneously in this resource plant

When growing L. punctata in distilled water, almost all “essential mineral nutrients” were deficient resulting in extreme nutrient starvation. To cope with this abiotic stress, L. punctata immediately increased expression of some transporters with the aim of increasing nutrient acquisition [25], but without success due to the absence of nutrients. The global physiological and metabolic status was altered and starch biosynthesis was enhanced, resulting in a high starch accumulation of 45% (dry weight) in 168 h [25]. These effects 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 (Fig. 6) [55]. As a class of important defense compounds, over-accumulation of flavonoids in plants can enhance stress tolerances by inhibiting the generation of ROS in plants [1315, 5661]. Manipulating flavonoid biosynthetic gene expression is an effective method to alter the accumulation of flavonoids in Arabidopsis and other plants [6265]. In nutrient starvation treated L. punctata, transcriptome analyses showed that most flavonoid involved transcripts were increased (Fig. 2, Additional file 2: Table S2), which was confirmed by iTRAQ based proteome results (Fig. 3, Additional file 3: Table S3). Metabolomic studies revealed a flavonoid accumulation from the original 4.51 to 5.56% (dry weight) after 168 h, with seven of the 17 detected flavonoid compounds having increased significantly [28], possibly due to the altered expression of flavonoid biosynthetic genes. Furthermore, purple color accumulation on the frond undersides correlated with the levels of flavonoids (Fig. 5). Overall, these integrated results from transcriptome, proteome, metabolome and morphology reveal a flavonoid based stress response in distilled water.

Fig. 6
figure 6

Growth status of L. punctata under different culture conditions. L. punctata 0202 monoclonal was cultivated in 1/6 × Hoagland, 1/6 × Hoagland spraying with 800 mg•L-1 uniconazole (Aoke Biotech Corp, Japan) solution on the surface, or distilled water for 12 days

Uniconazole, a plant growth retardant, has been extensively applied in plants to increase tolerance and improve quality by regulating endogenous hormone levels [6669]. Culturing duckweed in 1/6 × Hoagland medium and spraying with 800 mg · L-1 uniconazole is an optimized method to accumulate high starch content for bioethanol fermentation and biomass accumulation [26, 27]. The content of starch was increased from 3.16% to 48.01% in 240 h [27]. Different from that in distilled water treated L. punctata, the biomass of uniconazole treated L. punctata was almost equal to the control (1/6 × Hoagland) (Fig. 6), indicating that 1/6 × Hoagland and 800 mg · L-1 uniconazole did not create stress conditions and consequently did not trigger extra demand for flavonoids. As expected, flavonoid content was decreased from 2.83% to 0.94% at 168 h. Similarly, expression profiling results showed that more than half of the flavonoid involved genes were suppressed by this growth retardant. In our previous study, it was found that uniconazole increased the content of abscisic acid (ABA) and cytokinins (CK), and suppressed the synthesis of gibberellin (GA) [26, 27]. As ABA, CK and GA usually crosstalk with salicylic acid (SA), jasmonic acid (JA), other endogenous hormones and small molecule regulators, alterations ion the levels of these regulators should affect the regulatory network in L. punctata. Previous studies had revealed that different endogenous hormones usually play different roles in flavonoid biosynthesis. For example, sucrose can induce anthocyanin biosynthesis, but its effect is repressed by the addition of GA, whereas JA and ABA have a synergic effect with sucrose [70]. Accordingly, the decreased level of flavonoid may be a result of interference to the whole hormonal regulatory network in uniconazole treated L. punctata. However, whether the flavonoid decrease is primarily caused by the change of ABA, GA, CK, or uniconazole directly affecting the expression of some flavonoid related key genes, still requires further investigation in the future.

Full nutrient (1/6 × Hoagland) is an optimized culture medium, which usually cannot provide abiotic stress. So that the physiological status of L. punctata would not be altered in this study, starch and flavonoid content were kept at normal levels. Although the growth status of the nutrient starvation group was obviously suppressed (Fig. 6), the total starch weight was increased by 42 times in 7 days [25]. Comparatively, the biomass of uniconazole treated L. punctata was almost equal to the control (1/6 × Hoagland), with starch weight increasing by 46 times in 7 days [26]. Therefore, although starvation limits the accumulation of biomass, it can still have the same effect on starch accumulation, which is caused by a much lower dry matter rate in uniconazole treated L. punctata. Since the flavonoid content was increased by nutrient starvation, it can be surmised that nutrient starvation is the optimized method for obtaining high starch and high flavonoid content simultaneously in L. punctata, while uniconazole treatment can only produce high starch content.

Although some sampling time points of the RNA-Seq, iTRAQ and metabolome studies were inconsistent, the combined omics data reflect the changing trends of mRNAs, proteins and flavonoid compounds, as these high throughput technologies can characterize global gene expression patterns and metabolic status. In addition, a few discordant results appeared in the expression results in this study. As most enzymes were encoded by more than one transcript, the non-matching results may have been due to functional redundancy and spatio-temporal expression specificity of enzyme encoding transcripts.

Conclusions

Transcriptome and iTRAQ based expression profiling revealed that high flavonoid and low lignin content of L. punctata resulted primarily from phenylalanine metabolic flux directed towards the flavonoid biosynthetic pathway. Together with the metabolome assays, it was found that full nutrient medium generated high biomass with low starch and stable flavonoid content, uniconazole only induced starch accumulation accompanied by a decreased flavonoid content, while nutrient starvation triggered the accumulation of starch and flavonoids simultaneously. L. punctata has the potential to be developed as a resource plant for biofuel fermentation and flavonoid extraction.

Methods

Plant materials and treatments

Monoclonal L. punctata 0202 was cultivated in 1/6 × Hoagland nutrient solution (total N = 58.3 mg/L, P = 25.8 mg/L) for 14 days under a 16/8 h day/night cycle, with a light intensity of 130 μmol/m2/s, and a temperature of 25 °C/15 °C during the day/night. For the nutrient starvation group (NS), fresh fronds were transferred into distilled water for further cultivation over a period of two weeks. For the uniconazole treated group (UT), fronds were subsequently cultivated in 1/6 × Hoagland solution and sprayed on the surface with 800 mg · L-1 uniconazole (Aoke Biotech Corp, Japan) solution. The other groups were cultivated in 1/6 × Hoagland solution (FN). Different time points following the transfer of fronds into different media were selected for flavonoid analyses. For each time point, more than 3 g fresh fronds were collected from three culture flasks for each sample, corresponding to a total of >800 individuals.

RNA extraction and RNA-Sequencing analyses

For each frond sample, more than 1 g fronds was ground into powder in liquid nitrogen. Total RNA was extracted using OMEGA™ Plant DNA/RNA kit (OMEGA, USA), following the manufacturer’s instructions, and genomic DNA was removed by DNase I (Fermentas, USA). More than 20 μg total RNA was then submitted to Beijing Genomics Institute (BGI)-Shenzhen, Shenzhen, China [http://www.genomics.cn] for quality control. The purity, concentration and RNA integrity number (RIN) were measured by an Agilent 2100 Bioanalyzer or SMA3000. Qualified total RNAs were used for the following mRNA purification and 200 bp fragmented cDNA library construction, identical to that described in our previous study [25].

The validated fragmented cDNA library was submitted to the Illumina Hiseq 2000 platform at BGI for transcriptome sequencing. The 90 bp paired-end (PE) read sequence and base-calling quality values were obtained following the manufacturer’s instructions. The raw PE reads were qualified by removing the reads with adapter sequence or excessive unknown bases. The clean reads from the different samples were then pooled together and de novo assembled using Trinity (v2012-06-08) with the default parameters [37]. Length distribution was assessed by common perl scripts to generate the N50 number, average length, max length and contig number during different length intervals.

To profile the genome-wide expression patterns, all reads were aligned back to the assembly using perl scripts in the Trinity package (v2012-06-08) [37] for each RNA-Seq sample separately. The aligned read number was calculated and presented as digital expression levels for each contig. These values were then normalized for each RNA-Seq sample by RESM-based algorithm using perl scripts in the Trinity package (v2012-06-08) [37] to get FPKM values.

A BlastX sequence similarity search against the non-redundant protein database (NR) of NCBI [http://www.ncbi.nlm.nih.gov/] was conducted by a locally installed blast program (ncbi-blast-2.2.28+, ftp://ftp.ncbi.nlm.nih.gov/blast/executables/blast+/) to investigate functional annotation for each contig. BlastX results were uploaded into the Blast2GO platform [39, 40] for Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) annotation.

Protein extraction and iTRAQ based proteomic analyses

For each frozen sample, total protein extraction, qualification and digestion were performed as the method described in our previous study [28]. The digested peptides were labeled following the manufacturer’s protocol with iTRAQ® Reagent 8-plex Kit (AB SCIEX, USA) and subsequently used for LC-MS/MS analyses using an AB SCIEX TripleTOF™ 5600 mass spectrometer (AB SCIEX, USA), coupled with an LC-20AB HPLC Pump system (Shimadzu, Kyoto, Japan).

MS/MS data acquisition was performed with Analyst®QS2.0 software (AB SCIEX, USA), and processed by searching against the database generated from the annotated transcriptome using the Paragon™ Algorithm and the Mascot search engine (Matrix Science, London, UK; version 2.3.02). The relative abundance was analyzed by the report ion peak areas as previously described [71]. For protein quantitation, it was required that a protein contains at least two unique peptides.

Flavonoid content and classification

Flavonoid extraction and isolation were performed following the methods described in our previous study [29]. The flavonoid content of each frond sample was measured by spectrophotometry with a spectrophotometer (Varioskan Flash, Thermo Corp, USA) and HPLC (Thermo spectra system AS3000, Thermo Corp, USA)-UV (Thermo UV6000 Detector, USA) following the methods [72, 73]. HPLC/MS analyses of flavonoids were performed on an Agilent series 1100 HPLC instrument (Agilent, Waldbronn, Germany) coupled with a quadrupole time-of-flight (Q-TOF) mass spectrometry (micrOTOF-Q II; Bruker, Bremen, Germany) mainly in positive-ion mode. The ESI source conditions were set following the method of Yang [74]. The mass data were processed by Bruker Compass DataAnalysis 4.0 software.

Abbreviations

ABA:

Abscisic acid

CKs:

Cytokinins

DW:

Dry weight

EC:

Enzyme codes

FN:

Full nutrient

FPKM:

Fragments Per Kilobase of transcripts per Million mapped fragments

GA:

Gibberellins

iTRAQ:

Isobaric tags for relative and absolute quantitation

JA:

Jasmonic acid

log2FC:

log2 fold-change

NGS:

Next-generation sequencing

NR:

Non-redundant protein database

NS:

Nutrient starvation

PE:

Paired-end

RNA-Seq:

RNA-sequencing

SA:

Salicylic acid

UT:

Uniconazole treated

References

  1. Harborne JB, Williams CA. Advances in flavonoid research since 1992. Phytochemistry. 2000;55(6):481–504.

    Article  CAS  PubMed  Google Scholar 

  2. Winkel-Shirley B. Flavonoid biosynthesis. A colorful model for genetics, biochemistry, cell biology, and biotechnology. Plant Physiol. 2001;126(2):485–93.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Mazur WM, Duke JA, Wähälä K, Rasku S, Adlercreutz H. Isoflavonoids and lignans in legumes: nutritional and health aspects in humans. J Nutr Biochem. 1998;9(4):193–200.

    Article  CAS  Google Scholar 

  4. Crowden R, Jarman S. 3-deoxyanthocyanins from the fern Blechnum procerum. Phytochemistry. 1974;13(9):1947–8.

    Article  CAS  Google Scholar 

  5. Awika JM, Rooney LW, Waniska RD. Properties of 3-deoxyanthocyanins from sorghum. J Agr Food Chem. 2004;52(14):4388–94.

    Article  CAS  Google Scholar 

  6. Brown DE, Rashotte AM, Murphy AS, Normanly J, Tague BW, Peer WA, Taiz L, Muday GK. Flavonoids act as negative regulators of auxin transport in vivo in Arabidopsis. Plant Physiol. 2001;126(2):524–35.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Buer CS, Djordjevic MA. Architectural phenotypes in the transparent testa mutants of Arabidopsis thaliana. J Exp Bot. 2009;60(3):751–63.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Martin C, Gerats T. Control of pigment biosynthesis genes during petal development. Plant Cell. 1993;5(10):1253–64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Koes RE, Van Blokland R, Quattrocchio F, Van Tunen AJ, Mol J. Chalcone synthase promoters in petunia are active in pigmented and unpigmented cell types. Plant Cell. 1990;2(5):379–92.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Taylor LP, Grotewold E. Flavonoids as developmental regulators. Curr Opin Plant Biol. 2005;8(3):317–23.

    Article  CAS  PubMed  Google Scholar 

  11. Shirley BW. Flavonoid biosynthesis:‘new’functions for an ‘old’pathway. Trends Plant Sci. 1996;1(11):377–82.

    Google Scholar 

  12. Brunetti C, Di Ferdinando M, Fini A, Pollastri S, Tattini M. Flavonoids as antioxidants and developmental regulators: relative significance in plants and humans. Int J Mol Sci. 2013;14(2):3540–55.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Agati G, Azzarello E, Pollastri S, Tattini M. Flavonoids as antioxidants in plants: location and functional significance. Plant Sci. 2012;196:67–76.

    Article  CAS  PubMed  Google Scholar 

  14. Akhtar TA, Lees HA, Lampi MA, Enstone D, Brain RA, Greenberg BM. Photosynthetic redox imbalance influences flavonoid biosynthesis in Lemna gibba. Plant Cell Environ. 2010;33(7):1205–19.

    CAS  PubMed  Google Scholar 

  15. Chalker-Scott L. Environmental significance of anthocyanins in plant stress responses. Photochem Photobiol. 1999;70(1):1–9.

    Article  CAS  Google Scholar 

  16. Cushnie TP, Lamb AJ. Antimicrobial activity of flavonoids. Int J Antimicrob Agents. 2005;26(5):343–56.

    Article  CAS  PubMed  Google Scholar 

  17. Friedman M. Antibacterial, antiviral, and antifungal properties of wines and winery byproducts in relation to their flavonoid content. J Agric Food Chem. 2014;62(26):6025–42.

    Article  CAS  PubMed  Google Scholar 

  18. Landolt E. Biosystematic investigations in the family of duckweed (Lemnaceae), vol. 2. Stiftung Rubel, Zúrich, Suiza: Geobotanischen Inst ETH; 1986.

    Google Scholar 

  19. Hillman WS, Culley Jr DD. The uses of duckweed: The rapid growth, nutritional value, and high biomass productivity of these floating plants suggest their use in water treatment, as feed crops, and in energy-efficient farming. Am Sci. 1978;66:442–51.

    Google Scholar 

  20. Cheng J, Bergmann BA, Classen JJ, Stomp AM, Howard JW. Nutrient recovery from swine lagoon water by Spirodela punctata. Bio-Resour Technol. 2002;81(1):81–5.

    Article  CAS  Google Scholar 

  21. Cheng JJ, Stomp AM. Growing duckweed to recover nutrients from wastewaters and for production of fuel ethanol and animal feed. Clean-Soil Air Water. 2009;37(1):17–26.

    Article  CAS  Google Scholar 

  22. El-Shafai SA, El-Gohary FA, Nasr FA, Peter van der Steen N, Gijzen HJ. Nutrient recovery from domestic wastewater using a UASB-duckweed ponds system. Bioresour Technol. 2007;98(4):798–807.

    Article  CAS  PubMed  Google Scholar 

  23. Swanson SM, Rickard CP, Freemark KE, MacQuarrie P. Testing for pesticide toxicity to aquatic plants: recommendations for test species. Plants for toxicity assessment vol. 2. Philadelphia: American Society for Testing and Materials; 1991. pp. 77–97.

  24. Blazey EB, McClure JW. The distribution and taxonomic significance of lignin in the Lemnaceae. Am J Bot. 1968;55:1240–5.

    Article  Google Scholar 

  25. Tao X, Fang Y, Xiao Y, Jin YL, Ma XR, Zhao Y, He KZ, Zhao H, Wang HY. Comparative transcriptome analysis to investigate the high starch accumulation of duckweed (Landoltia punctata) under nutrient starvation. Biotechnol Biofuels. 2013;6(1):72.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Liu Y, Fang Y, Huang M, Jin Y, Sun J, Tao X, Zhang G, He K, Zhao Y, Zhao H. 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 

  27. Liu Y, Fang Y, Huang M, Jin Y, Sun J, Tao X, Zhang G, He K, Zhao Y, Zhao H. Uniconazole-induced starch accumulation in the bioenergy crop duckweed (Landoltia punctata) II: transcriptome alterations of pathways involved in carbohydrate metabolism and endogenous hormone crosstalk. Biotechnol Biofuels. 2015;8:64.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Huang MJ, Fang Y, Xiao Y, Sun JL, Jin YL, Tao X, Ma XR, He KZ, Zhao H. Proteomic analysis to investigate the high starch accumulation of duckweed (Landoltia punctata) under nutrient starvation. Ind Crop Prod. 2014;59:299–308.

    Article  CAS  Google Scholar 

  29. Wang N, Xu G, Fang Y, Yang T, Zhao H, Li G. New flavanol and cycloartane glucosides from Landoltia punctata. Molecules. 2014;19(5):6623–34.

    Article  PubMed  Google Scholar 

  30. Liu B, Zhu Y. Extraction of flavonoids from flavonoid-rich parts in tartary buckwheat and identification of the main flavonoids. J Food Eng. 2007;78(2):584–7.

    Article  CAS  Google Scholar 

  31. Xie J, Zhu L, Luo H, Zhou L, Li C, Xu X. Direct extraction of specific pharmacophoric flavonoids from gingko leaves using a molecularly imprinted polymer for quercetin. J Chromatogr A. 2001;934(1):1–11.

    Article  CAS  PubMed  Google Scholar 

  32. Wang J, Cao F, Su E, Wu C, Zhao L, Ying R. Improving flavonoid extraction from Ginkgo biloba leaves by prefermentation processing. J Agric Food Chem. 2013;61(24):5783–91.

    Article  CAS  PubMed  Google Scholar 

  33. Cui WH, Xu JL, Cheng JJ, Stomp AM. Starch accumulation in duckweed for bioethanol production. Biol Eng. 2011;3(4):187–97.

    Article  Google Scholar 

  34. Chen Q, Jin YL, Zhang GH, Fang Y, Xiao Y, Zhao H. Improving production of bioethanol from duckweed (Landoltia punctata) by pectinase pretreatment. Energies. 2012;5(8):3019–32.

    Article  CAS  Google Scholar 

  35. Anderson KE, Lowman Z, Stomp A-M, Chang J. Duckweed as a feed ingredient in laying hen diets and its effect on egg production and composition. Int J Poult Sci. 2011;10(1):4–7.

    Article  CAS  Google Scholar 

  36. Lam E, Appenroth KJ, Michael T, Mori K, Fakhoorian T. Duckweed in bloom: the 2nd International Conference on Duckweed Research and Applications heralds the return of a plant model for plant biology. Plant Mol Biol. 2014;84(6):737–42.

    Article  CAS  PubMed  Google Scholar 

  37. Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng Q, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29(7):644–52.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139–40.

    Article  CAS  PubMed  Google Scholar 

  39. Conesa A, Götz S, García-Gómez JM, Terol J, Talón M, Robles M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005;21(18):3674–6.

    Article  CAS  PubMed  Google Scholar 

  40. Conesa A, Gotz S. Blast2GO: a comprehensive suite for functional analysis in plant genomics. Int J Plant Genomics. 2008;2008:619832.

    Article  PubMed  Google Scholar 

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

    Article  CAS  Google Scholar 

  42. Van Dyke JM, Sutton DL. Digestion of duckweed (Lemna spp.) by the grass carp (Ctenopharyngodon idella). J Fish Biol. 1977;11(3):273–8.

    Article  Google Scholar 

  43. Oron G. Duckweed culture for wastewater renovation and biomass production. Agr Water Manage. 1994;26(1–2):27–40.

    Article  Google Scholar 

  44. Ge X, Zhang N, Phillips GC, Xu J. Growing Lemna minor in agricultural wastewater and converting the duckweed biomass to ethanol. Bioresour Technol. 2012;124:485–8.

    Article  CAS  PubMed  Google Scholar 

  45. Chaiprapat S, Cheng JJ, Classen JJ, Liehr SK. Role of internal nutrient storage in duckweed growth for swine wastewater treatment. T Asae. 2005;48(6):2247–58.

    Article  CAS  Google Scholar 

  46. El-Shafai SA, Abdel-Gawad F, Samhan F, Nasr FA. Resource recovery from septic tank effluent using duckweed-based tilapia aquaculture. Environ Technol. 2013;34(1-4):121–9.

    Article  CAS  PubMed  Google Scholar 

  47. Ozengin N, Elmaci A. Performance of duckweed (Lemna minor L.) on different types of wastewater treatment. J Environ Biol. 2007;28(2):307–14.

    CAS  PubMed  Google Scholar 

  48. Xu J, Cheng JJ, Stomp AM. Nutrient removal from swine wastewater by growing duckweed: a pilot study. T Asabe. 2012;55(1):181–5.

    Article  CAS  Google Scholar 

  49. Xu J, Cheng JJ, Stomp A-M. Growing Spirodela polyrrhiza in swine wastewater for the production of animal feed and fuel ethanol: a pilot study. CLEAN-Soil Air Water. 2012;40(7):760–5.

    Article  CAS  Google Scholar 

  50. Xu JL, Cui WH, Cheng JJ, Stomp AM. Production of high-starch duckweed and its conversion to bioethanol. Biosyst Eng. 2011;110(2):67–72.

    Article  Google Scholar 

  51. Wang W, Haberer G, Gundlach H, Glasser C, Nussbaumer T, Luo MC, Lomsadze A, Borodovsky M, Kerstetter RA, Shanklin J, et al. The Spirodela polyrhiza genome reveals insights into its neotenous reduction fast growth and aquatic lifestyle. Nat Commun. 2014;5:3311.

    CAS  PubMed  PubMed Central  Google Scholar 

  52. Seth PN, Venkataraman R, Maheshwari SC. Studies on the growth and flowering of a short-day plant. Wolffia Microscopica Planta. 1970;90(4):349–59.

    CAS  PubMed  Google Scholar 

  53. Huang J, Li J, Zheng S, Wu J, Zhang W, Sun T, Dewan SK, Kalionis B, Shen Z, Tai X, et al. Epimedium flavonoids counteract the side effects of glucocorticoids on hypothalamic-pituitary-adrenal axis. Evid Based Complement Altern Med. 2013;2013:938425.

    Google Scholar 

  54. Riemer DN. Duckweed aquaculture: a new aquatic farming system for developing countries. Soil Sci. 1994;157(3):200–1.

    Article  Google Scholar 

  55. Blanvillain R, Wei S, Wei P, Kim JH, Ow DW. Stress tolerance to stress escape in plants: role of the OXS2 zinc-finger transcription factor family. EMBO J. 2011;30(18):3812–22.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Isshiki R, Galis I, Tanakamaru S. Farinose flavonoids are associated with high freezing tolerance in fairy primrose (Primula malacoides) plants. J Integr Plant Biol. 2014;56(2):181–8.

    Article  CAS  PubMed  Google Scholar 

  57. Nakabayashi R, Yonekura-Sakakibara K, Urano K, Suzuki M, Yamada Y, Nishizawa T, Matsuda F, Kojima M, Sakakibara H, Shinozaki K, et al. Enhancement of oxidative and drought tolerance in Arabidopsis by overaccumulation of antioxidant flavonoids. Plant J. 2014;77(3):367–79.

    Article  CAS  PubMed  Google Scholar 

  58. Dixon RA, Paiva NL. Stress-induced phenylpropanoid metabolism. Plant Cell. 1995;7(7):1085–97.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Lu Y, Du J, Tang J, Wang F, Zhang J, Huang J, Liang W, Wang L. Environmental regulation of floral anthocyanin synthesis in Ipomoea purpurea. Mol Ecol. 2009;18(18):3857–71.

    Article  CAS  PubMed  Google Scholar 

  60. Rowan DD, Cao M, Lin-Wang K, Cooney JM, Jensen DJ, Austin PT, Hunt MB, Norling C, Hellens RP, Schaffer RJ, et al. Environmental regulation of leaf colour in red 35S:PAP1 Arabidopsis thaliana. New Phytol. 2009;182(1):102–15.

    Article  CAS  PubMed  Google Scholar 

  61. Stracke R, Favory JJ, Gruber H, Bartelniewoehner L, Bartels S, Binkert M, Funk M, Weisshaar B, Ulm R. The Arabidopsis bZIP transcription factor HY5 regulates expression of the PFG1/MYB12 gene in response to light and ultraviolet-B radiation. Plant Cell Environ. 2010;33(1):88–103.

    CAS  PubMed  Google Scholar 

  62. Shirley BW, Kubasek WL, Storz G, Bruggemann E, Koornneef M, Ausubel FM, Goodman HM. Analysis of Arabidopsis mutants deficient in flavonoid biosynthesis. Plant J. 1995;8(5):659–71.

    Article  CAS  PubMed  Google Scholar 

  63. Borevitz JO, Xia Y, Blount J, Dixon RA, Lamb C. Activation tagging identifies a conserved MYB regulator of phenylpropanoid biosynthesis. Plant Cell. 2000;12(12):2383–94.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  64. Butelli E, Titta L, Giorgio M, Mock HP, Matros A, Peterek S, Schijlen EG, Hall RD, Bovy AG, Luo J, et al. Enrichment of tomato fruit with health-promoting anthocyanins by expression of select transcription factors. Nat Biotechnol. 2008;26(11):1301–8.

    Article  CAS  PubMed  Google Scholar 

  65. Butelli E, Licciardello C, Zhang Y, Liu J, Mackay S, Bailey P, Reforgiato-Recupero G, Martin C. Retrotransposons control fruit-specific, cold-dependent accumulation of anthocyanins in blood oranges. Plant Cell. 2012;24(3):1242–55.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  66. Saito S, Okamoto M, Okamoto M, Kushiro T, Koshiba T, Kamiya Y, Hirai N, Todoroki Y, Sakata K, Nambara E, et al. A plant growth retardant, uniconazole, is a potent inhibitor of ABA catabolism in Arabidopsis. Biosci Biotech Bioch. 2006;70(7):1731–9.

    Article  CAS  Google Scholar 

  67. English PJ, Lycett GW, Roberts JA, Jackson MB. Increased 1-Aminocyclopropane-1-carboxylic acid oxidase activity in shoots of flooded tomato plants raises ethylene production to physiologically active levels. Plant Physiol. 1995;109(4):1435–40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  68. Izumi K, Nakagawa S, Kobayashi M, Oshio H, Sakurai A, Takahashi N. Levels of IAA, cytokinins, ABA and ethylene in rice plants as affected by a gibberellin biosynthesis inhibitor, uniconazole-P. Plant Cell Physiol. 1988;29(1):97–104.

    CAS  Google Scholar 

  69. Zhou WJ, Leul M. Uniconazole-induced alleviation of freezing injury in relation to changes in hormonal balance, enzyme activities and lipid peroxidation in winter rape. Plant Growth Regul. 1998;26(1):41–7.

    Article  CAS  Google Scholar 

  70. Loreti E, Povero G, Novi G, Solfanelli C, Alpi A, Perata P. Gibberellins, jasmonate and abscisic acid modulate the sucrose-induced expression of anthocyanin biosynthetic genes in Arabidopsis. New Phytol. 2008;179(4):1004–16.

    Article  CAS  PubMed  Google Scholar 

  71. Ford KL, Cassin A, Bacic A. Quantitative proteomic analysis of wheat cultivars with differing drought stress tolerance. Front Plant Sci. 2011;2:44.

    Article  PubMed  PubMed Central  Google Scholar 

  72. Qiao X, He W, Xiang C, Han J, Wu L, Guo D, Ye M. Qualitative and quantitative analyses of flavonoids in Spirodela polyrrhiza by high-performance liquid chromatography coupled with mass spectrometry. Phytochem Anal. 2011;22(6):475–83.

    Article  CAS  PubMed  Google Scholar 

  73. Gülçin İ, Kireçci E, Akkemik E, Topal F, Hisar O. Antioxidant, antibacterial, and anticandidal activities of an aquatic plant: duckweed (Lemna minor L. Lemnaceae). Turk J Biol. 2010;34:175–88.

    Google Scholar 

  74. Yang W, Fang D-M, He H-P, Hao X-J, Wu Z-J, Zhang G-L. Analysis of mexicanolide- and phragmalin-type limonoids from Heynea trijuga using high-performance liquid chromatography/electrospray tandem mass spectrometry. Rapid Commun Mass Sp. 2013;27(11):1203–12.

    Article  CAS  Google Scholar 

  75. Deng W, Wang Y, Liu Z, Cheng H, Xue Y. HemI: a toolkit for illustrating heatmaps. Plos One. 2014;9(11):e111988.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This work was supported by the National Natural Science Foundation of China (No. 31400218), the National Key Technology R&D Program of China (No. 2015BAD15B01) and the Projects of International Cooperation of Ministry of Science and Technology of China (No. 2014DFA30680).

Availability of data and materials

The transcriptome datasets supporting the conclusions of this article are available in the [NCBI’s Sequence Read Archive database (SRA) database] repository under the accession number of PRJNA185389 [unique persistent identifier and hyperlink to dataset(s) in http://www.ncbi.nlm.nih.gov/sra?term=PRJNA185389].

Authors’ contributions

XT conceived the study, carried out the data analysis, drafted and revised the manuscript. YF drafted and revised the manuscript. YX carried out the data analysis and revised the manuscript. MJH carried out the biochemical assays and HPLC experiment, drafted and revised the manuscript. YL carried out the biochemical assays and revised the manuscript. XRM conceived the study and revised the manuscript. HZ conceived the study and revised the manuscript. All authors read and approved the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Not applicable.

Ethics approval and consent to participate

Not applicable.

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Xin-Rong Ma or Hai Zhao.

Additional files

Additional file 1: Table S1.

Sequence annotation and expression profiling. (XLSX 20158 kb)

Additional file 2: Table S2.

Transcript levels of flavonoid metabolism related contigs. (XLSX 58 kb)

Additional file 3: Table S3.

Quantities of flavonoid metabolism related proteins. (XLSX 11 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Tao, X., Fang, Y., Huang, MJ. et al. High flavonoid accompanied with high starch accumulation triggered by nutrient starvation in bioenergy crop duckweed (Landoltia punctata). BMC Genomics 18, 166 (2017). https://doi.org/10.1186/s12864-017-3559-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-017-3559-z

Keywords