Skip to main content

Transcriptome profiling of the floating-leaved aquatic plant Nymphoides peltata in response to flooding stress



Waterlogging or flooding is one of the most challenging abiotic stresses experienced by plants. Unlike many flooding-tolerant plants, floating-leaved aquatic plants respond actively to flooding stress by fast growth and elongation of its petioles to make leaves re-floating. However, the molecular mechanisms of this plant group responding to flood have not been investigated before. Here, we investigated the genetic basis of this adaptive response by characterizing the petiole transcriptomes of a floating-leaved species Nymphoides peltata under normal and flooding conditions.


Clean reads under normal and flooding conditions with pooled sampling strategy were assembled into 124,302 unigenes. A total of 8883 unigenes were revealed to be differentially expressed between normal and flooding conditions. Among them, top ranked differentially expressed genes were mainly involved in antioxidant process, photosynthesis process and carbohydrate metabolism, including the glycolysis and a modified tricarboxylic acid cycle – alanine metabolism. Eight selected unigenes with significantly differentiated expression changes between normal and flooding conditions were validated by qRT-PCR.


Among these processes, antioxidant process and glycolysis are commonly induced by waterlogging or flooding environment in plants, whereas photosynthesis and alanine metabolism are rarely occurred in other flooding-tolerant plants, suggesting the significant contributions of the two processes in the active response of N. peltata to flooding stress. Our results provide a valuable genomic resource for future studies on N. peltata and deepen our understanding of the genetic basis underlying the response to flooding stress in aquatic plants.


Understanding how individuals respond to ever-changing environments is of fundamental importance in all organisms. Plants encounter various biotic and abiotic stresses throughout their life span. One of the most dramatic abiotic stresses is flooding [1]. Plants cannot actively escape a flooding environment due to their immobile character [2], and many plants are injured or killed by flooding events due to oxygen shortage in their cells [3, 4]. However, over the long period of adaptive evolution, plants have evolved the capacity to survive flooding habitats via escape phenotypes, e.g., the shoot elongation, the formation of aerenchyma, and the induction of gas films [5]. Recently, the molecular mechanisms of physiology and metabolic modulation behind these adaptive traits have been investigated in crops and wetland plants [610]. However, few such studies have been conducted on true aquatics, floating-leaved or submerged aquatic plants [11], which are specialized to life in water and likely have different responding mechanisms to flooding compared with other plants.

The fringed water lily Nymphoides peltata (S. G. Gmelin) Kuntze is a typical floating-leaved plant with a widespread distribution in temperate and subtropical regions of Eurasia [12]. Nymphoides peltata usually roots in the bottom mud and maintains its leaves afloat on the water surface with the connection of petioles. When subjected to flooding, the leaves of N. peltata can rapidly reach the water surface by rapid elongation of the petioles [13, 14], an ability that is also present in many other floating-leaved plants [15]. This trait makes N. peltata an ideal aquatic species for investigating the molecular mechanisms of plants to avoid submergence stress. Previous study showed that ethylene played a major role in the case where submergence promotes petiole elongation in N. peltata [13]. which is similar with many aquatic or flooding-tolerant species, e.g., Rumex palustris [6] and rice [10]. However, the molecular mechanisms of flood-adaption in N. peltata still remain unclear. As a non-model plant, genomic information of N. peltata is scarce, except for the development of several molecular markers [16, 17], which hinders the exploration of the underlying flood-adaption mechanisms in this plant. With the emergence of next-generation sequencing technologies, a new technology RNA-Seq (RNA sequencing) independent of genetic background has been developed [18, 19]. Recently, RNA-Seq has been utilized to elucidate the response of non-model plants to various environmental stresses, including flooding and waterlogging [2023].

In this study, we examined the global gene expression changes of N. peltata under both normal and flooding conditions using Illumina RNA-Seq technology. The results provide a comprehensive view of the complex molecular events involved in the response of floating-leaved plants to flooding stress and expand our understanding of response to flooding stress.


Illumina sequencing and assembly

In total, 78,037,588 and 103,266,542 clean reads were obtained from the untreated sample (US) and the treated sample (TS), respectively (Table 1). Assembly of those reads from US and TS separately generated 87,673 and 95,372 unigenes, respectively (Table 1). The strategy of pooling all clean reads from US and TS together generated 124,302 unigenes with a N50 length of 1449 bp after assembly (Table 1). The number of unigenes from the pooled strategy was larger than those from separate assemble because some unigenes with low expression levels generated with the pooled strategy cannot be found when using strategy of separate assemble due to their less reads. Among these 124,302 unigenes, the total number of unigenes longer than 500 bp was 56,943, accounting for 45.81% of all unigenes (Additional file 1). The detailed length distribution of 124,302 unigenes predicted from the pooled assembly can be also found in Additional file 1.

Table 1 Overview of transcriptome sequencing and unigene annotations

Gene annotation and functional classification

Unigenes generated using the pooled strategy were used for further analysis. Among these 124,302 unigenes, a total of 53,870 (43.34%) unigenes were annotated (Table 1). According to the NCBI non-redundant proteins (NR) annotation, 53,744 (43.24%) unigenes had homologous proteins in the NR protein database (Table 1). Meanwhile, unigenes were also matched with GO category, and 30,326 (24.40%) unigenes were assigned to one or more GO terms (Table 1). Using GO annotation, the functions of the unigenes were divided into three categories: biological process, cellular component, and molecular function. To obtain a detailed view of the GO classification, each GO category was further clustered to its parent term (Additional file 2). The results indicate that most of the sequenced genes were responsible for fundamental biological regulation and metabolism.

