ORIGINAL RESEARCH article
Microarray Meta-Analysis Focused on the Response of Genes Involved in Redox Homeostasis to Diverse Abiotic Stresses in Rice
- Abiotic Stress Tolerance, Department of Plant Nutrition, Instituts für Nutzpflanzenwissenschaften und Ressourcenschutz, Rheinische Friedrich-Wilhelms-Universität Bonn, Bonn, Germany
Plants are exposed to a wide range of abiotic stresses (AS), which often occur in combination. Because physiological investigations typically focus on one stress, our understanding of unspecific stress responses remains limited. The plant redox homeostasis, i.e., the production and removal of reactive oxygen species (ROS), may be involved in many environmental stress conditions. Therefore, this study intended to identify genes, which are activated in diverse AS, focusing on ROS-related pathways. We conducted a meta-analysis (MA) of microarray experiments, focusing on rice. Transcriptome data were mined from public databases and fellow researchers, which represented 36 different experiments and investigated diverse AS, including ozone stress, drought, heat, cold, salinity, and mineral deficiencies/toxicities. To overcome the inherent artifacts of different MA methods, data were processed using Fisher, rOP, REM, and product of rank (GeneSelector), and genes identified by most approaches were considered as shared differentially expressed genes (DEGs). Two MA strategies were adopted: first, datasets were separated into shoot, root, and seedling experiments, and these tissues were analyzed separately to identify shared DEGs. Second, shoot and seedling experiments were classed into oxidative stress (OS), i.e., ozone and hydrogen peroxide treatments directly producing ROS in plant tissue, and other AS, in which ROS production is indirect. In all tissues and stress conditions, genes a priori considered as ROS-related were overrepresented among the DEGs, as they represented 4% of all expressed genes but 7–10% of the DEGs. The combined MA approach was substantially more conservative than individual MA methods and identified 1001 shared DEGs in shoots, 837 shared DEGs in root, and 1172 shared DEGs in seedlings. Within the OS and AS groups, 990 and 1727 shared DEGs were identified, respectively. In total, 311 genes were shared between OS and AS, including many regulatory genes. Combined co-expression analysis identified among those a cluster of 42 genes, many involved in the photosynthetic apparatus and responsive to drought, iron deficiency, arsenic toxicity, and ozone. Our data demonstrate the importance of redox homeostasis in plant stress responses and the power of MA to identify candidate genes underlying unspecific signaling pathways.
While most animals can move and escape from harmful conditions, plants cannot. Drought, flood, salinity, extremes of temperature, nutrient deficiency, UV radiation, pollutants, herbicides, and pathogens are some of the factors a plant needs to cope with to survive and grow, usually with more than one at time. These factors limit crop yields and quality (Gill and Tuteja, 2010; Miller et al., 2010; Wang and Frei, 2011). Understanding, how a plant responds to such stresses and the mechanisms underlying stress tolerance can give us a better view of how to improve the global food production. Investigating the stress responses of rice (Oryza sativa L.) is particularly rewarding, because it is both a global staple food of great agronomical importance, and a well-studied model organism, for which many transcriptome profiling studies have been published. About 32% of the rice annual production fluctuations (corresponding to around three million tons) can be attributed abiotic stresses (AS), and specifically variations in climate including precipitation and temperature (Ray et al., 2015).
A common responses to different environmental stresses, both abiotic and biotic, is the excessive generation of reactive oxygen species (ROS) including superoxide (), perhydroxy radical (), hydrogen peroxide (H2O2), hydroxyl radical (OH⋅), singlet oxygen (1O2), and organic hydroperoxide (ROOH; Gill and Tuteja, 2010; Bhattacharjee, 2012). Plant cells continuously produce ROS as by-products of various metabolic processes including respiration and photosynthesis. However, these molecules can directly react with DNA, proteins and lipids causing severe damage to individual cells and whole organisms. For this reason, every aerobic organism tightly controls its ROS concentration by ROS-scavenging pathways composed of many enzymatic and non-enzymatic antioxidant components (Gill and Tuteja, 2010).
However, many studies have shown that ROS play important roles in plants’ stress signaling processes (Fujita et al., 2006; Miller et al., 2010; Mittler et al., 2011; Kim et al., 2012; Pucciariello et al., 2012). There are many advantages in the use of ROS as signaling molecules. The cell can rapidly produce and scavenge different forms of ROS in a simultaneous manner, enabling rapid and dynamic changes in ROS levels (caused by simply tilting the balance between cellular production and scavenging rates; Bhattacharjee, 2012). Each form of ROS has its own characteristics, such as mobility, process and location of origin and its reactivity with DNA, proteins or lipids (Bhattacharjee, 2012). Another advantage is the tight control over the subcellular localization of ROS signals in cells, with the regulation of enzymes specific of each of the cell compartments (Mittler et al., 2011). ROS may even act as “secondary messengers” modulating the activities of specific proteins or expression of genes by changing the redox balance of the cell. The network of redox signals orchestrates the metabolism for regulating energy production to utilization, interfering with primary signaling agents (hormones) to respond to changing environmental conditions at every stage of plant development (Bhattacharjee, 2012). Oxidative signaling is now considered to be a key in the responses to stress, involved not only in the defense to these stresses but also in the regulation of the plant growth and development (Noctor et al., 2014).
One limitation in the understanding of ROS homeostasis in AS is that most of the knowledge on molecular mechanisms of stress response was obtained from experiments under controlled laboratory conditions and focused on only one stress at a time. However, plants are often simultaneously exposed to many biotic and AS in their natural or agronomic habitats (Hazen et al., 2003; Rasmussen et al., 2013). By combining the data from different experiments it is possible to identify common and specific elements expressed in response to different stresses (Rabbani et al., 2003; Fujita et al., 2006; Yang et al., 2013).
With the rapid advances in biological high-throughput technology, a large and diverse set of genomic data has become publicly accessible. Combining information from multiple existing studies can increase the reliability and generalizability of results. The use of statistical techniques to combine results from independent but related studies is called “meta-analysis (MA).” Through MA, we can increase the statistical power to obtain a more precise estimate of gene expression differentials (Ramasamy et al., 2008; Tseng et al., 2012). The most common types of Microarray MA are a combination of p-values, a combination of effect sizes (fold change) or a combination of ranks, each method with its limitations and advantages (Tseng et al., 2012).
In the present study, instead of choosing one of those, we used an integrative approach to combine the forces of those different microarray MA methods and overcome possible biases, such as a “fishing for significance” effect for preferring one method over the others (Ioannidis, 2005; Ostlund and Sonnhammer, 2014). However, within these three categories of MA, there are many statistical methods to choose from, we selected them based on their evaluation in recent studies (Tseng et al., 2012; Chang et al., 2013; Ostlund and Sonnhammer, 2014). Among the methods of combination of p-values, we chose two distinct methods: Fisher’s (Rhodes et al., 2002), the most commonly implemented MA method, and rOP (rth ordered p-value) that is more restrictive but still flexible (Song and Tseng, 2014). Between the combined effect size methods, the random effect model, REM is the most adequate for a heterogeneous group of samples, as the one processed in the present study (Choi et al., 2003). Furthermore, the MA by ranks was computed by the program GeneSelector that produces a rank combining seven distinct statistic methods (Boulesteix and Slawski, 2009; Ostlund and Sonnhammer, 2014). Only genes elected by the majority of methods (at least three of those four) were further processed in our analyses. Although this strategy may be rather restrictive, it increases the power of our analysis, since each differentially expressed gene (DEG) was confirmed by at least three different statistic methods.
The goal of this MA was to detect DEGs involved in ROS homeostasis, which respond to AS treatments. The transcriptome data from rice plants subjected to different classes of AS were analyzed: drought, submergence, salinity, cold, heat, excess and/or deficiency of essential nutrients, such as phosphorus (P), zinc (Zn), and iron (Fe), and toxicity of heavy metals such as arsenic (As), cadmium (Cd), chrome (Cr), and lead (Pb). Also, to identify genes specifically involved in the redox homeostasis, data from experiments with ozone (O3) and hydrogen peroxide were included (Table 1), because these treatments are presumed to produce direct oxidative stimulus, unlike with the indirect forms of oxidative stress (OS) occurring in other environmental conditions.
More specifically, the following questions were addressed in this study:
(1) How does a MA integrating several of the MA approaches described above compare to individual MA approaches, and are the results more conservative?
(2) To what extent do stress responses differ between different stresses and plant tissues, and what is the role of ROS-related genes?
(3) Through MA, can we nominate possible key genes as major hubs of ROS-related stress signaling, and what are their putative functions?
Materials and Methods
For this study, expression data of rice plants exposed to diverse AS were combined. The raw expression data of different experiments were obtained from the Rice Oligonucleotide Array Database1 (Cao et al., 2012), the NCBI’s Gene Expression Omnibus repository2 (Barrett et al., 2013), and from fellow researches. The criteria for inclusion of those dataset were: relatively similar genetic background, plants must be from O. sativa indica or japonica subgroups; the RNA should have been extracted only from shoots, roots, or seedlings (whole plant), excluding other tissues such as flowers, seeds, or cali; the experiments must involve only against AS treatments; the data must originate from Affymetrix or Agilent microarray platforms, and the original study must follow the “Minimal Information About a Microarray Experiment” (MIAME) requirements (Brazma et al., 2001).
Individual Datasets Analysis
Microarray expression data from each source study was pre-processed separately as individual datasets. Agilent microarray data were processed in the R program using the packages LIMMA (Smyth, 2005), while the package Affy was used for data from the Affymetrix platforms. The raw data from both were treated using Robust Multi-array Average (RMA) background correction and quantile normalization. Non-informative and low-intensity probes were declared following the program standard settings, while duplicated probes were converted into their corresponding genomic locus. The ArrayQualityMetrics package (Kauffmann and Huber, 2013) was used to assess the quality of the normalized datasets. A sample was included or excluded in further analysis based on three different evaluations made by this program: distance between array, array intensity distributions, and individual array quality. The values normalized by RMA were used for the subsequent MA.
A contrast between each treatment and its control was estimated with the LIMMA package. In studies including several genotypes the genotype factor was not considered. After fitting the data into a linear model, the standard errors were corrected using a simple empirical Bayes model. Moderated t-statistic and log-odds of differential expression were computed for each contrast for each gene. Genes that showed significant P-value (FDR = 5%) were considered as DEG and log2-fold-change values of each experiment/dataset were saved for further analysis. Relative gene expression values corresponding to the same stress in different datasets were averaged for a simplified evaluation of the gene’s response to each condition.
Two different strategies were implemented in the present study. The first (MA.1) investigated the effect of the stresses in different tissues, for which the datasets were separated into shoots, roots and seedling microarrays and processed separately. In the second approach (MA.2), data from shoots and seedlings were combined, while the data set was separated into direct OS and others AS to compare the effects of direct and indirect OS on gene expression. Hydrogen peroxide and ozone stress were considered as OS, because these treatments directly lead to the production of ROS in plant tissue (Uchida et al., 2002; Kangasjärvi et al., 2005). For both approaches, the normalized expression values were used. Since the probe nomenclature differs between platforms, the MSU-ID was used, and when multiple probes matched to the same gene they were averaged.
The datasets were merged using the packages metaDE (Wang et al., 2012), and this merged dataset was once more filtered, excluding 20% of un-expressed genes (with small expression intensities) and 20% of non-informative genes (genes with small variation). Each independent study sample was processed with a modified two-sample t-statistics (modt) contrasting treated and control samples.
In an effort to overcome the inherent artifacts of each MA statistical method, we ran our data through the three common types of MA (i.e., by p-value, by effect size and by rank) to identify genes that are considered differentially expressed by distinct methods. Two methods of combination by p-value were used: Fisher and rOP. The classical Fisher’s method sums the log-transformed p-values obtained from individual studies and, under null hypothesis, follows a chi-squared distribution with 2K degrees of freedom, where K is the number of studies combined (Rhodes et al., 2002). However an extremely small p-value in only one study can be sufficient to cause statistical significance, even if the same gene are not significant in any other study. A more restrictive but flexible method is rOP, that combines Fisher with a generalized vote counting statistic. It uses the rth order statistic among sorted p-values of K combined studies, where r is a pre-determined minimum number of studies, in which a gene’s p-values must be small to be significant (Song and Tseng, 2014). In our analysis we implemented K/2 ≤r ≤ K, i.e., each gene must be differentially expressed in at least half of the combined studies to be significant. Taking into consideration the heterogeneity of our cluster of studies and the residual “noise” data derived from technical and biological differences between the studies, a more restrictive approach was not implemented (i.e., r = K).
A second way to combine expression data across different microarray studies and platforms is using effect size values. The REM method was implemented because it possesses a random effect element corresponding to the unknown heterogeneities between very distinct studies, such element is not present in the alternative method, the fixed effect model, FEM (Choi et al., 2003). The results of each of these analyses were corrected with the Benjamini and Hochberg procedure with a false discovery rate (FDR) of 5% as threshold.
For a MA based on the rank method, the normalized values were processed using GeneSelector (Boulesteix and Slawski, 2009). This package allows a ranking analysis of the data with seven distinct methods: ordinary t-test; Baldi and Long Bayesian t-test (Baldi and Long, 2001); Winconxon–Mann–Whitney U test; Fox and Dimmic Bayesian t-test (Fox and Dimmic, 2006); SAM statistics (Tusher et al., 2001); limma: moderated t-statistic based on a Bayesian hierarchical model which is estimated by an empirical Bayes approach (Smyth et al., 2003), and simple fold-change estimation (in log2). The obtained gene rankings were then aggregated by the mean value of the rank positions given by each method (AggMean), or aggregated on the basis of a Markov Chain model, AggMC (DeConde et al., 2006). A combined rank was produced by combining the top genes obtained by each of the seven ranking methods, together with the AggMean and AggMC lists. Genes of this combined list were compared with DEGs obtained by Fisher, rOP, and REM methods.
The DEGs shared by at least three of these four methods were further studied.
The genes elected by the combined meta-analyses MA.1 and MA.2, were further analyzed and characterized. First, the result lists were combined and compared with a list of genes described in the literature as involved in ROS scavenging and signaling processes. This list was curated based on recent studies (Frei et al., 2010a; Liu et al., 2010; Kim et al., 2012; Shaik et al., 2014). This list includes genes involved in the biosynthesis and recycling of enzymatic [such as superoxide dismutase (SOD), catalase (CAT), ascorbate peroxidase (APX), glutathione peroxidase (GPX), and glutathione reductase (GR)) and non-enzymatic antioxidants (such as ascorbic acid (AsA), reduced glutathione (GSH) and thioredoxin (Trx)], and also transcription factors such as zinc-fingers, MYB (myeloblastosis) and WRKY (that contains the WRKYGQK amino acid conserved sequence) families and other elements described as directly or indirectly in involved in signaling and response to OS. Co-occurrence of DEG with the resulting list was represented with Venn’s diagrams using Venny 2.0 (Oliveiros, 2015).
Gene Ontology (GO) Enrichment analysis of the DEGs obtained by MA.1 and MA.2 was conducted using the AgriGO platform (Du et al., 2010). A Singular Enrichment Analysis was performed using the Rice Gramene Locus set as reference (Jaiswal et al., 2002).
Genes differentially expressed in both OS and AS (MA.2 approach) were studied in depth. The STRING (Search Tool for the Retrieval of Inter-acting Genes/Proteins) database was used to detect functional association between those genes. This database constructs associations based on distinct lines of evidences: Experimental evidence from protein–protein interaction assays; co-expression data based on the expression data stored in the NCBI GEO database; co-occurrence of the genes in the same organisms; available information of other databases; recurring neighborhood of the genes in known genomes; events of fusion between those genes or orthologs; pathway annotation data in other databases such as Gene Ontology or Kyoto Encyclopedia of Genes and Genomes (KEGG), and automated text-mining based on Medline abstracts and a large collection of full-text articles (Szklarczyk et al., 2015). STRING computes a confidence score for those interactions based on the available evidences, from medium (score above 0.4) to highest (above 0.9).
Using the expression values obtained in the individual dataset analysis (a log2 fold change difference between control and treated expression values of each gene), an average value of the most relevant conditions was obtained and represented in a heatmap. The heatmap was made using the package gplots in R (Warnes et al., 2015).
The Percentage of ROS Related DEGs is Constant in the Response to Different Stresses
After search in the available databases and quality control analysis, raw microarray data from 36 independent experiments were selected. To obtain a global analysis, AS from different categories were chosen: drought, submergence, salinity, cold, heat, excess and/or deficiency of essential nutrients as Phosphorus (P), zinc (Zn), and iron (Fe), heavy metal toxicity (As, Cd, Cr, and Pb), and direct OS (O3 and H2O2; Table 1).
Some source data originated from experiments that tested different stressors (Jain et al., 2007; Zheng et al., 2009; Dubey et al., 2010; Mittal et al., 2012a,b); several time points after the exposure to the stressor (Cho et al., 2008; Mustroph et al., 2010; Dai et al., 2012; Ma et al., 2012; Mittal et al., 2012a); several concentrations of a stressor (Takehisa et al., 2015), and/or several genotypes, usually using contrasting lines (Pariasca-Tanaka et al., 2009; Frei et al., 2010a; Cotsaftis et al., 2011; Pandit et al., 2011; Norton et al., 2014; Supplementary Table S1). In these cases, every possible contrast of stress condition versus control was treated as a separate experiment. In total, plants from 21 different genotypes of domestic rice were used, 13 from the indica and eight from the japonica subspecies (Supplementary Table S1).
In parallel, a list of 1972 genes previously described as involved in ROS scavenging and signaling processes was made by combining information from the literature (Frei et al., 2010a; Liu et al., 2010; Kim et al., 2012; Shaik et al., 2014). These genes correspond to about 4% of the rice genes represented by the microarray probes included in this analysis (Supplementary Table S2) and represent a broad spectrum of functions, ranging from reductase and peroxidase enzymes to transcription factors.
In the different treatments and tissues, 7–10% of the DEGs were included in the list of ROS-related genes, indicating that this category of genes was overrepresented (Table 2). The proportion of DEG considered as ROS-related was similar in OS and AS experiments (Table 3).
TABLE 2. Differentially expressed genes (DEGs) in response to different abiotic stresses and the proportion of those involved in ROS scavenging or signaling mechanisms (ROS).
MA.1: DEGs in Response to Diverse Abiotic Stresses Differ Between Tissues
The microarray expression data was processed by two different MA approaches. In the first, MA.1, our objective was to identify DEGs in response to many AS in shoot, root, and seedling samples separately. Using combined MA methods, 1001 DEGs were identified in response to different stresses in shoots, 837 in roots and 1172 in seedlings, although only 14 were identified in all three tissues (Figure 1B, Table 3). From the 2691 genes elected by MA.1, 236 (9%) were included in the ROS-related list, including the coding genes of 42 zinc-finger signaling proteins, 18 MYB transcript factors, nine glutathione S-transferases, a Cu-Zn SOD, a copper chaperone for SOD, 19 peroxidase precursors, APX2, GPX2, and GPX3 (Figure 1B, Supplementary Table S5).
FIGURE 1. Venn Diagrams showing co-occurrence of differentially expressed genes (DEGs) in response to many abiotic stresses in: (A) Shoots samples processed by four different meta-analytical (MA) statistical methods: Fisher, rOP, REM, and GeneSelector. (B) MA.1 combined analysis of shoots, seedlings, and roots samples, compared with an a priori ROS-related gene list. (C) Oxidative stress (MA.2 OS) and other abiotic stresses (MA.2 AS), compared with MA.1 and ROS-related genes.
MA.2: DEGs Detected in Both Oxidative and Others Abiotic Stresses
In a second approach, expression data from seedling and shoot experiments were classed into two groups: OS (direct OS, i.e., O3 and H2O2) and AS (all the other AS). It was found that 679 DEGs were exclusive to OS and 1416 DEGs were exclusive to AS (Table 3, Figure 1C). In total 1212 genes were detected by both MA.2 and MA.1, of which 197 coincided with OS, 760 with AS, and 255 with both (Figure 1C). On the other hand, 50 ROS-related genes were only detected by the OS analysis, including APX7 and nine WRKY proteins, and 47 by the AS analysis (Supplementary Table S5).
We then focused on the 311 DEGs shared by both OS and AS. Among those, 33 were ROS-related genes, while others belonged to many distinct gene families associated with the response to stress and plant growth, such as Zn fingers, WRKY, and TIFY transcript factors, cytochromes, photosystem subunits, heat shock proteins, HIPPs (Heavy Metal Isoprenylated Plant Proteins), kinases and phosphatases (Supplementary Table S5).
A network analysis using the platform STRING 10.0 detected a cluster between 214 of those genes (Figure 2A). By increasing the stringency of this analysis (confidence score > 0.9), it was possible to isolate a cluster of 47 protein-coding genes (Figure 2B). Most hub genes of this network code proteins involved in the photosynthetic apparatus, e.g., ATP synthase (LOC_Os02g51470), oxygen evolving enhancer protein 3 (LOC_Os07g36080), photosystem I reaction center subunit III (LOC_Os03g56670), photosystem II core complex proteins psbY (LOC_Os08g02630), photosystem II reaction center W protein (LOC_Os05g43310) and many others chloroplastic protein (Supplementary Table S5). Many of these genes showed a similar expression pattern in the samples analyzed in this study (Supplementary Table S7). They were more highly expressed in response to iron deficiency and drought in roots, while suppressed in the samples of submitted to As toxic level, ozone, and submergence (Figure 3).
FIGURE 2. Network analysis of the DEGs identified both in oxidative (OS) and abiotic stress (AS) analyses of MA.2 approach. (A) Networks formed by STRING 10.0 database shows all connections of those genes with a confidence score > 0.4. (B) Even with a more conservative analysis (confidence score > 0.9), 47 of those genes formed a cluster (genes listed in Figure 3 and Supplementary Table S5). The connection colors represent the types of evidence for inferring association: recurring Neighborhood in different genomes (green line), events of Gene Fusion (red), Co-occurrence of those genes in the same organisms (dark blue), co-expression (black), Experimental protein–protein interaction data (pink), pathway described by other databases (light blue), literature text-mining (yellow), and homology (purple lines). Source: STRING 10.0 (Szklarczyk et al., 2015).
FIGURE 3. Relative expression of the genes shown in Figure 2B in response to different treatments and in different tissues (R, Roots; S, Shoots; and P, Seedlings) followed by each gene’s MSU locus, STRING node id and description. Each gene is represented in Figure 2B as a node identified by a STRING node number or the locus name (loc). Color scale shows relative expression of each gene in relation to all tested treatments.
Shared DEGs are Involved Many Metabolic, Response to Stimuli, and Regulatory Processes
Gene Ontology enrichment analysis was conducted to explore other possible functions of the DEGs detected in the different MA approaches (Supplementary Table S6). The most frequent and significant GO terms associated with DEGs in MA.1, MA.2, AS, and OS are represented in Figure 4. The biological processes terms suggest constitutive roles for those genes, as part of metabolic and biosynthetic processes, but also in the regulation of those processes on different levels (transcription, post-translational protein modification, macromolecule biosynthesis, phosphorylation, signal transduction, transport, and proteolysis). The enrichment of terms such as nucleic acid binding, transcription regulator, kinase activity, transmembrane transporter, and phosphatase activity indicates that many of those shared genes are also involved in signaling processes (Figure 4).
FIGURE 4. Gene Ontology (GO) – Frequency of the most representative Biological Process terms. In blue DEGs obtained by the first meta-analysis approach (MA.1), in red by oxidative stress (MA.2 OS) and in green by abiotic stress (MA.2 AS) of the second approach. The frequency of these terms in the reference background (BG), the Rice Gramene Locus set is shown in black. GO analysis made in the AgriGO platform (FDR = 5%).
Advantages of a Combined MA Approach
With the objective of identifying genes involved in the plant response to different AS and direct or indirect OS, the expression data of rice plants exposed to many stress conditions were combined. Instead of using one specific MA approach, we combined the results of different MA statistical methods to obtain the most relevant genes. The stringency of this methodology may exclude many important genes, but also excludes many false positives that could be obtained as artifacts of each statistical method (Tseng et al., 2012; Chang et al., 2013).
While methods such as Fisher’s identify a great number of genes, many of those were not identified as differentially expressed by size effect or ranking product methods (e.g., Figure 1A). Only a small fraction of the genes differentially expressed was shared between the independent methods. For example, while with the Fisher’s MA method the number of DEGs was between 40 and 98% of the number of DEGs detected in at least one experiment, our shared DEGs approach reached 13–25% (Table 3B).
MA.1: The Response to Stresses Varies Between Tissues
In our first MA approach, microarray data were separated into shoots, roots, and seedlings, and although those groups were composed of samples of similar size (Tables 1 and 2), they showed greatly distinct numbers of DEGs in response to stresses, and only 14 of those were shared between the three groups (Figure 2B). This is in agreement with the concept that the transcriptome and how it responds to a stress differs between tissues or organs of the same organism (Hazen et al., 2003).
ROS-Related DEGs are Overrepresented Among the DEGs Detected by MA.1 and MA.2
Although the DEGs in response to different stresses varied greatly, the proportion of those that were a priori considered as ROS-related was rather constant (Table 2). And while these ROS-related genes correspond to only 4% of the rice expressed genes, this category accounted for up to 10% the shared DEGs, highlighting the importance of ROS scavenging and signaling to the response to stress (Table 3C). Interestingly, the proportion was similar in OS and AS experiments. ROS imbalance or oxidative bursts that characterize OS are often indirect consequences of another environmental stresses (Guo et al., 2006). In the MA.2 approach, direct OS were separated and compared with other AS, such as salinity, drought, P deficiency, or Fe toxicity (AS). Ozone and hydrogen peroxide were used to generate the OS in the OS experiments. While hydrogen peroxide is a normal product of the plant metabolism, tropospheric ozone mostly originates from anthropogenic gas emissions (Baier et al., 2005; Wang and Frei, 2011). Most of the damage caused by this air pollutant is caused by its immediate degradation into ROS in the apoplastic space of plant cells, including hydrogen peroxide. Directly and indirectly O3 induces a cascade of active ROS production and signaling (Vaultier and Jolivet, 2015). Both OS treatments thus induce direct OS, though from opposite starting sites. While the O3 treatment involved fumigating of leaves (Frei et al., 2010a; Cho et al., 2013), the H2O2 treatment was performed by growing seedlings in a solution containing 10 mM H2O2, starting the oxidative reactions from the plant roots (Mittal et al., 2012a,b).
Despite the theoretical differences between OS and AS experiments, the percentage of DEGs considered to be ROS-related was similar (Tables 2 and 3). This data enforces the concept that most AS generate OS or ROS signaling to some degree and highlights the importance of ROS homeostasis in AS response. Diverse studies claimed associations of OS with the stresses included in this MA: drought (Noctor et al., 2014); salinity (Miller et al., 2010; Chawla et al., 2012); temperature extremes (Guo et al., 2006; Mittal et al., 2012a); zinc deficiency (Frei et al., 2010b; Höller et al., 2014); phosphorous deficiency (Hernandez and Munne-Bosch, 2015); potassium deficiency (Cakmak, 2005; Ma et al., 2012), and iron deficiency (Zheng et al., 2009); iron toxicity (Matthus et al., 2015); cadmium toxicity (Uraguchi et al., 2011; Ogo et al., 2014); arsenic toxicity (Azizur Rahman et al., 2007); and lead toxicity (Li et al., 2012). In agreement with these studies, our results suggest that ROS homeostasis plays a similarly important role in all of these stresses.
Putative Functions of Multi-Stress-Responsive Genes Detected by Combined MA Approaches
Differentially expressed genes detected by both approaches (i.e., MA.1 and MA.2) shared many GO terms, and even if the individual DEGs detected by each approach diverged (Figures 1 and 4), their functions were rather conserved. The DEGs identified by MA.1 and MA.2 were mostly related to metabolic and cellular processes, but also involved in the response to stimuli, regulation, transcription, and transport processes. Many of those genes can be associated with signaling pathways possessing functions such as catalytic, transferase, hydrolase, transport activity, DNA binding, and transcript regulator activity (Figure 3). Interestingly, the terms signal transduction process, metal ion binding activity, transcription factor activity and pyrophosphatase activity were found only among OS DEGs (Figure 4).
As the GO enrichment analysis demonstrated, the DEGs elected by these MA approaches represented different classes of proteins and are involved in many distinct processes. In MA.2, 311 genes were identified in both OS and AS. Between those genes, a great number of distinct transporters were present, such as the calcium transporter ATPases LOC_Os12g39660 and LOC_Os05g02940, the citrate transporter protein LOC_Os02g57620, the inorganic phosphate transporter LOC_Os02g38020, the sulfate transporter LOC_Os03g09930, the amino acid transporter LOC_Os06g36180, the aquaporin LOC_Os02g41860, and others (Supplementary Table S5). Among those were also four HIPPs, proteins that can act as cytoplasmic transporters of metallic ions and co-factor in the transcription of many stress related genes (de Abreu-Neto et al., 2013).
Using the network tool STRING 10.0, we identified a cluster of 36 DEGs that are mostly involved in the photosynthetic apparatus and its regulation (Figure 2). Photosynthesis is a main source of ROS in plants, which can be enhanced by AS, leading to oxidative damage if not controlled (Foyer and Shigeoka, 2011; Voss et al., 2013). Interestingly, only one of those DEGs is directly involved in the ROS scavenging pathway. The cytosolic monodehydroascorbate reductase 2 (MDAR2), coded by LOC_Os08g44340, is an enzyme that regenerates ascorbate back from its oxidized form (monodehydroascorbate; Noctor et al., 2014). Most genes identified by this “interactomic approach” did not belong to this group of well-studied antioxidants and antioxidant enzymes, e.g., AsA, GSH, Trx, CAT, APX, and GPX (Foyer and Shigeoka, 2011; Maruta et al., 2012). Instead, the elected genes were directly involved in the photosynthetic apparatus or regulatory elements, such as WRKY, MYB, and TYFF transcription factors (Supplementary Table S5). More than half of these genes showed a similar expression pattern, being highly expressed in response to iron deficiency and drought in roots, while suppressed in As toxicity, ozone, and submergence (Figure 3). Curiously, they were induced in roots and shoots under drought stress, but repressed in seedlings in the same treatments. This apparent contradiction could reflect differences in the age of the samples or the methods used to simulate the stress. In the experiments where roots and shoots were collected, the hydroponic solution in which the plants grew were slowly drained (Wang et al., 2011), while the whole seedlings were dried in tissue paper (Jain et al., 2007). MDAR2 did not represents the same expression pattern as this group of photosynthesis-related genes (Figure 3) and was positioned as a terminal node of the predicted cluster, connected only with LOC_Os06g0539, that codes a plastid gene of unknown function (Figure 2B). Although we cannot be certain about the involvement of MDAR2 with the others DEGs of this cluster, many studies have shown the importance of MDAR enzymes to the response and tolerance to AS (Sultana et al., 2012; El Airaj et al., 2013).
Among the DEGs shared between OS and AS, which did not fall into that cluster, other genes possessing hub roles (connecting distinct signaling and metabolic pathways) were also identified, for example OsSRO1c (Similar to Radical-Induced Cell Death-One 1c, LOC_Os03g12820). Radical-induced Cell Death1 (AtRCD1) received its name due to the ozone hypersensitive phenotype observed in plant knock-outs to this gene (Ahlfors et al., 2004; Miao et al., 2006). AtRCD1 activity is modulated through oxidation by a GPX (AtGPX3; Ahlfors et al., 2004; Miao et al., 2006). Recent studies have demonstrated that GPX proteins play important roles as redox sensors and connect ROS signaling with hormonal signaling pathways (Fourquet et al., 2008; Passaia et al., 2013; Passaia and Margis-pinheiro, 2015). One way this connection occurs is by SRO proteins (SRO), that were shown to interact with many different transcription factors (e.g., DREB2A and COL10) and are involved in transcription factor regulation and complex formation (Ahlfors et al., 2004; Jaspers et al., 2009). AtRCD1 plays a role in the plant development and response to stress, mutants experiments show the participation of this protein in ethylene, ROS, salicylic acid, abiscisic acid (ABA), and jasmonic acid (JA) signaling pathway (Ahlfors et al., 2004; Jaspers et al., 2009, 2010). A recent study have demonstrated that OsSRO1c is induced in response to multiple stresses and was show to improve drought and OS tolerance by promoting stomatal closure and H2O2 accumulation (You et al., 2013).
A MA approach integrating different statistical methods allowed us to narrow down shared DEGs to a relatively small number that should be further investigated in detail. The comparison of shared DEGs with a list of genes a priori considered to be ROS-related highlighted the importance of redox homeostasis in stress response and signaling. Among the shared DEGs identified in this study are interesting candidates such as OsSRO1c, which regulate a great number of other proteins and connect different signaling pathways.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The reviewer Jun You and handling Editor Zhulong Chan declared their shared affiliation, and the handling Editor states that, nevertheless, the process met the standards of a fair and objective review.
The authors wish to thank Lin-Bo Wu (INRES, University of Bonn) for providing unpublished microarray data. The study was jointly funded by a Capes-Humboldt Research Fellowship (JBAN) – a cooperation between the Brazilian CAPES (Coordenação de Aperfeiçoamento de Pessoal de Nível Superior) and the German Alexander von Humboldt Foundation – and the Deutsche Forschungsgemeinschaft (FR2952/1-1).
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fpls.2015.01260
Ahlfors, R., Lang, S., and Overmyer, K. (2004). Arabidopsis RADICAL-INDUCED CELL DEATH1 belongs to the WWE protein–protein interaction domain protein family and modulates abscisic acid, ethylene, and methyl jasmonate responses. Cell 16, 1925–1937. doi: 10.1105/tpc.021832.2
Azizur Rahman, M., Hasegawa, H., Mahfuzur Rahman, M., Nazrul Islam, M., Majid Miah, M. A., and Tasmen, A. (2007). Effect of arsenic on photosynthesis, growth and yield of five widely cultivated rice (Oryza sativa L.) varieties in Bangladesh. Chemosphere 67, 1072–1079. doi: 10.1016/j.chemosphere.2006.11.061
Baier, M., Kandlbinder, A., Golldack, D., and Dietz, K. J. (2005). Oxidative stress and ozone: perception, signalling and response. Plant Cell Environ. 28, 1012–1020. doi: 10.1111/j.1365-3040.2005.01326.x
Baldi, P., and Long, A. D. (2001). A Bayesian framework for the analysis of microarray expression data: regularized t -test and statistical inferences of gene changes. Bioinformatics 17, 509–519. doi: 10.1093/bioinformatics/17.6.509
Barrett, T., Wilhite, S. E., Ledoux, P., Evangelista, C., Kim, I. F., Tomashevsky, M., et al. (2013). NCBI GEO: archive for functional genomics data sets–update. Nucleic Acids Res. 41, D991–D995. doi: 10.1093/nar/gks1193
Brazma, A., Hingamp, P., Quackenbush, J., Sherlock, G., Spellman, P., Stoeckert, C., et al. (2001). Minimum information about a microarray experiment (MIAME)-toward standards for microarray data. Nat. Genet. 29, 365–371. doi: 10.1038/ng1201-365
Chang, L., Lin, H., Sibille, E., and Tseng, G. C. (2013). Meta-analysis methods for combining multiple expression profiles: comparisons, statistical characterization and an application guideline. BMC Bioinformatics 14:368. doi: 10.1186/1471-2105-14-368
Chawla, S., Jain, S., and Jain, V. (2012). Salinity induced oxidative stress and antioxidant system in salt-tolerant and salt-sensitive cultivars of rice (Oryza sativa L.). J. Plant Biochem. Biotechnol. 22, 27–34. doi: 10.1007/s13562-012-0107-4
Cho, K., Shibato, J., Agrawal, G. K., Jung, Y. H., Kubo, A., Jwa, N. S., et al. (2008). Integrated transcriptomics, proteomics, and metabolomics analyses to survey ozone responses in the leaves of rice seedling. J. Proteome Res. 7, 2980–2998. doi: 10.1021/pr800128q
Cho, K., Shibato, J., Kubo, A., Kohno, Y., Satoh, K., Kikuchi, S., et al. (2013). Genome-wide mapping of the ozone-responsive transcriptomes in rice panicle and seed tissues reveals novel insight into their regulatory events. Biotechnol. Lett. 35, 647–656. doi: 10.1007/s10529-012-1118-x
Cotsaftis, O., Plett, D., Johnson, A. A. T., Walia, H., Wilson, C., Ismail, A. M., et al. (2011). Root-specific transcript profiling of contrasting rice genotypes in response to salinity stress. Mol. Plant 4, 25–41. doi: 10.1093/mp/ssq056
Dai, X., Wang, Y., Yang, A., and Zhang, W.-H. (2012). OsMYB2P-1, an R2R3 MYB transcription factor, is involved in the regulation of phosphate-starvation responses and root architecture in rice. Plant Physiol. 159, 169–183. doi: 10.1104/pp.112.194217
de Abreu-Neto, J. B., Turchetto-Zolet, A. C., de Oliveira, L. F. V., Bodanese Zanettini, M. H., and Margis-Pinheiro, M. (2013). Heavy metal-associated isoprenylated plant protein (HIPP): characterization of a family of proteins exclusive to plants. FEBS J. 280, 1604–1616. doi: 10.1111/febs.12159
DeConde, R. P., Hawley, S., Falcon, S., Clegg, N., Knudsen, B., and Etzioni, R. (2006). Combining results of microarray experiments: a rank aggregation approach. Stat. Appl. Genet. Mol. Biol. 5, 15. doi: 10.2202/1544-6115.1204
Dubey, S., Misra, P., Dwivedi, S., Chatterjee, S., Bag, S. K., Mantri, S., et al. (2010). Transcriptomic and metabolomic shifts in rice roots in response to Cr (VI) stress. BMC Genomics 11:648. doi: 10.1186/1471-2164-11-648
El Airaj, H., Gest, N., Truffault, V., Garchery, C., Riqueau, G., Gouble, B., et al. (2013). Decreased monodehydroascorbate reductase activity reduces tolerance to cold storage in tomato and affects fruit antioxidant levels. Postharvest Biol. Technol. 86, 502–510. doi: 10.1016/j.postharvbio.2013.07.035
Fourquet, S., Huang, M.-E., D’Autreaux, B., and Toledano, M. B. (2008). The dual functions of thiol-based peroxidases in H 2 O 2 scavenging and signaling. Antioxid. Redox Signal. 10, 1565–1576. doi: 10.1089/ars.2008.2049
Frei, M., Tanaka, J. P., Chen, C. P., and Wissuwa, M. (2010a). Mechanisms of ozone tolerance in rice: characterization of two QTLs affecting leaf bronzing by gene expression profiling and biochemical analyses. J. Exp. Bot. 61, 1405–1417. doi: 10.1093/jxb/erq007
Frei, M., Wang, Y., Ismail, A. M., and Wissuwa, M. (2010b). Biochemical factors conferring shoot tolerance to oxidative stress in rice grown in low zinc soil. Funct. Plant Biol. 37, 74–84. doi: 10.1071/FP09079
Fujita, M., Fujita, Y., Noutoshi, Y., Takahashi, F., Narusaka, Y., Yamaguchi-Shinozaki, K., et al. (2006). Crosstalk between abiotic and biotic stress responses: a current view from the points of convergence in the stress signaling networks. Curr. Opin. Plant Biol. 9, 436–442. doi: 10.1016/j.pbi.2006.05.014
Gill, S. S., and Tuteja, N. (2010). Reactive oxygen species and antioxidant machinery in abiotic stress tolerance in crop plants. Plant Physiol. Biochem. 48, 909–930. doi: 10.1016/j.plaphy.2010.08.016
Guo, Z., Ou, W., Lu, S., and Zhong, Q. (2006). Differential responses of antioxidative system to chilling and drought in four rice cultivars differing in sensitivity. Plant Physiol. Biochem. 44, 828–836. doi: 10.1016/j.plaphy.2006.10.024
Höller, S., Meyer, A., and Frei, M. (2014). Zinc deficiency differentially affects redox homeostasis of rice genotypes contrasting in ascorbate level. J. Plant Physiol. 171, 1748–1756. doi: 10.1016/j.jplph.2014.08.012
Hu, W., Hu, G., and Han, B. (2009). Genome-wide survey and expression profiling of heat shock proteins and heat shock factors revealed overlapped and stress specific response under abiotic stresses in rice. Plant Sci. 176, 583–590. doi: 10.1016/j.plantsci.2009.01.016
Jain, M., Nijhawan, A., Arora, R., Agarwal, P., Ray, S., Sharma, P., et al. (2007). F-box proteins in rice. Genome-wide analysis, classification, temporal and spatial gene expression during panicle and seed development, and regulation by light and abiotic stress. Plant Physiol. 143, 1467–1483. doi: 10.1104/pp.106.091900
Jaiswal, P., Ware, D., Ni, J., Chang, K., Zhao, W., Schmidt, S., et al. (2002). Gramene: development and integration of trait and gene ontologies for rice. Comp. Funct. Genomics 3, 132–136. doi: 10.1002/cfg.156
Jaspers, P., Blomster, T., Brosché, M., Salojärvi, J., Ahlfors, R., Vainonen, J. P., et al. (2009). Unequally redundant RCD1 and SRO1 mediate stress and developmental responses and interact with transcription factors. Plant J. 60, 268–279. doi: 10.1111/j.1365-313X.2009.03951.x
Jaspers, P., Overmyer, K., Wrzaczek, M., Vainonen, J. P., Blomster, T., Salojärvi, J., et al. (2010). The RST and PARP-like domain containing SRO protein family: analysis of protein structure, function and conservation in land plants. BMC Genomics 11:170. doi: 10.1186/1471-2164-11-170
Kauffmann, A., and Huber, W. (2013). Introduction: Microarray Quality Assessment with ArrayQualityMetrics, 1–7. Available at: http://bioconductor.org/packages/arrayQualityMetrics/
Kim, S.-H., Song, M., Lee, K. J., Hwang, S.-G., Jang, C. S., Kim, J.-B., et al. (2012). Genome-wide transcriptome profiling of ROS scavenging and signal transduction pathways in rice (Oryza sativa L.) in response to different types of ionizing radiation. Mol. Biol. Rep. 39, 11231–11248. doi: 10.1007/s11033-012-2034-9
Li, X., Bu, N., Li, Y., Ma, L., Xin, S., and Zhang, L. (2012). Growth, photosynthesis and antioxidant responses of endophyte infected and non-infected rice under lead stress conditions. J. Hazard. Mater. 21, 55–61. doi: 10.1016/j.jhazmat.2012.01.052
Liu, F., Xu, W., Wei, Q., Zhang, Z., Xing, Z., Tan, L., et al. (2010). Gene Expression Profiles Deciphering Rice Phenotypic Variation between Nipponbare (Japonica) and 93-11 (Indica) during Oxidative Stress. PLoS ONE 5:e8632. doi: 10.1371/journal.pone.0008632
Maruta, T., Noshi, M., Tanouchi, A., Tamoi, M., Yabuta, Y., Yoshimura, K., et al. (2012). H2O2-triggered retrograde signaling from chloroplasts to nucleus plays specific role in response to stress. J. Biol. Chem. 287, 11717–11729. doi: 10.1074/jbc.M111.292847
Matthus, E., Wu, L.-B., Ueda, Y., Höller, S., Becker, M., and Frei, M. (2015). Loci, genes, and mechanisms associated with tolerance to ferrous iron toxicity in rice (Oryza sativa L.). Theor. Appl. Genet. 128, 2085–2098. doi: 10.1007/s00122-015-2569-y
Miao, Y., Lv, D., Wang, P., Wang, X.-C., Chen, J., Miao, C., et al. (2006). An Arabidopsis glutathione peroxidase functions as both a redox transducer and a scavenger in abscisic acid and drought stress responses. Plant Cell 18, 2749–2766. doi: 10.1105/tpc.106.044230
Miller, G., Suzuki, N., Ciftci-Yilmaz, S., and Mittler, R. (2010). Reactive oxygen species homeostasis and signalling during drought and salinity stresses. Plant Cell Environ. 33, 453–467. doi: 10.1111/j.1365-3040.2009.02041.x
Mittal, D., Madhyastha, D. A., and Grover, A. (2012a). Gene expression analysis in response to low and high temperature and oxidative stresses in rice: combination of stresses evokes different transcriptional changes as against stresses applied individually. Plant Sci. 197, 102–113. doi: 10.1016/j.plantsci.2012.09.008
Mittal, D., Madhyastha, D. A., and Grover, A. (2012b). Genome-wide transcriptional profiles during temperature and oxidative stress reveal coordinated expression patterns and overlapping regulons in rice. PLoS ONE 7:e40899. doi: 10.1371/journal.pone.0040899
Mustroph, A., Lee, S. C., Oosumi, T., Zanetti, M. E., Yang, H., Ma, K., et al. (2010). Cross-kingdom comparison of transcriptomic adjustments to low-oxygen stress highlights conserved and plant-specific responses. Plant Physiol. 152, 1484–1500. doi: 10.1104/pp.109.151845
Norton, G. J., Douglas, A., Lahner, B., Yakubova, E., Guerinot, M. L., Pinson, S. R., et al. (2014). Genome wide association mapping of grain arsenic, copper, molybdenum and zinc in rice (Oryza sativa L.) grown at four international field sites. PLoS ONE 9:e89685. doi: 10.1371/journal.pone.0089685
Ogo, Y., Kakei, Y., Itai, R. N., Kobayashi, T., Nakanishi, H., Takahashi, H., et al. (2014). Spatial transcriptomes of iron-deficient and cadmium-stressed rice. New Phytol. 201, 781–794. doi: 10.1111/nph.12577
Oliveiros, J. (2015). Venny. An Interactive Tool for Comparing Lists with Venn’s Diagrams. Available at: http://bioinfogp.cnb.csic.es/tools/venny/index.html [Accessed July 27, 2015].
Pandit, A., Rai, V., Sharma, T. R., Sharma, P. C., and Singh, N. K. (2011). Differentially expressed genes in sensitive and tolerant rice varieties in response to salt-stress. J. Plant Biochem. Biotechnol. 20, 149–154. doi: 10.1007/s13562-010-0022-5
Pariasca-Tanaka, J., Satoh, K., Rose, T., Mauleon, R., and Wissuwa, M. (2009). Stress response versus stress tolerance: a transcriptome analysis of two rice lines contrasting in tolerance to phosphorus deficiency. Rice 2, 167–185. doi: 10.1007/s12284-009-9032-0
Passaia, G., Fonini, L. S., Caverzan, A., Jardim-Messeder, D., Christoff, A. P., Gaeta, M. L., et al. (2013). The mitochondrial glutathione peroxidase GPX3 is essential for H2O2 homeostasis and root and shoot development in rice. Plant Sci. 208, 93–101. doi: 10.1016/j.plantsci.2013.03.017
Rabbani, M. A., Maruyama, K., Abe, H., Khan, M. A., Katsura, K., Ito, Y., et al. (2003). Monitoring expression profiles of rice genes under cold, drought, and high-salinity stresses and abscisic acid application using cDNA microarray and RNA gel-blot analyses. Plant Physiol. 133, 1755–1767. doi: 10.1104/pp.103.025742.genes
Rasmussen, S., Barah, P., Suarez-Rodriguez, M. C., Bressendorff, S., Friis, P., Costantino, P., et al. (2013). Transcriptome responses to combinations of stresses in Arabidopsis. Plant Physiol. 161, 1783–1794. doi: 10.1104/pp.112.210773
Rhodes, D. R., Barrette, T. R., Rubin, M. A., Ghosh, D., and Chinnaiyan, A. M. (2002). Meta-analysis of microarrays: interstudy validation of gene expression profiles reveals pathway dysregulation in prostate cancer. Cancer Res. 62, 4427–4433.
Shaik, R., Ramakrishna, W., and Open, R. C. W. (2014). Machine learning approaches distinguish multiple stress conditions using stress-responsive genes and identify candidate genes for broad resistance in rice. Plant Physiol. 164, 481–495. doi: 10.1104/pp.113.225862
Shankar, A., Singh, A., Kanwar, P., Srivastava, A. K., Pandey, A., Suprasanna, P., et al. (2013). Gene expression analysis of rice seedling under potassium deprivation reveals major changes in metabolism and signaling components. PLoS ONE 8:e70321. doi: 10.1371/journal.pone.0070321
Smyth, G. K. (2005). “limma: linear models for microarray data,” in Bioinformatics and Computational Biology Solutions Using R and Bioconductor Statistics for Biology and Health, eds R. Gentleman, V. Carey, W. Huber, R. Irizarry, and S. Dudoit (New York, NY: Springer-Verlag), 397–420.
Smyth, G. K., Yang, Y. H., and Speed, T. (2003). “Statistical Issues in cDNA microarray data analysis,” in Functional Genomics Methods in Molecular Biology, eds M. Brownstein and A. Khodursky (New Jersey, NY: Humana Press), 111–136.
Sultana, S., Khew, C. Y., Morshed, M. M., Namasivayam, P., Napis, S., and Ho, C. L. (2012). Overexpression of monodehydroascorbate reductase from a mangrove plant (AeMDHAR) confers salt tolerance on rice. J. Plant Physiol. 169, 311–318. doi: 10.1016/j.jplph.2011.09.004
Szklarczyk, D., Franceschini, A., Wyder, S., Forslund, K., Heller, D., Huerta-Cepas, J., et al. (2015). STRING v10: protein-protein interaction networks, integrated over the tree of life. Nucleic Acids Res. 43, D447–D452. doi: 10.1093/nar/gku1003
Tseng, G. C., Ghosh, D., and Feingold, E. (2012). Comprehensive literature review and statistical considerations for microarray meta-analysis. Nucleic Acids Res. 40, 3785–3799. doi: 10.1093/nar/gkr1265
Tusher, V. G., Tibshirani, R., and Chu, G. (2001). Significance analysis of microarrays applied to the ionizing radiation response. Proc. Natl. Acad. Sci. U.S.A. 98, 5116–5121. doi: 10.1073/pnas.091062498
Uchida, A., Jagendorf, A. T., Hibino, T., Takabe, T., and Takabe, T. (2002). Effects of hydrogen peroxide and nitric oxide on both salt and heat stress tolerance in rice. Plant Sci. 163, 515–523. doi: 10.1016/S0168-9452(02)00159-0
Uraguchi, S., Kamiya, T., Sakamoto, T., Kasai, K., Sato, Y., Nagamura, Y., et al. (2011). Low-affinity cation transporter (OsLCT1) regulates cadmium transport into rice grains. Proc. Natl. Acad. Sci. U.S.A. 108, 20959–20964. doi: 10.1073/pnas.1116531109
Voss, I., Sunil, B., Scheibe, R., and Raghavendra, A. S. (2013). Emerging concept for the role of photorespiration as an important part of abiotic stress response. Plant Biol. (Stuttg.) 15, 713–722. doi: 10.1111/j.1438-8677.2012.00710.x
Wang, D., Pan, Y., Zhao, X., Zhu, L., Fu, B., and Li, Z. (2011). Genome-wide temporal-spatial gene expression profiling of drought responsiveness in rice. BMC Genomics 12:149. doi: 10.1186/1471-2164-12-149
Wang, X., Kang, D. D., Shen, K., Song, C., Lu, S., Chang, L.-C., et al. (2012). An R package suite for microarray meta-analysis in quality control, differentially expressed gene analysis and pathway enrichment detection. Bioinformatics 28, 2534–2536. doi: 10.1093/bioinformatics/bts485
Warnes, A. G. R., Bolker, B., Bonebakker, L., Huber, W., Liaw, A., Lumley, T., et al. (2015). Package “gplots.” Available at: https://cran.r-project.org/web/packages/gplots/index.html
You, J., Zong, W., Li, X., Ning, J., Hu, H., Li, X., et al. (2013). The SNAC1-targeted gene OsSRO1c modulates stomatal closure and oxidative stress tolerance by regulating hydrogen peroxide in rice. J. Exp. Bot. 64, 569–583. doi: 10.1093/jxb/ers349
Keywords: microarray, meta-analysis, abiotic stress, rice, redox
Citation: de Abreu Neto JB and Frei M (2016) Microarray Meta-Analysis Focused on the Response of Genes Involved in Redox Homeostasis to Diverse Abiotic Stresses in Rice. Front. Plant Sci. 6:1260. doi: 10.3389/fpls.2015.01260
Received: 28 August 2015; Accepted: 24 December 2015;
Published: 12 January 2016.
Edited by:Zhulong Chan, Chinese Academy of Sciences, China
Reviewed by:Hans-Peter Piepho, University of Hohenheim, Germany
Jun You, Chinese Academy of Sciences, China
Gisele Passaia Prietsch, Carnegie Institution for Science, USA
Copyright © 2016 de Abreu Neto and Frei. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Michael Frei, email@example.com