The Histone Deacetylase Inhibitor Suberoylanilide Hydroxamic Acid Alleviates Salinity Stress in Cassava

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) 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.
Keywords: cassava, high salinity stress, epigenetics, histone modification, suberoylanilide hydroxamic acid (SAHA) INTRODUCTION 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 Ca 2+ 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.

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 CuSO 4 , and 3 g L −1 Gelrite (Wako) under a light intensity of 40-80 µmol photons m −2 s −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 weekincubation 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

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 23 Na and 39 K 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. HNO 3 was added to the cassava samples. The cassava samples were incubated with concentrated HNO 3 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.

RNA Extraction
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.

Microarray Analysis
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).

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. Genespecific 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.
To analyze whether SAHA treatment enhances tolerance to salinity stress in cassava, we counted the survival rate of SAHAtreated 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.

Concentration of Na + and K + in Stems and Leaves under High Salinity Stress
The concentration of sodium ( 23 Na) and potassium ( 39 K) 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.

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.

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

AHG1
In addition to ABA, the hormones JA and methyl jasmonate (MeJA) are also involved in salt stress response.
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.
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 SAHApretreated 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 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. 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.
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.

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 highsalinity 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 Ca 2+ -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 . 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, coexpression 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 Zn 2+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 Zn 2+ -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 environmentalstress 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.