Original Research ARTICLE
The Histone Deacetylase Inhibitor Suberoylanilide Hydroxamic Acid Alleviates Salinity Stress in Cassava
- 1Plant Biochemistry and Molecular Genetics Laboratory, Department of Biotechnology, Faculty of Science, Mahidol University, Bangkok, Thailand
- 2Plant Genomic Network Research Team, RIKEN Center for Sustainable Resource Science, Yokohama, Japan
- 3CREST, Japan Science and Technology Agency, Saitama, Japan
- 4Plant Productivity Systems Research Group, RIKEN Center for Sustainable Resource Science, Yokohama, Japan
- 5Chemical Genomics Research Group, RIKEN Center for Sustainable Resource Science, Saitama, Japan
- 6Plant Genomic Network Science Division, Kihara Institute for Biological Research, Yokohama City University, Yokohama, Japan
Cassava (Manihot esculenta Crantz) demand has been rising because of its various applications. High salinity stress is a major environmental factor that interferes with normal plant growth and limits crop productivity. As well as genetic engineering to enhance stress tolerance, the use of small molecules is considered as an alternative methodology to modify plants with desired traits. The effectiveness of histone deacetylase (HDAC) inhibitors for increasing tolerance to salinity stress has recently been reported. Here we use the HDAC inhibitor, suberoylanilide hydroxamic acid (SAHA), to enhance tolerance to high salinity in cassava. Immunoblotting analysis reveals that SAHA treatment induces strong hyper-acetylation of histones H3 and H4 in roots, suggesting that SAHA functions as the HDAC inhibitor in cassava. Consistent with increased tolerance to salt stress under SAHA treatment, reduced Na+ content and increased K+/Na+ ratio were detected in SAHA-treated plants. Transcriptome analysis to discover mechanisms underlying salinity stress tolerance mediated through SAHA treatment reveals that SAHA enhances the expression of 421 genes in roots under normal condition, and 745 genes at 2 h and 268 genes at 24 h under both SAHA and NaCl treatment. The mRNA expression of genes, involved in phytohormone [abscisic acid (ABA), jasmonic acid (JA), ethylene, and gibberellin] biosynthesis pathways, is up-regulated after high salinity treatment in SAHA-pretreated roots. Among them, an allene oxide cyclase (MeAOC4) involved in a crucial step of JA biosynthesis is strongly up-regulated by SAHA treatment under salinity stress conditions, implying that JA pathway might contribute to increasing salinity tolerance by SAHA treatment. Our results suggest that epigenetic manipulation might enhance tolerance to high salinity stress in cassava.
Cassava (Manihot esculenta Crantz) originated in South America and is an important root crop, of which worldwide cultivation has progressed throughout tropical and subtropical regions (Olsen and Schaal, 1999). This perennial crop grows a starchy root, with starch making up 70–90% of the total dry weight (El-Sharkawy, 2004; Nuwamanya et al., 2008). Over 500 million people use cassava root starch as a source of carbohydrates (FAO, 1998; El-Sharkawy, 2004). Moreover, cassava has multiple applications including as animal feed and as a raw material for biofuel production, which contributes to building a sustainable ecosystem (Fu et al., 2016).
Soil salinity is one of the leading factors that hinder crop production globally, and development of cassava plants that are more tolerant to salinity stress is required (Carretero et al., 2008). It is well known that soil salinity affects plant cells in two ways: water deficit caused by high concentrations of salt in soil leading to decreasing water uptake by roots (osmotic stress); and high accumulation of salt in the plant, which alters Na+/K+ ratios as well as leading to excessive Na+ and Cl− content (ion cytotoxicity; Munns and Tester, 2008; Julkowska and Testerink, 2015). Previous studies have revealed that several mechanisms such as maintenance of ion homeostasis, accumulation of compatible solutes, hormonal control, antioxidant systems, and Ca2+ signaling are essential for plants to survive under high salinity stress (Jia et al., 2015). Based on those findings, genetic engineering and conventional breeding have been widely used to develop salt-tolerant plants. To improve salt tolerance of transgenic plants, genes involved in several pathways against salt stress have been used as targets for expression modification. These include transporters for ion homeostasis such as NHX1 (Apse et al., 1999), SOS1/2/3 (Shi et al., 2003; Yang et al., 2009), and HKT1 (Mølle et al., 2009); and for the accumulation of osmolytes such as proline (Kishor et al., 1995) and glycine betaine (Sakamoto and Murata, 1998); late embryogenesis abundant (LEA) proteins (Xu et al., 1996); and enzymes for antioxidant synthesis such as GST/GPX (Roxas et al., 1997) and SOD (McKersie et al., 1999).
In addition to osmolytes, treatment with small molecules such as plant hormones has been used to enhance salt tolerance in plants. It is reported that the use of salicylic acid (SA) increases tolerance to drought and salt stress in wheat (Shakirova et al., 2003). In Arabidopsis, β-aminobutyric acid (BABA) functions as priming effect on abscisic acid (ABA) responses for drought and salinity stress, resulted in increasing these stress tolerance (Jakab et al., 2005). These results show that application of small molecules allows the improvement of plant traits responsible for stress tolerance, particularly for crops in which it is difficult to introduce traits by transformation or crossing. Furthermore, constitutive expression of stress responsive genes often induces growth inhibition. Treatment with small molecules during the optimized period may have the advantage of minimizing the growth inhibition, because stopping small-molecules application can release from the growth inhibition.
Histone deacetylase (HDAC) inhibitors are effective small molecules against environmental stresses. HDAC controls the level of histone acetylation with histone acetyltransferase (HAT; Seto and Yoshida, 2014; Verdin and Ott, 2015). Histone acetylation is one of the histone modifications involved in epigenetic regulation, and recent evidence has increasingly revealed that the balance of histone acetylation plays a pivotal role in response to salinity stress (Kim et al., 2015). For example, the transcriptional co-activator ADA2b is a component of several multiprotein HAT complexes that contain GCN5 as their catalytic subunit. The mutant ada2b shows hypersensitivity to salt. This result suggests that ADA2b is involved in the regulation of histone acetylation levels for salt-induced genes through the modulation of HAT activity (Kaldis et al., 2010). Conversely, the loss of an HDAC, HDA9, resulted in reduced sensitivity to salt and drought stresses in Arabidopsis (Zheng et al., 2016). These studies imply that modulation of induced histone acetylation is inextricably linked to salinity stress response, although other histone modifications such as methylation are also known to be involved (Shen et al., 2014; Kim et al., 2015; Asensi-Fabado et al., 2016). Interestingly, the HDAC inhibitor, Ky-2, can increase gene expression for an ion transporter (SOS1) and an enzyme involved in osmolyte accumulation (P5CS), leading to increased tolerance to salt stress (Sako et al., 2016). Therefore, there is the possibility that HDAC inhibitors could be used to manipulate crops and obtain desired characteristics.
In this study, we attempt to enhance tolerance to high salinity in cassava using the commercially available HDAC inhibitor suberoylanilide hydroxamic acid (SAHA, vorinostat; Dokmanovic et al., 2007). The survival rate shows that SAHA helps cassava become less sensitive to high salinity stress. Immunoblotting analysis reveals that SAHA induces hyperacetylation of histones H3 and H4 in roots. Transcriptome analysis using a microarray revealed that mRNA expression of genes, involved in phytohormone ABA and jasmonic acid (JA) biosynthesis pathways, is strongly induced under high salinity stress condition. In addition to above two phytophormone pathways, SAHA treatment is likely to enhance ethylene and gibberellin biosynthesis pathways. Among them, an allene oxide cyclase (MeAOC4) catalyzing a crucial step in JA biosynthesis is strongly up-regulated. We discuss which pathway plays a pivotal role in inducing tolerance to salinity stress by SAHA treatment in cassava. Taken together, our results provide fundamental information on the high salinity stress response and suggest the possibility that the epigenetic manipulation via HDAC inhibition might be applicable for increasing tolerance to high salinity stress in cassava.
Materials and Methods
Plant Materials and Growth Conditions
Cassava cultivar TMS60444 was asexually propagated in a plastic pot with Murashige and Skoog (MS) media (pH5.8, KOH) containing 4.4 g L−1 MS salts (Murashige and Skoog, 1962), 20 g L−1 sucrose, 2 μM CuSO4, and 3 g L−1 Gelrite (Wako) under a light intensity of 40–80 μmol photons m−2s−1 with a 12 h/12 h photoperiod at 28 ± 1°C. The 5-cm length shoot tops of 3-month-old cassava plantlets were cut and then transferred to liquid MS media without agar. After 1 week incubation, the plantlets showing root elongation were subjected to each analysis as follows.
Effect of SAHA Treatment on Survival Rate and Biomass under Salinity Stress Condition
To reveal to which extent SAHA treatment alleviates cassava plants from salinity stress, survival rate and biomass (fresh and dry weight) were measured as follows. At 1 week-incubation after the transfer of cut plantlets to liquid medium, the plantlets were incubated on liquid MS media supplemented with 100 μM SAHA (suberoylanilide hydroxamic acid: N-hydroxy-N′-phenyloctanediamide: vorinostat, from Tokyo Chemical Industry, code-H1388) for 24 h. SAHA was dissolved in dimethylsulfoxide (DMSO) and DMSO treatment served as the control for all analyses. Then, the plantlets were transferred to 200 mM NaCl containing liquid MS media for 10 days and subsequently transferred to normal liquid MS medium for 1 month to observe phenotype and count survival rate in each condition. At least 13 cassava plantlets were used for each condition. The survival rate was based on the percentage of surviving plants. Plants with regenerated leaves and green stems were counted as surviving plants, while the plants with white stems having no new leaves were counted as dead ones. The fresh weight (FW) of root and shoot samples was measured immediately after harvesting. To measure the dry weight (DW), cassava samples were incubated at 60°C for 1 week. FW and DW were measured from the regenerated plants that had 3 or more leaves. Three biological replicates were conducted for all analyses.
Measurement of Na+ and K+ Contents
The cassava samples were washed for 5 s three times using distilled water. Stems and leaves were cut separately using mineral-free carbon knives (Feather 5) and kept at −80°C until use. Analytical measurement of 23Na and 39K content was conducted using an inductively coupled plasma mass spectrometer (ICP-MS; NexION300, Perkin Elmer). Firstly, the cassava samples were dried using a freeze dryer for 3 days. The cassava sample weight was measured. Then 5 mL of conc. HNO3 was added to the cassava samples. The cassava samples were incubated with concentrated HNO3 overnight. After that, the cassava samples were digested using a microwave sample preparation system (MultiWave-3000, Perkin Elmer). All samples were completely digested prior to running on the ICP-MS. The digested samples were adjusted to a volume of 50 mL with Milli-Q water, then filtered through 5B filter paper (Advantec). The ICP-MS analysis was carried out as described previously (Itouga et al., 2014). Three independent biological replicates were generated for each condition. Na+ and K+ content were analyzed statistically with a t-test.
Protein Extraction and Immunoblotting Analysis
The leaf and root samples were ground using a Multi-beads shocker (Yasui Kikai) and 200 μL of 2 × Laemmli sample buffer [2% SDS, 50.4% (w/v) glycerol, 0.02 M Tris-HCl (pH 6.8) with 2% 2-mercaptoethenol] was immediately added to the powder. The mixture was heated at 95°C for 3 min then placed on ice for 5 min, and centrifuged at 4°C for 10 min at 14,000 × g. The protein supernatant was transferred to a new tube, and separated by 12.5% SDS-PAGE gel electrophoresis. After gel electrophoresis, the protein was blotted on equilibrated Immobilon-P PVDF membrane (Millipore). Then the blotted membrane was blocked using 5% skim milk at room temperature for 1 h. Next, the membrane was incubated at 4°C overnight for the primary antigen–antibody interaction. Secondary antibodies were incubated with the membrane for 1 h at room temperature. Finally, detection was carried out using an ImageQuant LAS 4000 (GE). Primary antibodies were H3 (Abcam, 1791), H4 (Abcam, 10158), acH3 (Merck Millipore, 06-599), and acH4 (Merck Millipore, 06-866), diluted 1:5000, 1:3000, 1:2000, and 1:5000, respectively. The secondary antibody was goat anti-rabbit IgG (GE, NA931), diluted 1:25,000.
Total RNA was extracted using Plant RNA Isolation Reagent (Invitrogen) following spin-column RNA precipitation using RNeasy Plant Mini Kit (Qiagen). Briefly, the sample tissues were ground into fine powder using a Multi-beads shocker. The tissue powder was used for RNA extraction using Plant RNA Isolation Reagent according to the manufacturer's instruction (Invitrogen). To clean the RNA solution, after adding 20 μL of RNase-free water to the RNA pellet, the RNA solution was purified using spin-column RNA precipitation followed by addition of RNase-Free DNase (Qiagen). In order to remove contaminated DNA completely, the process for RNA purification using spin-column RNA precipitation with DNase treatment was repeated twice. All samples were stored at −80°C until use.
For microarray analysis, the plantlets were treated with 100 μM SAHA for 24 h, followed by treatment with 200 mM NaCl for 2 or 24 h. The quality of total RNA was evaluated using a Bioanalyzer system (Agilent). Microarray analysis was performed according to the procedure of Utsumi et al. (2012). Briefly, total RNA was labeled with cyanin-3 (Cy3) using the Quick Amp Labeling kit (Agilent Technologies). Next, the Cy3-labeled cRNA was fragmented and hybridized to the Agilent microarray (GPL14139; Utsumi et al., 2016). After hybridization, the microarrays were washed and immediately scanned on an Agilent DNA Microarray Scanner (G2505B) using one color scan setting for 8 × 60 K array slides. Feature Extraction software (Ver. 9.1, Agilent Technologies) was used to process the digital images. Three independent experiments were performed for each condition. After normalization, microarray probes were detected based on one-way ANOVA with Benjamini–Hochberg correction (FDR) < 0.0001. An AGI gene code is shown in Tables 1–6, if the protein encoded in each gene (probe ID) has the sequence similarity with those in Arabidopsis at E ≤ 1 × 10−5. The E-value for each probe ID is available in Tables S1–S5. Gene ontology (GO) analyses were carried out using agriGO at http://bioinfo.cau.edu.cn/agriGO/. The information from the microarray data is available on the GEO website (GEO ID: GSE84715).
Table 1. List of top 40 cassava genes upregulated (log2 ratio > 1; FDR < 0.0001) in roots by treatment with 200 mM NaCl for 2 h in non-SAHA-pretreated plants.
Quantitative Real-Time RT-PCR (qRT-PCR) Analysis
First-strand cDNA was synthesized from 500 ng total RNA with random primers. ReverTra Ace (TOYOBO) was used for the reverse transcription reaction according to the manufacturer's instructions. Transcript levels were assayed using Fast SYBR Green Master Mix (Applied Biosystems) and a StepOnePlus Real-Time PCR System (Applied Biosystems) according to the manufacturer's protocols. Gene-specific primers were designed using the PrimerQuest tool (http://sg.idtdna.com/primerquest/Home/Index). Melting curve analysis was conducted to validate the specificity of the PCR amplification. For qRT-PCR, three biological replicates were conducted. Actin was used as reference gene to normalize data (cassava4.1_009807m). The relevant primers are listed in Table S6. Three independent biological replicates were generated for each condition. Changes in gene expression were analyzed statistically with a t-test.
Construction of Phylogenetic Tree
An evolutionary tree was constructed using the maximum likelihood method in MEGA7 (Kumar et al., 2016). The tree was evaluated with 1000 bootstrap replicates.
SAHA Enhances Salinity Stress Tolerance in Cassava
To confirm which NaCl concentration inhibits growth in cassava, plantlets were incubated in liquid medium containing different NaCl concentrations. A concentration of 200 mM NaCl clearly inhibited the growth of cassava plantlets (Figure 1A). After transfer to 200 mM NaCl salt medium for 1 week, the plantlets' older leaves started to turn yellow. At 2 weeks after transfer to salt medium, leaves showed chlorosis at both 150 and 200 mM NaCl. After 14-d incubation with 0, 150, and 200 mM NaCl, plants showed root growth of 7.4 ± 2.1, 3.6 ± 0.8, and 1 ± 0.2 cm, respectively (Figure 1B). The 200 mM NaCl concentration was selected for further experiments.
Figure 1. Plant phenotype after treating with various NaCl concentrations. (A) Morphological observation of cassava after 14-day salt stress treatment. NaCl concentrations of 0, 150, and 200 mM were used. (B) Inhibition of primary root elongation by salinity stress. Error bars represent the means ± SD (n = 5–6).
To analyze whether SAHA treatment enhances tolerance to salinity stress in cassava, we counted the survival rate of SAHA-treated plants (Figure 2A). SAHA-treated plants showed a 30.6% higher survival rate than non-treated plants under salinity stress conditions (Figure 2B). Fresh and dry weight measurement revealed that SAHA treatment increased biomass in roots (Figures 2C,D). These results suggest that SAHA treatment improves tolerance to high salinity stress in cassava.
Figure 2. Increase in survival rate and plant weight by SAHA treatment under high salinity stress conditions in cassava. (A) Experimental condition. Cassava plantlets were treated with 100 μM SAHA for 24 h then subjected to 200 mM NaCl medium for 10 days. The survival rate was counted after transferring cassava to normal medium for 1 month. (B) Survival rate under high salinity stress of non-treated plants and SAHA-treated plants. Three independent experiments were performed (at least 13 plants per experiment). (C,D) Comparison of plant fresh weight (C) and dry weight (D) between non-treated plants and SAHA-treated plants. Asterisks indicate significantly different means (*p < 0.05, **p < 0.005) as determined with a t-test. Error bars represent the mean ± SD.
Concentration of Na+ and K+ in Stems and Leaves under High Salinity Stress
The concentration of sodium (23Na) and potassium (39K) in stem and leaf samples was quantitatively analyzed by ICP-MS to reveal the extent of SAHA treatment's influence on ion homeostasis under salinity stress conditions. SAHA-treated plants had reduced Na+ concentrations in both stems (2.9 mg g−1 DW) and leaves (0.87 mg g−1 DW) compared with those in control plants (stems and leaves: 4.36 and 1.52 mg g−1 DW, respectively; Figure 3A). However, there was no obvious difference in K+ concentrations in SAHA-treated plants and control plants, resulting in higher K+/Na+ ratios in SAHA-treated plants than in control plants (Figures 3B,C). These results suggest that SAHA treatment allows cassava to control ion homeostasis under high salinity stress conditions.
Figure 3. Na+ and K+ content in stems and leaves under high salinity stress. (A,B) Na+ and K+ concentrations in the stems and leaves of plants exposed to 100 μM SAHA then subjected to 200 mM NaCl for 2 h. (C) the K+/Na+ ratio of stems and leaves of cassava plants. Asterisks indicate significantly different means (*p < 0.02, **p < 0.05) as determined with a t-test. Error bars represent the mean ± SD. Three independent biological replicates were performed for each condition.
Immunoblotting Analysis of Histones H3 and H4 Acetylation
It has so far been unclear whether SAHA functions as the HDAC inhibitor in plants. We performed immunoblotting analysis to determine whether SAHA has an effect on histone acetylation levels. Cassava plantlets treated with 100 μM SAHA for 12 or 24 h, and cassava plantlets treated with 100 μM SAHA for 24 h then subjected to 200 mM NaCl for 6 h, were used for detecting histone acetylation status. Hyperacetylation of histones H3 and H4 was detected in the SAHA-treated plants compared with the control plants, especially in root samples (Figures 4A,C). Levels of histone acetylation were significantly increased at 12 and 24 h of SAHA treatment and remained higher up to 30 h after additional NaCl treatment. In contrast to roots, the acetylation level of histone H3 in leaf samples was slightly induced at 12 h and returned to a level similar to that of non-SAHA treated plants at 24 h (Figure 4B). There was no difference in acetylation levels of histone H4 in leaves between SAHA-treated plants and control plants (Figure 4D). These results suggest that SAHA mainly induces hyperacetylation of histones H3 and H4 in the roots, which may lead to transcriptional changes that enhance tolerance to high salinity stress in cassava.
Figure 4. Changes in histones H3 and H4 acetylation levels during SAHA-treatment. Cassava plantlets were treated with 100 μM SAHA for 12 and 24 h. The 24 h SAHA treated-plants were the subjected to 200 mM NaCl for 6 h. (A,B) Detection of histone acetylation levels by western blotting (WB) with anti-AcH3 antibody in root (A) and leaf (B) samples. (C,D) The detection of histone acetylation levels by WB with anti-AcH4 antibody in root (C) and leaf (D) samples.
Transcriptome Analysis in Response to Salinity Stress
To understand responses to salt stress and unveil the mechanism underlying increased tolerance to salinity stress under SAHA treatment in cassava, the transcriptional changes in roots and leaves were analyzed. For root samples, there were 4389 genes with significant expression changes [one-way ANOVA with Benjamini–Hochberg correction (FDR) < 0.0001] in at least one condition (Figure 5). In contrast, only 53 genes showed significant differences in mRNA expression in leaves (Table S7). These data are consistent with the results of immunoblotting analysis.
Figure 5. Hierarchical cluster analysis of salt-responsive genes and SAHA-responsive genes in cassava roots. The expression profiles of cassava genes were obtained from cassava plantlets treated with 100 μM SAHA for 24 h then subjected with 200 mM NaCl for 2 and 24 h. Transcript data was generated from 3 replicates. The heat map represents the Z-score. The key shows the Z-score region from −2 to 2. Red represents up-regulated genes while blue represents down-regulated genes. The 4389 differentially expressed genes were detected using one-way ANOVA, BH FDR < 0.0001.
The heatmap shows significant differences among treatment conditions in root samples. According to their expression profiles, genes differentially expressed in roots can be divided into five classes (Figures 6A–E). GO analysis was used to gain an overview of the functions of genes in each class. Class 1 contains SAHA-unaffected genes with functions assigned to GO:0000003 (reproduction), GO:0009791 (post-embryonic development), and GO:0008152 (metabolic process; Figure 6A). Class 2 contains genes up-regulated by SAHA treatment. GO enrichment analysis revealed that the majority of the genes in this group belong to GO:0050896 (response to stimulus), GO:0044237 (cellular metabolic process), and GO:0044238 (primary metabolic process; Figure 6B). Classes 3 and 4 contain genes down-regulated by the effect of SAHA treatment at 2 h and 24 h NaCl treatment. Salt stress-unresponsive genes and transient salt stress-responsive genes were enriched in class 3 (Figure 6C) and class 4 (Figure 6D), respectively. GO enrichment analysis indicated that GO:0051179 (localization), GO:0006810 (transport), GO:0016043 (cellular component organization), and GO:0048856 (anatomical structure development) are included in class 3, whereas GO:0051179 (localization) and GO:0006810 (transport) are included in class 4. Class 5 contains the most genes (1103 genes). Class 5 is a group of salt-responsive genes that are up-regulated at 2 h under both NaCl and SAHA treatment. Several gene ontology (GO) categories were enriched among these genes, such as GO:0006950 (response to stress), GO:0009414 (response to water deprivation), GO:0009651 (response to salt stress), GO:0009737 (response to abscisic acid stimulus), and GO:0006082 (organic acid metabolic process; Figure 6E).
Figure 6. Clustering of genes in root transcriptome. Each cluster includes genes that responded to treatment conditions: (A) Class 1, SAHA-unaffected genes; (B) Class 2, genes up-regulated by SAHA treatment; (C) Class 3, enriched salt-stress unresponsive genes down-regulated at 2 h after NaCl addition in SAHA-pretreated plants (D) Class 4, enriched salt-stress responsive genes down-regulated at 2 h after NaCl addition in SAHA-pretreated plants, (E) Class 5, salt-responsive genes up-regulated at 2 h after NaCl addition in SAHA-pretreated plants.
Genes Related to Salinity Stress Response in Cassava
As GO ontology analysis revealed, class 5 contains high salinity stress-responsive genes involved in ABA biosynthesis and signal transduction, including homologs of Arabidopsis genes NCED3 (MeNCED3: RknMes02_025528: cassava4.1_026283m, Tables 1, 2), EDL3 (MeEDL3: RknMes02_033569: cassava4.1_030886m, Tables 1, 2), ABI1 (MeABI1: RknMes02_033475: cassava4.1_020355m, Tables S1, S2), ABI2 (MeABI2: RknMes02_010484: cassava4.1_010060m, Tables 1, 2), and other PP2Cs, such as AHG1 (MeAHG1: RknMes02_003939: cassava4.1_013309m, Tables 1, 2), AHG3 (MeAHG3: RknMes02_013552: cassava4.1_008067m, Tables S1, S2), HAI1 (MeHAI1: RknMes02_002573: cassava4.1_007998m, Table 1), HAI2 (MeHAI2: RknMes02_006068: cassava4.1_007913m, Tables 1, 2), HAB1 (MeHAB1: RknMes02_048414: cassava4.1_005959m, Tables S1, S2) (Ma et al., 2006). In Arabidopsis, the salt-induced upregulation of several transcription factors including ATAF1, ATHB12, NAP, AZF2, HSF2, RD26, and ATERF4, has been reported (Ma et al., 2006; Matsui et al., 2008). In cassava, salt stress induced the expression of their homologs MeATAF1 (RknMes02_016900: cassava4.1_013132m, Tables S1, S2), MeATHB12 (RknMes02_016407: cassava4.1_015049m, Table 2), MeNAP (RknMes02_036224: cassava4.1_013467m, Tables S1, S2), and MeRD26 (RknMes02_031261: cassava4.1_010999m, Tables 1, 2). LEA proteins, whose overexpression gives tolerance to salinity stress, are regulated by ABA in land plants and thought to function as molecular shields to prevent aggregation caused by dehydration (Shinde et al., 2012). MeLEA (RknMes02_006505: cassava4.1_025947m) showed the highest the highest up-regulated gene by salt stress (Table 1, Table S1). These results suggest that ABA synthesis and its signal transduction pathway coordinate the basis of the response to salinity stress in cassava. This is consistent with previous data in Arabidopsis (Matsui et al., 2008).
Table 2. List of top 40 cassava genes upregulated (log2 ratio > 1; FDR < 0.0001) in roots by treatment with 200 mM NaCl for 24 h in non-SAHA-pretreated plants.
In addition to ABA, the hormones JA and methyl jasmonate (MeJA) are also involved in salt stress response. Consistent with previous reports that MeJA induces salt responsive genes in roots (Ma et al., 2006), our transcriptome analysis from root sample detected the induced mRNA expression of genes for JA/MeJA biosynthesis (MeACX2: RknMes02_027394: cassava4.1_002855m, Table S1; MeACX3: RknMes02_012214: cassava4.1_002966m, Table S1; and MeJMT: RknMes02_028787, cassava4.1_010155m, Table 1), suggesting that the JA signaling pathway is also involved in response to salt stress in cassava.
The expression levels of osmoprotectant biosynthesis-related genes, such as Δ-1-pyrroline-5-carboxylate synthetase (P5CS) and raffinose synthase were induced under salinity stress condition. P5CS1 (MeP5CS1: RknMes02_003952: cassava4.1_002374m) catalyzes the rate-limiting step in the biosynthesis of proline. MeP5CS1 was up-regulated at 2 and 24 h NaCl treatment (Tables S1, S2). Soluble sugars of the raffinose family have been associated with plant response to abiotic stresses. Raffinose synthase (RS) is responsible for raffinose biosynthesis. It has been reported that high salinity stress increased RS5 transcription (Egert et al., 2013; ElSayed et al., 2014). The expression of MeRS5 (RknMes02_026274: cassava4.1_002019m) was upregulated at 2 h NaCl treatment (Table S1). Proline and raffinose are likely to be synthesized and function as osmoprotectants under salinity stress in cassava.
SAHA Pretreatment-Upregulated Genes Identified by Transcriptome Analysis
We identified 421 genes whose expression was upregulated by 24 h SAHA treatment (Table S3). After pre-treatment of SAHA, 745 and 268 genes were upregulated by 2 and 24 h salt treatment, respectively (Tables S4, S5).
Salt Overly Sensitive1 (SOS1) is a gene responsible for increased salinity tolerance in Arabidopsis plants treated with the HDAC inhibitor Ky-2. The HDAC inhibitor strongly induced the expression of AtSOS1 to 5-, 3.5-, and 1.67-fold during NaCl treatment for 0, 2 and 10 h, respectively (Sako et al., 2016). Overexpression of SOS1 enhances tolerance to salinity stress in Arabidopsis (Shi et al., 2003; Yang et al., 2009) and tobacco (Yue et al., 2012). SAHA treatment enhanced MeSOS1 gene expression to 1.35-fold compared with untreated plants (Figure 7A). Although, the upregulation of MeSOS1 gene by SAHA treatment was not so significant (Figure 7B), the slight enhanced expression of MeSOS1 might contribute to the increased salinity stress tolerance by SAHA pretreatment in cassava.
Figure 7. Expression profiles of MeSOS1 (Salt Overly Sensitive 1) and MeAOC4 (Allene Oxide Cyclase 4) genes using quantitative real-time RT-PCR (qRT-PCR) analysis. (A) Cassava plantlets were treated with 100 μM SAHA for 24 h then subjected to 200 mM salt medium for 2 h. (B) Cassava plantlets were treated with 100 μM SAHA for 24 h. Root samples were collected. Asterisks indicate significantly different means (*p < 0.05, **p < 0.005, ***p < 0.001) as determined with a t-test. Actin was used as reference gene. Error bars represent the means ± SD. Three independent biological replicates were performed for each condition.
Transcriptome analysis revealed that the mRNA expression of genes, involved in phytohormone [abscisic acid (ABA), jasmonic acid (JA), ethylene, and gibberellin] biosynthesis pathways, was up-regulated after NaCl treatment in SAHA-pretreated roots (Table S8). Among them, the expression of an allene oxide cyclase gene (MeAOC4: RknMes02_051874: cassava4.1_022180m) was strongly induced by SAHA treatment (Tables 3–5). SAHA treatment enhanced MeAOC4 expression to 32.8-fold compared with untreated plants (Figure 7A). Previous studies have revealed that overexpression of AOCs can confer salinity stress tolerance to several crops (Yamada et al., 2002; Pi et al., 2009; Zhao et al., 2014). According to the current cassava genome database, AOCs constitute a small gene family (MeAOC3-1: RknMes02_049533: cassava4.1_014582m; MeAOC3-2: RknMes02_054684: cassava4.1_026961m; and MeAOC4: RknMes02_051874: cassava4.1_022180m) in cassava (Figure S1). Our transcriptome analysis revealed that MeAOC3-2 is a high salinity stress-responsive gene and the expression of MeAOC4 is not upregulated by high-salinity stress (Figure S1). SAHA treatment increased the expression of MeAOC4 to 35-fold (Figure 7B), suggesting that activation of the JA pathway mediated by overexpression of MeAOC4 might play a pivotal role in alleviating high-salinity stress in cassava.
Table 3. List of top 40 cassava genes upregulated (log2 ratio > 1; FDR < 0.0001) in roots by SAHA treatment.
Table 4. List of top 40 cassava genes with higher expression (log2 ratio > 1; FDR < 0.0001) in roots of SAHA-pretreated plants compared with non-SAHA-pretreated plants in the presence of NaCl for 2 h.
Table 5. List of top 40 cassava genes with higher expression (log2 ratio > 1; FDR < 0.0001) in roots of SAHA-pretreated plants compared with non-SAHA-pretreated plants in the presence of NaCl for 24 h.
In the case of the HDAC inhibitor Ky-2, 72.9% of the HDAC inhibitor-inducible genes are salt-responsive genes in Arabidopsis (Sako et al., 2016). In contrast, 28.3% of SAHA-upregulated genes are high salinity stress-responsive (Figure S2A), and only 27 SAHA-upregulated genes were upregulated by high-salinity stress (2 and/or 24 h NaCl treatment; Table 6). In mangrove plants, lignin accumulation functions in blocking metal-ion influx with higher suberization, suggesting that lignification can prohibit ions from flowing inside (Cheng et al., 2014). We analyzed the following three genes that are believed to play critical roles in regulating lignin accumulation: L-phenylalanine ammonia lyase (PAL), cinnamyl alcohol dehydrogenase (CAD), and caffeic acid 3-Omethyltransferase (COMT) genes (Whetten and Sederoff, 1995; Ma and Xu, 2008; Nguyen et al., 2016). Among them, the expression of only COMT was significantly enhanced by SAHA treatment (Figure S3), suggesting that the lignin accumulation might be changed by overexpression of COMT induced by SAHA treatment and block ion uptake.
Table 6. List of 27 SAHA-upregulated and high-salinity-stress-upregulated genes (log2 ratio > 1; FDR < 0.0001) in roots.
SAHA treatment decreased the expression of 141 high salinity stress-responsive genes including cell wall-related genes such as beta-galactosidases (RknMes02_003559, cassava4.1_001503m; RknMes02_030412, cassava4.1_001733m; RknMes02_034218, cassava4.1_001885m), beta-glucosidase (RknMes02_058028, cassava4.1_032518m), pectin lyase (RknMes02_051310, cassava4.1_021247m), pectinesterase (RknMes02_057990, cassava4.1_032455m), and pectin methylesterase inhibitor 4 (RknMes02_055023, cassava4.1_027551m) (Figure S2, Table S9). The decreased expression of these genes might contribute to changes in cell-wall components caused by SAHA treatment (see Discussion).
Here we demonstrated that the HDAC inhibitor, SAHA, can decrease sodium ion content particularly in the stems, resulting in increased survival rates under high salinity stress conditions in cassava. Transcriptome analysis identified candidate genes, such as an allene oxide cyclase catalyzing an essential step in the biosynthesis of JA, whose expression is strongly upregulated by SAHA. This study shows evidence that the HDAC inhibitor is an effective small molecule for alleviating salinity stress in crops, and will improve understanding of the mechanisms by which histone acetylation regulates response to abiotic stress in cassava. In this study, we identified many salinity stress-upregulated genes in cassava. This study is the first report on the transcriptome analysis using microarray under high-salinity stress in cassava and the information will be useful for the development of high-salinity stress tolerant cassava plants.
This study showed that SAHA treatment reduced Na+ concentration in both stems and leaves. Plants are able to survive high salinity stress conditions by the maintenance of K+ and Na+ homeostasis using several transporters (Zhu, 2003; Ji et al., 2013; Julkowska and Testerink, 2015). Several transporters function in the alleviation of high-salinity stress. SOS1 encodes the plasma-membrane Na+/H+ antiporter. Its activity can be detected only during salt stress conditions and is mediated by the Ca2+-responsive SOS3-SOS2 protein kinase complex (Qiu et al., 2002). When plants are exposed to high salinity stress conditions, SOS1 functions in the control of Na+ efflux to maintain ion homeostasis. When the expression of SOS1 is constitutively driven in Arabidopsis, the overexpressors become considerably tolerant to salt stress (Shi et al., 2003). In the woody plant Populus, introduction of the constitutively active SOS2 isoform with induced expression of SOS1 can enhance tolerance to salinity stress (Zhou et al., 2014). Plants also prevent excessive Na+ in cells by Na+ compartmentation by Na+/H+ exchanger 1 (NHX1). NHX1 functions to compartmentalize Na+ in the vacuole, resulting in low Na+ concentration and thus adjusting osmotic pressure to maintain water uptake. In Arabidopsis, co-expression of SOS1 and NHX1 enhances tolerance to NaCl concentrations up to 250 mM, suggesting the crucial functions of these genes. Furthermore, high-affinity K+ transporters (HKT) play an important role in Na+ exclusion from leaves in monocots and dicots. The overexpression of AtHKT1;1 improves salt tolerance in Arabidopsis (Horie et al., 2009; Mølle et al., 2009). Our transcriptome analysis did not find significant alteration of their mRNA expression (data not shown).
Screening of HDAC inhibitors that may induce higher expression of genes for ion transporters (e.g., NHX1 and HKT) is of interest as it may help to improve salt tolerance in cassava. HDACs in plants can be separated into three distinct families. The largest family is type I (RPD3-like), which consists of Zn2+-dependent deacetylases and is classified into three classes (I, II, and IV; Yang and Seto, 2008; Seto and Yoshida, 2014). They are generally conserved in eukaryotes, and there are 12 and 6 putative members in Arabidopsis (Hollender and Liu, 2008) and cassava (Figure S4), respectively. The second group of HDACs is plant-specific, and consists of HD-tuins. Four HD-tuins (HDT1-4) and 3 putative HD-tuins have been identified in Arabidopsis and cassava, respectively (Figure S4). The last group consists of homologs to the yeast Sir2 protein, which is a NAD+-dependent deacetylase. Two sirtuin proteins, SRT1 and SRT2 have been found in Arabidopsis (Hollender and Liu, 2008). Two putative sirtuins have been identified in cassava (Figure S4). HDAC inhibitors including Trichostatin A, SAHA, and sodium butyrate, whose target is Zn2+-dependent deacetylase, do not inhibit sirtuin activity (Richon, 2006). SAHA and related hydroxamic acid-based HDAC inhibitors have inhibitory effects on class I (HDAC1, 2, 3 and 8), class II (HDAC4, 5, 6, 7, 9 and 10), and class IV (HDAC11) human HDACs (Bolden et al., 2006). SAHA is a non-class selective inhibitor. Different selective inhibitors might be valuable for increasing tolerance to salt stress in cassava, because different gene sets are activated based on class or type selectivity of HDAC inhibitors and transcript abundance is also changed depending on inhibitory efficiency of each HDAC.
SAHA treatment strongly induced the mRNA expression of MeAOC4. AOC regulates a crucial step in JA biosynthesis, and the JA derivative, MeJA, alleviates salt stress in soybean (Yoon et al., 2009). The constitutive expression of AOCs can confer salinity tolerance to plants such as tobacco cell lines (Yamada et al., 2002) and wheat (Zhao et al., 2014), and references for the involvement of JA in environmental stresses can be found in Riemann et al. (2015). In light of these findings, it is highly possible that overexpression of MeAOC4 contributes to increased tolerance to salt stress under SAHA treatment in cassava. AOC enzyme activity is altered by heteromerization with its own isomers, and the specific heteromerized pair AtAOC2 and AtAOC4 shows the highest activity of the four AOC proteins in Arabidopsis (Otto et al., 2016). As tissue-specific expression of AOC genes is observed in Arabidopsis (Stenzel et al., 2012), each AOC seems to be regulated in a tissue- or environmental-stress specific manner in cassava. Our transcriptomic data reveal that MeAOC4 is not salt-responsive under our growth conditions and might be specifically expressed in a tissue such as flower or tuber that we have not analyzed. SAHA treatment induced the expression of MeAOC4, which might allow heteromerization of artificially induced MeAOC4 with salt-induced MeAOC3-2, and play a pivotal role in increased tolerance under SAHA treatment.
OP, MU, and MS designed the experiments; OP, MU, MI, YK, YU, AM, MT, CU, and HS conducted the experiments; OP, MU, YU, and AM analyzed the data; OP, MU, MY, JN, and MS wrote the manuscript.
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.
We would like to show our appreciation to Ms. E. Moriya, Ms. Y. Okamoto, and Ms. K. Mizunashi for their kind experimental support. OP was supported by RIKEN as an International Program Associate (IPA). This work was supported by the grant to MU from RIKEN and the Japan Society for the Promotion of Sciences (KAKENHI grant No. 25650126), and grants to MS from RIKEN, the Japan Science and Technology Agency (JST) [Core Research for Evolutionary Science and Technology (CREST)], and Innovative Areas Grant No. 16H01476 and Challenging Exploratory Research Grant No. 16K14832 of the Ministry of Education Culture, Sports, and Technology of Japan.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fpls.2016.02039/full#supplementary-material
Apse, M. P., Aharon, G. S., Snedden, W. A., and Blumwald, E. (1999). Salt tolerance conferred by overexpression of a vacuolar Na+/H+ antiport in Arabidopsis. Science 285, 1256–1258. doi: 10.1126/science.285.5431.1256
Asensi-Fabado, M. A., Amtmann, A., and Perrella, G. (2016). Plant responses to abiotic stress: The chromatin context of transcriptional regulation. Biochim. Biophys. Acta. doi: 10.1016/j.bbagrm.2016.07.015. [Epub ahead of print].
Carretero, C. L., Cantos, M., García, J. L., Azcón, R., and Troncoso, A. (2008). Arbuscular-mycorrhizal contributes to alleviation of salt damage in cassava clones. J. Plant Nutr. 31, 959–971. doi: 10.1080/01904160802043296
Cheng, H., Jiang, Z. Y., Liu, Y., Ye, Z. H., Wu, M. L., Sun, C. C., et al. (2014). Metal (Pb, Zn and Cu) uptake and tolerance by mangroves in relation to root anatomy and lignification/suberization. Tree Physiol. 34, 646–656. doi: 10.1093/treephys/tpu042
Egert, A., Keller, F., and Peters, S. (2013). Abiotic stress-induced accumulation of raffinose in Arabidopsis leaves is mediated by a single raffinose synthase (RS5, At5g40390). BMC Plant Biol. 13:218. doi: 10.1186/1471-2229-13-218
ElSayed, A. I., Rafudeen, M. S., and Golldack, D. (2014). Physiological aspects of raffinose family oligosaccharides in plants: protection against abiotic stress. Plant Biol. (Stuttg). 16, 1–8. doi: 10.1111/plb.12053
FAO (1998). Storage and Processing of Roots and Tubers in the Tropics [Online]. Available online at: http://www.fao.org/docrep/x5415e/x5415e01.htm (Accessed 21 July, 2016).
Fu, L., Ding, Z., Han, B., Hu, W., Li, Y., and Zhang, J. (2016). Physiological investigation and transcriptome analysis of polyethylene glycol (PEG)-induced dehydration stress in cassava. Int. J. Mol. Sci. 17:283. doi: 10.3390/ijms17030283
Horie, T., Hauser, F., and Schroeder, J. I. (2009). HKT transporter-mediated salinity resistance mechanisms in Arabidopsis and monocot crop plants. Trends Plant Sci. 14, 660–668. doi: 10.1016/j.tplants.2009.08.009
Itouga, M., Kato, Y., and Sakakibara, H. (2014). Phenotypic plasticity and mineral nutrient uptake of the moss Polytrichum commune Hedw. (Polytrichaceae, Bryophyta) during acclimation to a change in light intensity. Hikobia 16, 459–466.
Jakab, G., Ton, J., Flors, V., Zimmerli, L., Metraux, J. P., and Mauch-Mani, B. (2005). Enhancing Arabidopsis salt and drought stress tolerance by chemical priming for its abscisic acid responses. Plant Physiol. 139, 267–274. doi: 10.1104/pp.105.065698
Ji, H., Pardo, J. M., Batelli, G., Van Oosten, M. J., Bressan, R. A., and Li, X. (2013). The salt overly sensitive (SOS) pathway: established and emerging roles. Mol. Plant 6, 275–286. doi: 10.1093/mp/sst017
Jia, H., Shao, M., He, Y., Guan, R., Chu, P., and Jiang, H. (2015). Proteome dynamics and physiological responses to short-term salt stress in Brassica napus leaves. PLoS ONE 10:e0144808. doi: 10.1371/journal.pone.0144808
Kaldis, A., Tsementzi, D., Tanriverdi, O., and Vlachonasios, K. E. (2010). Arabidopsis thaliana transcriptional co-activators ADA2b and SGF29a are implicated in salt stress responses. Planta 233, 749–762. doi: 10.1007/s00425-010-1337-0
Kim, J. M., Sasaki, T., Ueda, M., Sako, K., and Seki, M. (2015). Chromatin changes in response to drought, salinity, heat, and cold stresses in plants. Front. Plant Sci. 6:114. doi: 10.3389/fpls.2015.00114
Kishor, P., Hong, Z., Miao, G. H., Hu, C., and Verma, D. (1995). Overexpression of [delta]-pyrroline-5-carboxylate synthetase increases proline production and confers osmotolerance in transgenic plants. Plant Physiol. 108, 1387–1394. doi: 10.1104/pp.108.4.1387
Matsui, A., Ishida, J., Morosawa, T., Mochizuki, Y., Kaminuma, E., Endo, T. A., et al. (2008). Arabidopsis transcriptome analysis under drought, cold, high-salinity and ABA treatment conditions using a tiling array. Plant Cell Physiol. 49, 1135–1149. doi: 10.1093/pcp/pcn101
Mølle, I. S., Gilliham, M., Jha, D., Mayo, G. M., Roy, S. J., Coates, J. C., et al. (2009). Shoot Na+ exclusion and increased salinity tolerance engineered by cell type-specific alteration of Na+ transport in Arabidopsis. Plant Cell 21, 2163–2178. doi: 10.1105/tpc.108.064568
Nguyen, T. N., Son, S., Jordan, M. C., Levin, D. B., and Ayele, B. T. (2016). Lignin biosynthesis in wheat (Triticum aestivum L.): its response to waterlogging and association with hormonal levels. BMC Plant Biol. 16:28. doi: 10.1186/s12870-016-0717-4
Nuwamanya, E., Baguma, Y., Kawuki, R., and Rubaihayo, P. (2008). Quantification of starch physicochemical characteristics in a cassava segregating population. Afr. Crop Sci. J. 16, 191–202. doi: 10.4314/acsj.v16i3.54380
Otto, M., Naumann, C., Brandt, W., Wasternack, C., and Hause, B. (2016). Activity regulation by heteromerization of Arabidopsis allene oxide cyclase family members. Plants (Basel) 5:3. doi: 10.3390/plants5010003
Pi, Y., Jiang, K., Cao, Y., Wang, Q., Huang, Z., Li, L., et al. (2009). Allene oxide cyclase from Camptotheca acuminata improves tolerance against low temperature and salt stress in tobacco and bacteria. Mol. Biotechnol. 41, 115–122. doi: 10.1007/s12033-008-9106-z
Qiu, Q. S., Guo, Y., Dietrich, M. A., Schumaker, K. S., and Zhu, J. K. (2002). Regulation of SOS1, a plasma membrane Na+/H+ exchanger in Arabidopsis thaliana, by SOS2 and SOS3. Proc. Natl. Acad. Sci. U.S.A. 99, 8436–8441. doi: 10.1073/pnas.122224699
Richon, V. M. (2006). Cancer biology: mechanism of antitumour action of vorinostat (suberoylanilide hydroxamic acid), a novel histone deacetylase inhibitor. Br. J. Cancer 95, S2–S6. doi: 10.1038/sj.bjc.6603463
Riemann, M., Dhakarey, R., Hazman, M., Miro, B., Kohli, A., and Nick, P. (2015). Exploring jasmonates in the hormonal network of drought and salinity responses. Front Plant Sci. 6:1077. doi: 10.3389/fpls.2015.01077
Roxas, V. P., Smith, R. K. Jr., Allen, E. R., and Allen, R. D. (1997). Overexpression of glutathione S-transferase/glutathione peroxidase enhances the growth of transgenic tobacco seedlings during stress. Nat. Biotechnol. 15, 988–991. doi: 10.1038/nbt1097-988
Sakamoto, A., and Murata, A. N. (1998). Metabolic engineering of rice leading to biosynthesis of glycinebetaine and tolerance to salt and cold. Plant Mol. Biol. 38, 1011–1019. doi: 10.1023/A:1006095015717
Sako, K., Kim, J. M., Matsui, A., Nakamura, K., Tanaka, M., Kobayashi, M., et al. (2016). Ky-2, a histone deacetylase inhibitor, enhances high-salinity stress tolerance in Arabidopsis thaliana. Plant Cell Physiol. 57, 776–783. doi: 10.1093/pcp/pcv199
Shakirova, F. M., Sakhabutdinova, A. R., Bezrukova, M. V., Fatkhutdinova, R. A., and Fatkhutdinova, D. R. (2003). Changes in the hormonal status of wheat seedlings induced by salicylic acid and salinity. Plant Sci. 164, 317–322. doi: 10.1016/S0168-9452(02)00415-6
Shen, Y., Conde e Silva, N., Audonnet, L., Servet, C., Wei, W., and Zhou, D. X. (2014). Over-expression of histone H3K4 demethylase gene JMJ15 enhances salt tolerance in Arabidopsis. Front Plant Sci. 5:290. doi: 10.3389/fpls.2014.00290
Shi, H., Lee, B. H., Wu, S. J., and Zhu, J. K. (2003). Overexpression of a plasma membrane Na+/H+ antiporter gene improves salt tolerance in Arabidopsis thaliana. Nat. Biotechnol. 21, 81–85. doi: 10.1038/nbt766
Shinde, S., Nurul Islam, M., and Ng, C. K. (2012). Dehydration stress-induced oscillations in LEA protein transcripts involves abscisic acid in the moss, Physcomitrella patens. New Phytol. 195, 321–328. doi: 10.1111/j.1469-8137.2012.04193.x
Stenzel, I., Otto, M., Delker, C., Kirmse, N., Schmidt, D., Miersch, O., et al. (2012). ALLENE OXIDE CYCLASE (AOC) gene family members of Arabidopsis thaliana: tissue- and organ-specific promoter activities and in vivo heteromerization. J. Exp. Bot. 63, 6125–6138. doi: 10.1093/jxb/ers261
Utsumi, Y., Tanaka, M., Kurotani, A., Yoshida, T., Mochida, K., Matsui, A., et al. (2016). Cassava (Manihot esculenta) transcriptome analysis in response to infection by the fungus Colletotrichum gloeosporioides using an oligonucleotide-DNA microarray. J. Plant Res. 129, 711–726. doi: 10.1007/s10265-016-0828-x
Utsumi, Y., Tanaka, M., Morosawa, T., Kurotani, A., Yoshida, T., Mochida, K., et al. (2012). Transcriptome analysis using a high-density oligomicroarray under drought stress in various genotypes of cassava: an important tropical crop. DNA Res. 19, 335–345. doi: 10.1093/dnares/dss016
Xu, D., Duan, X., Wang, B., Hong, B., Ho, T., and Wu, R. (1996). Expression of a late embryogenesis abundant protein gene, HVA1, from barley confers tolerance to water deficit and salt stress in transgenic rice. Plant Physiol. 110, 249–257. doi: 10.1104/pp.110.1.249
Yamada, A., Saitoh, T., Mimura, T., and Ozeki, Y. (2002). Expression of mangrove allene oxide cyclase enhances salt tolerance in Escherichia coli, yeast, and tobacco cells. Plant Cell Physiol. 43, 903–910. doi: 10.1093/pcp/pcf108
Yang, Q., Chen, Z. Z., Zhou, X. F., Yin, H. B., Li, X., Xin, X. F., et al. (2009). Overexpression of SOS (Salt Overly Sensitive) genes increases salt tolerance in transgenic Arabidopsis. Mol. Plant 2, 22–31. doi: 10.1093/mp/ssn058
Yue, Y., Zhang, M., Zhang, J., Duan, L., and Li, Z. (2012). SOS1 gene overexpression increased salt tolerance in transgenic tobacco by maintaining a higher K+/Na+ ratio. J. Plant Physiol. 169, 255–261. doi: 10.1016/j.jplph.2011.10.007
Zhao, Y., Dong, W., Zhang, N., Ai, X., Wang, M., Huang, Z., et al. (2014). A wheat allene oxide cyclase gene enhances salinity tolerance via jasmonate signaling. Plant Physiol. 164, 1068–1076. doi: 10.1104/pp.113.227595
Zheng, Y., Ding, Y., Sun, X., Xie, S., Wang, D., Liu, X., et al. (2016). Histone deacetylase HDA9 negatively regulates salt and drought stress responsiveness in Arabidopsis. J. Exp. Bot. 67, 1703–1713. doi: 10.1093/jxb/erv562
Zhou, J., Wang, J., Bi, Y., Wang, L., Tang, L., Yu, X., et al. (2014). Overexpression of PtSOS2 enhances salt tolerance in transgenic poplars. Plant Mol. Biol. Rep. 32, 185–197. doi: 10.1007/s11105-013-0640-x
Keywords: cassava, high salinity stress, epigenetics, histone modification, suberoylanilide hydroxamic acid (SAHA)
Citation: Patanun O, Ueda M, Itouga M, Kato Y, Utsumi Y, Matsui A, Tanaka M, Utsumi C, Sakakibara H, Yoshida M, Narangajavana J and Seki M (2017) The Histone Deacetylase Inhibitor Suberoylanilide Hydroxamic Acid Alleviates Salinity Stress in Cassava. Front. Plant Sci. 7:2039. doi: 10.3389/fpls.2016.02039
Received: 21 August 2016; Accepted: 20 December 2016;
Published: 09 January 2017.
Edited by:Vijay Pratap Singh, Ramanuj Pratap Singhdev Post Graduate College, India
Reviewed by:Shikha Singh, University of Allahabad, India
Asad Jan, The University of Agriculture Peshawar, Pakistan
Copyright © 2017 Patanun, Ueda, Itouga, Kato, Utsumi, Matsui, Tanaka, Utsumi, Sakakibara, Yoshida, Narangajavana and Seki. 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: Motoaki Seki, email@example.com