Skip to main content

ORIGINAL RESEARCH article

Front. Plant Sci., 02 December 2014
Sec. Plant Physiology
This article is part of the Research Topic Environmental control of plant development View all 4 articles

Elucidation of the molecular responses to waterlogging in Jatropha roots by transcriptome profiling

\r\nPiyada JuntawongPiyada Juntawong1Anchalee SirikhachornkitAnchalee Sirikhachornkit1Rachaneeporn PimjanRachaneeporn Pimjan1Chutima SonthirodChutima Sonthirod2Duangjai SangsrakruDuangjai Sangsrakru2Thippawan YoochaThippawan Yoocha2Sithichoke TangphatsornruangSithichoke Tangphatsornruang2Peerasak Srinives*Peerasak Srinives3*
  • 1Special Research Unit in Microalgal Molecular Genetics and Functional Genomics, Department of Genetics, Faculty of Science, Kasetsart University, Bangkok, Thailand
  • 2National Center for Genetic Engineering and Biotechnology, Pathumthani, Thailand
  • 3Department of Agronomy, Faculty of Agriculture at Kamphaeng Saen, Kasetsart University, Nakhon Pathom, Thailand

Jatropha (Jatropha curcas) is a promising oil-seed crop for biodiesel production. However, the species is highly sensitive to waterlogging, which can result in stunted growth and yield loss. To date, the molecular mechanisms underlying the responses to waterlogging in Jatropha remain elusive. Here, the transcriptome adjustment of Jatropha roots to waterlogging was examined by high-throughput RNA-sequencing (RNA-seq). The results indicated that 24 h of waterlogging caused significant changes in mRNA abundance of 1968 genes. Comprehensive gene ontology and functional enrichment analysis of root transcriptome revealed that waterlogging promoted responses to hypoxia and anaerobic respiration. On the other hand, the stress inhibited carbohydrate synthesis, cell wall biogenesis, and growth. The results also highlighted the roles of ethylene, nitrate, and nitric oxide in waterlogging acclimation. In addition, transcriptome profiling identified 85 waterlogging-induced transcription factors including members of AP2/ERF, MYB, and WRKY families implying that reprogramming of gene expression is a vital mechanism for waterlogging acclimation. Comparative analysis of differentially regulated transcripts in response to waterlogging among Arabidopsis, gray poplar, Jatropha, and rice further revealed not only conserved but species-specific regulation. Our findings unraveled the molecular responses to waterlogging in Jatropha and provided new perspectives for developing a waterlogging tolerant cultivar in the future.

Introduction

Waterlogging is an adverse abiotic stress that can heavily damage crop production worldwide. The condition can be defined as the saturation of soils with water. Waterlogging occurs during the heavy rainy season in lowland areas, including Thailand. Due to the limited diffusion of gas under water, waterlogging creates low oxygen (hypoxia) environments in the root areas, causing a shortage of ATP from the inhibition of oxidative phosphorylation. Upon experiencing long-term waterlogging stress (WS), stomatal closure can lead to impaired root hydraulic conductivity, thereby reducing the photosynthetic rate and the nutrient and water uptake of the plant.

Jatropha curcas (common names: Jatropha, physic nut, and purging nut) is a promising crop for the generation of biodiesel. Jatropha seeds contain between 40 and 50% of high quality oil (Basha and Sujatha, 2009). Based on dominant marker analysis among worldwide populations, Jatropha possesses a narrow genetic base (Popluechai et al., 2009). Indigenous to South America, Jatropha is drought tolerant, however, the species is highly sensitive to waterlogging. Gimeno et al. (2012) reported that 10 days of waterlogging caused an approximately 30% reduction in Jatropha biomass. Moreover, a decrease in leaf photosynthetic rate and a decline in carbohydrate concentration in both leaves and roots were observed in waterlogged Jatropha (Gimeno et al., 2012). Despite these lines of evidence, the molecular mechanism underlying the waterlogging response in Jatropha remains unknown.

In recent years, studies in Arabidopsis (Arabidopsis thaliana), gray poplar (Populus × canescens), cotton (Gossypium hirsutum), and rice (Oryza sativa) have revealed that low oxygen can cause drastic changes in transcription, translation, and metabolite levels (Branco-Price et al., 2008; Kreuzwieser et al., 2009; Mustroph et al., 2009; Narsai et al., 2009; Christianson et al., 2010; Lee et al., 2011; Juntawong et al., 2014). The low oxygen condition results in limited ATP production and negatively affects cellular energy status. Therefore, a key feature for the acclimation to low oxygen environment is to activate genes encoding proteins and enzymes for anaerobic fermentation, glycolysis, transcription factors, and signaling pathways in order to allow biological and physiological adjustments to the low oxygen conditions (Bailey-Serres et al., 2012). While several published studies of model plants have provided some fundamental clues, quite a few studies were performed in hypoxic chambers, in which the entire plants were subjected to low oxygen conditions, grown on artificial media, or subjected to complete darkness. These conditions would not imitate the impact observed in plants exposed to soil waterlogging. Recently, high-throughput sequencing technology has been used as a powerful tool for genomic analysis in both model and non-model species. This method can provide quantitative descriptions of gene expression at the genome-scale level, accompanying high accuracy, low background noise, and large dynamic ranges. High-throughput transcriptome analysis could help provide a basic understanding of the molecular responses to waterlogging in Jatropha.

The goal of this study is to provide new insights into the molecular responses of Jatropha to waterlogging. Here, we profiled transcriptome changes in Jatropha roots subjected to 24 h of waterlogging, using high-throughput sequencing by the Ion Proton platform. Bioinformatic analysis of transcriptome data was performed to allow identification and functional annotation of differentially expressed genes (DEGs). Additionally, we comparatively analyzed waterlogging transcriptomes from Arabidopsis, Jatropha, gray poplar and rice to identify conserved and species-specific responses. Lastly, we discussed the use of genetic engineering to target some candidate genes for generation of waterlogging-tolerant Jatropha.

Materials and Methods

Plant Materials and Stress Conditions

Jatropha curcas seeds (cv. “Chai Nat”—a Thai local variety) were germinated in soil containing 50% (v/v) peat moss, 25% (v/v) perlite, and 25% (v/v) coconut fiber with a daily water supply. Plants were grown outdoors between July and August of 2013 and 2014 at Kasetsart University, Bang Khen campus. Thirty day-old, six-leaf-stage plants were used in the waterlogging treatment. In brief, plant pots were placed in plastic containers filled with tap water. The level of water was set at 3 cm above the soil. WS began at mid-day and continued for 24 h. For the control, non-treated plants were placed in a container with no water. Root tissue was harvested at the end of the treatment, immediately placed in liquid nitrogen, ground into a fine powder, and kept at −80°C.

RNA Extraction, Library Preparation, and Sequencing

Total RNA was extracted with TRIzol reagent (Invitrogen), according to the manufacturer's protocol. Total RNA samples were subjected to DNase treatment and RNA cleanup by a GF-1 RNA extraction kit (Vivantis). Poly(A)+ mRNAs were isolated from the total RNAs by using an Absolutely mRNA purification kit (Agilent). Jatropha mRNA (200 ng) was used to construct a sequencing library, by using a Ion Total RNA seq kit (Life Technologies). For each treatment, two independent biological replicates were sequenced on a Ion Proton sequencer (Life Technologies). FASTQ files were obtained with the base caller provided by the instrument and subjected to quality control filtering and 3′end trimming using the FASTX-toolkit (http://hannonlab.cshl.edu/fastx_toolkit/index.html) included in the Torrent Suite™ Software (Life Technologies). The preprocessed reads with a minimum of 35 nt and a quality value of 17, were deposited in the NCBI GEO database under the accession number GSE57428.

Transcriptome Analysis

FASTQ reads were mapped to the J. curcas genome release 4.5, downloaded from http://www.kazusa.or.jp/jatropha/ (Sato et al., 2011; Hirakawa et al., 2012) with TMAP mapping program (https://github.com/iontorrent/TMAP), using the default settings and “mapall” command. A binary format of sequence alignment files (BAM) was generated and subsequently used to create read count tables using HTseq library (Anders et al., 2014) for differential gene expression analysis. Statistical analysis of DEGs was calculated in the R environment (R Development Core Team, 2005), by using a generalized linear model (GLM) approach from the edgeR package (Robinson et al., 2010). Genes with the count per million (CPM) values ≥ 1 in at least two library samples were included in this analysis (Supplementary Table S1). In brief, the calcNormFactors function was applied for data normalization using the trimmed mean of M-values (TMM) method. Subsequently, the GLMTagwiseDisp function was used to generate the estimated tagwise dispersion. Next, the glmFit function was applied to fit the negative binomial GLM for each tag. Finally, the glmLTR function was applied to carry out the likelihood ratio test. Significant DEGs were filtered, based on false discovery values (FDR < 0.05).

