Genome-wide expression analysis of a rice mutant line under salt stress
Abstract
Salinity is a major environmental stress to plants. In this study, the ability of plants to tolerate salt was investigated by studying growth, physiological characteristics, and expression levels of genes related to the salt-stress response in the salt-tolerant rice mutant (Till-II-877), which was derived from γ-ray irradiation. Compared to plants grown under normal conditions, the height and root length of wild type (WT) were reduced by approximately 40 and 29% following exposure to salt stress for 3 weeks, whereas Till-II-877 line showed 29 and 23% reductions in plant height and root length, respectively. No significant changes were observed in total chlorophyll content, and the malondialdehyde content of the mutant increased less than that of the WT under salt treatment. Gene expression was compared between the WT and mutant lines using microarray analysis. An unbiased analysis of the gene expression datasets allowed us to identify the pathways involved in salt-stress responses. Among the most significantly affected pathways, changes in gene expression were observed in α-linolenic acid and linoleic acid metabolism (in lipid metabolism), fructose and mannose metabolism and glycolysis-gluconeogenesis (in carbohydrate metabolism), cysteine and methionine metabolism (in amino acid metabolism), and carbon fixation (in the energy metabolism of photosynthetic organisms) under salt stress. These results show that the differential response of plants subjected to salt stress was due to changes in multiple metabolic pathways. These findings increase our understanding of the effects of salt stress in rice and may aid in the development of salt-tolerant rice cultivars.
INTRODUCTION
Salinity is a major abiotic stress that affects plant growth and is a problem in many regions. The ability of plants to cope with salinity stress is an important determinant of crop distribution and productivity in many areas; thus, it is important to understand the mechanisms that confer tolerance to salty environments (Yu et al., 2016). Many parameters can be used to estimate salt tolerance in plant, including shoot and root growth, malondialdehyde (MDA) levels, and chlorophyll content (Hala et al., 2005; Kim et al., 2010). Salinity stress induces ionic perturbations, generation of reactive oxygen species (ROS), and cell damage. In particular, ROS generated by salt stress can cause cell damage (Hasegawa et al., 2000). MDA is a reactive and naturally occurring compound in plants, which acts as a marker of oxidative stress and these are used as makers of salt stress responses. Therefore, many studies have been carried out to the elimination of MDA content for selected salt tolerance mutants. Salinity can also cause a progressive loss of chlorophyll. Additionally, enhanced tolerance to salt stress can lead to changes in the chloroplast membrane, which mainly comprises the saturation of fatty acid (linoleic and alpha-linolenic acid) in many plants (Simopoulos, 2004).
Rice (Oryza sativa L.) is an important crop that feeds approximately half of the world’s population (Goff et al., 2002). Additionally, rice is a significant resource for genetic research on the development of animal and plant species. In Korea, rice is grown on an area of about 80 million hectares, producing 432 million tonnes with an average yield of 5420 kg per hectare. Since the year 2000, 5132 varieties of rice have been collected, and more than 300 useful cultivars have been created for national plot testing each year (Qi et al., 2009). Despite significant output, there are many restrictions on tillage. To overcome the problems of salinity, various improvement strategies have been employed. However, these strategies have proven to be impractical due to their high costs and the low return on investments. The most practical approach is to develop tolerant varieties; thus, many researchers have considered breeding for salt tolerance (Qi et al., 2009; Kim et al., 2010). This trait requires genetic variability, screening parameters, and an understanding of the genetic and physiological mechanisms of tolerance (Kikuchi et al., 2003; Moradi et al., 2003; Liu et al., 2015; Silveira et al., 2015; Yu et al., 2016).
Studies on rice genomic are very advanced because the rice genome (430 Mb) is the smallest among members of Poaceae, and rice is used as a model crop for cereal genomics research (Kondou et al., 2006). Recently, investigations into the genomes of many species and the assembly of cDNA collections have been actively promoted. In 2002, sequencing of the entire rice genome was completed, and more than 20,000 full-length cDNA clones were annotated (Goff et al., 2002). Thus, detailed analyses of genetic information have become possible using public databases, and the information obtained from cDNA is particularly useful for investigating the structure and function of genes (Kikuchi et al., 2003; Liu et al., 2015; Silveira et al., 2015; Yu et al., 2016).
Although rice is generally considered salt-sensitive, the development of cultivars with increased salt tolerance is required (Moradi et al., 2003; Kim et al., 2010). Utilization of genetic variation to improve the salt tolerance of rice may benefit from an improved understanding of the response of rice to salinity stress, and radiation-induced mutations have become an important method used to induce genetic variability (Kim et al., 2010). One aspect of the stress response occurs at the transcriptional level and involves alterations in gene expression (Moradi et al., 2003; Kondou et al., 2006). Owing to the large number of genes involved in the response to various abiotic stresses, microarrays are increasingly being used to monitor global changes in gene expressions in Arabidopsis thaliana, rice, and in other living organisms (Katsuhara et al., 2005; Kondou et al., 2006; Li et al., 2016; Luo et al., 2016).
To explain the response to salt stress, it is necessary to identify the genes involved in various metabolic pathways that undergo altered expression. The identification of differentially expressed genes (DEGs) in response to salt treatment could help us to elucidate the response to salt stress in rice. In our study on the changes that occur during the salt-stress response, we performed a microarray analysis to identify DEGs using the rice mutant line Till-II-877, which was induced by γ-ray irradiation.
MATERIAL AND METHODS
Plant materials and growth conditions
O. sativa L. Dongan (wild type, WT) were irradiated with γ-rays using a [60Co] γ-irradiator (150 TBq capacity; ACEL, Ontario, Canada) at the Korea Atomic Energy Research Institute. Among the 2961 M2 mutants derived from cv. Dongan (WT), a salt-tolerant γ-ray-irradiated mutant line, Till-II-877 (877), was selected based on the results of our previous study (Kim et al., 2010). Genetically fixed mutant lines (M6 generation) with excellent agricultural characteristics were selected. Dongan (WT) and Till-II-877 lines were planted in trays in triplicate, with 15 seeds per pot for each of the lines. The trays were filled with soil, which was commercially available for use in rice culture (Sunghwa Inc., Korea). During this period, the plant culture was achieved when the plants were grown in a controlled greenhouse at 27 ± 1°C with a light/dark cycle of 16 h light and 8 h darkness. To compare the growth characteristics of WT and 877 lines, the plants were treated with 171 mM (1%, EC = 16 mS) NaCl for 3 weeks in triplicate. Following salt treatment, whole plants were collected and stored at -80°C until the RNA was extracted for microarray analysis (Figure 1).
Comparison of mutant (Till-II-877) and wild type (WT) plants (Dongan) grown under salt stress conditions.

