- Research article
- Open Access
Transcriptome profiling of the floating-leaved aquatic plant Nymphoides peltata in response to flooding stress
BMC Genomics volume 18, Article number: 119 (2017)
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 . Plants cannot actively escape a flooding environment due to their immobile character , 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 . Recently, the molecular mechanisms of physiology and metabolic modulation behind these adaptive traits have been investigated in crops and wetland plants [6–10]. However, few such studies have been conducted on true aquatics, floating-leaved or submerged aquatic plants , 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 . 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 . 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 . which is similar with many aquatic or flooding-tolerant species, e.g., Rumex palustris  and rice . 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 [20–23].
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.
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 (http://www.phytozome.net/) 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.
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).
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).
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).
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).
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).
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 (http://www.mobot.org/MOBOT/research/APweb/), 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  (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  (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).
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 . 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 . 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 [29–31]. 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, 32–34]. 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 . 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 , 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 . 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 . 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 . The yield produced from the metabolic reprogramming associated with alanine metabolism (4 ATP) doubled energy produced from glycolysis (2 ATP) . 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 . 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 . 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 .
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 , 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 . 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 . Moreover, the unigene expression was normalized using the fragments per kilo bases per million reads (FPKM) method described by Mortazavi . Subsequently, the differential gene expression between US and TS were analyzed using the edgeR software , 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 https://genomevolution.org/CoGe (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 . 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 . 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 . 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.
Cell division cycle protein 48 homolog
Cell division cycle 7
Cluster of orthologous groups
Differentially expressed genes
False discovery rate
Fragments per kilo bases per million reads
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
Reactive oxygen species
Sequence read archive
Swiss-prot protein sequence database
Perata P, Armstrong W, Voesenek LACJ. Plants and flooding stress. New Phytol. 2011;190(2):269–73.
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.
Fukao T, Bailey-Serres J. Plant responses to hypoxia – is survival a balancing act? Trends Plant Sci. 2004;9(9):449–56.
Bailey-Serres J, Voesenek LACJ. Flooding stress: acclimations and genetic diversity. Annu Rev Plant Biol. 2008;59(59):313–39.
Voesenek LACJ, Bailey-Serres J. Flood adaptive traits and processes: an overview. New Phytol. 2015;206(1):57–73.
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.
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.
Nishiuchi S, Yamauchi T, Takahashi H, Kotula L, Nakazono M. Mechanisms for coping with submergence and waterlogging in rice. Rice. 2012;5(1):2.
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.
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.
Barrett SC, Eckert CG, Husband BC. Evolutionary processes in aquatic plant populations. Aquat Bot. 1993;44:105–45.
Ho TNO, Ornduff RM. Menyanthaceae in Flora of China, vol 16. Beijing and St. Louis: Science Press and Missouri Botanical Garden Press; 1995.
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.
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.
Cookson C, Osborne D. The stimulation of cell extension by ethylene and auxin in aquatic plants. Planta. 1978;144(1):39–47.
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.
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.
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.
Wang Z, Gerstein M, Snyder M. RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet. 2009;10:57–63.
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.
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.
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.
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.
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.
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.
Bowler C, Montagu MV, Inze D. Superoxide dismutase and stress tolerance. Annu Rev Plant Bio. 1992;43:83–116.
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.
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.
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.
Blokhina O, Virolainen E, Fagerstedt KV. Antioxidants, oxidative damage and oxygen deprivation stress: a review. Ann Bot. 2003;91(2):179–94.
Halliwell B. Reactive species and antioxidants. Redox biology is a fundamental theme of aerobic life. Plant phisiol. 2006;141(2):312–22.
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.
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.
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.
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.
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.
Gibbs J, Greenway H. Mechanisms of anoxia tolerance in plants. I. Growth, survival and anaerobic catabolism. Funct Plant Biol. 2003;30:1–47.
Maberly SC, Madsen TV. Freshwater angiosperm carbon concentrating mechanisms: processes and patterns. Funct Plant Biol. 2002;29:393–405.
Prins HBA, Elzenga JTM. Bicarbonate utilization: function and mechanism. Aquat Bot. 1989;34:59–83.
Keeley JE. CAM photosynthesis in submerged aquatic plants. Bot Rev. 1998;64:121–75.
Keeley JE, Rundel PW. Evolution of CAM and C4 Carbon-concentrating mechanisms. Int J Plant Sci. 2003;164:S55–77.
Colmer TD, Winkel A, Pedersen O. A perspective on underwater photosynthesis in submerged terrestrial wetland plants. Ann Bot. 2011;62:567–88.
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.
Sand-Jensen K, Frost-Christensen H. Plant growth and photosynthesis in the transition zone between land and stream. Aquat Bot. 1999;63:23–35.
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.
Tadege M, Dupuis I, Kuhlemeier C. Ethanolic fermentation: new functions for an old pathway. Trends Plant Sci. 1999;4(8):320–5.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008;5:621–8.
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.
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.
Castresana J. Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. Mol Bio Evol. 2000;17(4):540–52.
Yang Z. PAML 4: phylogenetic analysis by maximum likelihood. Mol Bio Evol. 2007;24(8):1586–91.
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.
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.
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.
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.
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.
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.
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.
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.
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.
The distribution of sequence lengths for the unigenes predicted from the pooling transcriptome assembly of Nymphoides peltata. (PDF 122 kb)
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)
Swissprot annotation for the total unigenes generated from the pooling transcriptome assembly. (XLS 4955 kb)
The deduced protein sequences of all predicted genes in the pooling transcriptome assembly of Nymphoides peltata. (TXT 21193 kb)
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)
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)
Evolutionary analysis of molecular adaptation of aquatic plants compared with non-aquatic plants. (PDF 158 kb)
Primer sequences used in the qRT-PCR experiment. (PDF 90 kb)
About this article
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). https://doi.org/10.1186/s12864-017-3515-y
- Nymphoides peltata
- Flooding stress