Assignment of the DEGs into functional bins and the creation of a mapping file from the Jatropha genome were performed using the Mercator pipeline for automated sequence annotation (Lohse et al., 2013). The inputs for the Mercator were Jatropha protein sequences supplied by http://www.kazusa.or.jp/jatropha/. This analysis was performed with a blast score cutoff > 80. For visualization, the MapMan (Thimm et al., 2004) and the PageMan (Usadel et al., 2006) programs were used.

Gene ontology (GO) enrichment analysis was performed in the R environment. Firstly, the GO terms for each gene were derived from homology searching of plant protein databases using GOANNA with an e-value cutoff < E−10 (McCarthy et al., 2006). Subsequently, a gene annotation file was generated. GO term enrichment analysis in up or down regulated DEGs was performed by the GOHyperGALL function (Horan et al., 2008). Significant GO terms were filtered based on the adjusted p-value < 0.05.

Identification of Jatropha AP2/ERF Genes

For mining of AP2/ERF genes from the Jatropha genome, the hidden Markov model (HMM) profile of the AP2/ERF superfamily was extracted from the Pfam and used to search Jatropha protein sequences.

To perform phylogenetic analysis comparing Jatropha and Arabidopsis AP2/ERF genes, the amino acid sequences of the Arabidopsis AP2/ERF genes were downloaded from the TAIR database (http://www.arabidopsis.org/). A multiple alignment was performed using the MUSCLE alignment in the MEGA5 software (Tamura et al., 2011). A phylogenetic tree of AP2/ERF genes was constructed by the Neighbor Joining method, with a bootstrap number set to 1000 replicates. A list of predicted Jatropha AP2/ERF genes can be found in Supplementary Table S2.

Quantitative Real-Time PCR

cDNA was synthesized from 1.2 μg of total RNAs using oligo(dT)20 and superscript III reverse transcriptase (Invitrogen), according to the manufacturer's protocol. Real-time PCR was assayed in a 20 μL reaction, containing KAPA SYBR FAST qPCR master mix (KAPABIOSYSTEMS) using the Stratagene Mx3000P real-time PCR system (Agilent Technologies). Amplification specificity was validated by melt-curve analysis at the end of each PCR experiment. Relative gene expression was calculated using the comparative ΔΔcT method (Livak and Schmittgen, 2001). Ubiquitin protein (UBC; Jcr4S00238.120) was used as a normalization control. Primer sequences and annealing temperatures were listed in Supplementary Table S3.

Analysis of Publicly Available Microarray Data

Publicly available microarrays were downloaded from the Gene Expression Omnibus (GEO) database (details listed in Supplementary Table S4). Microarray data were normalized using the Robust Multi-chip Average (RMA) method. Differential expression analysis was performed on RMA-normalized data. False discovery rates (FDR) for significant differences between genes in each comparison were generated using p-value distributions (Smyth, 2004). For each comparison, DEGs were selected based on criteria of FDR < 0.05 and the absolute values of log2 fold change (logFC) ≥ 1.

Ortholog Identification

Orthologs among Jatropha, Arabidopsis, rice, and poplar were identified based on protein sequence homology using ORTHOMCL (OMCL) version 1.4 (Li et al., 2003). The inflate number used to perform this analysis was set to 1.2. Protein sequences used in this step were downloaded from Jatropha, Arabidopsis, rice, and poplar genome databases (ftp://ftp.kazusa.or.jp/pub/jatropha/JAT_r4.5.aa.gz, ftp://ftp.arabidopsis.org/home/tair/Sequences/blast_datasets/TAIR10_blastsets/TAIR10_pep_20110103_representative_gene_model_updated, ftp://ftp.jgi-psf.org/pub/compgen/phytozome/v9.0/Osativa/annotation/Osativa_204_protein_primaryTranscriptOnly.fa.gz, and ftp://ftp.jgi-psf.org/pub/compgen/phytozome/v9.0/Ptrichocarpa/annotation/Ptrichocarpa_210_protein_primaryTranscriptOnly.fa.gz, respectively). A list of OMCL clusters can be found in Supplementary Table S5.

Analysis of Leaf Chlorophyll Content

Chlorophyll content was measured using the atLEAF+ chlorophyll meter (FT Green LLC, Wilmington, DE). The youngest fully expanded leaves of approximately 30 day-old plants were measured three times and the averages were used in subsequence analysis. Six plants were analyzed for each time point. The total chlorophyll content of the leaves was obtained by converting the atLEAF+ values in SPAD using an online webtool: http://www.atleaf.com/SPAD.aspx.

Carbohydrate Assay

One hundred mg of frozen root tissue was used to quantify the total carbohydrate content using a method described by Sadasivam and Manickam (2007). Non-structural carbohydrates were extracted and hydrolyzed by adding 5 mL of 2.5 N HCl and incubated in a boiling water bath for 3 h. The extract was neutralized by adding 0.75 g of Na2CO3. The anthrone method was used to determine total carbohydrate content relative to a standard series of glucose. In brief, the extract (300 μL) and distilled water (700 μL) were mixed with 4 mL of 0.14% (w/v) anthrone solution in 95% H2SO4, incubated in a boiling water bath for 8 min, and rapidly cooled on ice. The absorbance was quantified at 630 nm.

Results

Waterlogging Stress Triggered Transcriptome Readjustment in Jatropha Roots

In our observation, we found that long-term waterlogging in young Jatropha seedlings resulted in leaf chlorosis (Supplementary Figure S1A) and reduction of total root carbohydrate content (Supplementary Figure S1B). To determine the molecular responses of Jatropha to waterlogging, six-leaf-stage plants (30-days old) were tested by applying 24 h of waterlogging. First, we examined the response of Jatropha leaves and roots to waterlogging by using the semi-quantitative reverse-transcription PCR. We found waterlogging induced the expression of low oxygen responsive marker genes, alcohol dehydrogenases (ADHs) and pyruvate decarboxylase (PDC) in roots, but not in leaves (Supplementary Figure S2). Therefore, our transcriptome analysis was focused only on the response in root organs.

To compare the effects of waterlogging on the Jatropha roots, we quantitatively profiled transcriptome from WS and non-stress (NS) samples. Over 4 million reads were obtained for each of the two independent biological replicates. The majority of reads (93–96%) mapped to the Jatropha genome (Supplementary Figure S3A). Biological replicate transcriptome data were highly correlated as shown by Pearson's correlation coefficient of CPM values (r = 0.93 and 0.85 from NS and WS libraries, respectively; Supplementary Figure S3B). A multi-dimension scaling plot of the gene expression data demonstrated that samples were clearly separated by treatments (Supplementary Figure S3C).

Transcriptome analysis identified 1968 DEGs with significant changes in expression evaluated by the false discovery rate (FDR < 0.05) (Supplementary Table S1). Of these, 931 genes (47%) were up-regulated and 1037 genes (53%) were down-regulated. Assignment of the DEGs to functional bins, using a well-established Mercator pipeline (Lohse et al., 2013), categorized the DEGs into 35 functional bins (Figure 1; Supplementary Table S1). Approximately 28% of the DEGs fell into a “not-assigned” bin. Most of the genes in this bin were annotated as “short,” “partial,” and “transposable element (TE)” by Sato et al. (2011) and Hirakawa et al. (2012). Nevertheless, the genes that related to cell wall, secondary metabolism, hormone metabolism, stress, RNA, protein, signaling, and transport were abundant in the DEG dataset.

FIGURE 1
www.frontiersin.org

Figure 1. Functional classification of waterlogging-responsive DEGs. DEGs (1968 genes) were assigned into 35 functional bins using the Mercator annotation pipeline.

To further evaluate the response to waterlogging, we considered the biological function of the DEGs that were co-regulated in response to waterlogging. GO enrichment analysis was performed with the up-regulated or down-regulated DEGs. Co-regulation of DEGs with similar functions was observed (Figure 2A; Supplementary Table S1). The up-regulated DEGs were involved in the response to stress (adjusted p-value: 1.45E-08), response to hypoxia (adjusted p-value: 2.37E-03), and response to ethylene (adjusted p-value: 7.34E-03). In addition, we also observed enrichment of genes encoding transcription factors (adjusted p-value: 4.08E-02) following WS. The down-regulated DEGs were enriched for cell wall organization or biogenesis (adjusted p-value: 9.03E-10), cellular carbohydrate biosynthetic process (adjusted p-value: 6.21E-04), secondary metabolite biosynthetic process (adjusted p-value: 1.20E-03), and growth (adjusted p-value: P1.96E-03).

FIGURE 2
www.frontiersin.org

Figure 2. Waterlogging reprograming transcriptome of Jatropha roots. (A) Selected gene ontology (GO) categories (p-values calculated by GOHyperG) from up-regulated (Up) and down-regulated (Down) groups are shown. (B) PAGEMAN enrichment analysis with significant up/down-regulation (Wilcoxon rank sum test; p-value < 0.05).

To confirm the results from the GO enrichment analysis, we performed functional enrichment analysis using PAGEMAN tools (Usadel et al., 2006) (Figure 2B). The PAGEMAN results were highly correlated with those derived from the GO enrichment analysis. In fact, results from PAGEMAN further demonstrated that DNA synthesis/chromatin structure related genes were down-regulated (p-value: 3.15E-03). On the other hand, the expression of the genes involved in protein degradation, calcium signaling, and sugar transport were enhanced (p-value: 4.68E-02, 9.59E-03, and 4.07E-03, respectively). Altogether, the GO enrichment and PAGEMAN results confirmed that waterlogging caused global gene expression reprogramming in Jatropha roots.

Differential Regulation of Genes Associated with Anaerobic Fermentation and Glycolysis

Because the low oxygen conditions could promote glycolysis and fermentation, we further examined regulation of the DEGs encoding central carbon metabolism enzymes (Figure 3A; Supplementary Table S1). We found that waterlogging promoted the expression of the genes involved in sugar and starch cleavage [e.g., alpha-amylase and sucrose synthase (SUSY)]. In contrast, waterlogging inhibited the accumulation of mRNAs encoding for enzymes controlling starch synthesis [e.g., ADP-glucose pyrophosphorylase (AGPase) and starch synthase]. Moreover, the expression of the genes associated with glycolysis was up-regulated by waterlogging [e.g., phosphofructokinase (PFK), pyruvate kinase (PK), and glyceraldehyde phosphate dehydrogenase (GAPDH)]. As expected, an increase in the expression of fermentative genes (PDC and ADH) was evident. Furthermore, the expression of alanine transaminase (AlaAT), an enzyme responsible for conversion of pyruvate to alanine, was found up-regulated in our analysis. Overall the results revealed that waterlogging promoted catabolism of carbohydrates and a switch from oxidative to anaerobic respiration in Jatropha roots.

FIGURE 3
www.frontiersin.org

Figure 3. Waterlogging caused differential expression of transcripts encoding proteins involved in energy, nitrate, and ethylene metabolisms. (A) Starch metabolism (Starch), glycolysis, and fermentation (B) Nitrate metabolism (C) Ethylene production and perception. A plus sign represents induction and a minus sign represents reduction. Numbers in brackets represent numbers of DEGs found in this analysis. Data can be found in Supplementary Table S1.

Differential Regulation of Genes Involved in Nitrate Metabolism and Nitric Oxide Production

Changes in nitrate metabolism were observed upon waterlogging in several plant species (Horchani et al., 2010; Oliveira et al., 2013; Oliveira and Sodek, 2013). Therefore, we evaluated the regulation of the genes involved in nitrate metabolism and found that a set of genes was up-regulated in response to waterlogging in Jatropha roots (Figure 3B; Supplementary Table S1). These included genes encoding nitrate reductase (NR), a key enzyme responsible for conversion of nitrate (NO3) to nitrite (NO2) and nitric oxide (NO), nitrite reductase (NIR), an enzyme responsible for reduction of NO2 to ammonium and conversion of NO2 to NO, and NO3 transporter. Besides this, waterlogging induced the expression of non-symbiotic hemoglobins (nsHbs; Supplementary Table S1), which have been found to be associated with NO removal (Heckmann et al., 2006). These results suggest the regulation of NO3 metabolism and that the modulation of endogenous NO levels might be important for waterlogging acclimation in Jatropha.

Differential Regulation of Genes Related to Ethylene Synthesis and Perception

Accumulation of a plant stress hormone, ethylene, in response to low oxygen has been reported (Voesenek et al., 1993; Vriezen et al., 1999; Rzewuski et al., 2007; Hinz et al., 2010). Our analysis demonstrated the differential regulation of DEGs encoding for proteins functioning in ethylene synthesis and response (Figure 3C; Supplementary Table S1). The expression of two key-genes controlling ethylene synthesis, 1-aminocyclopropane-1-carboxylate (ACC) synthase (ACS) and ACC oxidase (ACO), was up-regulated by waterlogging. In addition, we found that waterlogging induced accumulation of transcripts encoding for ethylene receptor (ETR). Interestingly, several genes encoding for ERF (ethylene responsive factors) transcription factors (TFs) were up-regulated by waterlogging (Figure 4A). These marks highlighted the roles of ethylene in acclimation to waterlogging in Jatropha.

FIGURE 4
www.frontiersin.org

Figure 4. The group-VII ERFs of Jatropha. (A) Graphical representation of waterlogging-regulated transcription factors based on their assigned protein families. “Up” and “Down” represent up-regulation and down-regulation in this analysis. (B) Phylogenetic tree of group-VII ERF proteins. The full-length protein sequences were analyzed with a Neighbor-joining method. Numbers above branches represent the bootstrapped value from 1000 replicates. (C) Multiple sequence analysis of N-terminal sequences of group-VII ERFs. Bold letters indicate a conserved motif at the N-terminus initiated with Met1-Cys2, as identified by multiple sequence alignment. Asterisks and semi-colons indicate identical and conserved substitution, respectively.

Differential Regulation of Genes Encoding TFs

Since low oxygen induced TFs could mediate the expression of anaerobic responsive genes, we analyzed for TFs that were differentially regulated by waterlogging. One-hundred and fifteen waterlogging-regulated TFs were classified based on their assigned protein families (Figure 4A; Supplementary Table S1). Of these, 85 TFs were up-regulated and 30 TFs were down-regulated. Remarkably, TFs accounted for nearly six percent of the DEGs analyzed in this study. ERFs (21 genes), MYBs (9 genes), and WRKYs (8 genes) were mostly pronounced from the up-regulated group, whereas Homeobox TFs (HBs; 5 genes) were mostly noticeable from the down-regulated group. Recent studies on Arabidopsis demonstrated that group-VII ERFs are modulators of the anaerobic response under low oxygen conditions (Gibbs et al., 2011; Licausi et al., 2011). Therefore, we examined the expression of the group-VII ERFs in Jatropha. Firstly, we searched the Jatropha genome using the HMM profile of the Apetalla2 (AP2)/ERF domain in order to identify all AP2/ERF genes. As a result, 133 putative AP2/ERF genes were identified. Secondly, we performed sequence alignment and phylogenetic analyses to compare Arabidopsis and Jatropha ERFs. Finally, 14 groups of AP2/ERF genes, including 11 ERF subgroups, AP2, RAV, and soloist were successfully identified (Supplementary Table S2). This analysis revealed that the Jatropha group-VII ERFs comprised three members, all possessing a conserved N-terminal domain [NH2-MCGGAII(A/S)D] (Figures 4B,C). Remarkably, two of them (Jcr4S00982.160 and Jcr4S01651.60; JcERFVII-2 and JcERFVII-3, respectively) were up-regulated in our transcriptome analysis (Supplementary Table S1). Together, these data implicated the essence of transcriptional regulation in response to waterlogging.

Validation of Differentially Regulated Genes by Quantitative Real-Time PCR

To verify the transcriptome data, 10 genes were selected for analysis by quantitative real-time PCR to determine their relative expression in response to WS. These selected genes included 6 waterlogging-induced DEGs [ADHa, ADHb, JcERFVII-2, JcERFVII-3, NR, and trihelix TF (trihelix)], two non-DEGs with limited changes in response to waterlogging (actin and JcERFVII-1), one waterlogging-repressed DEG (proline-rich), and one control for fold change analysis (UBC). Real-time PCR results of the same samples used for RNA-seq clearly confirmed that waterlogging induced the accumulation of ADHa, ADHb, JcERFVII-2, JcERFVII-3, NR, and trihelix but not actin and JcERFVII-1 (Figure 5). By contrast, mRNA accumulation of proline-rich, a gene encoding a cell-wall protein, was clearly reduced upon waterlogging (Figure 5). These data provided additional support that transcriptome analysis by RNA-seq could be used to quantify genome-scale gene expression in the roots of waterlogged Jatropha.

FIGURE 5
www.frontiersin.org

Figure 5. Quantitative real-time PCR validation of transcriptome data for selected genes. Relative expression was normalized to the abundance of Ubiquitin (UBC; Jcr4S00238.120). Data represent mean ± SE (n = 3). *p-value < 0.01 (t-test).

Comparative Analysis of Transcriptome Response to Waterlogging in Arabidopsis, Jatropha, Gray Poplar and Rice

The variation in tolerance to the low oxygen conditions among plant species has been reported (Bailey-Serres and Voesenek, 2008). In order to identify common or specific changes in transcriptional regulation between waterlogging intolerant Arabidopsis and Jatropha and tolerant gray poplar and rice, we first analyzed public microarray data to identify DEGs in response to waterlogging or submergence in Arabidopsis, gray poplar and rice (Supplementary Table S4). Additionally, we performed a PAGEMAN over-representation analysis (ORA) using Fisher's exact test and a cut-off value of two to obtain an overview of comparative comparison. PAGEMAN ORA analysis allowed comparison of genes with similar functions. We observed an over-represented induction of genes encoding for anaerobic fermentation in Arabidopsis, Jatropha, and rice (Figure 6; Supplementary Table S5). In addition, AP2/ERF DNA binding TFs were up-regulated in all species. The down-regulation of genes involved in cell wall processes and secondary metabolism was common in both tolerant and intolerant species (Figure 6; Supplementary Table S5).

FIGURE 6
www.frontiersin.org

Figure 6. Overview of transcriptome response for selected functional categories to waterlogging or submergence in Arabidopsis (ATH), Jatropha (JCA), gray poplar (PTR), and rice (OSA). PAGEMAN analysis of the gene expression data (|log2 fold change| > 1; FDR < 0.05). Statistical analysis of over-represented functional bins was carried out using Fisher method. Z-scores indicate over/under representation (Numbers indicate z-scores; Green, over-represented; Red, under-represented). Data used to generate this figure can be found in Supplementary Table S5.

Additionally, we attempted to compare changes in expression of orthologs in the genomes of Arabidopsis, Jatropha, gray poplar and rice. OrthoMCL (OMCL) cluster analysis using a Markov Cluster algorithm was applied in order to group orthologs (Li et al., 2003). This analysis identified 19,103 OMCL clusters (Supplementary Table S5). Over 7000 OMCL clusters were differentially regulated upon waterlogging (Supplementary Figure S4). We found that using a list of 49 core hypoxia genes reported to be up-regulated in response to low oxygen stress in shoots and roots of Arabidopsis seedlings (Mustroph et al., 2009), we were able to identify 37 OMCL clusters among Arabidopsis, poplar, Jatropha, and rice (Supplementary Figure S5). Of these, 30 OMCL clusters were induced in at least two species. OMCL236 (wound-responsive proteins), OMCL442 (pyruvate decarboxylase), OMCL450 (alcohol dehydrogenase), and OMCL1101 (ethylene receptor) were induced in all species. We further compared the expression of OMCL clusters involved in carbohydrate cleavage, glycolysis, and fermentation and, in general, found up-regulation of OMCL clusters in these categories (Supplementary Figure S6). All together, the results suggested that, regardless of the variation in waterlogging tolerance among plant species, common responses could still be identified.

To examine whether there is a link between tolerance levels and specific molecular mechanisms, we further inspected the PAGEMAN ORA results. We found that multiple processes were specifically regulated in Jatropha. For example, the genes involved in non-mevalonate pathway biosynthesis of isoprenoids, response to stress including heat and touch/wounding, and calcium signaling were up-regulated, while lipid metabolism, cold stress response, and cell organization were down-regulated (Figure 6). In rice, we found several protein kinases genes were specifically up-regulated (Figure 6). Interestingly, we also found that the down-regulation of genes controlling NO3 metabolism, including NR, was specific to rice (Figure 6). An OMCL cluster analysis also confirmed that, while NRs were specifically down-regulated in rice, this was not the case in Arabidopsis and Jatropha (Figure 7). Additionally, the analysis also showed that up-regulation of NIR was specific to Jatropha. Evidently, the up-regulation of class I nsHbs (OMCL3239), putative NO scavengers, was found in Arabidopsis, Jatropha, and rice. However, the induction levels were much higher in Jatropha than in Arabidopsis and rice (Figure 7).

FIGURE 7
www.frontiersin.org

Figure 7. Response of OMCL clusters with genes involved in nitrogen metabolism and NO production to waterlogging. Yellow, Up-regulated DEGs; Blue, Down-regulated DEGs; White, Orthologs not differentially expressed. Numbers represent mean log2 fold change values.

Comparative analysis of TFs in response to waterlogging revealed that the expression pattern of TFs was generally conserved in all plant species examined (Supplementary Table S5). OMCL clusters of ERFs, bZIPs, NACs and WRKYs were induced by waterlogging. Interestingly, OMCL clusters of trihelix TFs were found to be up-regulated in Arabidopsis, Jatropha, and gray poplar but not in rice (Supplementary Table S5). For protein kinases, OMCL clusters of leucine rich repeat kinases XI and XII and wall-associated kinases were induced in much greater quantities in rice (Supplementary Table S5). Examination of OMCL clusters of heat stress and redox regulation revealed that the induction of alternative oxidase was found only in Jatropha. Up-regulation of heat-shock proteins (HSPs), were commonly found in both intolerant and tolerant species (Supplementary Table S5). Interestingly, OMCL clusters of peroxidases were greatly induced in rice (Supplementary Table S5).

Discussion

In this study, high-throughput RNA-sequencing technology was employed to compare differential gene expression profiles of Jatropha roots subjected to 24 h WS. Our study provides a new insight into understanding the molecular mechanisms underlying the response to waterlogging in Jatropha. WS reprogrammed expression of 1968 DEGs; seven of these were further confirmed to be differentially expressed by quantitative real-time PCR (Figure 5). By applying bioinformatic analyses for non-model species, our study demonstrated that a number of cellular processes in Jatropha roots were affected by WS. These included anaerobic fermentation, carbohydrate metabolism, glycolysis, ethylene synthesis and perception, NO3 metabolism, NO production, regulation of transcription, protein degradation, transport, and signaling, suggesting that diverse physiological processes were affected by WS.

Plants subjected to low oxygen conditions shift their metabolism from oxidative phosphorylation to anaerobic fermentation to maintain ATP production (Bailey-Serres et al., 2012). Our results supported a longstanding notion that waterlogging promotes anaerobic respiration as observed by the up-regulation of DEGs encoding enzymes in glycolysis and fermentation (Figure 3; Supplementary Table S1). Remarkably, the expression of Alternative oxidase (AOX) was specifically enhanced in roots of waterlogged Jatropha (Jcr4S02312.80, Supplementary Table S1) implying that AOX may function as an alternative to cytochrome oxidase under low oxygen conditions. The increase in AOX gene expression could prevent reactive oxygen species (ROS) formation from the over-reduction of the ubiquinone pool. Plants can initiate additional responses to low oxygen conditions, including down-regulation of energy consuming processes such as storage metabolism (Geigenberger et al., 2000) and switching from invertase to sucrose synthase for hydrolysis of sucrose (e.g., Guglielminetti et al., 1995, 1997; Biemelt et al., 1999; Zeng et al., 1999; Albrecht and Mustroph, 2003; Bologa et al., 2003; Bieniawska et al., 2007). Here, we observed in the roots of waterlogged Jatropha, the down-regulation of genes involved in the synthesis of the cell wall, DNA, secondary metabolites, and starch (Figure 2) and the up-regulation of sucrose synthase (Figure 3A). These results, together with the finding that waterlogging resulted in the decrease in total root carbohydrate (Supplementary Figure 1B) provided a strong conclusion that Jatropha roots responded to waterlogging via the regulation of energy consumption and production.

Low oxygen conditions in plants promote the utilization of NO3 and the production of NO to facilitate anaerobic survival (Horchani et al., 2010; Oliveira et al., 2013; Oliveira and Sodek, 2013). Alternative to anaerobic fermentation, NO3 could be used as NADH acceptors. This allows NAD+ to be reused in glycolysis. It has been observed that longer periods of low oxygen survival can be achieved when NO3 is provided (Allegre et al., 2004; Horchani et al., 2010). In waterlogged Jatropha roots, the expression of NR, NiR, and nitrate transporter was induced (Figure 4B; Supplementary Table S1), emphasizing the roles of NO3 in waterlogging acclimation. Future studies should focus more on investigating the significant of NR, NiR, and NO3 on tolerance to waterlogging in Jatropha.

Accumulation of the gas hormone ethylene is extremely important for the induction of plant responses to low oxygen. Since waterlogging inhibits gas diffusion, ethylene can be trapped in the waterlogged plant organs. Our transcriptome data revealed that root waterlogging promoted mRNA accumulation of two key enzymes in ethylene biosynthesis, ACS and ACO (Figure 3C; Supplementary Table S1). Concomitantly, the expression of ETR, was enhanced in waterlogged Jatropha roots. These findings suggested that waterlogging might increase ethylene synthesis and perception in Jatropha roots. Our study recognized three group-VII ERFs from the Jatropha genome. Two of these, designated JcERFVII-2 and JcERFVII-3, were noticeably induced in response to waterlogging (Figures 4, 5; Supplementary Table S1). Previous studies identified the two group-VII ERF genes, Snorkel and Sub1, as key players in the submergence response in deepwater rice and lowland rice varieties, respectively, (Xu et al., 2006; Hattori et al., 2009). These two rice varieties utilize two contrasting strategies that allow adaptive growth responses to submergence. Snorkel promotes the internode elongation in the deepwater rice, whereas, Sub1 restricts the shoot elongation in the lowland rice. In Arabidopsis, studies have also shown that the group-VII ERFs function in the low-oxygen response (Hinz et al., 2010; Gibbs et al., 2011; Licausi et al., 2011). Taken together, we propose that the group-VII ERFs are promising candidates for engineering of waterlogging tolerant Jatropha. In support of this idea, overexpression of the Arabidopsis group-VII ERFs (HRE1, HRE2, RAP2.2, and RAP2.12) significantly improved low oxygen survival by promoting expression of the genes involved in low oxygen adaptation (Hinz et al., 2010; Gibbs et al., 2011; Licausi et al., 2011).

Enhanced NO production by low oxygen conditions has been documented in plants (Rockel et al., 2002; Dordas et al., 2003). NO is believed to be a signaling molecule. As previously mentioned, NO can be synthesized from NO3 via NR and NiR (Stohr et al., 2001). In Arabidopsis, it was reported that accumulation of NO could induce AOX expression via inhibition of aconitase resulting in accumulation of citrate and a shift of metabolism toward nitrogen assimilation under hypoxia (Gupta et al., 2012). A recent study revealed that NO could control the stability of Arabidopsis group-VII ERFs via proteolytic modulation (Gibbs et al., 2014). We speculate that enhanced NO accumulation in Jatropha roots under low oxygen conditions could destabilize the group-VII ERF proteins and, therefore, decrease low oxygen tolerance. Modulation of NO accumulation via nsHbs in roots under low oxygen conditions has been reported (Dordas et al., 2003, 2004; Perazzolli et al., 2004; Hebelstrup et al., 2006). Our study demonstrated that waterlogging could induce the expression of nsHbs in Jatropha roots (Figure 7; Supplementary Table S1). The function of nsHbs possibly involves fine-tuning the accumulation of NO in waterlogged roots of Jatropha. Future research to improve waterlogging tolerance should emphasize genetic engineering to modulate NO levels.

Nearly six percent of DEGs encoding for TFs was identified in this study, implying that transcriptional regulation plays an important role in the waterlogging response. In agreement with the finding that ethylene synthesis and perception were activated, this study recognized AP2/ERF TFs as the most pronounced TFs in response to waterlogging (Figure 4A). Waterlogging induced AP2/ERFs included two of group-VII ERFs and an additional 19 ERF, whose function might be associated with waterlogging acclimation (Supplementary Table S1). The ERF IX group was the most represented waterlogging-induced ERFs in this study (Supplementary Table S1). The function of the group-IX ERFs was associated with hormonal response, such as the responses to ethylene (Wang et al., 2007) or ethylene in combination with jasmonate (Champion et al., 2009). The induction of JcERFX-6 (Jcr4S03895.40; OMCL 2117), an ortholog of Arabidopsis RAP2.6L, was evident in this analysis (Supplementary Table S1). Overexpression of RAP2.6L in Arabidopsis improved the waterlogging response by delaying waterlogging-induced premature senescence (Liu et al., 2012). In this study, the expression of trihelix TF (Jcr4S02762.30: OMCL11085), an ortholog of Arabidopsis HYPOXIA RESPONSE ATTENUATOR1 (HRA1), was found to be strongly up-regulated by waterlogging (Supplementary Table S1). Recently, Giuntoli et al. (2014) demonstrated that HRA1 negatively regulates RAP2.12 through protein–protein interaction and proposed that the modulation of the aerobic response by group-VII ERFs and HRA1 is important to control plant responses to low oxygen stress. In addition, our analysis demonstrated that waterlogging enhanced the expression of several members of MYBs and WRKYs. An in-depth functional analysis of waterlogging-induced TFs is required to further identify candidate genes for growth improvement of waterlogged Jatropha.

Previous studies reported the use of comparative transcriptome comparison for identification of conserve and species-specific responses to low oxygen stress (Mustroph et al., 2010; Narsai et al., 2011; Narsai and Whelan, 2013). In this study, comparative transcriptome analysis between Jatropha and other plant species identified a common response to waterlogging including the up-regulation of carbohydrate cleavage, glycolysis and fermentative genes and the down-regulation of genes involved in cell-wall and secondary metabolite biosynthesis (Figure 6; Supplementary Table S5). These findings also confirmed previous reports that the alteration of carbohydrate metabolism appeared to be a common response to low oxygen stress in plants. Additionally, Arabidopsis, gray poplar and Jatropha showed the up-regulation of genes involved in ethylene production and signaling (Figure 6; Supplementary Table S5), which were commonly found in response to low oxygen in several plants. The up-regulation of AP2/ERF genes by waterlogging was commonly identified in all species examined (Figure 6; Supplementary Table S5). These results strongly emphasized the roles of AP2/ERF genes in low oxygen adaptation. While comparative analysis exposed the common response to waterlogging in plants, the analysis showed that various processes were specifically regulated in each species. Contrary to what has been observed in waterlogged Jatropha, the expression of NR was down-regulated in submerged rice (Figure 7) suggesting the down-regulation of the NR gene may be a characteristic of waterlogging-tolerant species.

Global warming results in extreme climates, such as drought, flooding and heat stress. Improvement of combined stress tolerance could benefit crop production. It is widely known that stresses including drought, flooding, and high temperatures, could trigger accumulation of ROS, which causes oxidative damage of cellular components. On the other hand, ROS serves as a critical signaling molecule in the oxidative stress response. Pucciariello et al. (2012) demonstrated a connection in Arabidopsis between redox sensing and low oxygen response of seedlings. Evidence from the study of Fukao et al. (2011) showed that submergence-tolerant rice overexpressing Sub1A gene exhibited drought-tolerance phenotypes through oxidative stress regulation. Additionally, Loreti et al. (2005) found that several HSPs were induced by low oxygen stress and heat pretreatment enhanced anoxia tolerance of seedlings. Later on, Banti et al. (2008) demonstrated a correlation between the up-regulation of HSP mRNAs and enhanced low oxygen tolerance in Arabidopsis. These findings highlighted the roles of HSPs and ROSs in plant low oxygen acclimation. In this study, we found that waterlogging induced several OMCL clusters of HSPs and ROS enzymes in the roots of Jatropha (Supplementary Table S5). Taken together, it is possible that genetically engineered, waterlogging-tolerant crops might confer multiple stress tolerance phenotypes.

In summary, Jatropha orchestrates a complex transcriptional adjustment in response to WS that yields changes in its root physiological response to WS. This study highlights several possibilities for future investigation, including the roles of specific pathways and genes involved in the waterlogging response. Thus, a detailed characterization of individual genes must be carried out as the first step to understand their specific functions. Due to the high similarity among global J. curcas accessions, as reported by Popluechai et al. (2009), the improvement of waterlogging-tolerant Jatropha through traditional breeding, can be difficult to achieve. With biotechnology advances, genetic engineering and gene transformation into Jatropha are now feasible (Mukherjee et al., 2011). A recently developed protocol by Jaganath et al. (2014) enables Agrobacterium-mediated in planta transformation into Jatropha that does not require in vitro-based multiplication of transformed plants. Clearly, this method has the potential to facilitate the genetic modification of Jatropha. Successfully validated candidate genes could be targeted for the engineering of waterlogging-tolerant Jatropha.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We thank Dr. Kittisak Yokthongwattana for the equipment and Dr. Anongpat Suttangkakul, Dr. Supachai Vuttipongchaikij, and Ms. Pattralak Songserm for technical assistance. This work was funded by the Coordinating Center for Thai Government Science and Technology Scholarship Students/Ministry of Science and Technology (SCH-NR2012-112), Faculty of Science, Kasetsart University, and Thailand Research Fund (MRG5680032) to Piyada Juntawong and NSTDA Chair Professor Grant (P-11-00599) to Peerasak Srinives.

Supplementary Material

The Supplementary Material for this article can be found online at: http://www.frontiersin.org/journal/10.3389/fpls.2014.00658/abstract

Supplementary Figure S1. Physiological responses of Jatropha seedlings to waterlogging. (A) Total leaf chlorophyll content. Data represent mean ± SE (n = 6). (B) Total root carbohydrate content. Data represent mean ± SE (n = 3). Letters represent significant differences calculated by One-Way ANOVA (p-value < 0.05).

Supplementary Figure S2. Change in the expression of selected fermentative genes in response to waterlogging in leaves and roots of Jatropha. Semi-quantitative PCR analysis of fermentative genes (ADHs and PDC) from 24-h waterlogged (WS) and non-stressed (NS) Jatropha roots and leaves.

Supplementary Figure S3. Statistic and quality control of RNA-seq data. (A) Mapping statistics (B) Pearson's correlation derived from CPM values. (C) A multi-dimension scaling plot of gene expression data. NSR, non-stressed root sample; WSR, waterlogged root sample. Numbers represent numbers of replicates.

Supplementary Figure S4. Differentially regulated OMCL clusters grouped by the number of organisms per cluster.

Supplementary Figure S5. Waterlogging responsive OMCL clusters containing “49 core hypoxia” genes (Mustroph et al., 2009, 2010). Arabidopsis (ATH), Jatropha (JCA), gray poplar (PTR), and rice (OSA). Numbers indicated numbers of induced genes found in clusters. Yellow, induced orthologs; White, no induced orthologs; Gray (NA), no orthologs identified.

Supplementary Figure S6. OMCL clusters with genes involved in glycolysis and fermentation. OMCL clusters with change in at least one species were included. Up, Up-regulated DEGs; Down, Down-regulated DEGs. Bright yellow, Cluster has >1 induced DEG; Light yellow, Cluster has 1 induced DEG; Bright blue, Cluster includes >1 reduced DEG in; Light blue, Cluster includes 1 reduced DEG; White, No DEG ortholog identified; Gray, ortholog not identified.

Supplementary Table S1. Analysis of differentially regulated genes.

Supplementary Table S2. A list of Jatropha AP2/ERF transcription factors.

Supplementary Table S3. Real-time PCR primer information.

Supplementary Table S4. Analysis of waterlogging transcriptome from Jatropha, Arabidopsis, gray poplar, and rice.

Supplementary Table S5. Comparative analysis of waterlogging transcriptome.

References

Albrecht, G., and Mustroph, A. (2003). Localization of sucrose synthase in wheat roots: increased in situ activity of sucrose synthase correlates with cell wall thickening by cellulose deposition under hypoxia. Planta 217, 252–260. doi: 10.1007/s00425-003-0995-6

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Allegre, A., Silvestre, J., Morard, P., Kallerhoff, J., and Pinelli, E. (2004). Nitrate reductase regulation in tomato roots by exogenous nitrate: a possible role in tolerance to long-term root anoxia. J. Exp. Bot. 55, 2625–2634. doi: 10.1093/Jxb/Erh258

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Anders, S., Pyl, P. T., and Huber, W. (2014). HTSeq; A Python framework to work with high-throughput sequencing data. bioRxiv. doi: 10.1101/002824

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Bailey-Serres, J., Fukao, T., Gibbs, D. J., Holdsworth, M. J., Lee, S. C., Licausi, F., et al. (2012). Making sense of low oxygen sensing. Trends Plant Sci. 17, 129–138. doi: 10.1016/j.tplants.2011.12.004

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Bailey-Serres, J., and Voesenek, L. A. (2008). Flooding stress: acclimations and genetic diversity. Annu. Rev. Plant Biol. 59, 313–339. doi: 10.1146/annurev.arplant.59.032607.092752

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Banti, V., Loreti, E., Novi, G., Santaniello, A., Alpi, A., and Perata, P. (2008). Heat acclimation and cross-tolerance against anoxia in Arabidopsis. Plant Cell Environ. 31, 1029–1037. doi: 10.1111/j.1365-3040.2008.01816.x

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Basha, S. D., and Sujatha, M. (2009). Genetic analysis of Jatropha species and interspecific hybrids of Jatropha curcas using nuclear and organelle specific markers. Euphytica 168, 197–214. doi: 10.1007/s10681-009-9900-0

CrossRef Full Text | Google Scholar

Biemelt, S., Hajirezaei, M. R., Melzer, M., Albrecht, G., and Sonnewald, U. (1999). Sucrose synthase activity does not restrict glycolysis in roots of transgenic potato plants under hypoxic conditions. Planta 210, 41–49. doi: 10.1007/s004250050652

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Bieniawska, Z., Paul Barratt, D. H., Garlick, A. P., Thole, V., Kruger, N. J., Martin, C., et al. (2007). Analysis of the sucrose synthase gene family in Arabidopsis. Plant J. 49, 810–828. doi: 10.1111/j.1365-313X.2006.03011.x

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Bologa, K. L., Fernie, A. R., Leisse, A., Loureiro, M. E., and Geigenberger, P. (2003). A bypass of sucrose synthase leads to low internal oxygen and impaired metabolic performance in growing potato tubers. Plant Physiol. 132, 2058–2072. doi: 10.1104/pp.103.022236

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Branco-Price, C., Kaiser, K. A., Jang, C. J., Larive, C. K., and Bailey-Serres, J. (2008). Selective mRNA translation coordinates energetic and metabolic adjustments to cellular oxygen deprivation and reoxygenation in Arabidopsis thaliana. Plant J. 56, 743–755. doi: 10.1111/j.1365-313X.2008.03642.x

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Champion, A., Hebrard, E., Parra, B., Bournaud, C., Marmey, P., Tranchant, C., et al. (2009). Molecular diversity and gene expression of cotton ERF transcription factors reveal that group IXa members are responsive to jasmonate, ethylene and Xanthomonas. Mol. Plant Pathol. 10, 471–485. doi: 10.1111/j.1364-3703.2009.00549.x

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Christianson, J. A., Llewellyn, D. J., Dennis, E. S., and Wilson, I. W. (2010). Global gene expression responses to waterlogging in roots and leaves of cotton (Gossypium hirsutum L.). Plant Cell Physiol. 51, 21–37. doi: 10.1093/pcp/pcp163

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Dordas, C., Hasinoff, B. B., Igamberdiev, A. U., Manac'h, N., Rivoal, J., and Hill, R. D. (2003). Expression of a stress-induced hemoglobin affects NO levels produced by alfalfa root cultures under hypoxic stress. Plant J. 35, 763–770. doi: 10.1046/j.1365-313X.2003.01846.x

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Dordas, C., Hasinoff, B. B., Rivoal, J., and Hill, R. D. (2004). Class-1 hemoglobins, nitrate and NO levels in anoxic maize cell-suspension cultures. Planta 219, 66–72. doi: 10.1007/s00425-004-1212-y

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Fukao, T., Yeung, E., and Bailey-Serres, J. (2011). The submergence tolerance regulator SUB1A mediates crosstalk between submergence and drought tolerance in rice. Plant Cell 23, 412–427. doi: 10.1105/tpc.110.080325

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Geigenberger, P., Fernie, A. R., Gibon, Y., Christ, M., and Stitt, M. (2000). Metabolic activity decreases as an adaptive response to low internal oxygen in growing potato tubers. Biol. Chem. 381, 723–740. doi: 10.1515/BC.2000.093

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Gibbs, D. J., Lee, S. C., Isa, N. M., Gramuglia, S., Fukao, T., Bassel, G. W., et al. (2011). Homeostatic response to hypoxia is regulated by the N-end rule pathway in plants. Nature 479, 415–418. doi: 10.1038/nature10534

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Gibbs, D. J., Md Isa, N., Movahedi, M., Lozano-Juste, J., Mendiondo, G. M., Berckhan, S., et al. (2014). Nitric oxide sensing in plants is mediated by proteolytic control of group VII ERF transcription factors. Mol. Cell 53, 369–379. doi: 10.1016/j.molcel.2013.12.020

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Gimeno, V., Syvertsen, J. P., Simon, I., Nieves, M., Diaz-Lopez, L., Martinez, V., et al. (2012). Physiological and morphological responses to flooding with fresh or saline water in Jatropha curcas. Environ. Exp. Bot. 78, 47–55. doi: 10.1016/j.envexpbot.2011.12.014

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Giuntoli, B., Lee, S. C., Licausi, F., Kosmacz, M., Oosumi, T., Van Dongen, J. T., et al. (2014). A trihelix DNA binding protein counterbalances hypoxia-responsive transcriptional activation in Arabidopsis. PLoS Biol. 12:e1001950. doi: 10.1371/journal.pbio.1001950

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Guglielminetti, L., Perata, P., and Alpi, A. (1995). Effect of anoxia on carbohydrate metabolism in rice seedlings. Plant Physiol. 108, 735–741.

Pubmed Abstract | Pubmed Full Text | Google Scholar

Guglielminetti, L., Wu, Y., Boschi, E., Yamaguchi, J., Favati, A., Vergara, M., et al. (1997). Effects of anoxia on sucrose degrading enzymes in cereal seeds. J. Plant Physiol. 150, 251–258. doi: 10.1016/S0176-1617(97)80116-2

CrossRef Full Text | Google Scholar

Gupta, K. J., Shah, J. K., Brotman, Y., Jahnke, K., Willmitzer, L., Kaiser, W. M., et al. (2012). Inhibition of aconitase by nitric oxide leads to induction of the alternative oxidase and to a shift of metabolism towards biosynthesis of amino acids. J. Exp. Bot. 63, 1773–1784. doi: 10.1093/jxb/ers053

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Hattori, Y., Nagai, K., Furukawa, S., Song, X. J., Kawano, R., Sakakibara, H., et al. (2009). The ethylene response factors SNORKEL1 and SNORKEL2 allow rice to adapt to deep water. Nature 460, 1026–1030. doi: 10.1038/nature08258

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Hebelstrup, K. H., Hunt, P., Dennis, E., Jensen, S. B., and Jensen, E. O. (2006). Hemoglobin is essential for normal growth of Arabidopsis organs. Physiol. Plant. 127, 157–166. doi: 10.1111/j.1399-3054.2006.00653.x

CrossRef Full Text | Google Scholar

Heckmann, A. B., Hebelstrup, K. H., Larsen, K., Micaelo, N. M., and Jensen, E. O. (2006). A single hemoglobin gene in Myrica gale retains both symbiotic and non-symbiotic specificity. Plant Mol. Biol. 61, 769–779. doi: 10.1007/s11103-006-0048-1

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Hinz, M., Wilson, I. W., Yang, J., Buerstenbinder, K., Llewellyn, D., Dennis, E. S., et al. (2010). Arabidopsis RAP2.2: an ethylene response transcription factor that is important for hypoxia survival. Plant Physiol. 153, 757–772. doi: 10.1104/pp.110.155077

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Hirakawa, H., Tsuchimoto, S., Sakai, H., Nakayama, S., Fujishiro, T., Kishida, Y., et al. (2012). Upgraded genomic information of Jatropha curcas L. Plant Biotechnol. 29, 123–130. doi: 10.5511/plantbiotechnology.12.0515a

CrossRef Full Text | Google Scholar

Horan, K., Jang, C., Bailey-Serres, J., Mittler, R., Shelton, C., Harper, J. F., et al. (2008). Annotating genes of known and unknown function by large-scale coexpression analysis. Plant Physiol. 147, 41–57. doi: 10.1104/pp.108.117366

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Horchani, F., Aschi-Smiti, S., and Brouquisse, R. (2010). Involvement of nitrate reduction in the tolerance of tomato (Solanum lycopersicum L.) plants to prolonged root hypoxia. Acta Physiol. Plant. 32, 1113–1123. doi: 10.1007/s11738-010-0503-0

CrossRef Full Text | Google Scholar

Jaganath, B., Subramanyam, K., Mayavan, S., Karthik, S., Elayaraja, D., Udayakumar, R., et al. (2014). An efficient in planta transformation of Jatropha curcas (L.) and multiplication of transformed plants through in vivo grafting. Protoplasma 251, 591–601. doi: 10.1007/s00709-013-0558-z

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Juntawong, P., Girke, T., Bazin, J., and Bailey-Serres, J. (2014). Translational dynamics revealed by genome-wide profiling of ribosome footprints in Arabidopsis. Proc. Natl. Acad. Sci. U.S.A. 111, E203–E212. doi: 10.1073/pnas.1317811111

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Kreuzwieser, J., Hauberg, J., Howell, K. A., Carroll, A., Rennenberg, H., Millar, A. H., et al. (2009). Differential response of gray poplar leaves and roots underpins stress adaptation during hypoxia. Plant Physiol. 149, 461–473. doi: 10.1104/pp.108.125989

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Lee, S. C., Mustroph, A., Sasidharan, R., Vashisht, D., Pedersen, O., Oosumi, T., et al. (2011). Molecular characterization of the submergence response of the Arabidopsis thaliana ecotype Columbia. New. Phytol. 190, 457–471. doi: 10.1111/j.1469-8137.2010.03590.x

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Li, L., Stoeckert, C. J. Jr., and Roos, D. S. (2003). OrthoMCL: identification of ortholog groups for eukaryotic genomes. Genome Res. 13, 2178–2189. doi: 10.1101/gr.1224503

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Licausi, F., Kosmacz, M., Weits, D. A., Giuntoli, B., Giorgi, F. M., Voesenek, L., et al. (2011). Oxygen sensing in plants is mediated by an N-end rule pathway for protein destabilization. Nature 479, 419–422. doi: 10.1038/Nature10536

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Liu, P., Sun, F., Gao, R., and Dong, H. (2012). RAP2.6L overexpression delays waterlogging induced premature senescence by increasing stomatal closure more than antioxidant enzyme activity. Plant Mol. Biol. 79, 609–622. doi: 10.1007/s11103-012-9936-8

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Livak, K. J., and Schmittgen, T. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 25, 402–408. doi: 10.1006/meth.2001.1262

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Lohse, M., Nagel, A., Herter, T., May, P., Schroda, M., Zrenner, R., et al. (2013). Mercator: a fast and simple web server for genome scale functional annotation of plant sequence data. Plant Cell Environ. 37, 1250–1258. doi: 10.1111/pce.12231

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Loreti, E., Poggi, A., Novi, G., Alpi, A., and Perata, P. (2005). A genome-wide analysis of the effects of sucrose on gene expression in Arabidopsis seedlings under anoxia. Plant Physiol. 137, 1130–1138. doi: 10.1104/pp.104.057299

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

McCarthy, F. M., Wang, N., Magee, G. B., Nanduri, B., Lawrence, M. L., Camon, E. B., et al. (2006). AgBase: a functional genomics resource for agriculture. BMC Genomics 7:229. doi: 10.1186/1471-2164-7-229

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Mukherjee, P., Varshney, A., Johnson, T. S., and Jha, T. B. (2011). Jatropha curcas: a review on biotechnological status and challenges. Plant Biotechnol. Rep. 5, 197–215. doi: 10.1007/s11816-011-0175-2

CrossRef Full Text | Google Scholar

Mustroph, A., Lee, S. C., Oosumi, T., Zanetti, M. E., Yang, H., Ma, K., et al. (2010). Cross-kingdom comparison of transcriptomic adjustments to low oxygen stress highlights conserved and plant-specific responses. Plant Physiol. 152, 1484–1500. doi: 10.1104/pp.109.151845

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Mustroph, A., Zanetti, M. E., Jang, C. J., Holtan, H. E., Repetti, P. P., Galbraith, D. W., et al. (2009). Profiling translatomes of discrete cell populations resolves altered cellular priorities during hypoxia in Arabidopsis. Proc. Natl. Acad. Sci. U.S.A. 106, 18843–18848. doi: 10.1073/pnas.0906131106

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Narsai, R., Howell, K. A., Carroll, A., Ivanova, A., Millar, A. H., and Whelan, J. (2009). Defining core metabolic and transcriptomic responses to oxygen availability in rice embryos and young seedlings. Plant Physiol. 151, 306–322. doi: 10.1104/pp.109.142026

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Narsai, R., Rocha, M., Geigenberger, P., Whelan, J., and van Dongen, J. T. (2011). Comparative analysis between plant species of transcriptional and metabolic responses to hypoxia. New Phytol. 190, 472–487. doi: 10.1111/j.1469-8137.2010.03589.x

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Narsai, R., and Whelan, J. (2013). How unique is the low oxygen response? An analysis of the anaerobic response during germination and comparison with abiotic stress in rice and Arabidopsis. Front. Plant Sci. 4:349. doi: 10.3389/fpls.2013.00349

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Oliveira, H. C., Freschi, L., and Sodek, L. (2013). Nitrogen metabolism and translocation in soybean plants subjected to root oxygen deficiency. Plant Physiol. Biochem. 66, 141–149. doi: 10.1016/j.plaphy.2013.02.015

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Oliveira, H. C., and Sodek, L. (2013). Effect of oxygen deficiency on nitrogen assimilation and amino acid metabolism of soybean root segments. Amino Acids 44, 743–755. doi: 10.1007/s00726-012-1399-3

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Perazzolli, M., Dominici, P., Romero-Puertas, M. C., Zago, E., Zeier, A., Sonoda, M., et al. (2004). Arabidopsis nonsymbiotic hemoglobin AHb1 modulates nitric oxide bioactivity. Plant Cell 16, 2785–2794. doi: 10.1105/tpc.104.025379

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Popluechai, S., Breviario, D., Mulpuri, S., Makkar, H., Raorane, M., Reddy, A., et al. (2009). Narrow genetic and apparent phenetic diversity in Jatropha curcas: initial success with generating low phorbol ester interspecific hybrids. Nat. Preced.

Pucciariello, C., Parlanti, S., Banti, V., Novi, G., and Perata, P. (2012). Reactive oxygen species-driven transcription in Arabidopsis under oxygen deprivation. Plant Physiol. 159, 184–196. doi: 10.1104/pp.111.191122

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

R Development Core Team. (2005). R: A Language and Environment for Statistical Computing [Online]. R Foundation for Statistical Computing (Vienna). Available online at: http://www.r-project.org/

Robinson, M. D., Mccarthy, D. J., and Smyth, G. K. (2010). edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26, 139–140. doi: 10.1093/bioinformatics/btp616

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Rockel, P., Strube, F., Rockel, A., Wildt, J., and Kaiser, W. M. (2002). Regulation of nitric oxide (NO) production by plant nitrate reductase in vivo and in vitro. J. Exp. Bot. 53, 103–110. doi: 10.1093/jexbot/53.366.103

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Rzewuski, G., Cornell, K. A., Rooney, L., Burstenbinder, K., Wirtz, M., Hell, R., et al. (2007). OsMTN encodes a 5′5′-methylthioadenosine nucleosidase that is up-regulated during submergence-induced ethylene synthesis in rice (Oryza sativa L.). J. Exp. Bot. 58, 1505–1514. doi: 10.1093/jxb/erm014

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Sadasivam, S., and Manickam, A. (2007). Biochemical Methods. New Delhi: New Age International (P) Ltd.

Google Scholar

Sato, S., Hirakawa, H., Isobe, S., Fukai, E., Watanabe, A., Kato, M., et al. (2011). Sequence analysis of the genome of an oil-bearing tree, Jatropha curcas L. DNA Res. 18, 65–76. doi: 10.1093/dnares/dsq030

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Smyth, G. K. (2004). Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat. Appl. Genet. Mol. Biol. 3, Article3. doi: 10.2202/1544-6115.1027

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Stohr, C., Strube, F., Marx, G., Ullrich, W. R., and Rockel, P. (2001). A plasma membrane-bound enzyme of tobacco roots catalyses the formation of nitric oxide from nitrite. Planta 212, 835–841. doi: 10.1007/s004250000447

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Tamura, K., Peterson, D., Peterson, N., Stecher, G., Nei, M., and Kumar, S. (2011). MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol. Biol. Evol. 28, 2731–2739. doi: 10.1093/molbev/msr121

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Thimm, O., Blasing, O., Gibon, Y., Nagel, A., Meyer, S., Kruger, P., et al. (2004). MAPMAN: a user-driven tool to display genomics data sets onto diagrams of metabolic pathways and other biological processes. Plant J. 37, 914–939. doi: 10.1111/j.1365-313X.2004.02016.x

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Usadel, B., Nagel, A., Steinhauser, D., Gibon, Y., Blasing, O. E., Redestig, H., et al. (2006). PageMan: an interactive ontology tool to generate, display, and annotate overview graphs for profiling experiments. BMC Bioinformatics 7:535. doi: 10.1186/1471-2105-7-535

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Voesenek, L., Banga, M., Thier, R. H., Mudde, C. M., Harren, F., Barendse, G., et al. (1993). Submergence-induced ethylene synthesis, entrapment, and growth in two plant species with contrasting flooding resistances. Plant Physiol. 103, 783–791.

Pubmed Abstract | Pubmed Full Text | Google Scholar

Vriezen, W. H., Hulzink, R., Mariani, C., and Voesenek, L. A. (1999). 1-aminocyclopropane-1-carboxylate oxidase activity limits ethylene biosynthesis in Rumex palustris during submergence. Plant Physiol. 121, 189–196. doi: 10.1104/pp.121.1.189

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Wang, A., Tan, D., Takahashi, A., Li, T. Z., and Harada, T. (2007). MdERFs, two ethylene-response factors involved in apple fruit ripening. J. Exp. Bot. 58, 3743–3748. doi: 10.1093/jxb/erm224

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Xu, K., Xu, X., Fukao, T., Canlas, P., Maghirang-Rodriguez, R., Heuer, S., et al. (2006). Sub1A is an ethylene-response-factor-like gene that confers submergence tolerance to rice. Nature 442, 705–708. doi: 10.1038/nature04920

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Zeng, Y., Wu, Y., Avigne, W. T., and Koch, K. E. (1999). Rapid repression of maize invertases by low oxygen. Invertase/sucrose synthase balance, sugar signaling potential, and seedling survival. Plant Physiol. 121, 599–608. doi: 10.1104/pp.121.2.599

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text | Google Scholar

Keywords: Jatropha curcas, genomics, low oxygen, ERFs, transcriptional control, RNA-seq

Citation: Juntawong P, Sirikhachornkit A, Pimjan R, Sonthirod C, Sangsrakru D, Yoocha T, Tangphatsornruang S and Srinives P (2014) Elucidation of the molecular responses to waterlogging in Jatropha roots by transcriptome profiling. Front. Plant Sci. 5:658. doi: 10.3389/fpls.2014.00658

Received: 20 May 2014; Accepted: 04 November 2014;
Published online: 02 December 2014.

Edited by:

Franck Anicet Ditengou, University of Freiburg, Germany

Reviewed by:

Angelika Mustroph, University Bayreuth, Germany
Veronique Storme, VIB-University of Ghent, Belgium

Copyright © 2014 Juntawong, Sirikhachornkit, Pimjan, Sonthirod, Sangsrakru, Yoocha, Tangphatsornruang and Srinives. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Peerasak Srinives, Department of Agronomy, Faculty of Agriculture at Kamphaeng Saen, Kasetsart University, Kamphaeng Saen, Nakhon Pathom 73140, Thailand e-mail: agrpss@yahoo.com

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.