Biochemical measurements
To evaluate the physiological response of the rice mutant to salt stress, chlorophyll and MDA contents were analyzed in triplicate. Chlorophyll was extracted according to the method described by Lichtenaler (1987), and chlorophyll content was determined with the formulae described by Arnon (1949). Fresh leaf samples (0.5 g) were ground in 5 mL 100% acetone. The samples were centrifuged at 10,000 g at 4°C for 10 min. The supernatant fraction (1 mL) was removed and re-centrifuged for 5 min at 10,000 rpm. For pigment determination, the supernatant was diluted in 100% acetone. The absorbance of the mixture was measured at 470, 644.8, and 661.6 nm with a spectrophotometer (UVIKON 923, Bio-Tek instruments, USA).
Levels of lipid peroxidation products were estimated in triplicate, according to the method described by El-Moshaty et al. (1993). Briefly, the concentration of MDA (as an end-product of lipid peroxidation) was measured via a reaction with thiobarbituric acid (TBA). Fresh leaf samples (0.5 g) were ground in 5 mL 0.1% trichloroacetic acid (TCA). The samples were centrifuged at 10,000 g for 15 min. The supernatant fraction (1 mL) was collected and mixed with 4 mL 20% (w/v) TCA containing 0.5% (w/v) TBA. The mixture was incubated at 95°C for 30 min and quickly cooled in an ice bath for 15 min. After centrifugation at 10,000 g for 10 min, the absorbance of the extract was read at 532 nm. A correction for non-specific turbidity was made by subtracting the absorbance value obtained at 600 nm. Due to the limited specificity of this method, the concentration of TBA-reactive species (TBARS) was calculated using an extinction coefficient of 155 mM/cm, and the results were expressed as nanomoles TBARS per gram fresh weight (FW).
Microarray analysis
Total RNA was extracted from frozen material using the RNeasy Plant Mini Kit (QIAGEN, Germany), according to the manufacturer protocol. Total RNA was labeled and used for microarray hybridization, according to the manufacturer protocol (Agilent Technologies, USA). cDNA was synthesized from total RNA and labeled with cyamine 3- or cyamine 5-CTP. This cDNA was then used for hybridization, which was performed at 60°C for 17 h in a hybridization chamber (Agilent Technologies). Following previously reported methods (Song et al., 2012), the products were washed with 6X sodium chloride/sodium citrate containing 0.005% Triton X-120 for 5 min after hybridization. The microarray experiments were performed with an Axon 4000B scanner (Axon Instruments, Union City, CA, USA). A feature extraction software program (Agilent Technologies) was used to determine the signal intensities. To establish DEGs in response to salt stress in rice mutants, we performed a microarray analysis using the Affymetrix Gene Chip Rice Genome Array (Santa Clara, CA, USA).
Functional annotations
Functional clustering was performed using the Database for Annotation, Visualization, and Integrated Discovery (DAVID) Bioinformatics Resource software (v.6.7;
RESULTS
Physiological responses of Till-II-877 and Dongan (WT) rice to salt stress
Height and root length of the 877 and the WT under salt- and non-salt conditions for 3 weeks are given in Table 1. The height of the WT and 877 line was 37.7 and 37.5 cm under non-salt conditions and 22.8 and 26.5 cm under salt conditions, respectively. The height of the WT under the salt condition was reduced by approximately 40% compared to that of the WT under the non-salt condition, while the height of 877 line was reduced by 29%. The root length of the WT grown under the salt condition was reduced by approximately 29% when compared to that of the WT under the non-salt condition, while the 877 line exhibited a 23% reduction, respectively. *NSF: Non-salt treatment, **SF: salt treatment (171 mM) ***a, b, ab, c: significant differences at the 5% level (Duncan’s multiple range tests; N = 3).Growth characteristics of wild type (WT) and 877 line grown under salt and non-salt conditions for 21 days.
Treatment
Plant height (cm) ***
Root length (cm)***
Wild type
NSF*
37.7a
8.3a
SF**
22.8c
6.0c
Till-II-877
NSF
37.5a
8.9a
SF
26.5b
6.9b
Changes in chlorophyll a (

Changes in malondialdehyde (MDA) contents following 3-weeks salt stress. White and gray boxes represent non-treatment and salt treatment, respectively. WT, wild type; 877, mutant line. a, b, ab, c: significant differences at the 5% level (Duncan’s multiple range tests; N = 3).

Overview and comparison of global gene expression among rice mutant lines and the original cultivar under salt stress
DEG analysis between mutated line (877) and WT under salt and non-salt conditions were developed based on log2 ratio (Figure 4). Among 57,381 genes represented on the Affymetrix Gene Chip, 630 were up- or down-regulated 2-fold under non-treatment conditions and 1563 genes were up- or down-regulated under the salt stress condition. We found that 436 genes were up-regulated and 194 were down-regulated under non-treatment conditions, whereas 772 genes were up-regulated and 791 were down-regulated under the salt stress condition in 877. Among the DEGs, 40 up-regulated and 11 down-regulated genes were common in plants subjected to non-treatment and salt stress conditions. The expression patterns of 110 genes were altered under salt stress. Sixteen genes switched from being down-regulated to being up-regulated, whereas 94 genes switched from being up-regulated to being down-regulated in a comparison between WT and 877.
Number of genes that were differentially expressed between WT and 877 line during salt treatment.

GO analysis revealed several classes of genes that are affected by salt stress
Figure 5 is a ReviGO representation summarizing GO biological processes that were detected by DAVID as overrepresented (P ≤ 0.01) in the set of specifically expressed genes, and those that showed altered expression under salt stress. Functional categories that were enriched in specifically expressed genes included ‘fatty acid biosynthesis’ (fatty acid metabolism, fatty acid biosynthesis, oxylipin biosynthesis, oxylipin metabolism, organic acid biosynthesis, and carboxylic acid biosynthesis), ‘carbohydrate catabolism’ (carbohydrate catabolism, aromatic compound catabolism, hemicellulose metabolism, cell-wall polysaccharide metabolism, xylan metabolism, phenylpropanoid metabolism, L-phenylalanine metabolism, L-phenylalanine catabolism, and aromatic amino acid family catabolism), ‘M phase’ (cell-cycle phase, cell-cycle process, M phase of the mitotic cell cycle, M phase, mitotic cell cycle, nuclear division, and mitosis), ‘lipid transport’ (lipid localization and lipid transport), ‘oxidation-reduction process’ and ‘cell cycle’ (Figure 5A). Genes with altered expression represented the functional categories ‘oligopeptide transport’ (peptide transport and oligopeptide transport), ‘oxidation-reduction process’, and ‘glutamine family amino acid metabolism’ (Figure 5B).
Representation of the gene ontology (GO) biological processes significantly enriched with genes specifically expressed (

MapMan analysis of salt-affected genes
To permit a visual comparison of the WT and mutant lines under salt stress, specifically expressed genes and genes with altered expression were imported into MapMan. The ‘metabolism overview’ and ‘cellular response overview’ are shown in Figure 6. The metabolism overview of specifically expressed genes revealed that genes from the 877 line were involved in primary metabolism that significantly up-regulated compared with the WT line (Figure 6A). Specifically, genes involved in photosynthetic light reactions and the Calvin cycle of the photosynthesis system, starch synthesis, and sugar degradation in major carbohydrate (CHO) metabolism, glycolysis, fermentation, cellulose synthesis, pectinesterases, and cell-wall proteins involved in cell-wall metabolism, and genes involved in fatty acid (FA) synthesis, and FA elongation in lipid metabolism were more highly expressed in the mutant compared with the WT. Genes related to waxes showed increased expression, although other genes involved in secondary metabolism were down-regulated under salt stress. In the cellular response overview, genes involved in drought/salt abiotic stresses, heme in redox reactions, cell division, and the cell cycle were found to be more highly expressed in the mutants than in the WT (Figure 6B). Metabolic pathways affected by the transcriptional changes that occur during salt stress.
Differential gene expression in metabolic pathways under salt stress
To identify expression changes in genes involved in metabolic pathways, genes that were specifically expressed under salt stress were added to the gene list of the KEGG pathway detected by DAVID functional annotation (Table 2). The specifically expressed genes were involved in α-linolenic acid and linoleic acid metabolism (in lipid metabolism), fructose and mannose metabolism and glycolysis-gluconeogenesis (in carbohydrate metabolism), cysteine and methionine metabolism (in amino acid metabolism), carbon fixation in photosynthetic organisms (in energy metabolism), biosynthesis of alkaloids, biosynthesis of plant hormones, and DNA replication. The expression of lipoxygenase, which is involved in lipid metabolism, was decreased in the 877 line as compared to that in the WT. Genes involved in carbohydrate metabolism and DNA replication were more highly expressed in the 877 line. Most genes involved in cysteine and methionine metabolism were up-regulated, with the exception of Os01g0720700 and Os01g0578000 compared to WT.
Genes involved in the salt-stress response with Kyoto encyclopedia of genes and genomes (KEGG) pathway information.
Probe ID | Gene symbol | Gene name | DEG |
---|---|---|---|
α-Linolenic acid metabolism | |||
Os.8266.1.A1_at | Os03g0767000 | Allene oxide synthase 1, chloroplastic | -1.49 |
Os.15829.1.S1_at | Os02g0194700 | Lipoxygenase | 1.66 |
Os.53604.1.S1_at | Os03g0738600 | Lipoxygenase 2; Lipoxygenase | -1.38 |
Os.4416.1.S1_at | Os08g0508800 | Lipoxygenase 7, chloroplastic; Lipoxygenase | -2.12 |
OsAffx.26361.1.S1_x_at | Os04g0447100 | -1.86 | |
OsAffx.26361.1.S1_at | Os04g0447100 | -1.67 | |
Os.7513.1.S1_at | Os11g0605500 | -1.05 | |
Os.2371.1.S1_at | Os12g0448900 | -1.87 | |
Linoleic acid metabolism | |||
Os.15829.1.S1_at | Os02g0194700 | Lipoxygenase | 1.66 |
Os.53604.1.S1_at | Os03g0738600 | Lipoxygenase 2; Lipoxygenase | -1.38 |
Os.4416.1.S1_at | Os08g0508800 | Lipoxygenase 7 chloroplastic; Lipoxygenase | -2.12 |
OsAffx.26361.1.S1_x_at | Os04g0447100 | -1.86 | |
OsAffx.26361.1.S1_at | Os04g0447100 | -1.67 | |
Biosynthesis of alkaloids | |||
Os.4167.1.S1_at | Os02g0236000 | Aspartate aminotransferase | 1.40 |
Os.27606.1.S1_at | Os03g0850400 | Aspartokinase | 1.52 |
Os.5725.1.S1_at | Os11g0171300 | Fructose-bisphosphate aldolase, chloroplastic; Fructose-bisphosphate aldolase | 1.04 |
Os.47303.1.S1_at | Os01g0743500 | Malic enzyme | -1.42 |
Os.47303.1.S1_s_at | Os01g0743500 | Malic enzyme | -1.48 |
Os.57475.1.S1_x_at | Os02g0601300 | 1.42 | |
Os.23560.1.S1_at | Os04g0543900 | -1.35 | |
Os.52150.1.S1_x_at | Os04g0543900 | -1.47 | |
Os.12014.1.S1_s_at | Os05g0524400 | 2.12 | |
OsAffx.27312.1.S1_at | Os05g0524400 | 1.91 | |
Os.3350.1.S1_a_at | Os10g0390500 | 1.21 | |
Os.3350.1.S1_at | Os10g0390500 | 1.53 | |
Os.52536.1.S1_at | Os04g0518400 | Phenylalanine ammonia-lyase | -1.69 |
Os.57474.1.S1_x_at | Os04g0677500 | Pyruvate kinase | 1.22 |
Biosynthesis of plant hormones | |||
Os.8266.1.A1_at | Os03g0767000, | Allene oxide synthase 1, chloroplastic | -1.49 |
Os.4167.1.S1_at | Os02g0236000, | Aspartate aminotransferase | 1.40 |
Os.27606.1.S1_at | Os03g0850400, | Aspartokinase | 1.52 |
Os.5725.1.S1_at | Os11g0171300, | Fructose-bisphosphate aldolase, chloroplastic; Fructose-bisphosphate aldolase | 1.04 |
Os.15829.1.S1_at | Os02g0194700, | Lipoxygenase | 1.66 |
Os.53604.1.S1_at | Os03g0738600, | Lipoxygenase 2; Lipoxygenase | -1.38 |
Os.4416.1.S1_at | Os08g0508800, | Lipoxygenase 7, chloroplastic; Lipoxygenase | -2.12 |
Os.10543.1.S1_at | Os02g0109100, | 1.02 | |
Os.57475.1.S1_x_at | Os02g0601300, | 1.42 | |
Os.28127.2.S1_at | Os03g0231700, | 2.07 | |
Os.28127.1.S1_at | Os03g0231700, | 1.90 | |
Os.53086.1.S1_at | Os03g0810100, | 2.86 | |
OsAffx.26361.1.S1_x_at | Os04g0447100, | -1.86 | |
OsAffx.26361.1.S1_at | Os04g0447100, | -1.67 | |
Os.49507.1.S1_at | Os04g0578000, | -2.01 | |
Os.12014.1.S1_s_at | Os05g0524400, | 2.12 | |
OsAffx.27312.1.S1_at | Os05g0524400, | 1.91 | |
Os.56273.1.S1_at | Os06g0225800, | -1.22 | |
Os.46633.1.A1_at | Os10g0533500, | -1.74 | |
Os.12743.1.S1_at | Os11g0525200, | 1.05 | |
Os.7513.1.S1_at | Os11g0605500, | -1.05 | |
Os.52536.1.S1_at | Os04g0518400, | Phenylalanine ammonia-lyase | -1.69 |
Os.57474.1.S1_x_at | Os04g0677500, | Pyruvate kinase | 1.22 |
Carbon fixation in photosynthetic organisms | |||
Os.4167.1.S1_at | Os02g0236000 | Aspartate aminotransferase | 1.40 |
Os.5725.1.S1_at | Os11g0171300 | Fructose-bisphosphate aldolase, chloroplastic; Fructose-bisphosphate aldolase | 1.04 |
Os.47303.1.S1_at | Os01g0743500 | Malic enzyme | -1.42 |
Os.47303.1.S1_s_at | Os01g0743500 | Malic enzyme | -1.48 |
Os.3350.1.S1_a_at | Os10g0390500 | 1.21 | |
Os.3350.1.S1_at | Os10g0390500 | 1.53 | |
Os.57474.1.S1_x_at | Os04g0677500 | Pyruvate kinase | 1.22 |
Os.10959.1.S2_at | Os12g0292400 | Ribulose bisphosphate carboxylase small chain | 1.37 |
Os.12423.1.S1_at | Os01g0841600 | Triosephosphate isomerase | 1.99 |
Cysteine and methionine metabolism | |||
Os.4167.1.S1_at | Os02g0236000 | Aspartate aminotransferase | 1.40 |
Os.27606.1.S1_at | Os03g0850400 | Aspartokinase | 1.52 |
Os.10215.1.S1_at | Os02g0105400 | L-lactate dehydrogenase | 1.54 |
Os.31913.1.S1_at | Os06g0104900 | L-lactate dehydrogenase | 1.68 |
Os.49507.1.S1_at | Os04g0578000 | Os04g0578000 | -2.01 |
Os.57351.1.S1_at | Os01g0720700 | Probable serine acetyltransferase 1 | -1.15 |
DNA replication | |||
Os.7052.1.S1_at | Os02g0829100 | 1.44 | |
Os.54130.1.S1_at | Os04g0476400 | 1.14 | |
Os.53549.1.S1_at | Os05g0160800 | 1.04 | |
Os.54267.1.S1_at | Os05g0235800 | 1.14 | |
Os.12974.1.S1_at | Os11g0484300 | 1.14 | |
Fructose and mannose metabolism | |||
Os.5725.1.S1_at | Os11g0171300 | Fructose-bisphosphate aldolase chloroplastic; Fructose-bisphosphate aldolase | 1.04 |
Os.12866.1.S1_at | Os02g0714200 | 1.05 | |
Os.12014.1.S1_s_at | Os05g0524400 | 2.12 | |
OsAffx.27312.1.S1_at | Os05g0524400 | 1.91 | |
Os.10991.1.S1_at | Os06g0247500 | 1.66 | |
Os.12423.1.S1_at | Os01g0841600 | Triosephosphate isomerase | 1.99 |
Glycolysis_gluconeogenesis | |||
Os.57464.1.S1_x_at | Os11g0210300 | Alcohol dehydrogenase 1 | 1.49 |
Os.5725.1.S1_at | Os11g0171300 | Fructose-bisphosphate aldolase chloroplastic; Fructose-bisphosphate aldolase | 1.04 |
Os.10215.1.S1_at | Os02g0105400 | L-lactate dehydrogenase | 1.54 |
Os.31913.1.S1_at | Os06g0104900 | L-lactate dehydrogenase | 1.68 |
Os.57475.1.S1_x_at | Os02g0601300 | 1.42 | |
Os.313.2.S2_at | Os02g0730000 | 1.43 | |
Os.12014.1.S1_s_at | Os05g0524400 | 2.12 | |
OsAffx.27312.1.S1_at | Os05g0524400 | 1.91 | |
Os.5319.1.S1_at | Os10g0189100 | 1.00 | |
Os.57474.1.S1_x_at | Os04g0677500 | Pyruvate kinase | 1.22 |
Os.12423.1.S1_at | Os01g0841600 | Triosephosphate isomerase | 1.99 |
DISCUSSION
Rice plants are highly sensitive to salinity during the early growth stage (Kim et al., 2010; Yu et al., 2016); therefore, we investigated the salt tolerance of 877 and WT. The plant height and root length of 877 line were much higher than those of WT. Chlorophyll and carotenoid levels are considered to be indicators of salt stress in plants. Because the levels of the main photosynthetic pigments, chlorophyll a, chlorophyll b, and carotenoids, are significantly affected by salt stress, these are used as parameters to select tolerant and sensitive cultivars of crop plants (Eryilmaz, 2007; Rahdari et al., 2012). Besides, the level of MDA is a measure of membrane lipid peroxidation damage caused by abiotic stresses and it is an important indication of oxidative damage under conditions of salt stress (Liang et al., 2003). MDA is a final product of lipid peroxidation in bio-membranes, and the MDA content reflects the level of lipid peroxidation while indirectly reflecting the extent of membrane injury (Wang et al., 2009). In the present study, the chlorophyll content of the rice mutant line Till-II-877 was altered compared with that of the WT line. The total chlorophyll content was not significantly altered in 877 line subjected to salt stress for 3 weeks, although WT showed a 4.2% decrease in total chlorophyll levels. Specifically, chlorophyll b levels were markedly decreased (15.9%) in WT, whereas no significant change was observed in 877 line. The accumulation of MDA in seedlings increased with prolonged salinity stress in all plants including those of the 877 and WT lines; however, a greater increase was observed in WT (128.3%) compared with 877 (18.8%) line. Based on the changes in chlorophyll and MDA contents, we suggest that mutant rice is more tolerant to salt stress than WT rice.
Ion cytotoxicity and osmotic stress adversely affect plant growth under salt stress, and lead to metabolic imbalances and oxidative stress (Zhu, 2002; Hussain et al., 2008). Thus, under stress conditions, lipid metabolism and carbohydrate metabolism occur, and the accumulation of ROS (such as superoxide, hydrogen peroxide, and hydroxyl radicals) is induced (Hasegawa et al., 2000). In the present study, GO analysis showed that the GO terms ‘fatty acid biosynthesis’, ‘carbohydrate catabolism’, ‘M phase’, ‘oxidation-reduction process’, and ‘cell cycle’ were significantly enriched in specifically expressed genes, and ‘oligopeptide transport’, ‘oxidation-reduction process’, and ‘glutamine family amino acid metabolism’ were significantly enriched in genes exhibiting altered expression under salt stress (Figure 5). Plants of the mutant line 877, showed similar metabolic responses under salt stress as previously described. ‘Oligopeptide transport’ and ‘glutamine family amino acid metabolism’ in genes with altered expression were also reported to be related to salt stress in previous studies (Jung et al., 2010).
There were differences in the expression of genes involved in α-linolenic acid metabolism and linoleic acid metabolism between WT and 877 line in the present study (Table 1 and
Romero et al. (2001) suggested that cysteine biosynthesis is induced by salt stress in plants in order to protect against damage caused by high Na+ concentrations. Cysteine determines the cellular glutathione (GSH) concentrations, and GSH is an important factor produced by plants in response to environmental stress, including oxidative stress (Marrs, 1996; Ruiz and Blumwald, 2002). The results of the present study showed there were changes in the expression levels of genes related to cysteine and methionine metabolism (Table 1 and
Photosynthesis is modulated by the enzyme RuBisCO activase (Tezara et al., 1999; Mott and Woodrow, 2000). Under salt stress, RuBisCO activity is decreased, which leads to a reduction in the photosynthetic rate (Tezara et al., 1999; Wingler and Roitsch, 2008; Wang et al., 2009). In our study, there was differential expression of genes related to carbon fixation in photosynthetic organisms between WT and 877 line. Among these, Os12g0292400 (ribulose bisphosphate carboxylase small chain) was more highly expressed in the 877 line than in the WT (log2 ratio, 1.37) (Table 1;
Six genes involved in fructose and mannose metabolism were more highly expressed in 877 than in WT. Altered sugar content is a common response to environmental stress (Wingler and Roitsch, 2008). Sugar plays important roles as an osmoprotectant and as an immediate energy source in plants (Wingler and Roitsch, 2008). Many studies have reported increased sugar contents in plants under salt stress (Brosché et al., 2005; Gong et al., 2005; Sanchez et al., 2008; Li et al., 2016; Luo et al., 2016).
All genes involved in glycolysis were more highly expressed in 877 line compared with WT in this study (Table 1 and
Salt stress causes huge losses in agricultural productivity worldwide (Moradi et al., 2003; Yu et al., 2016). Therefore, strategies aimed at overcoming severe environmental stresses need to be quickly and fully implemented. In this study, a mutant rice line was examined for its increased salt tolerance using microarray and gene-network analyses, and we found that its differential response under salt stress was due to changes in the expression of components of various metabolic pathways. This salt-tolerant mutant rice line (Till-II-877) will further our understanding of the effects of salt stress in rice and may aid the development of salt-tolerant rice cultivars.