Furthermore, the possible functions of 124,302 unigenes were predicted using searches against the Cluster of orthologous groups (COG) database as well as Swiss-Prot Protein Sequence (Swissprot) database. A total of 48,156 unigenes were matched with the COG database (Table 1) and classified into 25 specific categories (Additional file 2). The “General function prediction only” was the largest group (20.08%), followed by “Posttranslational modification, protein turnover, chaperones” (9.99%) and “Signal transduction mechanisms” (9.21%) (Additional file 2). We also obtained 37,556 hits when searched against Swissprot protein database with an E-value of 1.0E-5, covering 31.62% of all unigenes (Additional file 3). All the deduced protein sequences corresponding to the predicted coding DNA sequences (CDS) were listed in Additional file 4. Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis was also conducted to predict the function of the unigenes during the development process. In total, 11,409 unigenes were found to be involved in 244 pathways (Table 1).

We also used the phytozome database ( to annotate unigenes. Similar as the NR database, the percent of annotated unigenes was about 43% (Additional file 5), which is much lower than a normal annotation percentage 70% in other plants. To examine whether these unannotated unigenes were non-coding RNA genes or lineage-specific genes, we predicted sequences of the unannotated unigenes and found that most of them have predicted CDS. We further used the unannotated unigenes as queries to search against the plant noncoding RNAs database and found that only a small proportion of unigenes matched the sequences in the noncoding RNAs database (Additional file 5). Therefore, most of the unannotated unigenes were inferred to be Nymphoides peltata lineage-specific genes.

Comparison between treated and untreated samples

Functional annotation and COG classification for the separately assembled unigenes for samples TS and US was also conducted after annotation of unigenes using the pooled strategy, the results were shown in Fig. 1. The number of unigenes in each category for GO and COG annotation showed parallel proportion between these two samples, indicating that transcriptomes of sample TS and US were comparable.

Fig. 1
figure 1

GO (Gene ontology) and COG (Cluster of orthologous groups) annotation. a GO annotation of the assembled unigenes for the treated sample (TS) and the untreated sample (US) separately, the functions of the unigenes were divided into three categories. b Information of COG classification for the treated sample (TS) and the untreated sample (US) separately, the unigenes were mainly clustered into 25 components

Furthermore, we listed the top 20 ranked unigenes with high expression level in sample TS and the corresponding unigenes in sample US (Table 2). Among these 20 unigenes, 10 were DEGs, which encoded proteins mainly involved in osmoregulation (e.g., aquaporin and osmotin-like protein) and proline-rich protein (Table 2). The former is very helpful in the regulation of water homeostasis and water transport, and the latter is a cell wall protein of plant regulating plant wound and defense responses. As for the other 10 unigenes, they were mainly involved in photosynthesis process (e.g., ribulose bisphosphate carboxylase and chlorophyll a-b binding protein) and antioxidant process (e.g., peroxidase and glutamine synthetase cytosolic isozyme) (Table 2).

Table 2 The top 20 ranked unigenes (based on expression level) in the TS sample and the corresponding unigenes in the US sample

Analysis of differentially expressed genes

After calculating the unigene expression level, with an FDR (false discovery rate) of 0.05 and |log2Fold Change| ≥ 1 as a cutoff, a total of 8883 (6401 up-regulated and 2482 down-regulated) unigenes were revealed to be significantly differentially expressed between the treated and untreated samples. The top 20 ranked differentially expressed unigenes were identified. Among them, most genes were involved in basic cell component (e.g., proline-rich protein and ribosomal protein), photosynthesis (e.g., chlorophyll a-b binding protein, chloroplastic glyceraldehyde 3-phosphate dehydrogenase, and ribulose bisphosphate carboxylase) and antioxidant process (e.g., L-ascorbate oxidase homolog, ferric reduction oxidase, and peroxidase) (Table 3).

Table 3 The top 20 ranked most differentially expressed unigenes between normal and flooding conditions

GO enrichment analysis of DEGs indicated that 153 GO terms were significantly enriched with the criteria of FDR < 0.01. Of them, 71 terms were enriched in the category of biological process, 52 in molecular function, and 30 in cellular component (Additional file 6). Further analysis showed that GO terms related to three aspects (energy, antioxidant, photosynthesis) were significantly important (Fig. 2).

Fig. 2
figure 2

GO enrichment analysis of putatively selected GO terms of DEG functions. DEG differentially expressed gene. Gene number is shown next to each GO term

KEGG pathway analysis of the DEGs indicated that various genes were working together to execute functions. Overall, the most significant pathways were Ribosome pathway (ko03010) with 315 DEGs enriched, followed by pathways of Photosynthesis (ko00195, 99) and Photosynthesis - antenna protein (ko00196, 91) (Fig. 3). Further analysis showed that pathways involved in energy metabolism, antioxidant process, and photosynthesis process were also present (Fig. 3).

Fig. 3
figure 3

Significantly enriched KEGG pathways of DEG functions. DEG differentially expressed gene

Transcriptome changes of Nymphoides peltata and five other plants in response to flooding

To gain a broad picture of plant responses to flooding, we listed transcriptome changes of Nymphoides peltata and five other plants to cope with flooding in Table 4. It appears that common responses to flooding in transcriptome were involved in glycolysis, the antioxidant process, the fermentation and the present of group VII ethylene response factor transcription factors (GVIIERFs) (Table 4). By contrast, expression changes of the alanine metabolism under submergence were only observed in N. peltata and Lotus japonicas (Table 4). As a species of true aquatic plant, N. peltata remained active photosynthesis under submergence, which has not been reported in other plants (Table 4).

Table 4 Transcriptome changes of Nymphoides peltata (this study) and five other plants in response to flooding

Molecular adaptation of aquatic plants compared with non-aquatic plants

To identify molecular adaptation in aquatic plants and help understand flooding response in plants, we applied a likelihood approach to detect molecular adaptation by estimating the ratio (ω) of nonsynonymous to synonymous substitution rates. Our dataset contained three aquatic plants and five non-aquatic plants. Within the aquatic plants, one (Utricularia gibba) has a draft genome, and the remaining two (Ranunculus bungei, and Nymphoides peltata) have transcriptome data. By contrast, all five non-aquatic plants (Arabidopsis thaliana, Oryza sativa, Solanum lycopersicum, Daucus carota, and Cucumis sativus) possess available genome sequences. A total of 5319 one-to-one orthologous genes were identified in our dataset. Using an established species tree from the Angiosperm Phylogeny Website (, we tested the possibility of differential selection between aquatic and non-aquatic plants, by comparing a one-ratio model with a two-ratio model. The one-ratio model assumes a uniform ω ratio across all plants, while the two-ratio model assumes two ω ratios: all non-aquatic plants have ω0, whereas all aquatic plants have ω1. Our analyses yielded 18 genes with significantly higher evolutionary rates in aquatic plants (Additional file 7), indicating that these genes may play important roles in adaptation to aquatic environments. Notably, two of these genes were outstanding for their antioxidant role. One is the gene encoding NAT2 (Nucleobase-ascorbate transporter), which is involved in the ascorbate transport and the ascorbate play an important role in antioxidant [24] (Additional file 7). The other gene is CAT2 encoding an enzyme Catalase-2, which is an anti-oxidant enzyme to protect cells from the toxic effects of hydrogen peroxide [25] (Additional file 7).

Validation of the RNA-Seq data

To verify our RNA-Seq data, the expression levels of eight unigenes were examined by quantitative real-time PCR (qRT-PCR) experiment. Of the eight genes, two genes are known to encode Glutathione S-transferases (GST) and Glutathione peroxidases (GPX), which are antioxidant enzymes (Table 5); two genes encode Photosystem I subunit XI (PsaL) and Photosynthetic electron transport ferredoxin (PetF), which are associated with photosynthesis process (Table 5); four genes encode Alcohol dehydrogenase (ADH), Lactate dehydrogenase (LDH), Alanine aminotransferase (AlaAT) and Aspartate transaminase (AspAT), which are involved in carbohydrate metabolism, including fermentative pathway and a modified tricarboxylic acid cycle mode - alanine metabolism (Table 5). All eight genes exhibited higher expression level in response to flooding (Table 5), indicating that these genes are indeed stimulated by flooding stress. In addition, their expression changes measured from qRT-PCR were consistent with those estimated from the RNA-Seq data (Table 5).

Table 5 Validation of the DEGs with qRT-PCR experiment


In this study, using the RNA-Seq technology, transcriptome changes of N. peltata were obtained under both normal and flooding conditions. In total, 78,037,588 and 103,266,542 high-quality reads were acquired for US and TS sample, respectively. Gene annotation towards various databases (NR, COG, Swissprot, GO, and KEGG) was conducted after reads assembly. After calculation of gene expression for each unigene and a strict criterion, a total of 8883 unigenes were defined as DEGs, which were either induced or depressed by flooding stress. Besides, the top 20 ranked differentially expressed unigenes were mainly involved in antioxidant, photosynthesis and protein-related activity. In addition, GO enrichment and KEGG enrichment analysis of DEGs also presented processes involved in carbohydrate source (photosynthesis), self-protection area (antioxidants) and energy supplier (glycolysis). Moreover, the molecular adaptation of the aquatic plants compared with non-aquatic plants suggested that genes involved in anti-oxidant processes having significantly higher evolutionary rates in aquatic plants. These findings indicate that processes mentioned above play important roles in response of N. peltata to flooding.

When plants under flooding pressure, many physiological responses will be induced. Ethylene release was one of these responses. Previous studies showed that the shoot elongation strategy was mainly driven by ethylene, such as rice and Rumex [10, 22]. In N. peltata, the shoot elongation was also driven by ethylene and has been physiologically validated [13]. In addition, the cell number was increased to a large number, which contributed to the rapid elongation. Our transcriptome results identified some ethylene response factors and some cell cycle related genes such as CDC7 and CAFP, which proved the fact that ethylene release was induced by flooding as a physiological response. Moreover, flooding often leads plants to oxygen shortage situation, which can induce ROS generation, and anti-oxidant enzymes were activated to protect the cell membrane [26, 27]. In N. peltata, though we were not able to detect enzyme activity, we can infer that the activity of antioxidant enzymes was induced under flooding stress. Several antioxidant genes such as GST and GPX and two fast evolving genes NAT2 and CAT2 involved in anti-oxidant were identified, indicating an activated ROS elimination process. Taken together, these results validated the physiological response of N. peltata under flooding stress.

Plants are challenged by various abiotic stresses. The increased accumulation of reactive oxygen species (ROS), including singlet oxygen, superoxide radicals, and hydrogen peroxide, is a key signature of abiotic stress at the molecular level [28]. To resist this harmful effect, plants employ a system that catalyzes the elimination of ROS and fights with oxidative damage via the formation of antioxidant enzymes, such as GST and GPX [2931]. In the present study, the expression of DEGs encoding antioxidant enzymes GST and GPX was greatly larger than that under normal conditions, which has been validated by the qRT-PCR experiment (Table 5). Consistently, many antioxidant enzymes have been shown to be essential for plant survival during adaptive responses to waterlogging or flooding stress [22, 3234]. Therefore, these antioxidant enzymes may be induced by flooding stress and critical for the survival of N. peltata in a submerged environment.

When plants live in a submerged environment, both light and carbon supplies are limited due to the slower diffusion rates in water [35, 36], which decreases plant photosynthesis performance. However, photosynthesis is important for plant survival, which makes the continuation of aerobic respiration through the elevated oxygen concentrations possible. The aerobic respiration is more efficient compared with anaerobic metabolism [37]. For N. peltata, transcriptome analysis revealed that many photosynthesis processes were predominantly enriched in GO terms (Fig. 2) and KEGG pathways (Fig. 3), suggesting the continuation of photosynthesis under flooding stress. Previous studies showed that submerged plants have physiological adaptation mechanisms termed carbon-concentrating mechanisms (CCMs) to produce more CO2 [38], including use of HCO3 and crassulacean acid metabolism (CAM) [39, 40]. Moreover, the use of HCO3 as a CO2 substrate is a common way for most true aquatic plants [41, 42]. Therefore, the continuation of photosynthesis of N. peltata under submergence may be attributed to the mechanism of using HCO3. In addition, species that have leaf gas films or that can produce new leaves under waterlogging often have higher CO2 affinity and higher CO2 concentrations [43, 44]. Therefore, the similar ability of N. peltata to produce new acclimated leaves during submergence might also contribute to the continuation of photosynthesis underwater. Although the activation of photosynthesis was indicated by our petiole transcriptional data, further studies should be added to draw a reliable conclusion on photosynthesis of N. peltata under flooding stress.

Large amounts of energy and carbohydrate are required for the rapid submergence-induced petiole elongation of N. peltata. The initial response to a submergence environment is the induction of anaerobic metabolism [3]. For N. peltata, many DEGs involved in the fermentative pathway were identified and also the glycolysis process and pyruvate metabolism (ko00710) was predominantly in the GO and KEGG enrichment analysis (Figs. 2 and 3). The similar performance of these two enzymes stimulated by waterlogging has also been reported in other plants [23, 45], indicating that the fermentative pathway was likely activated to provide essential energy. However, this universal anoxia metabolism has a side effect: low efficiency [46]. Considering the active and rapid response of N. peltata to flooding stress, there should be alternative metabolic forms with high efficiency to supply the energy required. Moreover, in order to keep the efficiency in glycolysis under oxygen deficiency, it is important to remove accumulated pyruvate. The enzyme AlaAT (Alanine aminotransferase) can transform pyruvate into alanine, and increased alanine accumulation, which is correlated with enhanced activity of AlaAT, under anoxic conditions has been reported in other plants [47, 48]. This process is also likely present in N. peltata because DEGs encoding AlaAT were identified and their up-regulation expression level was validated by the qRT-PCR experiment (Table 5). In addition, a modified tricarboxylic acid (TCA) cycle mode-alanine metabolism was found to be induced by waterlogging in Lotus japonicus [49]. The yield produced from the metabolic reprogramming associated with alanine metabolism (4 ATP) doubled energy produced from glycolysis (2 ATP) [49]. Another DEGs encoding critical enzyme in alanine metabolism, AspAT (Aspartate transaminase), whose expression level was validated by the qRT-PCR, indicating that alanine metabolism was likely induced as energy source by flooding in N. peltata.

It is well known that ethylene is one of the main drivers for depth adaptation in flooding-tolerant plants [13]. Furthermore, group VII ethylene response factors (GVIIERFs) were identified to activate the expression of hypoxia-related genes by an N-end rule pathway under low oxygen conditions [50]. GVIIERF proteins were found in many plants such as Arabidopsis thaliana, Oryza sativa and Rumex palustris (Table 4). In N. peltata, we identified 3 GVIIERFs with the characteristic MCGGAIL amino-terminus in our transcriptome assembly. Indeed, GVIIERFs were found in other aquatic plants as Ranunculus bungei and Utricularia gibba, even can be found in the marine angiosperm Zostera marina by searching its genome sequence. These findings indicate that the oxygen sensing mechanism via GVIIERFs may be conserved in N. peltata and in other higher plants. Notably, we didn’t find GVIIERFs in the differential gene sets of N. peltata, possibly because of our long-duration sampling. Specifically, in this study we were mainly interested in the expression changes of unigenes at the time point of the seventh day, which may be too late to detect the expression changes of GVIIERFs. Indeed, the release of ethylene is an earlier signal in response to flooding [22].


In the present study, comprehensive and valuable genomic resources were built by comparative transcriptome of petioles under normal and flooding conditions of Nymphoides peltata. Our data suggests that two processes rarely occurred in other flooding-tolerant plants, active photosynthesis and alanine metabolism, are likely contributed to the active response of N. peltata to flooding stress. These results deepen our understanding of the genetic basis underlying the response to flooding stress in aquatic plants. The response of plants to abiotic stresses is a complex network functioning with the regulation of stress-related genes [51], therefore, further investigations are still essential to detail the active responses of aquatic plants to flooding stress.


Plant growth and flooding treatment

Young, healthy N. peltata plants were identified by Professor Dan Yu and collected from Liangzi Lake (30°15′29″N, 114°33′30″E) and cultured in glass tanks in a greenhouse at Wuhan University, China. After a week of culturation, the experiment was conducted in May 2014. We chose 20 plants with a height of about 15 cm and transplanted them into two glass tanks, 10 in each tank. We added water in the two tanks to reach water depth of 15 cm and acclimatized plants two days. One tank was used as the untreated group keeping the water depth of 15 cm, the other served as the treated group increasing the water depth to 100 cm. The submerged leaves in the treated tank reached the water surface in the seventh day by elongation of their petioles. Then plant petioles were collected. Petioles from the control tank were used as untreated sample (US) and those from the treated tank were considered as treated sample (TS). Each sample was the mixture of petioles from 5 plants. Samples were frozen in liquid nitrogen and stored at −80 °C prior to RNA extraction. Total RNA of each sample was extracted using HiPure Plant RNA Kits (Magen, China) following the manufacturer’s instruction.

Library construction and transcriptome sequencing

Sequencing was conducted commercially following the manufacturer’s instructions after checking the quality and concentration of RNAs. The procedure was as follows: fragmenting the mRNA after purification, synthesis of the first and second strand cDNA, and adding specific sequence adaptors. After that, cDNA fragments of ~200 bp were chosen to conduct with PCR amplification. The original image data were transferred into raw reads and saved as “fastq” files. The raw reads generated in this study have been deposited in NCBI database under accession number SRA259910.

The adapter sequences and low quality base calls were removed. Firstly, the pooled strategy was employed, the left files from both samples were mixed into left.fq file, and right files from both samples were also mixed into the right.fq file. Transcriptome assembly was accomplished based on the left.fq and right.fq using Trinity with default parameters with the pooled strategy [52]. Then separate assembly of sample US and sample TS were made using Trinity in order to prove our assembly quality. For both strategies, redundancy and over-representation were reduced after Trinity assemble by finding similar sequences using TIGR Gene Indices clustering program (TGICL) and Cluster database at high identity with tolerance (CD-HIT) with minimum 90 and 95% similarity cut off respectively [53, 54]. The TransDecoder was then used to identify the possible coding sequence (CDS) from the assembled sequences.

Functional annotation of pooled unigenes

To gain a better comprehension of the transcriptom information, the pooled assembly unigene sequences was used since some unigenes with low expression levels generated with the pooled strategy cannot be found when using separate sample strategy due to their less reads in sample TS and sample US, which is an advantage of pooled strategy over separate assemble. These unigene sequences were used towards the Phytozome, NCBI NR, COG and Swissprot database, respectively, with a cutoff E-value of 1.0E-5. Only the best match gene ID was assigned to each unigene. Functional annotation by gene ontology analysis was analyzed by Blast2GO software with an E-value ≤1.0E-5. In addition, KEGG pathway analysis was conducted using the KEGG Automatic Annotation Server (KAAS).

Identification of differentially expressed genes

To examine the expression level of each unigene in both samples, the expression of each unigene generated with the pooled strategy was calculated using the Cufflinks program [55]. Moreover, the unigene expression was normalized using the fragments per kilo bases per million reads (FPKM) method described by Mortazavi [56]. Subsequently, the differential gene expression between US and TS were analyzed using the edgeR software [57], with an FDR of 0.05 and |logFC| ≥ 1 as the threshold. DEGs were conducted GO enrichment analysis and KEGG enrichment analysis using R based on hypergeometric distribution. Significantly enriched GO terms and KEGG pathways were identified based on the corrected P-value (P < 0.01 and P < 0.05, respectively).

Transcriptome changes of Nymphoides peltata and five other plants in response to flooding

To have a better understanding of the plants responses under flooding stress and to find the similarities/differences of waterlogging responses in Nymphoides peltata with other land plants, transcriptome comparisons among Nymphoides peltata and five other plants (Arabidopsis thaliana, Oryza sativa, Rumex palustris, Lotus japonicas, Taxodium mucronatum × T. distichum) were conducted. The transcriptome changes of other five plants were from previously published paper. In this method, candidate genes and biological processes involved in flooding response were compared across Nymphoides peltata and five other plants.

Molecular adaptation of aquatic plants compared with non-aquatic plants

A comparative analysis between aquatic plants and non-aquatic plants was carried out to identify commonalities in aquatic plants. The genome of aquatic species Utricularia gibba was downloaded from (v4.1, ID 19475) and transcriptome reads of aquatic Ranunculus bungei were from SRR1822529 under the NCBI. Genomes of five non-aquatic plants including Solanum lycopersicum (assembly SL2.50), Daucus carota (assembly ASM162521v1), Cucumis sativus (assembly ASM407v2), Oryza sativa (assembly Build 4.0) were downloaded from NCBI and the genome sequences of Arabidopsis thaliana were from The Arabidopsis Information Resource (TAIR10). To identify one-to-one orthologous genes, the reciprocal BLAST approach was performed with each species’ CDS sequences and Arabidopsis proteins as query. E-value of 1e-5 was applied and the best hit was retained. The protein IDs of Arabidopsis were used as reference, a total of 5319 one-to-one orthologous in all species were identified and CDS sequences were extracted using a perl script. Each orthologous gene set was aligned used the PRANK program [58]. All of the genes were aligned at the codon level with the following settings: -shortnames +F -termgap -codon -f = fasta. Following alignment, Gblocks program was employed to identify the conserved regions at the codon level [59]. To detect fast evolving genes in aquatic plant group, we estimated a two-ratio branch model allowing different ω values (the ratios of nonsynonymous to synonymous substitution rates) in aquatic and non-aquatic plants and one-ratio model assuming a uniform ω value in all plants, using codeml in the PAML4.8 package [60]. In our analysis, we set the aquatic plants as the foreground, the other branches were set as the background. If a given gene was estimated to have a significantly higher in the foreground branches (ω1) than in the background branches (ω0) (corrected p < 0.05, FDR method), the gene would be considered as a candidate undergoing molecular adaptation.

Validation of quantitative real-time PCR (qRT-PCR)

To validate the Illumina sequencing results, eight unigenes involved in those processes that may be responsible for the adaptation were selected for analysis using qRT-PCR. Gene-specific primers were designed with the Primer Premier 5.0 software, and these primer sequences can be found in Additional file 8. RT-PCR was performed as follows: 95 °C for 3 min, 35 cycles at 95 °C for 30 s, 54–64 °C for 30 s and 72 °C for 20 s; and final extension at 72 °C for 3 min. qRT-PCR was conducted using the SuperReal PreMix Plus (SYBR Green) (Tiangen Biotech, Beijing, China) and a CFX Real-Time PCR System (Bio-Rad, CA, USA). The actin like gene was used as internal reference controls to standardize the results. Statistical analysis was performed using the 2- ΔΔCT method. The final values were presented as means of three independent biological trials.



Alanine aminotransferase


Aspartate transaminase


Cell division cycle protein 48 homolog




Carbon-concentrating mechanisms


Cell division cycle 7


Cluster of orthologous groups


Differentially expressed genes


False discovery rate


Fragments per kilo bases per million reads


Gene ontology






Group VII ethylene response factors


Kyoto encyclopedia of genes and genomes


Nucleobase-ascorbate transporter 2


NCBI non-redundant proteins


Polymerase chain reaction


Photosynthetic electron transport ferredoxin


Photosystem I subunit XI


Quantitative real-time PCR


RNA sequencing


Reactive oxygen species


Sequence read archive


Swiss-prot protein sequence database


Treated sample


Untreated sample


  1. Perata P, Armstrong W, Voesenek LACJ. Plants and flooding stress. New Phytol. 2011;190(2):269–73.

    Article  PubMed  Google Scholar 

  2. Peeters AJM, Cox MCH, Benshop JJ, Vreeburg RAM, Bou J, Voesenek LACJ. Submergence research using Rumex palustris as a model: looking back and going forward. J Exp Bot. 2002;53(368):391–8.

    Article  CAS  PubMed  Google Scholar 

  3. Fukao T, Bailey-Serres J. Plant responses to hypoxia – is survival a balancing act? Trends Plant Sci. 2004;9(9):449–56.

    Article  CAS  PubMed  Google Scholar 

  4. Bailey-Serres J, Voesenek LACJ. Flooding stress: acclimations and genetic diversity. Annu Rev Plant Biol. 2008;59(59):313–39.

    Article  CAS  PubMed  Google Scholar 

  5. Voesenek LACJ, Bailey-Serres J. Flood adaptive traits and processes: an overview. New Phytol. 2015;206(1):57–73.

    Article  CAS  PubMed  Google Scholar 

  6. Vreeburg RAM, Benschop JJ, Peeters AJM, Colmer TD, Ammerlaan AHM, Staal M, Elzenga TM, Staals RHJ, Darley CP, McQueen-Mason SJ, et al. Ethylene regulates fast apoplastic acidification and expansin A transcription during submergence-induced petiole elongation in Rumex palustris. Plant J. 2005;43(4):597–610.

    Article  CAS  PubMed  Google Scholar 

  7. Hossain MA, Uddin SN. Mechanisms of waterlogging tolerance in wheat: morphological and metabolic adaptations under hypoxia or anoxia. Aust J Crop Sci. 2011;5(9):1094–101.

    CAS  Google Scholar 

  8. Nishiuchi S, Yamauchi T, Takahashi H, Kotula L, Nakazono M. Mechanisms for coping with submergence and waterlogging in rice. Rice. 2012;5(1):2.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Yu F, Han X, Geng C, Zhao Y, Zhang Z, Qiu F. Comparative proteomic analysis revealing the complex network associated with waterlogging stress in maize (Zea mays. L) seedling root cells. Proteomics. 2014;15(1):135–47.

    Article  PubMed  Google Scholar 

  10. Hattori Y, Nagai K, Furukawa S, Song XJ, Kawano R, Sakakibara H, Wu J, Matsumoto T, Yoshimura A, Kitano H, et al. The ethylene response factors SNORKEL1 and SNORKEL2 allow rice to adapt to deep water. Nature. 2009;460(7258):1026–30.

    Article  CAS  PubMed  Google Scholar 

  11. Barrett SC, Eckert CG, Husband BC. Evolutionary processes in aquatic plant populations. Aquat Bot. 1993;44:105–45.

    Article  Google Scholar 

  12. Ho TNO, Ornduff RM. Menyanthaceae in Flora of China, vol 16. Beijing and St. Louis: Science Press and Missouri Botanical Garden Press; 1995.

    Google Scholar 

  13. Ridge I, Amarasinghe I. Ethylene and growth control in the fringed waterlily (Nymphoides peltata): Stimulation of cell division and interaction with buoyant tension in petioles. Plant Growth Regul. 1984;2(3):235–49.

    Article  CAS  Google Scholar 

  14. Yu L, Yu D. Differential responses of the floating-leaved aquatic plant Nymphoides peltata to gradual versus rapid increases in water levels. Aquat Bot. 2011;94(2):71–6.

    Article  Google Scholar 

  15. Cookson C, Osborne D. The stimulation of cell extension by ethylene and auxin in aquatic plants. Planta. 1978;144(1):39–47.

    Article  CAS  PubMed  Google Scholar 

  16. Uesugi R, Tani N, Goka K, Nishihiro J, Tsumura Y, Washitani I. Isolation and characterization of highly polymorphic microsatellites in the aquatic plant, Nymphoides peltata (Menyanthaceae). Mol Ecol Notes. 2005;5(2):343–5.

    Article  CAS  Google Scholar 

  17. Liao YY, Yue XL, Guo YH, Gituru WR, Wang QF, Chen JM. Genotypic diversity and genetic structure of populations of the distylous aquatic plant Nymphoides peltata (Menyanthaceae) in China. J Syst Evol. 2013;51(5):536–44.

    Article  Google Scholar 

  18. Marioni JC, Mason CE, Mane SM, Stephens M, Gilad Y. RNA-seq: an assessment of technical reproducibility and comparison with gene expression arrays. Genome Res. 2008;18:1509–17.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Wang Z, Gerstein M, Snyder M. RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet. 2009;10:57–63.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Dong C-H, Li C, Yan X-H, Huang S-M, Huang J-Y, Wang L-J, Guo R-X, Lu G-Y, Zhang X-K, Fang X-P, et al. Gene expression profiling of Sinapis alba leaves under drought stress and rewatering growth conditions with Illumina deep sequencing. Mol Biol Rep. 2012;39(5):5851–7.

    Article  CAS  PubMed  Google Scholar 

  21. Pang T, Ye CY, Xia X, Yin W. De novo sequencing and transcriptome analysis of the desert shrub, Ammopiptanthus mongolicus, during cold acclimation using Illumina/Solexa. BMC Genomics. 2013;14(488):1471–2164.

    Google Scholar 

  22. van Veen H, Mustroph A, Barding GA, Vergeer-van Eijk M, Welschen-Evertman RA, Pedersen O, Visser E, Larive C, Pierik R, Bailey-Serres J, et al. Two Rumex species from contrasting hydrological niches regulate flooding tolerance through distinct mechanisms. Plant Cell. 2013;25(11):4691–707.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Qi B, Yang Y, Yin Y, Xu M, Li H. De novo sequencing, assembly, and analysis of the Taxodium ‘Zhongshansa’ roots and shoots transcriptome in response to short-term waterlogging. BMC Plant Biol. 2014;14(1):201.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Maurino VG, Grube E, Zielinski J, Schild A, Fischer K, Flügge UI. Identification and expression analysis of twelve members of the nucleobase–ascorbate transporter (NAT) gene family in Arabidopsis thaliana. Plant Cell Physiol. 2006;47(10):1381–93.

    Article  CAS  PubMed  Google Scholar 

  25. Satoh M, Tokaji Y, Nagano AJ, Hara-Nishimura I, Hayashi M, Nishimura M, Ohta H, Masuda S. Arabidopsis mutants affecting oxylipin signaling in photo-oxidative stress responses. Plant Physiol Bioch. 2014;81:90–5.

    Article  CAS  Google Scholar 

  26. Bowler C, Montagu MV, Inze D. Superoxide dismutase and stress tolerance. Annu Rev Plant Bio. 1992;43:83–116.

    Article  CAS  Google Scholar 

  27. Gautam P, Lal B, Raja R, Baig MJ, Haldar D, Rath L, Shahid M, Tripathi R, Mohanty S, Bhattacharyya P, et al. Post–flood nitrogen and basal phosphorus management affects survival, metabolic changes and anti-oxidant enzyme activities of submerged rice (Oryza sativa). Funct Plant Biol. 2014;41(12):1284–94.

    Article  CAS  Google Scholar 

  28. Hasanuzzaman M, Hossain M, da Silva J, Fujita M, editors. Plant response and tolerance to abiotic oxidative stress: antioxidant defense is a key factor. In: Crop stress and its management: perspectives and strategies. Berlin: Springer; 2012.

  29. Alscher RG, Erturk N, Heath LS. Role of superoxide dismutases (SODs) in controlling oxidative stress in plants. J Exp Bot. 2002;53(372):1331–41.

    Article  CAS  PubMed  Google Scholar 

  30. Blokhina O, Virolainen E, Fagerstedt KV. Antioxidants, oxidative damage and oxygen deprivation stress: a review. Ann Bot. 2003;91(2):179–94.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Halliwell B. Reactive species and antioxidants. Redox biology is a fundamental theme of aerobic life. Plant phisiol. 2006;141(2):312–22.

    Article  CAS  Google Scholar 

  32. Pourabdal L, Heidary R, Farboodnia T. The effects of flooding stress on induction of oxidative stress and antioxidant enzymes activity in Zea mays L. seedlings. Res J Biol Sci. 2008;3:391–4.

    Google Scholar 

  33. Pucciariello C, Parlanti S, Banti V, Novi G, Perata P. Reactive oxygen species-driven transcription in arabidopsis under oxygen deprivation. Plant Physiol. 2012;159(1):184–96.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Zou X, Tan X, Hu C, Zeng L, Lu G, Fu G, Cheng Y, Zhang X. The transcriptome of Brassica napus L. roots under waterlogging at the seedling stage. In J Mol Sci. 2013;14(2):2637–51.

    Article  CAS  Google Scholar 

  35. Madsen TV, Sand-Jensen K. The interactive effects of light and inorganic carbon on aquatic plant growth. Plant Cell Environ. 1994;17(8):955–62.

    Article  CAS  Google Scholar 

  36. Colmer TD. Long-distance transport of gases in plants: a perspective on internal aeration and radial oxygen loss from roots. Plant Cell Environ. 2003;26(1):17–36.

    Article  CAS  Google Scholar 

  37. Gibbs J, Greenway H. Mechanisms of anoxia tolerance in plants. I. Growth, survival and anaerobic catabolism. Funct Plant Biol. 2003;30:1–47.

    Article  CAS  Google Scholar 

  38. Maberly SC, Madsen TV. Freshwater angiosperm carbon concentrating mechanisms: processes and patterns. Funct Plant Biol. 2002;29:393–405.

    Article  CAS  Google Scholar 

  39. Prins HBA, Elzenga JTM. Bicarbonate utilization: function and mechanism. Aquat Bot. 1989;34:59–83.

    Article  CAS  Google Scholar 

  40. Keeley JE. CAM photosynthesis in submerged aquatic plants. Bot Rev. 1998;64:121–75.

    Article  Google Scholar 

  41. Keeley JE, Rundel PW. Evolution of CAM and C4 Carbon-concentrating mechanisms. Int J Plant Sci. 2003;164:S55–77.

    Article  CAS  Google Scholar 

  42. Colmer TD, Winkel A, Pedersen O. A perspective on underwater photosynthesis in submerged terrestrial wetland plants. Ann Bot. 2011;62:567–88.

    Google Scholar 

  43. Sand-Jensen K, Pedersen MF, Nielsen SL. Photosynthetic use of inorganic carbon among primary and secondary water plants in streams. Freshwater Biol. 1992;27:283–93.

    Article  Google Scholar 

  44. Sand-Jensen K, Frost-Christensen H. Plant growth and photosynthesis in the transition zone between land and stream. Aquat Bot. 1999;63:23–35.

    Article  CAS  Google Scholar 

  45. Christianson JA, Llewellyn DJ, Dennis ES, Wilson IW. Global gene expression responses to waterlogging in roots andleaves of cotton (Gossypium hirsutum L.). Plant Cell Physiol. 2010;51(1):21–37.

    Article  CAS  PubMed  Google Scholar 

  46. Tadege M, Dupuis I, Kuhlemeier C. Ethanolic fermentation: new functions for an old pathway. Trends Plant Sci. 1999;4(8):320–5.

    Article  CAS  PubMed  Google Scholar 

  47. Miyashita Y, Dolferus R, Ismond KP, Good AG. Alanine aminotransferase catalyses the breakdown of alanine after hypoxia in Arabidopsis thaliana. Plant J. 2007;49(6):1108–21.

    Article  CAS  PubMed  Google Scholar 

  48. De Sousa C, Sodek L. Alanine metabolism and alanine aminotransferase activity in soybean (Glycine max) during hypoxia of the root system and subsequent return to normoxia. Environ Exp Bot. 2003;50(1):1–8.

    Article  Google Scholar 

  49. Rocha M, Licausi F, Araújo WL, Nunes-Nesi A, Sodek L, Fernie AR, van Dongen JT. Glycolysis and the tricarboxylic acid cycle are linked by alanine aminotransferase during hypoxia induced by waterlogging of Lotus japonicus. Plant Physiol. 2010;152(3):1501–13.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Gibbs DJ, Lee SC, Isa NM, Gramuglia S, Fukao T, Bassel GW, Holdsworth MJ. Homeostatic response to hypoxia is regulated by the N-end rule pathway in plants. Nature. 2011;479(7373):415–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Mazzucotelli E, Mastrangelo AM, Crosatti C, Guerra D, Stanca AM, Cattivelli L. Abiotic stress response in plants: when post-transcriptional and post-translational regulations control transcription. Plant Sci. 2008;174(4):420–31.

    Article  CAS  Google Scholar 

  52. 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 Biotechno. 2011;29(7):644–52.

    Article  CAS  Google Scholar 

  53. Li W, Godzik A. Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics. 2006;22(13):1658–9.

    Article  CAS  PubMed  Google Scholar 

  54. Pertea G, Huang X, Liang F, Antonescu V, Sultana R, Karamycheva S, Lee Y, White J, Cheung F, Parvizi B, et al. TIGR Gene Indices clustering tools (TGICL): a software system for fast clustering of large EST datasets. Bioinformatics. 2003;19(5):651–2.

    Article  CAS  PubMed  Google Scholar 

  55. Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, Pimentel H, Salzberg SL, Rinn JL, Pachter L. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat Protoc. 2012;7(3):562–78.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008;5:621–8.

    Article  CAS  PubMed  Google Scholar 

  57. 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 

  58. Löytynoja A, Goldman N. An algorithm for progressive multiple alignment of sequences with insertions. Proc Natl Acad Sci U S A. 2005;102(30):10557–62.

    Article  PubMed  PubMed Central  Google Scholar 

  59. Castresana J. Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. Mol Bio Evol. 2000;17(4):540–52.

    Article  CAS  Google Scholar 

  60. Yang Z. PAML 4: phylogenetic analysis by maximum likelihood. Mol Bio Evol. 2007;24(8):1586–91.

    Article  CAS  Google Scholar 

  61. Licausi F, Van Dongen JT, Giuntoli B, Novi G, Santaniello A, Geigenberger P, Perata P. HRE1 and HRE2, two hypoxia- inducible ethylene response factors, affect anaerobic responses in Arabidopsis thaliana. Plant J. 2010;62(2):302–15.

    Article  CAS  PubMed  Google Scholar 

  62. Bailey-Serres J, Fukao T, Gibbs DJ, Holdsworth MJ, Lee SC, Licausi F, van Dongen JT. Making sense of low oxygen sensing. Trends Plant Sci. 2012;17(3):129–38.

    Article  CAS  PubMed  Google Scholar 

  63. Van Der Knaap E, Jagoueix S, Kende H. Expression of an ortholog of replication protein A1 (RPA1) is induced by gibberellin in deepwater rice. Proc Natl Acad Sci U S A. 1997;94(18):9979–83.

    Article  PubMed  PubMed Central  Google Scholar 

  64. Fukao T, Bailey-Serres J. Submergence tolerance conferred by Sub1A is mediated by SLR1 and SLRL1 restriction of gibberellin responses in rice. Proc Natl Acad Sci U S A. 2008;105(43):16814–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  65. Xiong H, Li Y, Yang J, Li Y. Comparative transcriptional profiling of two rice genotypes carrying SUB1A-1 but exhibiting differential tolerance to submergence. Funct Plant Biol. 2012;39(6):449–61.

    Article  CAS  Google Scholar 

  66. Striker GG, Casas C, Manzur ME, Ploschuk RA, Casal JJ. Phenomic networks reveal largely independent root and shoot adjustment in waterlogged plants of Lotus japonicus. Plant Cell Environ. 2014;37(10):2278–93.

    CAS  PubMed  Google Scholar 

Download references


The authors thank Gangzheng Wang for help in qRT-PCR experiment, and Kai Wang and Wencheng Liu for their valuable comments.


This study was supported by grants from the National Natural Science Foundation of China (31270265).

Availability of data and material

The reads available in the NCBI Sequence Read Archive (SRA) database under the PRJNA282005. This Transcriptome Shotgun Assembly project has been deposited at GenBank under the accession GFBT00000000. The version described in this paper is the first version, GFBT01000000.

Authors’ contributions

JW performed the experiment; JW and HZ analyzed the data and wrote the manuscript; HZ, DY and XX designed the project 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

All authors have read and approved the final manuscript.

Ethics approval and consent to participate

No specific permits were required, and no endangered or protected species were involved.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Xinwei Xu.

Additional files

Additional file 1:

The distribution of sequence lengths for the unigenes predicted from the pooling transcriptome assembly of Nymphoides peltata. (PDF 122 kb)

Additional file 2:

GO (Gene ontology) and COG (Cluster of orthologous groups) annotation for the pooling transcriptome assembly of Nymphoides peltata. (a) Functional annotation of the assembled unigenes based on GO category, the functions of the unigenes were divided into three categories. (b) Information of clusters of orthologous groups (COG) classification, the unigenes were mainly clustered into 25 components. (PDF 295 kb)

Additional file 3:

Swissprot annotation for the total unigenes generated from the pooling transcriptome assembly. (XLS 4955 kb)

Additional file 4:

The deduced protein sequences of all predicted genes in the pooling transcriptome assembly of Nymphoides peltata. (TXT 21193 kb)

Additional file 5:

Unigene predictions from the Nymphoides peltata transcriptome assemblies based on three databases: Phytozome, NR (NCBI non-redundant proteins), and Non-coding RNA. (PDF 251 kb)

Additional file 6:

This table provides all differentially expressed genes meeting the threshold (FDR ≤ 0.01) and the GO terms that the differentially expressed genes were enriched. (XLS 207 kb)

Additional file 7:

Evolutionary analysis of molecular adaptation of aquatic plants compared with non-aquatic plants. (PDF 158 kb)

Additional file 8:

Primer sequences used in the qRT-PCR experiment. (PDF 90 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, 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 ( 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

Wu, J., Zhao, HB., Yu, D. et al. Transcriptome profiling of the floating-leaved aquatic plant Nymphoides peltata in response to flooding stress. BMC Genomics 18, 119 (2017).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: