Figures
Abstract
Waterlogging is a common adverse environmental condition that limits plant growth. Sesame (Sesamum indicum) is considered a drought-tolerant oil crop but is typically susceptible to harmful effects from waterlogging. The present study used comparative analysis to explore the waterlogging stress response associated with two sesame genotypes. The RNA-seq dataset generated during a time course of 0, 3, 9 and 15 h of waterlogging as well as 20 h post-drainage indicated that stress gradually suppressed the expression of sesame genes, with 9 h as the critical time point for the response of sesame to waterlogging stress. Of the 19,316 genes expressed during waterlogging, 72.1% were affected significantly. Sesame of both tolerant and susceptible genotypes showed decreased numbers of upregulated differentially expressed genes (DEGs) but increased numbers of downregulated DEGs at the onset of waterlogging. However, the tolerant-genotype sesame exhibited 25.5% more upregulated DEGs and 29.7% fewer downregulated DEGs than those of the susceptible-genotype strain between 3 and 15 h. The results indicated that the tolerant sesame displayed a more positive gene response to waterlogging. A total of 1,379 genes were significantly induced and commonly expressed in sesame under waterlogging conditions from 3 to 15 h regardless of tolerance level; of these genes, 98 are known homologous stress responsive genes, while the remaining 1,281 are newly reported here. This gene set may represent the core genes that function in response to waterlogging, including those related mainly to energy metabolism and phenylpropanoid biosynthesis. Furthermore, a set of 3,016 genes functioning in energy supply and cell repair or formation was activated in sesame recovery from waterlogging stress. A comparative analysis between sesame of the tolerant and susceptible genotypes revealed 66 genes that may be candidates for improving sesame tolerance to waterlogging. This study provided a comprehensive picture of the sesame gene expression pattern in response to waterlogging stress. These results will help dissect the mechanism of the sesame response to waterlogging and identify candidate genes to improve its tolerance.
Citation: Wang L, Li D, Zhang Y, Gao Y, Yu J, Wei X, et al. (2016) Tolerant and Susceptible Sesame Genotypes Reveal Waterlogging Stress Response Patterns. PLoS ONE 11(3): e0149912. https://doi.org/10.1371/journal.pone.0149912
Editor: Xianlong Zhang, National Key Laboratory of Crop Genetic Improvement, CHINA
Received: November 14, 2015; Accepted: February 6, 2016; Published: March 2, 2016
Copyright: © 2016 Wang et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: The RNA-seq data sets supporting the results of this article were submitted to GeneBank with SRA accession numbers SRR2886791 (WS) and SRR2886790 (WT) under experiment SRX1406650.
Funding: This work was supported by the National Natural Science Foundation of China (no. 31101182), China Agriculture Research System (no. CARS−15), Agricultural Science and Technology Innovation Project of the Chinese Academy of Agricultural Sciences (CAAS-ASTIP-2013-OCRI), and the National Crop Germplasm Resources Infrastructure (NB2012-2130135).
Competing interests: The authors have declared that no competing interests exist.
Abbreviations: DEG, differentially expressed gene; FDR, false discovery rate; GO, gene ontology; IGPD, imidazoleglycerol-phosphate dehydratase; KEGG, kyoto encyclopedia of genes and genomes; RPKM, reads per kilo base of exon model per million mapped reads; WT, waterlogging-tolerant sesame; WS, waterlogging-susceptible sesame
Introduction
Crops have developed a variety of adaptations to natural and agronomic habitats that have contributed to their survival and regeneration over the course of evolution. However, abnormal or altered climates tend to overwhelm the endurance capacity of a crop, imposing severe negative influences on the productivity of arable farmland as the vast majority of crops are ill-suited to grow under stress conditions. Waterlogging is a common adverse environmental condition that limits plant growth. Waterlogging reduces gas exchange between plant tissues and the atmosphere, resulting in an imbalance between slow diffusion and rapid consumption of oxygen in the rhizosphere that drastically reduces the oxygen supply and induces anoxia in plants [1].
Sesame, a crop with high oil content, has the potential capacity to combat nutritional deficiencies in developing regions and countries. Most current cultivars contain 50–60% oil and 18–24% protein in their seeds [2, 3] (http://ndb.nal.usda.gov/). In particular, greater than 80% of its oil is in the form of unsaturated fatty acids, which are more beneficial for human health than are saturated fatty acids. In addition, the antioxidant properties of sesame lignans, primarily sesamin and sesamolin, are used for therapeutic and cosmetic applications [4, 5]. Sesame is typically considered drought-tolerant but susceptible to waterlogging, a property that can be ascribed to its suspected origin in Africa or India and its subsequent dispersal to tropical or semitropical regions [3, 6]. According to the Food and Agricultural Organization, the average sesame yield was alarmingly low at only 617 kg/ha worldwide in 2011 and ranked second to last among 22 oil crops between 2007 and 2011 (http://faostat.fao.org). This low yield may be attributed to several reasons, but waterlogging is a primary factor that has a severe effect in countries such as China and Korea due to changing climate.
To understand the effects of abiotic stress in an effort to maintain a stable food supply, a number of studies have investigated the responses of model plants and crops to stresses [7]. These studies have revealed that plant responses to different stresses are coordinated by complex and often interconnected signaling pathways that regulate numerous metabolic networks [8, 9]. At the protein level, low oxygen selectively induces the synthesis of anaerobic proteins, especially enzymes involved in sugar metabolism, glycolysis and fermentation [10, 11]. The vast majority of these proteins have been investigated in waterlogging-susceptible or -tolerant strains of Arabidopsis or rice [8, 12, 13]. Despite an increased understanding of the adaptive mechanisms and molecular regulation at play in other crops, the mechanisms underlying the sesame response to waterlogging need further elucidation [14, 15].
With the recent completion of the de novo assembly of the sesame genome [16, 17], it is now possible to profile gene expression in response to waterlogging by sequencing [15, 18]. Thus, the present study explored the molecular mechanisms of the sesame response to waterlogging stress by performing a comparative RNA-seq-based analysis between waterlogging-tolerant and -susceptible genotypes.
Materials and Methods
Materials
The waterlogging-tolerant cultivar Zhongzhi No. 13 (WT) and the waterlogging-susceptible strain ZZM0563 (WS) were used in this study. The two plants were potted under the same growth and experimental conditions and treated with waterlogging simultaneously for 15 h during the flowering stage. For the sesame of both genotypes, a total of 10 samples were collected at equivalent stages: 0, 3, 9, and 15 h under waterlogging stress as well as 20 h post-drainage (S1 Table). For each sample, the root tips of five plants were obtained and mixed. All samples were subjected to RNA-seq analysis.
RNA extraction and library preparation
RNA extraction and sequencing were performed according to the procedure described by Wei et al. [19]. Total RNA was extracted from 10 samples using TRIzol reagent (Invitrogen), and its concentration and quality were tested using an ultraviolet spectrophotometer and denaturing agarose gels (1.0%). Next, the RNA was treated with DNase I and magnetic oligo (dT) beads. Divalent cations and heat were used to fragment the purified poly-(A) mRNA before library construction. The RNA fragments were transcribed into double-stranded cDNA, which were subjected to end repair and index adapter ligation using T4 DNA polymerase. The adapter-ligated cDNA fragments of the desired size range (200 ± 25 bp) were excised from an agarose gel. The cDNA fragments were then enriched selectively via polymerase chain reaction (PCR).
Data generation and quality assessment
The 10 cDNA libraries generated from RNA samples isolated from the roots of WT and WS at 0, 3, 9 and 15 h during stress and 20 h post stress (P20h) were paired-end sequenced using the Illumina platform. We used FastQC (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/) to determine the base quality of the RNA-seq reads and removed the paired-end reads containing more than 5% ambiguous residues (Ns) and those containing > 10% bases with a Phred quality score < 20. The remaining reads were referred to as “clean reads”[20]. After cleaning and quality assessment, approximately 25.6–26.8 million high-quality reads of 90-bp lengths remained in each sample (S1 Fig and S2 Table). The high-quality reads were mapped to a reference genome using SOAPaligner/SOAP2 [16, 21, 22], allowing no more than one mismatch in the alignment. Approximately 65.1–81.4% of the clean reads were uniquely mapped to the reference genome [16], with 50.3–64.2% of them uniquely mapped to genic regions. Most gene models were covered by reads over 80% (S2 Fig), and the reads were evenly distributed on these genes (S3 Fig).
Statistical analysis of gene expression
Gene expression levels were calculated based on the number of unique matched reads to the sesame genome [16] and were normalized to reads per kilobase of transcript per million mapped reads (RPKM) using Cufflinks 2.0 software [23]. Differentially expressed genes (DEGs) were identified as described by Chen et al. [24] and Wang et al. [15] based on a Poisson distribution [25]. In addition to P-values, false discovery rates (FDR) were used to determine the threshold P-value in multiple tests. FDR ≤ 0.001 and the absolute value of log2Ratio ≥ 1 were used as thresholds to determine the significance of the DEGs [26].
Gene annotation and enrichment analysis
Target genes were annotated by gene ontology (GO; http://www.geneontology.org) and metabolic pathway. GO is an international standardized gene function classification system that offers a dynamically updated, controlled vocabulary and a strictly defined conceptual framework for the comprehensive description of gene properties and products in any organism. GO term and metabolic pathway enrichment analyses were performed according to Wang et al. [15]. GO terms were summarized using the online tool GOSlimAuto [27].
Homologous stress-responsive genes and hormone-related genes
The 3,150 biotic and abiotic stress-responsive genes of Arabidopsis thaliana were downloaded from the Stress Responsive Transcription Factor Database (STIFDB V2.0) [28]. The 1,318 Arabidopsis hormone-related genes were downloaded from Arabidopsis Hormone Database 2.0 (AHD2.0) [29]. Then, their homologous genes in sesame were predicted using the Reciprocal Best Blast Hit (RBH) method [30, 31].
Real-time quantitative PCR (qRT-PCR)
RT-PCR analyses of target gene expression were performed according to Wang et al. [15] using the iQ™5 Real-Time PCR Detection System (Bio-Rad, San Diego, CA, USA). Samples were run in triplicate on the same plate with a negative control that lacked cDNA. The positive control for each sample was the sesame gene ubiquitin-conjugating enzyme 9 (UBC9). The PCR efficiency was determined by a series of 10-fold dilutions of cDNA. The calculated efficiency of all primers was 0.9 to 1.0. The relative expression levels of the genes were calculated using the 2-ΔΔCT method [32].
Results and Discussion
WT and WS showed similar expressed gene numbers under waterlogging stress
After mapping the reads to the reference sequence [16], gene expression levels were calculated as RPKMs. A range of 15,499 to 17,875 genes was expressed at each time point, with RPKMs ranging from 1 to 23,367 and averaging 48.6 (Fig 1a and 1b). Comparative numbers of genes were expressed in the sesame of both genotypes (18,714 in WT and 19,008 in WS)) (Fig 1c and 1d), resulting in a total of 19,316 genes expressed. Of these genes, 14,528 and 14,463 were expressed commonly at each time point in WT and WS plants, respectively, with 13,959 shared by the two genotypes (Figs 1b, 1d and S4), accounting for approximately 51.4% of sesame genes. The activation of these genes before and during waterlogging and post-drainage suggests their importance in survival of sesame. GO enrichment of the common gene set indicated that most genes were involved in metabolic processes and were strongly associated with initiating functions, such as protein binding, small molecule binding, nucleotide binding and nucleic acid binding (S5 Fig). In addition to the genes commonly expressed between WT and WS, there were 569 genes uniquely expressed in WT and 504 in WS (RPKM ≥ 1) (S4 Fig). The 569 genes in WT were enriched in the molecular function term imidazoleglycerol-phosphate dehydratase (IGPD) activity, but no GO terms were enriched among the 504 genes in WS. IGPD catalyzes the sixth step in the histidine biosynthesis pathway, an essential process in plants and microorganisms [33]. We speculated that IGPD genes might help WT endure waterlogging stress.
(a) and (b) indicate the gene expression levels (RPKM ≥ 1, RPKM: reads per kilobase of transcript per million mapped reads) in WT and WS plants at each time point, respectively. (A), (B), (C), (D), and (E) correspond to the samples at 0, 3, 9, and 15 h during waterlogging and 20 h post drainage, respectively. (F) represents the gene density (mRNA, 500-kb window) in the sesame linkage groups, and (G) represents the 16 assembled pseudomolecules of sesame. (c) and (d) indicate the shared and uniquely expressed gene numbers during a time-point assay of waterlogging. WT: waterlogging-tolerant sesame genotype, WS: waterlogging-susceptible sesame genotype.
Gradual suppression of expressed genes by waterlogging
The statistical graph showed that the global transcriptional responses in root tissue were suppressed gradually by waterlogging, with the rate of expressed genes decreasing from 62.8 to 58.5% in WT, and from 65.0 to 57% in WS in the time-course assay under stress, from 3 to 9 h of waterlogging. When the water was removed (P20), many genes in the sesame of both genotypes resumed expression, but the numbers of active genes were lower in the tolerant plants compared with the susceptible plants (Fig 2a and S6 Fig).
(a) Percentages of the expressed genes during waterlogging from 0 to 15 h as well as 20 h post drainage. (b) The number of genes significantly affected by the stress. (c) Upregulated genes using 0 h as a control. (d) Downregulated genes using 0 h as a control. (e) Upregulated genes using the former sampling point as a control, illustrating the dynamic changes in differentially expressed genes (DEGs) in sesame upon waterlogging exposure. (f) Downregulated genes using the former sampling point as a control.
Next, we determined how many of these genes were affected significantly by stress. Using the 0 h time point as a control, the DEGs at each time point were identified for the two genotypes. Of the 19,316 expressed genes in sesame under waterlogging stress, 72.1% were differentially expressed in the samples during waterlogging stress or post-drainage in both WT and WS, compared with 4–15% of the assayed genes in Arabidopsis [34], rape [18], cotton [35] and poplar [36]. This result suggested that most of these genes are affected by low rhizosphere oxygen or are involved in recovery. The greater preponderance of affected genes in sesame underscores its sensitivity to waterlogging. However, the tolerant strain exhibited a smoother trend in fluctuation than did the susceptible plants, suggesting its enhanced ability to stabilize metabolism to combat abiotic stress (Fig 2a and 2b).
For the remaining 5,397 genes that were not significantly affected by waterlogging stress, enrichment analysis revealed that most were associated with a few GO terms (P < 0.01) (S3 Table). Using GO-SLIM, these enriched insensitive genes were further categorized into fundamental biological processes, including cellular metabolic, primary metabolic, and biosynthetic and transport processes, together with binding functions, catalytic activity, oxidoreductase activity, kinase activity and transferase activity (S7 Fig). Thus, these insensitive genes may be primarily housekeeping genes.
Sesame endurance to waterlogging showed a critical time point of 9 h
Over the time course of waterlogging stress, both WT and WS showed gradually increasing DEG numbers from 3 to 15 h (using the 0 h samples as controls), which decreased after the stress was removed. In comparison with the control (0 h), the number of upregulated DEGs decreased from 0 to 9 h in both genotypes and then increased, with 25.5 to 29.7% more upregulated DEGs in WT than WS over the treatment duration. However, the downregulated DEGs exhibited a different pattern, in that WT expressed 36.0% fewer downregulated DEGs than those of WS at 9 h and 15 h (Fig 2c and 2d), with no significant differences noted at the onset of waterlogging or post-drainage. By comparing each sample (except for 0 h) against those collected at the previous time-points, the immediate changes in DEG numbers in both WT and WS were examined. The two genotypes had similar numbers of upregulated genes at each time point but exhibited differences in downregulated gene numbers, in that more genes were downregulated in WS than WT at 9 h. These results demonstrated that 9 h was a critical time point for endurance of waterlogging stress in sesame (Fig 2e and 2f).
The core gene set in sesame responsible for the response to waterlogging stress
WT and WS showed different DEG numbers over the course of waterlogging, using 0 h as control; 2,401 genes were commonly and differentially expressed in the tolerant strain from 3 to 15 h and 2,593 in the susceptible plant. Of these DEGs, 1,379 were shared by the two genotypes (Fig 3a and S4 Table). The shared genes may represent a core gene set in sesame that function in response to waterlogging stress, despite tolerance level. These 1,379 genes were further clustered into 18 groups, and various expression patterns of the genes in each group were observed (S8 Fig). Of these genes, 326 were upregulated, representing positive responses to waterlogging stress, whereas 839 were downregulated in both WT and WS strains.
(a) Venn diagrams of the special and unique DEGs at time points between 3 and 15 h under waterlogging conditions. (b) Enriched gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway categories among the upregulated (up channel) and downregulated (down channel) genes in the core set.
Anoxic or hypoxic conditions of the rhizosphere affect the plant’s maintenance of numerous pathways [37]. To generate energy under conditions of low O2 availability, plants must regulate their molecular pathways to survive stress. Functional annotation of the 326 upregulated genes in response to stress revealed that sesame enlists the molecular functions of ribulose-bisphosphate carboxylase activity, carbon-carbon lyase activity, acyl-[acyl-carrier-protein] desaturase activity and lyase activity by activating various biological processes, such as those involved in the response to biotic stimuli (GO:0009607), response to stress (GO:0006950), defense response (GO:0006952), carbon fixation (GO:0015977), glycolysis (GO:0006096), glucose metabolism (GO:0006006) and generation of precursor metabolites and energy (GO:0006091).
Energy metabolism plays crucial roles in plant physiological adaptation to hypoxia [10]. Many studies have demonstrated that plants switch from respiration to fermentation as a rescue mechanism that regenerates NAD+ to maintain glycolysis and produce ATP under anaerobic conditions [11, 35, 36]. A high rate of fermentation increases the demand for carbohydrates. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that carbon fixation in photosynthetic organisms, as well as glycolysis and gluconeogenesis, were enriched pathways among the upregulated genes (Figs 3b, S9 and S10). Carbon fixation is a process used by plants to generate energy under abiotic stress [38, 39] and occurs during the second half of the photosynthetic process of converting light to energy. The upregulation of genes related to carbon fixation and glycolysis suggests that sesame will invoke or enhance energy generation to maintain growth during long periods of root hypoxia [36].
In addition to upregulating genes involved in carbon fixation and fermentative metabolism [34], an alternative abiotic stress survival mechanism relies on decreased energy consumption. Protein degradation plays a negative regulatory role in response to hypoxia [40], resulting in large-scale energy savings under hypoxic conditions [41]. We found that waterlogging stress in sesame resulted not only in the downregulation of genes enriched in nitrogen metabolism and cyanoamino acid metabolism, but also in secondary metabolite and phenylpropanoid biosynthesis (Fig 3c). Coupled with these processes, oxidation-reduction processes usually slow down to reduce the release of biological energy [42, 43]. Phenylpropanoid biosynthesis is regulated by biotic and abiotic stimuli [44–46], and the phenylpropanoid-based polymers such as lignin, suberin and tannin contribute substantially to the stability and robustness of plants in the face of mechanical or environmental damage. Some of the homologous genes (SIN_1026962, SIN_1023248, SIN_1022230, SIN_1013457, SIN_1012086, SIN_1009800, SIN_1007000) involved in phenylpropanoid biosynthesis had been reported function in response to biotic or abiotic stress in Arabidopsis [47, 48]. Their downregulation in expression suggested the harmful effect of waterlogging. In plants, biotic and abiotic stressors tend to intersect. The presence of abiotic stress can result in the reduced enhancement of its susceptibility to a biotic pest or pathogen, and vice versa [49]. For example, a change in temperature or humidity may induce the expression of defense genes in response to certain pathogens. These results indicate that energy metabolism and phenylpropanoid biosynthesis play crucial roles in sesame endurance against waterlogging stress.
The genes involved in recovery from waterlogging
After drainage, some sesame genes may be selectively employed in the recovery and reconstruction of metabolic processes. We observed that the number of DEGs differed between the two genotypes, as more genes were expressed in the tolerant strain (5,722) than the susceptible strain (4,585); 3,016 shared genes were expressed in both strains (S11 Fig). Further analysis subdivided the shared genes into 22 expression patterns in the WT and WS plants (S12 Fig). The enrichment of GO terms suggested that these genes significantly enhance biological processes related to energy metabolism, such as carbohydrate metabolism, cellular glucan metabolism and oxidation-reduction processes. In addition, such enhanced genes could mediate processes involved in molecular complex generation, including protein polymerization, cellular macromolecular complex assembly, cellulose biosynthesis, and polysaccharide metabolism. These results indicated that energy supply and cell reparation were crucial for sesame recovery following waterlogging stress (S13 Fig). Furthermore, we found enrichment of genes involved in microtubule-based processes and movement, DNA replication and the response to biotic stimuli. Microtubules are critical for a number of cellular processes, as they participate in fundamental aspects of cell division, growth, and differentiation and maintain the structure of the cell by forming the cytoskeleton in concert with microtubule-associated proteins [50]. Microtubules are essential for guard-cell function in Vicia and Arabidopsis [51]. Guard cells can regulate leaf gas exchange by altering their shape, thereby regulating the stomatal aperture. Thus, microtubules may affect photosynthesis and transpiration in sesame after drainage. The 3,016 shared genes are thought to participate in sesame recovery from waterlogging stress, functioning mainly in energy supply and cell repair.
Waterlogging response of the homologous stress-responsive and hormone-related genes in sesame
Based on STIFDB V2.0 (Stress Responsive Transcription Factor Database) [28], 1,740 homologous biotic and abiotic stress-responsive genes in sesame compared with Arabidopsis thaliana were identified using the RBH method [30, 31]. In total, 968 of these genes were expressed differentially during at least one waterlogging time point, and 1000 DEGs were observed post drainage between WT and WS plants (S14 Fig). Thus, these genes may also function in sesame in response to waterlogging stress. Among these genes, 217 and 187 were detected in the roots of WT and WS, respectively, during waterlogging between 3 and 15 h, with 98 shared genes identified (S5–S7 Tables). Clustering analysis of the 98 shared genes just distinguished the time points under waterlogging (3, 9, and 15 h) from others (0, P20 h) (S15 Fig). The 98 genes were also expected to be included in the 1,379 core waterlogging response genes. This finding suggested the presence of 1,281 newly discovered core responsive genes (S4 Table). Similarly, 276 homologous stress responsive genes were included in the 3,016 genes that expected to participate in recovery from waterlogging stress (S8 Table).
Hormones, such as abscisic acid, auxin, ethylene, gibberellin, salicylic acid, brassinosteroid, cytokinin and jasmonic acid, play essential roles throughout the lifespan of plants and can mediate many pathways required for both biotic resistance and abiotic tolerance [52]. Many genes involved in hormone production and signaling are affected in both the roots and shoots of waterlogged plants. For example, ethylene is involved in the responses to hypoxia and is thought to contribute to adventitious root production and aerenchyma formation [35, 53]. Gibberellin (GA) helps protect Oryza sativa L from drowning in deep water by increasing the level of replication protein A1 (RPA1), a heterotrimeric protein involved in the growth of the intercalary meristem of the internode [54]. According to the Arabidopsis Hormone Database 2.0 [29], 601 hormone homologous genes were identified in sesame, and 39 were also members of the aforementioned 1,379 shared DEGs, including abscisic acid, auxin, ethylene, GA, salicylic acid, brassinosteroid, cytokinin and jasmonic acid. Among them, abscisic acid was predominant under waterlogging stress, as 10 of its pathway genes were included (S9 Table). Among the 3,016 genes functioning in recovery, 94 homologous genes representing all eight hormones were detected, and the top three classes included abscisic acid (22), salicylic acid (18), and auxin (16) (S10 Table). These homologous stress-responsive hormone genes may be especially helpful for increasing sesame endurance to waterlogging stress.
Candidate waterlogging-tolerant genes in sesame
The susceptibility of sesame to waterlogging stress is generally attributed to its genotype. At the onset of stress, some sesame genes may begin, increase, or terminate expression to adapt to abnormal nutrient contents in the rhizosphere. In contrast to susceptible sesame, the tolerant genotype may contain genes that are expressed or function differently to improve its survivability in an abiotic environment. In the present study, the number of DEGs between WT and WS plants varied along the time course before and after initiation of waterlogging. Before stress, WT differentially expressed 2,528 genes compared with WS, and 73.3% of the genes were downregulated. After the onset of stress, the number of DEGs decreased until 9 h but consisted mainly of downregulated genes (Fig 4a). Approximately 9 h of waterlogging represented a critical point for sesame, as not only the DEG number between WT and WS decreased to a minimum, but also the number of upregulated genes was more than two-fold the number of downregulated genes. After 9 h of waterlogging, the DEGs between WT and WS increased and tended to equilibrate to a normal status until P20, with the numbers and rates of up- and downregulated genes comparable to those prior to stress. These results reflected that WT responded were more positively than did WS to waterlogging stress.
(a) DEGs in WT versus WS plants at various waterlogging time points revealed that waterlogging has a severe effect on gene expression at 9 h. (b) Expression patterns of the 66 common DEGs between WT and WS under waterlogging conditions. The RPKM values were log2 transformed. (c) A Venn diagram depicting the common and unique DEGs between 3 and 15 h of waterlogging in WT and WS plants.
Annotation for the DEGs between WT and WS at each time point revealed these genes were enriched mainly in pathways involving circadian rhythm, ascorbate and aldarate metabolism, indole alkaloid biosynthesis, flavonoid biosynthesis, phenylpropanoid biosynthesis, plant hormone signal transduction and plant-pathogen interaction. Notably, the enriched categories of phenylpropanoid biosynthesis and plant hormone signal transduction predominated at more than one time point, indicating their particular importance in the WT genotype. The expression profiles of genes in the two pathways indicated that while most genes were expressed at lower levels in WT plants at the beginning of waterlogging, some became more active than those in WS plants with ongoing stress (Figs 5 and S16).
Red indicates the genes upregulated in WT compared with WS plants, and green indicates the downregulated genes.
Between 3 and 15 h of waterlogging, 90 DEGs common to WT and WS were detected at all three time points (Fig 4b and 4c). Of these genes, 66 were expressed (RPKM > 10) during at least one time point and may function to improve the endurance of sesame in the face of waterlogging stress. In addition, seven of these genes were differentially expressed prior to waterlogging and 18 expressed post-drainage, suggesting they are long-acting genes involved in guarding against and recovery from waterlogging stress. A hierarchical cluster analysis based on Pearson correlation was used to divide the 66 genes into five groups (Fig 4b), suggesting these genes function in various patterns to enhance the tolerance of the WT genotype. We found five homologous stress-responsive genes (SIN_1024017, SIN_1022545, SIN_1021853, SIN_1021706, and SIN_1012279) among the 66 genes, and their expression patterns were validated by RT-PCR (S17 Fig). Notably, SIN_1021706 encodes a gene homologous to a hormone with abscisic acid 8'-hydroxylase activity involved in abscisic acid catabolism. Mutant analyses demonstrated that disruption of this gene results in increased drought tolerance [55, 56]. The expression of SIN_1021706 was increased in WT plants but downregulated in WS plants during waterlogging. In addition, the 66 active genes were highlighted in the functional categories of chitin catabolic processing, cell wall macromolecule catabolic processing with chitin binding, chitinase activity, and carbohydrate binding (S11 Table). Chitin is a characteristic component that contributes significantly to the mechanical strength of the cell wall, and chitin synthesis is essential for the plant’s response to environmental stress [57–59]. Therefore, these genes are associated with strengthening the endurance of sesame to waterlogging stress. Therefore, these genes are candidate waterlogging-tolerance genes that should be studied further.
Conclusions
Sesame is affected easily and negatively by waterlogging. The present study provides comprehensive insight into the genetic response of sesame to waterlogging stress over time. Waterlogging stress gradually suppresses the expression of sesame genes, with approximately half being significantly induced or suppressed. We found that ~9 h after the onset of waterlogging stress was a critical time point for sesame, providing useful information for farmers in planning how to ameliorate or prevent sesame water damage. The fact that more genes in sesame than in other species were affected by waterlogging demonstrated the high susceptibility of sesame. However, sesame of different genotypes exhibited differential response patterns under waterlogging stress; the tolerant species remained more stable throughout the waterlogging. The core gene set involved in the sesame response to waterlogging stress—and in particular, the newly discovered responsive genes—will be useful in expanding our knowledge of the mechanism underlying waterlogging stress in sesame. These data also revealed that energy metabolism and phenylpropanoid biosynthesis play crucial roles in regulating sesame endurance against waterlogging stress. Based on the comparison between tolerant and susceptible genotypes, a list of candidate genes for improving sesame tolerance to waterlogging stress was generated, which will serve as a guideline for future research. However, these genes require further verification in additional experiments. Future studies utilizing independent approaches—such as transcriptomics, proteomics, physiology, biochemistry, metabolite profiling, and genetic analysis—will lead to a better understanding of the response of sesame to waterlogging stress and will aid the effort to develop waterlogging-tolerant sesame.
Supporting Information
S1 Fig. Clean data for 10 sesame transcriptomes acquired by next generation sequencing technology.
https://doi.org/10.1371/journal.pone.0149912.s001
(TIF)
S3 Fig. Distributions of reads among gene regions.
https://doi.org/10.1371/journal.pone.0149912.s003
(TIF)
S4 Fig. Venn diagram depicting the genes commonly expressed in WT and WS plants at each time point.
https://doi.org/10.1371/journal.pone.0149912.s004
(TIF)
S5 Fig. GO-SLIM categories for the genes commonly expressed in both WT and WS plants.
BP: biological processes; MF: molecular function; CC: cellular component.
https://doi.org/10.1371/journal.pone.0149912.s005
(TIF)
S6 Fig. Statistics of the expressed genes in each sample during waterlogging and post-drainage.
https://doi.org/10.1371/journal.pone.0149912.s006
(TIF)
S7 Fig. GO-SLIM categories of the insensitive genes commonly expressed in both WT and WS plants.
BP: biological processes; MF: molecular function; CC: cellular component.
https://doi.org/10.1371/journal.pone.0149912.s007
(TIF)
S8 Fig. Expression patterns of the 18 clusters of the 1,379 core genes responsive to waterlogging.
https://doi.org/10.1371/journal.pone.0149912.s008
(TIF)
S9 Fig. The positions of the core responsive genes in the carbon fixation pathway in photosynthetic organisms.
https://doi.org/10.1371/journal.pone.0149912.s009
(TIF)
S10 Fig. The positions of the core responsive genes in the glycolysis and gluconeogenesis pathway.
https://doi.org/10.1371/journal.pone.0149912.s010
(TIF)
S11 Fig. Venn diagrams of the unique and common DEGs in the WT and WS genotypes at 20 h post-waterlogging.
https://doi.org/10.1371/journal.pone.0149912.s011
(TIF)
S12 Fig. Expression patterns of the 3,016 core genes responsible for recovery from waterlogging in sesame.
https://doi.org/10.1371/journal.pone.0149912.s012
(TIF)
S13 Fig. GO categories of the enriched core gene set involved in recovery after drainage.
BP: biological processes; MF: molecular function; CC: cellular component.
https://doi.org/10.1371/journal.pone.0149912.s013
(TIF)
S14 Fig. Expression profiles of sesame homologous stress-responsive genes during and post-waterlogging.
https://doi.org/10.1371/journal.pone.0149912.s014
(TIF)
S15 Fig. Expression patterns of the 98 homologous stress-responsive genes.
https://doi.org/10.1371/journal.pone.0149912.s015
(TIF)
S16 Fig. DEGs between WT and WS plants in the phenylpropanoid biosynthesis pathway.
https://doi.org/10.1371/journal.pone.0149912.s016
(TIF)
S17 Fig. RT-PCR validation of the genes in WT and WS.
https://doi.org/10.1371/journal.pone.0149912.s017
(TIF)
S1 Table. The samples and their appointed names used in the present study.
https://doi.org/10.1371/journal.pone.0149912.s018
(XLSX)
S3 Table. GO enrichment of the 5,397 genes insensitive to waterlogging.
https://doi.org/10.1371/journal.pone.0149912.s020
(XLSX)
S4 Table. The 1,379 core genes involved in sesame response to waterlogging.
https://doi.org/10.1371/journal.pone.0149912.s021
(XLSX)
S5 Table. 217 homologous stress responsive genes differently expressed in WT plants during waterlogging from 3 to 15 h.
https://doi.org/10.1371/journal.pone.0149912.s022
(XLSX)
S6 Table. 187 homologous stress-responsive genes differentially expressed in WS during waterlogging between 3 and 15 h.
https://doi.org/10.1371/journal.pone.0149912.s023
(XLSX)
S7 Table. 98 shared homologous stress-responsive DEGs in both WT and WS plants during waterlogging between 3 and 15 h.
https://doi.org/10.1371/journal.pone.0149912.s024
(XLSX)
S8 Table. 276 homologous stress-responsive genes differentially expressed in both WT and WS plants post-drainage.
https://doi.org/10.1371/journal.pone.0149912.s025
(XLSX)
S9 Table. List of the 39 core hormone homologous genes responsive to waterlogging in sesame.
https://doi.org/10.1371/journal.pone.0149912.s026
(XLSX)
S10 Table. List of the 94 hormone homologous genes functioning post-drainage in sesame.
https://doi.org/10.1371/journal.pone.0149912.s027
(XLSX)
S11 Table. GO enrichment of the 66 common DEGs between WT and WS under waterlogging from 3 to 15h.
https://doi.org/10.1371/journal.pone.0149912.s028
(XLSX)
Acknowledgments
The authors wish to acknowledge the contributions by the staff from BGI-Shenzhen, who assisted with the data generation and analysis.
Author Contributions
Conceived and designed the experiments: XZ LW. Performed the experiments: YZ DL YG. Analyzed the data: LW JY XW. Wrote the paper: LW XZ.
References
- 1. Sachs MM, Freeling M, Okimoto R. The anaerobic proteins of maize. Cell. 1980;20(3):761–7. pmid:7418006.
- 2. Mondal N, Bhat KV, Srivastava PS. Variation in Fatty Acid Composition in Indian Germplasm of Sesame. J Am Oil Chem Soc. 2010;87(11):1263–9. ISI:000283120500004.
- 3.
Ram R, Catlin D, Romero J, Cowley C. Sesame: New approaches for crop improvement. In: Simon JJaJE, editor. Advances in new crops. Portland, OR.: Timber Press; 1990. p. 225–8.
- 4. Nakano D, Kwak CJ, Fujii K, Ikemura K, Satake A, Ohkita M, et al. Sesamin metabolites induce an endothelial nitric oxide-dependent vasorelaxation through their antioxidative property-independent mechanisms: possible involvement of the metabolites in the antihypertensive effect of sesamin. J Pharmacol Exp Ther. 2006;318(1):328–35. pmid:16597711.
- 5. Miyahara Y, Komiya T, Katsuzaki H, Imai K, Nakagawa M, Ishi Y, et al. Sesamin and episesamin induce apoptosis in human lymphoid leukemia Molt 4B cells. Int J Mol Med. 2000;6(1):43–6. pmid:10851264.
- 6. Bedigian D. History and lore of sesame in Southwest Asia. Econ Bot. 2004;58(3):329–53.
- 7. Rasmussen S, Barah P, Suarez-Rodriguez MC, Bressendorff S, Friis P, Costantino P, et al. Transcriptome responses to combinations of stresses in Arabidopsis. Plant Physiol. 2013;161(4):1783–94. pmid:23447525.
- 8. Nakashima K, Ito Y, Yamaguchi-Shinozaki K. Transcriptional regulatory networks in response to abiotic stresses in Arabidopsis and grasses. Plant Physiol. 2009;149(1):88–95. pmid:19126699.
- 9. Miro B, Ismail AM. Tolerance of anaerobic conditions caused by flooding during germination and early growth in rice (Oryza sativa L.). Frontiers in plant science. 2013;4:269. pmid:23888162.
- 10. Huang S, Greenway H, Colmer TD, Millar AH. Protein synthesis by rice coleoptiles during prolonged anoxia: implications for glycolysis, growth and energy utilization. Annals of botany. 2005;96(4):703–15. pmid:16027131.
- 11. Komatsu S, Yamamoto R, Nanjo Y, Mikami Y, Yunokawa H, Sakata K. A comprehensive analysis of the soybean genes and proteins expressed under flooding stress using transcriptome and proteome techniques. J Proteome Res. 2009;8(10):4766–78. pmid:19658438.
- 12. Bahmanyar MA. The influence of continuous rice cultivation and different waterlogging periods on morphology, clay mineralogy, Eh, pH and K in paddy soils. Pakistan journal of biological sciences: PJBS. 2007;10(17):2844–9. pmid:19090186.
- 13. Atkinson NJ, Lilley CJ, Urwin PE. Identification of genes involved in the response of Arabidopsis to simultaneous biotic and abiotic stresses. Plant Physiol. 2013;162(4):2028–41. pmid:23800991.
- 14. Wei W, Li D, Wang L, Ding X, Zhang Y, Gao Y, et al. Morpho-anatomical and physiological responses to waterlogging of sesame (Sesamum indicum L.). Plant science: an international journal of experimental plant biology. 2013;208:102–11. pmid:23683935.
- 15. Wang L, Zhang Y, Qi X, Li D, Wei W, Zhang X. Global gene expression responses to waterlogging in roots of sesame (Sesamum indicum L.). Acta Physiol Plant. 2012;34(6):2241–9.
- 16. Wang L, Yu S, Tong C, Zhao Y, Liu Y, Song C, et al. Genome sequencing of the high oil crop sesame provides insight into oil biosynthesis. Genome Biol. 2014;15(2):R39. pmid:24576357.
- 17. Wang L, Yu J, Li D, Zhang X. Sinbase: An Integrated Database to Study Genomics, Genetics and Comparative Genomics in Sesamum indicum. Plant Cell Physiol. 2014. pmid:25480115.
- 18. Lee YH, Kim KS, Jang YS, Hwang JH, Lee DH, Choi IH. Global gene expression responses to waterlogging in leaves of rape seedlings. Plant Cell Rep. 2014;33(2):289–99. ISI:000330994700009. pmid:24384821
- 19. Wei W, Qi X, Wang L, Zhang Y, Hua W, Li D, et al. Characterization of the sesame (Sesamum indicum L.) global transcriptome using Illumina paired-end sequencing and development of EST-SSR markers. BMC Genomics. 2011;12:451. pmid:21929789.
- 20. Yang F, Li WS, Jorgensen HJL. Transcriptional Reprogramming of Wheat and the Hemibiotrophic Pathogen Septoria tritici during Two Phases of the Compatible Interaction. PLOS ONE. 2013;8(11). ISI:000327546400073.
- 21. Li R, Zhu H, Ruan J, Qian W, Fang X, Shi Z, et al. De novo assembly of human genomes with massively parallel short read sequencing. Genome research. 2010;20(2):265–72. pmid:20019144.
- 22. Wang X, Wang H, Wang J, Sun R, Wu J, Liu S, et al. The genome of the mesopolyploid crop species Brassica rapa. Nature genetics. 2011;43(10):1035–9. pmid:21873998.
- 23. Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, et al. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010;28(5):511–5. pmid:20436464.
- 24. Chen S, Yang P, Jiang F, Wei Y, Ma Z, Kang L. De novo analysis of transcriptome dynamics in the migratory locust during the development of phase traits. PLOS ONE. 2010;5(12):e15633. pmid:21209894.
- 25. Audic S, Claverie JM. The significance of digital gene expression profiles. Genome Res. 1997;7(10):986–95. pmid:9331369.
- 26. Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Meth. 2008;5(7):621–8. pmid:18516045.
- 27. McCarthy FM, Gresham CR, Buza TJ, Chouvarine P, Pillai LR, Kumar R, et al. AgBase: supporting functional modeling in agricultural organisms. Nucleic Acids Res. 2011;39(Database issue):D497–506.
- 28. Naika M, Shameer K, Mathew OK, Gowda R, Sowdhamini R. STIFDB2: an updated version of plant stress-responsive transcription factor database with additional stress signals, stress-responsive transcription factor binding sites and stress-responsive genes in Arabidopsis and rice. Plant Cell Physiol. 2013;54(2):e8. pmid:23314754.
- 29. Jiang Z, Liu X, Peng Z, Wan Y, Ji Y, He W, et al. AHD2.0: an update version of Arabidopsis Hormone Database for plant systematic studies. Nucleic Acids Res. 2011;39(Database issue):D1123–9. pmid:21045062.
- 30. Moreno-Hagelsieb G, Latimer K. Choosing BLAST options for better detection of orthologs as reciprocal best hits. Bioinformatics. 2008;24(3):319–24. pmid:18042555.
- 31. Koonin EV. Orthologs, paralogs, and evolutionary genomics. Annu Rev Genet. 2005;39:309–38. pmid:16285863.
- 32. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001;25(4):402–8. pmid:11846609.
- 33. Glynn SE, Baker PJ, Sedelnikova SE, Davies CL, Eadsforth TC, Levy CW, et al. Structure and mechanism of imidazoleglycerol-phosphate dehydratase. Structure. 2005;13(12):1809–17. pmid:16338409.
- 34. Liu FL, VanToai T, Moy LP, Bock G, Linford LD, Quackenbush J. Global transcription profiling reveals comprehensive insights into hypoxic response in Arabidopsis. Plant Physiol. 2005;137(3):1115–29. ISI:000227957200033. pmid:15734912
- 35. Christianson JA, Llewellyn DJ, Dennis ES, Wilson IW. Global gene expression responses to waterlogging in roots and leaves of cotton (Gossypium hirsutum L.). Plant Cell Physiol. 2010;51(1):21–37. pmid:19923201.
- 36. Kreuzwieser J, Hauberg J, Howell KA, Carroll A, Rennenberg H, Millar AH, et al. Differential Response of Gray Poplar Leaves and Roots Underpins Stress Adaptation during Hypoxia. Plant Physiol. 2009;149(1):461–73. ISI:000262261500048. pmid:19005089
- 37. Drew MC. Oxygen deficiency and root metabolism: Injury and acclimation under hypoxia and anoxia. Annu Rev Plant Physiol Plant Mol Biol. 1997;48:223–50. ISI:A1997XD14000010. pmid:15012263
- 38. Forsthoefel NR, Vernon DM, Cushman JC. A salinity-induced gene from the halophyte M. crystallinum encodes a glycolytic enzyme, cofactor-independent phosphoglyceromutase. Plant Mol Biol. 1995;29(2):213–26. pmid:7579174.
- 39. Deng XF, Fu FL, Ni N, Li WC. Differential Gene Expression in Response to Drought Stress in Maize Seedling. Agricultural Sciences in China. 2009;8(7):767–76. ISI:000278076500001.
- 40. Zou X, Jiang Y, Liu L, Zhang Z, Zheng Y. Identification of transcriptome induced in roots of maize seedlings at the late stage of waterlogging. BMC Plant Biol. 2010;10:189. pmid:20738849.
- 41. Zou X, Tan X, Hu C, Zeng L, Lu G, Fu G, et al. The Transcriptome of Brassica napus L. Roots under Waterlogging at the Seedling Stage. Int J Mol Sci. 2013;14(2):2637–51. pmid:23358252.
- 42. Schopfer LM, Massey V, Ghisla S, Thorpe C. Oxidation-reduction of general acyl-CoA dehydrogenase by the butyryl-CoA/crotonyl-CoA couple. A new investigation of the rapid reaction kinetics. Biochemistry. 1988;27(17):6599–611. pmid:3219356.
- 43. Kotliar AB. [Activation of complex I in the reaction of NADH oxidation and delta mu H+-dependent NAD+ reduction by succinate]. Biokhimiia. 1990;55(2):195–200. pmid:2111181.
- 44. La Camera S, Gouzerh G, Dhondt S, Hoffmann L, Fritig B, Legrand M, et al. Metabolic reprogramming in plant innate immunity: the contributions of phenylpropanoid and oxylipin pathways. Immunol Rev. 2004;198:267–84. ISI:000220069200019. pmid:15199968
- 45. Wang LH, Zhang YX, Qi XQ, Li DH, Wei WL, Zhang XR. Global gene expression responses to waterlogging in roots of sesame (Sesamum indicum L.). Acta Physiol Plant. 2012;34(6):2241–9. ISI:000310165900018.
- 46. Fini A, Guidi L, Ferrini F, Brunetti C, Di Ferdinando M, Biricolti S, et al. Drought stress has contrasting effects on antioxidant enzymes activity and phenylpropanoid biosynthesis in Fraxinus ornus leaves: An excess light stress affair? J Plant Physiol. 2012;169(10):929–39. ISI:000307413000001. pmid:22537713
- 47. Jiang Y, Deyholos MK. Comprehensive transcriptional profiling of NaCl-stressed Arabidopsis roots reveals novel classes of responsive genes. BMC Plant Biol. 2006;6:25. pmid:17038189.
- 48. Oono Y, Seki M, Satou M, Iida K, Akiyama K, Sakurai T, et al. Monitoring expression profiles of Arabidopsis genes during cold acclimation and deacclimation using DNA microarrays. Funct Integr Genomics. 2006;6(3):212–34. pmid:16463051.
- 49. Atkinson NJ, Urwin PE. The interaction of plant biotic and abiotic stresses: from genes to the field. J Exp Bot. 2012;63(10):3523–43. pmid:22467407.
- 50. Hepler PK, Hush JM. Behavior of Microtubules in Living Plant Cells. Plant Physiol. 1996;112(2):455–61. pmid:12226402.
- 51. Eisinger W, Ehrhardt D, Briggs W. Microtubules are essential for guard-cell function in Vicia and Arabidopsis. Mol Plant. 2012;5(3):601–10. pmid:22402260.
- 52. Hsu FC, Chou MY, Chou SJ, Li YR, Peng HP, Shih MC. Submergence confers immunity mediated by the WRKY22 transcription factor in Arabidopsis. Plant Cell. 2013;25(7):2699–713. pmid:23897923.
- 53. Sairam RK, Kumutha D, Ezhilmathi K, Deshmukh PS, Srivastava GC. Physiology and biochemistry of waterlogging tolerance in plants. Biol Plant. 2008;52(3):401–12.
- 54. 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. pmid:9275237.
- 55. Umezawa T, Okamoto M, Kushiro T, Nambara E, Oono Y, Seki M, et al. CYP707A3, a major ABA 8'-hydroxylase involved in dehydration and rehydration response in Arabidopsis thaliana. Plant J. 2006;46(2):171–82. pmid:16623881.
- 56. Saito S, Hirai N, Matsumoto C, Ohigashi H, Ohta D, Sakata K, et al. Arabidopsis CYP707As encode (+)-abscisic acid 8'-hydroxylase, a key enzyme in the oxidative catabolism of abscisic acid. Plant Physiol. 2004;134(4):1439–49. pmid:15064374.
- 57. Osmond BC, Specht CA, Robbins PW. Chitin synthase III: synthetic lethal mutants and "stress related" chitin synthesis that bypasses the CSD3/CHS6 localization pathway. Proc Natl Acad Sci U S A. 1999;96(20):11206–10. pmid:10500155.
- 58. Marisco G, Saito ST, Ganda IS, Brendel M, Pungartnik C. Low ergosterol content in yeast adh1 mutant enhances chitin maldistribution and sensitivity to paraquat-induced oxidative stress. Yeast. 2011;28(5):363–73. pmid:21360751.
- 59. Chang S, Puryear J, Funkhouser EA, Newton RJ, Cairney J. Cloning of a cDNA for a chitinase homologue which lacks chitin-binding sites and is down-regulated by water stress and wounding. Plant Mol Biol. 1996;31(3):693–9. pmid:8790302.