- Research article
- Open Access
Unraveling the molecular basis of oxidative stress management in a drought tolerant rice genotype Nagina 22
BMC Genomics volume 17, Article number: 774 (2016)
Drought stress tolerance for crop improvement is an important goal worldwide. Drought is a complex trait, and it is vital to understand the complex physiological, biochemical, and molecular mechanisms of drought tolerance to tackle it effectively. Osmotic adjustment, oxidative stress management (OSM), and cell membrane stability (CMS) are major components of cellular tolerance under drought stress. In the current study, we explored the molecular basis of OSM in the drought tolerant rice variety, Nagina 22 and compared it with the popular drought sensitive rice variety, IR 64, under drought imposed at the reproductive stage, to understand how the parental polymorphisms correlate with the superiority of Nagina 22 and tolerant bulk populations under drought.
We generated recombinant inbred lines (RIL) from contrasting parents Nagina 22 and IR 64 and focussed on spikelet fertility (SF), in terms of its correlation with OSM, which is an important component of drought tolerance in Nagina 22. Based on SF under drought stress and its correlations with other yield related traits, we used superoxide dismutase (SOD), glutathione reductase (GR), and ascorbate peroxidase (APX) activity assays to establish the relationship between SF and OSM genes in the tolerant and sensitive lines. Among the OSM enzymes studied, GR had a significant and positive correlation with single plant yield (SPY) under drought stress. GR was also positively correlated with APX but negatively so with SOD. Interestingly, none of the enzyme-morphology correlations were significant under irrigated control (IC). Through genome-wide SNP analysis of the 21 genes encoding for OSM enzymes, we identified the functional polymorphisms between the parents and identified superior alleles. By using network analysis of OSM genes in rice, we identified the genes that are central to the OSM network.
From the biochemical and morphological data and the SNP analysis, the superiority of Nagina 22 in spikelet fertility under drought stress is because of its superior alleles for SOD (SOD2, SODCC1, SODA) and GR (GRCP2) rather than for APX, for which IR 64 had the superior allele (APX8). Nagina 22 can bypass APX8 by directly interacting with SODA. For nine of the 11 genes present in the central network, Nagina 22 had the superior alleles. We propose that Nagina 22 tolerance could mainly be because of SODA which is a reactive oxygen scavenger in mitochondria which is directly associated with spikelet fertility.
Plant photosynthesis occurs in specialized organelles known as chloroplast where visible light is absorbed mainly by the pigment molecules chlorophyll a, b and carotenoids . The reaction centres of PSI and PSII present in chloroplast thylakoids are the major sites of reactive oxygen species (ROS), which are partially reduced forms of molecular oxygen; e.g. superoxide radical (O2.-), hydroxyl free radical (.OH), hydrogen peroxide (H2O2) and singlet oxygen, formed during electron transport activities in the chloroplasts, mitochondria and plasma membrane . Photoreduction of oxygen to hydrogen peroxide (H2O2) in PSI was discovered in 1951 , and later, superoxide anion (O2 .) was identified as the reduced product of molecular oxygen and superoxide dismutase (SOD), a key enzyme for the conversion of O2 .- to H2O2 and O2 . H2O2 cannot accumulate in chloroplasts since it interferes with photosynthesis by rapidly oxidizing thiol-regulated enzymes , and therefore it is reduced to water by ascorbate peroxidases (APX) with the electrons derived from water splitting in PSII . Glutathione reductase (GR), an NAD(P)H-dependent enzyme, effectively maintains a reduced pool of glutathione (GSH) levels inside the chloroplasts to prevent oxidative damage by ROS [6, 7]. Furthermore, in the ascorbate-glutathione pathway, APX and GR work synergistically to maintain ROS levels below the threshold. Catalases (CAT) convert H2O2 to water and molecular oxygen. These enzymes have extremely high maximum catalytic rates but low substrate affinities, and the reaction requires the simultaneous access of two H2O2 molecules in the active site [8, 9].
Disequilibrium between ROS production and antioxidant protection causes oxidative stress, and the accumulating ROS causes peroxidation of lipids, denaturation of proteins, mutation in DNA sequence and various other types of cellular oxidative damage [10–13]. Enhanced generation of ROS has been found in plants under a variety of abiotic stresses, such as salt, drought, chilling, heat, and metal or metalloid stresses [14–18].
Rice is one of the most widely consumed food crops, and is the most sensitive to drought stress during reproductive development [19–21], as rice spikelet fertility is extremely sensitive to water stress. In the absence of drought, yield potential determines grain production, but as drought becomes severe, drought escape or drought tolerance becomes important . Spikelet fertility is directly related to yield and is known to be sensitive to water and heat stress in rice [23–25]. The upland rice variety, Nagina 22, is an internationally recognised hardy aus genotype frequently used in drought and heat tolerance studies [23, 24, 26]. It is reported to have comparatively less spikelet sterility under drought stress because of a more efficient antioxidant defence system for scavenging of H2O2 in panicles [27, 28].
Rice, besides being a staple food crop, is also a model crop for genomics, particularly in monocots, owing to its relatively smaller genome size and the availability of high quality genome sequences  and a wealth of genomic resources [30–32]. Numerous transcription factors and genes are implicated to be differentially expressed in rice under drought stress , and some of these genes been demonstrated to increase drought tolerance [34–36]. There have been multiple efforts to identify quantitative trait loci (QTLs) related to yield under drought tolerance and their subsequent introgression in major rice varieties [37–39]. Although complex traits such as yield are routinely dissected into their component traits namely grain size, test weight, and number of productive tillers per plant in rice , and functional markers have been developed [41, 42] for plant breeders, the same is not true in drought stress research.
Drought resistance can occur through either through drought avoidance (better root system, waxy leaves, etc.,) or other drought tolerance mechanisms such as better water use efficiency and cellular tolerance. Osmotic adjustment (OA), oxidative stress management (OSM), and cell membrane stability (CMS) are the major components of cellular drought tolerance. So far, only limited literature is available in rice in the detection of QTLs and their fine mapping and cloning of the causal genes associated with these physiological/biochemical traits related to drought tolerance. Some researchers have identified QTLs for CMS , osmotic adjustment [44–46], and water use efficiency  and all of these were coarse or preliminary mapping studies. However, the components associated with drought tolerance (OA, OSM and CMS) remain as associations for which neither the genetic basis nor the molecular basis has been unequivocally established. Since the enzymes governing OSM are well known, we decided to explore their role in spikelet sterility under drought stress in a mapping population of rice. This will validate any associations by using a genetical genomics approach.
In the present work, we have used a recombinant inbred population (RIL) derived from the genetic cross between the drought tolerant Nagina 22 and the rice mega variety IR 64, which is sensitive to drought stress at the reproductive stage. We have identified the extreme bulks for spikelet fertility under drought stress from this population and used them to investigate the mechanism of OSM.
Plant material and growth conditions
Rice genotypes, Nagina 22 (drought tolerant), IR 64 (drought sensitive) and their 281 RILs, (F9:10 generation) were phenotyped for drought tolerance at the Experimental Farm, Indian Agricultural Research Institute (IARI), New Delhi (28°38′N and 77°12′E, 293 m elevation above sea level) during the 2011 and 2014 Kharif seasons. Nagina 22 was chosen as the tolerant parent because of its international importance and positive attributes such as better drought and heat tolerance . IR 64 which is sensitive to drought stress at reproductive stage, was chosen as the other parent because it is a mega variety and regularly used by rice breeders for its improvement through marker- assisted backcross breeding . Staggered nursery sowing at 5-days interval was done to synchronize the flowering time of RILs as closely as possible. This was done based on the flowering data of RILs recorded in the previous two generations (data not shown). The experiment was done in an augmented block design with check lines repeated after every 16 RILs. The check lines used were the rice cultivars Nagina 22 and IR 64 (parents) in both seasons. Lines were planted in plots of 2 m × 0.6 m with a spacing of 20 cm × 20 cm between hills and rows, both in irrigated and water stress conditions. In 2014, a subset of RILs (200) was planted in an augmented design under a rainout shelter at IARI in 3 m × 0.20 m plot size with 15 cm × 20 cm spacing. Standard cropping practices were adopted for proper nutrient and weed management. All plots were maintained at field capacity and water stress was imposed by withholding irrigation at the initiation of the booting stage to 10 days after anthesis in the majority of the lines. Phenotyping for physiological traits was done after the RILs showed leaf rolling and drying, whereas yield related parameters were measured at maturity.
On the basis of phenotype data recorded on spikelet sterility (SS), and days to 50 % flowering (DTF) under drought stress in the 2011 Kharif, 41 RILs were selected for validation of their SS phenotype. In 2014, these 41 RILs were carefully phenotyped in the field for SS and yield related traits under stress. From these data, two extreme groups termed ‘tolerant bulk’ (TB) and ‘susceptible bulk’ (SB) were selected. A third group, termed ‘random bulk’ (RB) was created from the rest of the RIL population using a random number method. Each group comprised of eight RILs and their flag leaves were collected and stored at −80 °C from both the irrigated and drought treatment plots. All the observations were individually measured in the lines constituting the bulks with adequate sampling (five plants per genotype) and averaged to obtain the trait means of the respective bulks.
Observation of physiological and agronomic parameters
Chlorophyll content of flag leaves was measured using the SPAD-502, a portable chlorophyll meter, which measures the greenness or the relative chlorophyll concentration of leaves. Plant height, panicle number, tiller number, number of spikelets per panicle (SPP), and single plant yield (SPY) were measured from five plants per RIL. SS was calculated as ratio of unfilled spikelets to total spikelets from at least three main panicles and was expressed as percent (%). To assess yield reduction under stress compared to irrigated environments, the Stress Susceptibility index (SSI) was calculated using the formula, SSI = 1-[(Yi)s/(Yi)ns]/1-[(Yps)/(Ypns)], where Yi, Yp, s, and Yns represent individual yield, population mean yield, stress condition, and non-stress conditions, respectively .
Preparation of enzyme extracts
Enzyme extracts for superoxide dismutase (SOD; EC 126.96.36.199), Glutathione reductase (GR; EC 188.8.131.52), Ascorbate peroxidase (APX; EC 184.108.40.206) and Catalase (EC 220.127.116.11) was prepared from the flag leaf tissue samples that were pre-frozen in liquid nitrogen to prevent proteolytic activity, and ground in 3 ml extraction buffer containing 0.1 M phosphate buffer (pH 7.5) and 0.5 mM EDTA. Extracts were centrifuged for 20 min at 15,000 g and the supernatant was used in the enzyme assay . The Bradford method was used to estimate protein concentrations of all the samples before performing enzyme assays. Each enzyme assay was scaled down to a 200 μl volume and all the spectrophotometric measurements were done in a 96 well plate reader (Varioskan™, Thermo Scientific,USA). Enzyme assays were performed using three biological and three technical replicates. For each enzyme assay, at least eight readings at an interval of 30 s were collected and simple regression analysis was performed on optical density (O.D) vs. time to calculate the change in absorbance per minute. Means of biological and technical replicates were collected to calculate final change in OD per minute. SOD activity in the samples was estimated by recording the decrease in the OD of nitro-blue tetrazolium dye by the enzyme . A volume of 200 μl reaction mixture contained 13 mM methionine, 25 mM nitroblue tetrazolium chloride (NBT), 0.1 mM EDTA, 50 mM phosphate buffer (pH 7.8), 50 mM sodium carbonate and 10 μl enzyme extract. The reaction was started by adding 2 mM riboflavin and placing the tubes under two 15 W fluorescent lamps for 15 min. The reaction was stopped by switching off the light and keeping the tubes in the dark. A complete reaction mixture without enzyme, which gave the maximal colour, served as a positive control. A non-irradiated complete reaction mixture served as a blank. Separate controls that lacked enzymes were used for total SOD and inhibitor studies. The absorbance was recorded at 560 nm, and one unit of enzyme activity was taken as the amount of enzyme that reduced the absorbance reading to 50 % in comparison to control tubes. APX was assayed by recording the decrease in optical density because of ascorbic acid at 290 nm . The 200 μl reaction mixture contained 50 mM potassium phosphate buffer (pH 7.0), 0.5 mM Ascorbic Acid, 0.1 mM EDTA, 1.5 mM H2O2 and 10 μl enzyme extract. The reaction was initiated with the addition of H2O2. The molar extinction coefficient of 2.8 mM−1 cm−1 was used for the calculation of enzyme activity. Absorbance was measured at 290 nm in a UV-visible spectrophotometer. GR was assayed by recording the increase in absorbance in the presence of oxidized glutathione and DTNB (5,5-dithiobis-2-nitrobenzoic acid) . The 200 μl reaction mixture contained 200 mM potassium phosphate buffer (pH 7.5) containing 1 mM EDTA, 1.5 mM DTNB (5,5-dithiobis-2-nitrobenzoic acid) in 0.01 M potassium phosphate buffer (pH 7.5), 0.2 mM NADPH, 10 μl enzyme extract and distilled water to make up the final volume. The reaction was initiated by adding 0.2 mM GSSG (oxidized glutathione or glutathione disulphide). The increase in absorbance at 412 nm was recorded and 6.22 mM−1cm-1 was used as the molar extinction coefficient for calculation of enzyme activity. Catalase was assayed by measuring the disappearance of H2O2 . Reaction mixture (200 μl) consisted of 10 μl of dilute enzyme extract and 100 μl of 0.1 M phosphate buffer (pH 7) and 60 μl of water. The reaction was initiated by adding 30 μl of 75 mM H2O2. A decrease in absorbance at 240 nm was observed every 30 s for over 3 min with UV- visible spectrophotometer.
Sequence analysis of Nagina 22 and IR 64 SOD, APX, GR and CAT genes
The SOD, APX, GR, and CAT gene CDS and protein sequences from Nagina 22 and IR 64 were retrieved from the Manually Curated Database of Rice Proteins . Nagina 22 sequences were aligned to IR 64 sequences using CLUSTALW and non-synonymous and synonymous single nucleotide polymorphisms (SNPs) were identified using BioEdit v7.2.5 software . I-mutant3  server was used to predict the effect of identified SNPs on protein stability. Network analysis in RiceNet v2  was performed using ‘gene prioritization based on network direct neighborhood’ option, and viewed using Cytoscape v3.3.0  for SOD, APX, GR and CAT genes.
Data analysis and visualization
Majority of the data analysis was done in an open source R package  for statistical computing. Multivariate exploratory data analysis was done in FactomineR , ggbiplot , corrplot  and, Hmisc . PerformanceAnalytics  was also used for data analysis and visualization. ANOVA for the enzyme assay was performed using the OP STAT web server (http://18.104.22.168/opstat/default.asp).
Population performance under irrigated conditions and drought stress conditions at the reproductive stage
Summary statistics for the agronomical and physiological attributes measured in parents, the entire RIL population, extreme RILs identified, and the three bulks namely, the tolerant, sensitive and random bulks under both the irrigated condition (IC) and reproductive stage drought stress (RS) condition during the two cropping seasons in 2011 and 2014 are presented in Tables 1 and 2. Significant differences were observed between the treatments across the traits, especially for plant height, number of spikelets/panicle, spikelet sterility and single plant yield in the parents and in the entire population. Chlorophyll content (Table 2) measured as SPAD value, was better in IR64 under both the treatment conditions compared to Nagina 22. Tiller number (TN) showed variation in the tolerant parent and the tolerant bulk between IC and RS. Conversely, neither the sensitive parent nor the SB showed any difference in TN between IC and RS. There was an increase in the coefficient of variation (CV) for plant height, panicle length and single plant yield (SPY) under RS compared to IC in both years. The mean of the 41 extreme RILs identified was nearly identical to the population mean (Table 2), across all of the traits; moreover, there was an increase in the variation, suggesting that the identification of the extreme lines was robust. Although the extreme lines were identified based only on SS, the other traits also followed the same trend presumably because of the in-built trait correlations. This trend was more evident in single plant yield and SSI (Tables 1 and 2), perhaps because the trait correlations for these traits are more robust. Interestingly, the random bulk identified, consisting of just eight lines, had a mean identical to the population of 281 lines across all the eight traits, except for number of spikelets/panicle. Overall, the morphological observations collectively established the differential response of the parents and RILs in the mapping population for drought stress, reflecting their diverse genetic background and the robust relationship between the target trait and the stress. The differential performance of bulks as well as the parents with respect to SS under the IC and RS suggested that they could be used as a background to study genetical genomics.
The enzyme assays for all four enzymes (SOD, APX, GR, and CAT) were performed in the parents as well as in the constituent RILs of all the three bulks. However, CAT results were found to be error-prone despite repeating the experiment three independent times. Hence, the results of only three enzyme assays are presented in Fig. 1.
Although, a general increase in SOD activity was observed in all the bulks and in the parental genotype under drought stress (Fig. 1a), the activity was more pronounced in TB, SB and Nagina 22 as compared to RB and IR 64 under stress situation. Thus, though distinct differences in SOD activity could be found between the parents, the same was not reflected in the extreme bulks.
Ascorbate peroxidase activity significantly decreased (Fig. 1b) across all the three bulks selected from the RIL population as well as in the drought tolerant parent, Nagina 22, under RS compared to IC. The quantum of decrease in activity was the least in Nagina 22 compared to SB and RB. However, in the other parent, IR 64, APX activity significantly increased under RS compared to IC. This could be because either APX is not a key player in OSM under drought stress or because IR 64 may have the better allele(s) for one or more of the APX genes.
Although GR activity was found to be enhanced in TB and in both of the parents under RS (Fig. 1c), the quantum of increase was much higher in the tolerant parent followed by TB. The other extreme bulk, SB showed a decreased activity under stress, whereas in RB, the difference in activity under IC and RS was negligible. Thus, GR was the only enzyme that had the expected activity in the contrasting parents and their extreme bulks.
The overall results of the enzyme assays though showed that different parents had better activity for different enzymes under stress, by comparing the quantum of ROS activity, Nagina 22 was superior to IR 64, at least for GR and SOD. The direction and quantum of expression in the extreme bulks established that there is a clear and robust correlation between the target trait SS and the GR activity in oxidative stress management under drought stress.
Principal component analysis and linear correlations
Since enzyme assay are the cumulative effect of multiple genes (as well as isoenzyme forms) that encode the enzymes in the genome, we undertook multivariate analysis. The multiple morphological and biochemical observations analysed for their possible correlated attributes by principal component analysis (PCA) and correlation analysis were expected to clarify these relationships because the observation of individual plants could be considered instead of the mean of the bulked RILs. Dimensionality of the data was revealed by an individual factor map (Fig. 2) which showed that the parents, TB, SB and RB had differentiated into distinct groups under RS compared to IC. Under RS, there was absolutely no overlap between the extreme bulks, whereas the random bulk was distributed across the extreme bulks thus, validating the suitability of the material under study. To capture the variability and interactions present in the morphological, physiological, and biochemical traits under different treatments, a variable factor map (Fig. 3) and correlogram (Fig. 4) were constructed. The first two principal components captured 24.1 % and 18.2 % of the variation under IC and 29.3 % and 17.8 % under RS, respectively. Under IC, TN, SPY, and the enzymes GR and SOD were found to interact and were highly associated to PC-1 while SPP and SS were associated with PC-2. This is consistent with our observation that the performance of TB and Nagina 22 was similar for TN, GR, and SOD (Tables 1, 2 and Fig. 1a and c). Under drought stress, the variability in APX, SPP, TN and SPAD were primarily captured by PC-1, whereas PC-2 explained the interactions for SOD and SSI. This suggests that APX is also an important component under drought stress, although Nagina 22, TB, and SB did not show any difference between the treatments (Fig. 1b). Thus, the multivariate analyses more clearly resolved the interactions between the enzyme assays and morphological traits under IC and RS.
Linear trait correlations including the enzyme assay from the parents and the RILs constituting the three bulks (Fig. 4) and the entire population were analysed individually (Additional file 1: Figure S1) under the two treatments. Since only the extreme RIL set, consisting of 26 individuals including the parents, had enzyme data, enzyme-morphological correlations were calculated using this data set instead of the entire population under IC and RS. Correlations among OSM enzymes revealed that all three had significant and positive correlations under IC (GR-APX was not statistically significant), while under RS, GR had a significant and positive correlation with APX (Fig. 4a and b and Additional file 2: Table S1). However, SOD was negatively correlated with both APX and GR under RS. This could be the reason why the extreme bulks had the expected differences for GR and APX enzyme assay but not for SOD (Fig. 1). Among the enzyme-morphological correlations, only the correlation between SPY and GR was significant and positive under RS; and GR was significantly correlated with APX and SOD, albeit in the opposite direction. APX and SPAD readings were positively correlated under RS while negatively so under IC. Interestingly, under IC, none of the 21 pairs of correlations was significant between the OSM enzymes and the morphological traits. This clearly shows that the OSM pathway is significant only under drought stress and not under well-managed conditions. The target trait, SS also showed a robust correlation with the major productivity traits, SPY and SSI (Fig. 4a and b). Correlation analysis in IC and RS done in the population clearly defined the associations among the agronomic traits (Additional file 1: Figures S1a and S1b). For instance, under irrigated control conditions, PH was positively correlated to PL, SS, and SPP, and inversely related to SPY. Similarly, TN and SS were positively correlated under IC. However the hierarchy and the direction of correlation completely changed under RS. A positive relationship was observed between SPY, PH, SS and SPAD, negative correlation was observed between SPY and SSI under RS (Additional file 1: Figure S1b).
SNPs and corresponding changes in proteins encoded by OSM genes in the parents
There was a total of 21 OSM genes in the rice genome, of which 7, 8, 3 and 3 code for SOD, APX, GR, and CAT, respectively. To investigate allelic differences in the multigenic ROS enzymes (allozymes) between the parents, SOD, APX, GR, and CAT gene sequences (18 of 21 genes; for two SOD and one CAT encoding genes, the sequences were not available) for the parents were retrieved. Both the CDS and protein sequences of SOD, APX, GR, and CAT were compared. Comparison of the CDS revealed that one-third of the ROS genes (7/21; two each from SOD, APX, and GR, and one from CAT) had no SNPs. Interestingly, most of the SNPs found resulted in non-synonymous amino acid substitutions in the proteins, especially in SOD and GR genes, which accounted for eight out of nine changes. In case of APX encoding genes, 45 % of SNPs (5/11) resulted in synonymous substitutions. In APX5, SNP substitution resulted in a missense mutation; still the resulting amino acid substitution (isoleucine to aspargine) did not alter the stability of the protein. In another APX gene (APX4), there was a long deletion (61 bp) as well as a frame shift mutation caused by a single nucleotide insertion that resulted in a completely different protein in IR 64. Overall, these changes led to a more thermodynamically stable SOD, APX, GR, and CAT proteins in Nagina 22 (nine out of eleven cases; Table 3 and Fig. 5). In only one gene APX8, the IR 64 allele had a more stable APX protein. The unfavorable allele from Nagina 22 for APX8 was probably inherited by most of the plants that constituted the RB, SB, and TB. That is why APX activity in all the three bulks was similar to the Nagina 22 parent (Fig. 1a). An interesting feature observed was the occurrence of two point mutations in the same codon at c.169, 171 positions in a SOD coding gene (SODA) in IR64, which is rare. This resulted in the acidic to basic amino acid substitution -p.D57Q mutation that, in turn, resulted in a decrease in protein stability in the drought sensitive parent IR 64. Furthermore, the mitochondrial target of this gene makes it an important candidate in this spikelet sterility study. A SNP found between the parents in GRCP2 but none of the other GR genes, corresponded perfectly with their enzyme assay results (Table 3 and Fig. 1). Although we did not have the CAT enzyme data, no SNPs were identified in CATB whereas CATA had three SNPs, all giving rise to the superior Nagina 22 allele.
Network analysis of genes involved in OSM in rice
Our enzyme studies and multi-factor analyses suggested that all of the three enzymes are important players with respect to SS under drought stress, and all the SNPs between the parents in OSM genes, except one, led to more stable proteins in the drought tolerant parent Nagina 22; therefore we carried out a functional network analysis of the genes involved in OSM in rice. An interaction analysis showed that while all of the SOD, GR, and CAT genes were present in the network, for APX only three out of the eight genes were present (Fig. 5 and Additional file 3: Table S2). APX3, APX4, and APX8 were found to interact with SODCP, CATC, SOD2, and GRCP1, respectively. APX8 had more stability in IR 64 compared to the drought tolerant parent Nagina 22. Whatever limited drought tolerance exhibited by IR 64 could be because of this APX8 gene. In Nagina 22, SODA could bypass the less stable form APX8 and directly interact with four more SOD genes, and the two GR genes present in the central network. Thus the lesser stable form APX8, does not negatively affect the drought tolerant parent in terms of its ultimate phenotype, SS or yield under RS. Of the 11 genes present in the central network, Nagina 22 had more stable proteins encoded by GRCP2, SODA and SODCC1, APX4, and CATA than IR 64; whereas, both parents shared the same allele for four more genes (one SOD, two GR, and one CAT; Fig. 5 and Table 3). For the two other genes present in the network (one SOD and CAT each), although we did not have the SNP data, the network analysis revealed that the CAT gene interacted with only APX4, which had the functional allele in Nagina 22 but not in IR 64 because of a frame-shift mutation. Thus, the differences in the enzyme activity in the parents and extreme bulks under drought stress (Fig. 1) could be explained by the SNP changes and resultant protein stability in the parents, as suggested by the network analysis of the genes involved in the OSM pathway, and the multivariate analyses.
Breeding for drought tolerance in rice through introgression of QTLs that govern yield under drought [65, 66] is still an open area of research, as the causal genes are not yet known. Hence, one of the major challenges that researchers face is to unravel the complex mechanisms of drought tolerance through more intensive and integrated studies in order to find the key players or machineries that can be effectively exploited for developing drought tolerant crops . Oxidative stress management (OSM) is one such key machinery or component of cellular tolerance under drought stress in different crop plants [68, 69]. Since rice has abundant genomic resources, it is possible to design studies without the need to invest huge financial resources, just by effectively combining forward and reverse genetic approaches. In this study, we made use of a mapping population we developed for the identification of QTLs of components for drought tolerance in rice, as well as utilized the genomic resources available in the public domain. One of the parents, Nagina 22, is tolerant to drought and known to have superior OSM  and higher spikelet fertility [23, 24] under drought.
We chose to use enzyme activity assays rather than real time transcript expression analysis for three reasons: (1) enzyme assays are more cost effective; (2) enzyme activity is a cumulative measure of activity of multiple genes; and (3) enzyme assays are closer to the phenotype than the transcript profiles to the phenotype, and our hypothesis for this work was that the biochemical or physiological component of phenotypes for complex traits should be identified along with morphological phenotypes through a genomic route that would allow for the identification of functional polymorphisms, and establish genotype-phenotype associations unequivocally in a Mendelian population.
The demonstration of the effect of drought stress on yield component traits and OSM enzyme activity under field conditions is not well documented in the available literature, especially in mapping populations. To our knowledge, this is the first report wherein enzyme activities have been measured under drought conditions imposed in the field and in well phenotyped extreme individuals following the principle of bulked segregant analysis in a synthetic Mendelian population [70, 71]. We decided to carry out the enzyme assays in extreme bulks for two reasons: (1) it would be difficult to do the enzyme assay in parents and all of the 281 RILs from the mapping population because biochemical assays need to be completed in sufficient biological and technical replicate: and (2) our target trait is well defined. We added a third bulk, called a random bulk, to reflect the scenario in the rest of the population.
Water stress caused significant reduction in PH, chlorophyll content, SF, and SPY in the RIL population including Nagina 22, IR 64, and the bulks (Tables 1 and 2) compared to the irrigated control. Although there was an overall decrease in productivity related traits in the population, Nagina 22 and TB showed more tolerance than did IR 64 and SB under drought. A stress susceptibility index that is a measure of drought resistance  showed Nagina 22 and TB to be better performers under drought stress than IR 64 and SB. RB showed an intermediate tolerance to drought as expected.
Water deficit causes a decline in photosynthetic rate and a limitation in carbon dioxide fixation, resulting in the exposure of chloroplasts to excess excitation energy and leading to the enhanced production of ROS, especially O2− and H2O2 . Drought stress caused increased SOD activity across all genotypes (Fig. 1a) compared to the irrigated control, indicating that high levels of H2O2 production and successful imposition of drought stress in field conditions. Lower APX activities in all bulks and Nagina 22 but not IR 64 (Fig. 1b) was observed under drought compared to the control, clearly showing an exhausted ascorbate pool because of drought stress. Increased GR activity in Nagina 22, IR 64, and TB (Fig. 1c) indicated an increased demand for reduced glutathione to combat the oxidative stress induced by drought. Both APX and GR work together to maintain the ascorbate-glutathione cycle, and the differential regulation of APX and GR enzymes in Nagina 22, IR 64, and their bulks under drought and control conditions is evidence that when one of the components of antioxidant defence is limiting it can be compensated for by the regulation of other components , suggesting a common regulatory mechanism of OSM genes  in rice.
Despite the compensatory mechanism of OSM, the enzyme activity patterns matched well with SNPs found in the OSM enzyme coding genes and their protein stability analyses (Fig. 1 and Table 3). For SOD and GR, the tolerant parent had higher activity under stress when compared to the sensitive parent. This is consistent with the observation that for two of the three GR enzymes, the parents had the same allele and for third one, the tolerant parent had the superior allele. Thus, the parents were monogenic with respect to GR activity. This explains why the direction of enzyme activity was the same while the quantum was different between the parents. Interestingly, the TB and SB behaved as though they have accumulated either all the favourable alleles or unfavourable alleles, respectively (Fig. 1c). This may account for the positive and significant GR-trait correlations (Additional file 2: Table S1). Such behaviour of extreme bulks could not be explained by simply accounting for differences in GR, but are explained when considering other OSM loci. In case of SOD, for two of the seven genes, the parents had identical alleles, and for the rest of the SOD genes, the tolerant parent had the superior allele. Although, the parents behaved in accordance with this allelic composition, the extreme bulks were identical to each other (Fig. 1a) and this explains why the SOD-trait associations were not significant (Additional file 2: Table S1). Such a complex system is bound to arise when multiple genes govern a target trait. Still, our network analysis along with the SNP analysis clearly showed the cross-talk among SOD, APX, GR, and CAT enzymes and for most of the key players (9/11) that are central to this interaction, the tolerant parent had the favourable alleles. Further expression analyses of all OSM enzyme encoding genes during RS and IC with leaf and root samples might suggest other clues on the difference in drought stress response between Nagina 22 and IR64. However, improving OSM is only one of the multiple factors that govern drought tolerance, especially SF under drought. Recently, the OsCPK9 gene in rice, which belongs to the group III-b CDPK family, has been demonstrated to be involved in controlling SF under drought stress . In RNA interference studies, OsHXK10 has also been implicated in spikelet fertility by its role in pollen germination and anther dehiscence .
For most of the work on drought tolerance breeding in rice, a Nagina 22 x IR 64 population is available, and now markers can be designed based on the SNP information from our study, and used for marker-assisted selection of superior alleles of OSM genes. An ideal drought tolerant cultivar shall possess all of the superior alleles known thus far from the literature and efforts must continue to identify the yet unknown partners.
The present study was carried out to unravel the oxidative stress management in rice under moisture deficit stress and associate the morphological, biochemical and genomic information to yield-related traits using a RIL mapping population of Nagina 22 and IR 64. Nagina 22 had superior alleles for SOD (SOD2, SODCC1, SODA) and GR (GRCP2) rather than for APX, for which IR 64 had the superior allele (APX8) in terms of predicted protein stability. Gene network and SNP analysis revealed the superiority of Nagina 22 in terms of drought tolerance at reproductive stage.
Oxidative stress management
Recombinant inbred lines
Reactive oxygen species
Single plant yield
Stress Susceptibility index
Whitmarsh J, Govindjee. The Photosynthetic Process. In: Concepts in Photobiology: Photosynthesis and Photomorphogenesis, ed. GS Singhal, G Renger, SK Sopory, K-D Irrgang, Govindjee. New Delhi: Narosa, Dordrecht: Kluwer Academic; 1999. p. 11–51.
Asada K. Univalent reduction of molecular oxygen by spinach chloroplasts on illumination. J Biol Chem. 1974;249:2175–82.
Mehler AH. Studies on reactions of illuminated chloroplasts: I. Mechanism of the reduction of oxygen and other hill reagents. Arch Biochem Biophys. 1951;33:65–77.
Kaiser WM. Reversible inhibition of the calvin cycle and activation of oxidative pentose phosphate cycle in isolated intact chloroplasts by hydrogen peroxide. Planta. 1979;145:377–82.
Asada K. Production and scavenging of reactive oxygen species in chloroplasts and their functions 1. Plant Physiol. 2006;141:391–6.
Shao HB, Liang ZS, Shao MA, Sun Q. Dynamic changes of anti-oxidative enzymes of 10 wheat genotypes at soil water deficits. Colloids Surfaces B Biointerfaces. 2005;42:187–95.
Gill SS, Anjum NA, Hasanuzzaman M, Gill R, Trivedi DK, Ahmad I, et al. Glutathione and glutathione reductase: A boon in disguise for plant abiotic stress defense operations. Plant Physiol Biochem. 2013;70:204–12.
Willekens H, Chamnongpol S, Davey M, Schraudner M, Langebartels C, Van Montagu M, et al. Catalase is a sink for H 2 O 2 and is indispensable for stress defence in C 3 plants. EMBO J. 1997;16:4806–16.
Dixit S, Grondin A, Lee C-R, Henry A, Olds T-M, Kumar A. Understanding rice adaptation to varying agro-ecosystems: trait interactions and quantitative trait loci. BMC Genet. 2015;16:86.
Evans MD, Dizdaroglu M, Cooke MS. Oxidative DNA damage and disease: induction, repair and significance. Mutat Res. 2004;567:1–61.
Foyer CH. Redox homeostasis and antioxidant signaling: a metabolic interface between stress perception and physiological responses. Plant Cell Online. 2005;17:1866–75.
Møller IM, Jensen PE, Hansson A. Oxidative modifications to cellular components in plants. Annu Rev Plant Biol. 2007;58:459–81.
Sharma P, Jha AB, Dubey RS, Pessarakli M. Reactive oxygen species, oxidative damage, and antioxidative defense mechanism in plants under stressful conditions. J Bot. 2012;2012:1–26.
Hossain MA, Mostofa MG, Fujita M. Heat-shock positively modulates oxidative protection of salt and drought-stressed mustard (Brassica campestris L.) seedlings. J Plant Sci Mol Breed. 2013;2:2.
Das K, Roychoudhury A. Reactive oxygen species (ROS) and response of antioxidants as ROS-scavengers during environmental stress in plants. Front Environ Sci. 2014;2:1–13.
Hazman M, Hause B, Eiche E, Nick P, Riemann M. Increased tolerance to salt stress in OPDA-deficient rice ALLENE OXIDE CYCLASE mutants is linked to an increased ROS-scavenging activity. J Exp Bot. 2015;66:3339–52.
Hossain MA, Bhattacharjee S, Armin SM, Qian PP, Xin W, Li HY, et al. Hydrogen peroxide priming modulates abiotic oxidative stress tolerance: insights from ROS detoxification and scavenging. Front Plant Sci. 2015;6:19.
Nahar K, Hasanuzzaman M, Alam MM, Fujita M. Exogenous glutathione confers high temperature stress tolerance in mung bean (Vigna radiata L.) by modulating antioxidant defense and methylglyoxal detoxification system. Environ Exp. 2015;112:44–54.
Cruz RT, O’Toole JC. Dryland rice response to an irrigation gradient at flowering Stage1. Agron J. 1984;76:178.
Boonjung H, Fukai S. Effects of soil water deficit at different growth stages on rice growth and yield under upland conditions. 2. Phenology, biomass production and yield. F Crop Res. 1996;48:47–55.
Pantuwan G, Fukai S, Cooper M, Rajataserreekul S, O’Toole JC. Yield response of rice (Oryza sativa L.) genotypes to different types of drought under rainfed lowlands Part 1. Grain yield and yield components. F Crop Res. 2002;73:153–68.
Fischer KS, Fukai S, Kumar A, Leung H, Jongdee B. Field phenotyping strategies and breeding for adaptation of rice to drought. Front Physiol. 2012;3:1–21.
Jagadish SVK, Muthurajan R, Oane R, Wheeler TR, Heuer S, Bennett J, et al. Physiological and proteomic approaches to address heat tolerance during anthesis in rice (Oryza sativa L.). J Exp Bot. 2010;61:143–56.
Rang ZW, Jagadish SVK, Zhou QM, Craufurd PQ, Heuer S. Effect of high temperature and water stress on pollen germination and spikelet fertility in rice. Environ Exp Bot. 2011;70:58–65.
Wei S, Hu W, Deng X, Zhang Y, Liu X, Zhao X, et al. A rice calcium-dependent protein kinase OsCPK9 positively regulates drought stress tolerance and spikelet fertility. BMC Plant Biol. 2014;14:133.
Mithra SVA, Kar MK, Mohapatra T, Robin S, Sarla N, Seshashayee M, et al. DBT propelled national effort in creating mutant resource for functional genomics in rice. Curr Sci. 2016;110:543–48.
Srivalli B, Sharma G, Khanna-Chopra R. Antioxidative defense system in an upland rice cultivar subjected to increasing intensity of water stress followed by recovery. Physiol Plant. 2003;119:503–12.
Selote DS, Khanna-Chopra R. Drought-induced spikelet sterility is associated with an inefficient antioxidant defence in rice panicles. Physiol Plant. 2004;121:462–71.
Matsumoto T, et al. International rice genome sequencing project. The map-based sequence of the rice genome. Nature. 2005;436:793–800.
Zhao K, Tung C-W, Eizenga GC, Wright MH, Ali ML, Price AH, et al. Genome-wide association mapping reveals a rich genetic architecture of complex traits in Oryza sativa. Nat Commun. 2011;2:467.
Kumar V, Singh A, Mithra SVA, Krishnamurthy SL, Parida SK, Jain S, et al. Genome-wide association mapping of salinity tolerance in rice (Oryza sativa). DNA Res. 2015;22:133–45.
The 3 000 rice genomes project. The 3,000 rice genomes project. Gigascience. 2014;3:7.
Bihani P, Char B, Bhargava S. Transgenic expression of sorghum DREB2 in rice improves tolerance and yield under water limitation. J Agric Sci. 2011;149:95–101.
Hu H, Dai M, Yao J, Xiao B, Li X, Zhang Q, et al. Overexpressing a NAM, ATAF, and CUC (NAC) transcription factor enhances drought resistance and salt tolerance in rice. Proc Natl Acad Sci U S A. 2006;103:12987–92.
Karaba A, Dixit S, Greco R, Aharoni A, Trijatmiko KR, Marsch-Martinez N, et al. Improvement of water use efficiency in rice by expression of HARDY, an Arabidopsis drought and salt tolerance gene. Proc Natl Acad Sci U S A. 2007;104:15270–5.
Jeong JS, Kim YS, Baek KH, Jung H, Ha S-H, Do Choi Y, et al. Root-specific expression of OsNAC10 improves drought tolerance and grain yield in rice under field drought conditions. Plant Physiol. 2010;153:185–97.
Dixit S, Singh A, Sta Cruz M, Maturan PT, Amante M, Kumar A. Multiple major QTL lead to stable yield performance of rice cultivars across varying drought intensities. BMC Genet. 2014;15:16.
Bimpong IK, Manneh B, Sock M, Diaw F, Amoah NKA, Ismail AM, et al. Improving salt tolerance of lowland rice cultivar “Rassi” through marker-aided backcross breeding in West Africa. Plant Sci. 2016;242:288–99.
Singh R, Singh Y, Xalaxo S, Verulkar S, Yadav N, Singh S, et al. From QTL to variety-harnessing the benefits of QTLs for drought, flood and salt tolerance in mega rice varieties of India through a multi-institutional network. Plant Sci. 2016;242:278–87.
Shanmugavadivel PS, Mithra SVA, Dokku P, Kumar KAR, Rao GJN, Singh VP, et al. Mapping quantitative trait loci (QTL) for grain size in rice using a RIL population from Basmati × indica cross showing high segregation distortion. Euphytica. 2013;194:401–16.
Wang C, Chen S, Yu S. Functional markers developed from multiple loci in GS3 for fine marker-assisted selection of grain length in rice. Theor Appl Genet. 2011;122:905–13.
Wang S, Li S, Liu Q, Wu K, Zhang J, Wang S, et al. The OsSPL16-GW7 regulatory module determines grain shape and simultaneously improves rice yield and grain quality. Nat Genet. 2015;47:949–54.
Tripathy JN, Zhang J, Robin S, Nguyen TT, Nguyen HT. QTLs for cell-membrane stability mapped in rice (Oryza sativa L.) under drought stress. Theor Appl Genet. 2000;100:1197–202.
Babu RC, Nguyen BD, Chamarerk V, Shanmugasundaram P, Chezhian P, Jeyaprakash P, et al. Genetic analysis of drought resistance in rice by molecular markers. Crop Sci. 2003;43:1457.
Robin S, Pathan MS, Courtois B, Lafitte R, Carandang S, Lanceras S, et al. Mapping osmotic adjustment in an advanced back-cross inbred population of rice. Theor Appl Genet. 2003;107:1288–96.
Zhang Y, Verhoeff NI, Chen Z, Chen S, Wang M, Zhu Z, et al. Functions of OsDof25 in regulation of OsC4PPDK. Plant Mol Biol. 2015;89:229–42.
Xu Y, This D, Pausch RC, Vonhof WM, Coburn JR, Comstock JP, et al. Leaf-level water use efficiency determined by carbon isotope discrimination in rice seedlings: Genetic variation associated with population structure and QTL mapping. Theor Appl Genet. 2009;118:1065–81.
Septiningsih EM, Pamplona AM, Sanchez DL, Neeraja CN, Vergara GV, Heuer S, et al. Development of submergence-tolerant rice cultivars: The Sub1 locus and beyond. Ann Bot. 2009;103:151–60.
Fischer RA, Maurer R. Drought resistance in spring wheat cultivars, 1. Grain yield responses. Aust J Agric Res. 1978;29:897–912.
Dhindsa RS, Plumb-dhindsa P, Thorpe TA. Leaf senescence: Correlated with increased levels of membrane permeability and lipid peroxidation, and decreased levels of superoxide dismutase and catalase. J Exp Bot. 1981;32:93–101.
Nakano Y, Asada K. Hydrogen peroxide is scavenged by ascorbate-specific peroxidase in spinach chloroplasts. Plant Cell Physiol. 1981;22:867–80.
Smith IK, Vierheller TL, Thorne CA. Assay of glutathione reductase in crude tissue homogenates using 5,5-dithiobis(2-nitrobenzoic acid). Anal Biochem. 1988;175:408–13.
Aebi H. Catalase in vitro. In: Packer L, editor. Methods in Enzymology 105. Orlando: Academic; 1984. p. 121–6.
Gour P, Garg P, Jain R, Joseph SV, Tyagi AK, Raghuvanshi S. Manually curated database of rice proteins. Nucleic Acids Res. 2014;42(1):D1214-D1221.
Hall T. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp Ser. 1999;41:95–8.
Capriotti E, Fariselli P, Casadio R. I-Mutant2.0: Predicting stability changes upon mutation from the protein sequence or structure. Nucleic Acids Res. 2005;33:W306-W310. doi:10.1093/nar/gki375.
Lee T, Oh T, Yang S, Shin J, Hwang S, Kim CY, et al. RiceNet v2: an improved network prioritization server for rice genes. Nucleic Acids Res. 2015;43:W122–7.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13:2498–504.
Team RC. R: A Language and Environment for Statistical Computing. 2015.
Lê S, Josse J, Husson F. FactoMineR: An R package for multivariate analysis. J Stat Softw. 2008;25:1–18.
Vincent Q. ggbiplot: A ggplot2 based biplot. 2011.
Taiyun W. corrplot: Visualization of a correlation matrix. 2013.
Harrell Jr FE, with contributions from Charles Dupont and many others. Hmisc: Harrell Miscellaneous.2015. R package version 3.17-4. https://CRAN.R-project.org/package=Hmisc.
Peterson BG, Peter C. PerformanceAnalytics: Econometric tools for performance and risk analysis. 2014.
Siangliw JL, Jongdee B, Pantuwan G, Toojinda T. Developing KDML105 backcross introgression lines using marker-assisted selection for QTLs associated with drought tolerance in rice. ScienceAsia. 2007;33:207–14.
Noraziyah AAS, Mallikarjuna Swamy BP, Ratnam WM, Raman A, Kumar A. Marker assisted pyramiding of drought yield QTLs into a popular Malaysian rice cultivar, MR 219. BMC Genet. 2016. doi:10.1186/s12863-016-0334-0.
Hu H, Xiong L. Genetic engineering and breeding of drought-resistant crops. Annu Rev Plant Biol. 2014;65:715–41.
Farooq M, Wahid A, Kobayashi N, Fujita D, Basra SMA. Plant drought stress: effects, mechanisms and management. Agron Sustain Dev. 2009;29:185–212.
Chugh V, Kaur N, Grewal MS, Gupta AK. Differential antioxidative response of tolerant and sensitive maize (Zea mays L.) genotypes to drought stress at reproductive stage. Indian J Biochem Biophys. 2013;50:150–8.
Michelmore RW, Paran I, Kesseli RV. Identification of markers linked to disease-resistance genes by bulked segregant analysis: a rapid method to detect markers in specific genomic regions by using segregating populations. Proc Natl Acad Sci U S A. 1991;88:9828–32.
Tiwari S, Krishnamurthy SL, Kumar V, Singh B, Rao A, Mithra SVA, et al. Mapping QTLs for salt tolerance in rice (Oryza sativa L.) by bulked segregant analysis of recombinant inbred lines using 50K SNP chip. Yadav RS, editor. PLoS One. 2016;11:e0153610.
Raman A, Verulkar S, Mandal N, Variar M, Shukla V, Dwivedi J, et al. Drought yield index to select high yielding rice lines under different drought stress severities. Rice. 2012;5:31.
Smirnoff N. The role of active oxygen in the response of plants to water deficit and desiccation. New Phytol. 1993;125:27–58.
Tsukamoto S, Morita S, Hirano E, Yokoi H, Masumura T, Tanaka K. A novel cis-element that is responsive to oxidative stress regulates three antioxidant defense genes in Rice1. Plant Physiol. 2005;137:317–27.
Xu FQ, Li XR, Ruan YL. RNAi-mediated suppression of hexokinase gene OsHXK10 in rice leads to non-dehiscent anther and reduction of pollen germination. Plant Sci. 2008;175:674–84.
The first author is grateful to Council of Scientific and Industrial Research, Government of India for the award of a senior research fellowship, ICAR-IARI and ICAR-NRCPB, New Delhi for providing fund and facilities. We sincerely thank Dr. P.K. Mandal, Principal Scientist, ICAR-NRCPB for advice and suggestions in the enzyme assay. Dr. Viswanathan Satheesh’s contribution in guiding the analyses using R package is also acknowledged.
Availability of supporting data
Only the raw data of the enzyme assay and phenotyping are not submitted. The latter one is too huge. This will be made available upon request.
CP carried out the field phenotyping, multivariate and bioinformatics analysis. CP and PKS carried out enzyme assay work. SVAM conceived the experiment, compiled and interpreted the results and coordinated the entire work. SVAM and CP drafted the manuscript. TM developed the mapping population and guided phenotyping under drought. NKS edited and approved the manuscript and supervised the work. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
NA; Sequence data already available in the public domain have been made use of.
Correlogram of the morpho-physio traits of the entire RIL population (Kharif 2014) under a Irrigated control (IC), SPY is positively associated with SPAD and SPP and b Reproductive stage drought stress (RS), SPY is negatively interacting with SS and SSI while positively with SPAD. (TIF 3276 kb)
Linear correlation r (Pearson) among morpho- physiological and enzyme results in parents and bulks under a IC: Irrigated control b RS: Reproductive stage drought stress (2014 Kharif data). (DOCX 17 kb)
Evidence for interaction network (.xlx). The file contains all the necessary scores and loci interaction details generated using RiceNet v2. (XLSX 14 kb)
About this article
Cite this article
Prakash, C., Mithra, S.V.A., Singh, P.K. et al. Unraveling the molecular basis of oxidative stress management in a drought tolerant rice genotype Nagina 22. BMC Genomics 17, 774 (2016). https://doi.org/10.1186/s12864-016-3131-2
- Drought tolerance
- Functional polymorphism
- Oxidative stress management
- Nagina 22
- Spikelet fertility