Transcriptomic complexity in young maize primary roots in response to low water potentials
© Opitz et al.; licensee BioMed Central Ltd. 2014
Received: 14 May 2014
Accepted: 22 August 2014
Published: 29 August 2014
Widespread and more frequently occurring drought conditions are a consequence of global warming and increase the demand for tolerant crop varieties to feed the growing world population. A better understanding of the molecular mechanisms underlying the water deficit response of crops will enable targeted breeding strategies to develop robust cultivars.
In the present study, the transcriptional response of maize (Zea mays L.) primary roots to low water potentials was monitored by RNA sequencing (RNA-Seq) experiments. After 6 h and 24 h of mild (-0.2 MPa) and severe (-0.8 MPa) water deficit conditions, the primary root transcriptomes of seedlings grown under water deficit and control conditions were compared. The number of responsive genes was dependent on and increased with intensification of water deficit treatment. After short-term mild and severe water deficit 249 and 3,000 genes were differentially expressed, respectively. After a 24 h treatment the number of affected genes increased to 7,267 and 12,838 for mild and severe water deficit, respectively, including more than 80% of the short-term responsive genes. About half of the differentially expressed genes were up-regulated and maximal fold-changes increased with treatment intensity to more than 300-fold. A consensus set of 53 genes was differentially regulated independently of the nature of deficit treatment. Characterization revealed an overrepresentation of the Gene Ontology (GO) categories “oxidoreductase activity” and “heme binding” among regulated genes connecting the water deficit response to ROS metabolism.
This study gives a comprehensive insight in water deficit responsive genes in young maize primary roots and provides a set of candidate genes that merit further genetic analyses in the future.
Population growth and global warming are major challenges for global food security. It is estimated that the demand for agricultural products will increase by ~50% until 2030 . This requires historically unprecedented annual production growth rates . Worldwide, about 70% of the food production is provided by cereals . Maize (Zea mays L.) outcompeted all other cereals with an estimated global yield of 863 million tons in 2012/2013 . While cereal production needs to be significantly increased, climate change adversely affects global maize production with an estimated loss of ~4% relative to what could have been achieved without the climate trends . Poor soil moisture is widespread among arable land and as a consequence of global warming more areas are affected by drought conditions each year . Since water availability is the most critical environmental factor for plant growth , drought can limit crop productivity more than any other abiotic stress. Furthermore, variations in water availability within fields can result in uneven crop stands that cause yield losses .
Under drought conditions, when water loss through transpiration is high, it is essential that roots maintain the capacity to acquire soil water and nutrients. This is reflected by the ability of roots to continue elongation even under severe water deficit conditions albeit at a slower rate . From a physiological perspective, root growth maintenance is predominantly regulated by the plant hormone abscisic acid (ABA). Accumulation of ABA suppresses excessive ethylene production and thereby prevents growth inhibition. ABA is further involved in the processes leading to osmotic adjustment as it promotes the transport of proline to the root apex. At the more basal regions of the root, hexoses are the predominant solutes providing osmotic adjustment and maintaining turgor pressure [reviewed in . From a cellular viewpoint, the processes related to the water deficit response begin with stress perception, followed by signal transduction, and a change in gene expression that finally confers the complex metabolic and physiological alterations necessary to gain stress tolerance [11, 12]. On the molecular level, genes regulated by water deficit can be grouped into two categories. The first group of genes encodes proteins providing direct stress tolerance such as chaperones, transporters, osmolytic and detoxifying proteins, and repair-enzymes . The second category includes proteins involved in stress response by regulating signal transduction and gene expression for instance transcription factors, protein kinases and phosphatases, and other signaling molecules . The high quantity of genes regulated upon water deficit reflects the complexity of the stress response . Nevertheless, details of the translation of environmental changes to metabolic responses i.e. the adjustment of transcriptional and post-transcriptional modifications of metabolic enzymes still remains unclear .
In the past, microarray chip hybridization experiments monitored gene expression profiles of maize leaves and roots to elucidate the transcriptional changes upon water deficit [14–19]. Recently developed next-generation sequencing approaches such as RNA sequencing (RNA-Seq) allow fully quantitative gene expression analyses  of all 39,656 (FGSv2; , release 5b.60) high-confidence maize gene models currently annotated . The digital nature of the method enables the detection of a large dynamic range of expression levels with absolute values and the capture of even subtle gene expression changes [23, 24].
In the present study, a tightly controlled, reproducible experimental setup was applied to expose maize seedlings to water deficit conditions. Young seedlings were grown in paper rolls soaked with polyethylene glycol (PEG) solutions. Previously, it was demonstrated that PEG treatment is an effective way to simulate drought stress conditions occurring in drying soil . Water potentials of -0.2 MPa and -0.8 MPa were used for mild and severe water deficit treatments, respectively. At -0.2 MPa root growth is not or only slightly affected while shoot growth is reduced by half . The lower water potential of -0.8 MPa completely inhibits elongation of maize shoots and leaves while roots continue to elongate consistently [9, 25, 26]. This maintenance of root growth in drying soil is beneficial to plants as they can reach deeper water resources. As top soil layers are prone to drying, it is particularly important for seedlings to adapt to low water potentials . Although maize is most susceptible to drought stress during the flowering period , drought conditions during the seedling stage can negatively affect its yield .
To gain a better understanding of the early molecular responses to water deficit we utilized an RNA-Seq approach and compared the root transcriptomes of stressed and control maize seedlings. Seedlings were subjected to mild and severe water deficit conditions for 6 h and 24 h. The overall goal of this study was to identify a set of genes involved in initial water deficit responses in maize primary roots. Furthermore, application of mild and severe water deficit conditions at two time points aimed at detecting specifically and commonly regulated genes across treatment intensities and time. This data set will be a resource for future genetic analyses of candidate genes involved in water deficit response in young maize primary roots.
Characterization of water deficit treatment and phenotypic responses
To analyze phenotypic stress responses, primary root length was measured before and after treatment (Figure 1B). Short-term water deficit did not affect root elongation significantly. On average roots elongated 0.75 cm in 6 h. Similarly, after 24 h of water deficit treatment no differences of root growth between control and mildly stressed seedlings were observed (2.5 cm). However, severe water deficit significantly reduced root elongation (1.8 cm in 24 h) by almost 30% (p <0.001).
RNA sequencing and mapping of maize primary root transcriptomes
To identify genes responsive to water deficit in young maize primary roots, global gene expression was surveyed by Illumina RNA-Seq. On average, the RNA-Seq experiments yielded between 20 and 34 million reads per sample. The raw sequencing data has been deposited at the NCBI sequencing read archive (SRA, , Accession: SRP032921). Among all reads, 83 to 90% mapped to unique positions in the maize reference genome (ZmB73_RefGen_v2; Additional file 1). After removal of redundant reads sharing the same start and end coordinate, sequencing direction, and sequence (“stacked reads”), 76 to 78% of the remaining reads mapped uniquely to the “filtered gene set” (FGSv2, release 5b.60; Additional file 1), a set of 39,656 high confidence gene models predicted by a combination of evidence-based and ab initio approaches followed by stringent TE filtering [21, 22]. A gene was declared expressed if a minimum of five reads mapped in all four replicates of a sample. As a result, 25,570 genes (64%) of the FGSv2 were expressed in at least one of the experimental conditions. A complete list of expressed genes with normalized expression values is provided in Additional file 2.
Exploration of differentially expressed genes in response to water deficit
Besides the pairwise comparisons a two-way analysis of variance was performed to determine treatment main effects, time main effects, and treatment by time interactions. This analysis yields genes differentially regulated in response to both water deficit levels, both treatment periods, and to the combination of treatment intensity and duration, respectively. The latter gene set includes genes whose expression does not change initially but changes over time and treatment intensity . As a result of this survey, 6,266 and 15,575 genes displayed treatment main effects and time main effects, respectively. Moreover, 23 genes displayed treatment by time interactions (Additional file 4). Between the three gene sets some overlap was identified: 15 genes were present in all data sets representing 65% (15/23) of the genes that displayed treatment by time interactions. Furthermore, 75% of the genes that showed treatment main effects were among the genes displaying time main effects (Figure 3B).
Functional categorization of stress responsive genes
GO enrichment analyses were also performed for treatment and time main effects and treatment by time interaction specific genes. For 8,609 (55%) of 15,575 time main effect, 3,826 (58%) of 6,266 treatment main effect genes, and 14 (61%) of 23 treatment by time interaction genes at least one GO term was assigned (Additional file 4). In contrast to the analysis of pairwise comparisons, comparison of these gene sets with the set of all expressed genes by SEA yielded many overrepresented cellular component categories. Among the treatment by time interaction genes no enriched term was found. Among treatment main effect and time main effect genes 10 and 16 categories were overrepresented, respectively (Additional file 7). Each of the terms enriched among treatment main effect genes was included in the enriched time main effect terms. Four terms were connected to protein synthesis (“translation” GO:0006412, “structural constituent of ribosome” GO:0003735, “ribonucleoprotein complex” GO:0030529, and “ribosome” GO:0005840). Further categories were “structural molecule activity” (GO:0005198), the cellular components organelles (“intracellular non-membrane-bounded organelle” GO:0043232, “non-membrane-bounded organelle” GO:0043228), cytosol (“cytoplasmic part” GO:0044444, “cytoplasm” GO:0005737), and “macromolecular complex” (GO:0032991). Furthermore, the time main effect terms included five transport-related terms (“hydrogen ion transmembrane transporter activity” GO:0015078, “monovalent inorganic cation transmembrane transporter activity” GO:0015077, “inorganic cation transmembrane transporter activity” GO:0022890, “proton-transporting ATP synthase complex” GO:0045259, “proton-transporting two-sector ATPase complex GO:0016469) and the nucleosome (GO:0000786) as additional cellular compartment.
Validation of differentially expressed genes by qRT-PCR
To independently confirm the RNA-Seq results a subset of nine genes differentially regulated in response to all treatments was selected for quantitative real-time PCR (qRT-PCR) analyses (Additional file 8). As a result, 30 of 36 RNA-Seq expression patterns were confirmed i.e. trends of gene expression were identical in the sequencing and qRT-PCR approach. For nine of the 30 data points the analysis did not have enough statistical power to detect significant differences while for 21/30 data points statistically significant differences (Student’s t-test: p-value <0.05) were measured.
To promote plant survival, roots are often able to continue growth at low water potentials that completely inhibit shoot elongation [9, 25, 26]. As a result, roots might escape dry soil layers and reach deeper water resources. When water availability is so low to be life threatening, yield will be very low even if plants survive. Thus, from an agronomic viewpoint, mere survival of grain crops is irrelevant . Therefore rather moderate water deficits were chosen for the present study. Water potentials of -0.2 MPa and -0.8 MPa were used to simulate mild and severe water deficit conditions, respectively. Short-term (6 h) water deficit treatment did not affect root elongation in the present study. This was expected as marked differences in steady maize root elongation rates between control conditions and low water potentials were observed only after 10 h of treatment . Similarly, 24 h of mild water deficit did not impair primary root growth, while severe water deficit significantly reduced primary root elongation. These findings are comparable to results by Sharp and co-workers [9, 25] whereas Westgate and Boyer  measured considerable growth reduction only at water potentials lower than -0.8 MPa. This discrepancy might at least in part be explained by the different maize genotypes that were analyzed in these studies. Since 6 h of water deficit treatment did not affect root growth this experimental period was used to monitor gene expression changes at the initial phase of the water deficit response. Long-term treatment with mild water deficit had likewise no primary root growth effect suggesting that the root is able to tolerate this mild drought. In contrast, severe water deficit decelerated primary root elongation implying that the root had adapted by metabolic alterations.
RNA-Seq technology was utilized to compare the primary root transcriptomes of seedlings treated with low water potentials to those grown under control conditions. The seedlings used for the experiment were selected by their root length rather than by age minimizing variability of gene expression due to developmental differences . In total, 25,570 genes were expressed under at least one of the time/treatment combinations. In the present study, 74 and 669 genes were identified that responded with expression changes to 6 h of mild and severe water deficit, respectively (|Fc| ≥2). After 24 h, 1,981 and 3,006 genes were differentially expressed (|Fc| ≥2) in response to mild and severe water deficit, respectively, reflecting complex acclimatization processes. A similar pattern of an increased number of responsive genes following more intense stress treatment has been reported for maize (Han21) seedling roots by Zheng, et al. , who considered 190 uniquely expressed sequence tags, and also for the whole maize (SRG-200) root system . Comparative analyses in this study revealed that most (≥80%) of the short-term water deficit responsive genes were a subset of the long-term water deficit responsive genes. Similarly, a large overlap (≥90%) between mild and severe water deficit responsive genes was detected. An analogous pattern was observed for drought-stressed maize (SRG-200) and barley root systems although stress treatment was more intense in these studies [15, 36]. However, neither maize nor barley leaves showed this pattern [15, 19, 36]. This indicates a root-specific response that is already established after the short period and low intensity of water deficit used in this study. Besides, there is further indication that the drought stress response is organ-specific. A microarray analysis of maize root, leaf, and shoot gene expression changes in response to drought stress revealed only a rather limited overlap between stress regulated genes . In line with that observation the functional categories affected by the transcriptional changes were very distinct across tissues . Similarly, a comparison of differentially expressed genes identified in the present study to genes identified by Kakumanu, et al.  via Illumina sequencing of RNAs isolated from basal leaf meristem tissues of maize plants grown under drought stress yielded very little overlap both on gene and functional category level. However, plant developmental stages (seedlings vs. adult plants, reproductive stage) and treatment period (24 h vs. 4 d) were quite different.
Among the water deficit responsive genes identified in this study several GO categories were overrepresented in comparison to all expressed genes. Cross-comparison of these categories resulted in one commonly enriched term, “oxidoreductase activity”. Its child term “monooxygenase activity” was also enriched among three of the four gene sets. Several previous studies detected differential expression of oxidoreductase and monooxygenase genes in response to water deficit in maize, barley, and soybean [14, 18, 36, 38–40]. However, the complexity of these studies was limited to a few hundred or thousand entities and only individual genes or proteins of the oxidoreductase category were identified along with other differentially expressed genes. The reactions catalyzed by oxidoreductases can result in scavenging as well as in the generation of reactive oxygen species (ROS). Differential expression, including both up and down-regulation, of ROS scavengers in response to drought stress has also been shown in Arabidopsis [reviewed in . While ROS can cause cellular damage , they also play an important role as signaling molecules [reviewed by  and are associated with growth and development [reviewed in . For instance, the expression and protein abundance of an oxalate oxidase that produces hydrogen peroxide increased in stressed maize roots [18, 40] and is probably involved in root growth regulation upon water deficit . Other commonly enriched GO terms among water deficit responsive genes in the present study were “heme binding”, its ancestor term “tetrapyrrole binding” and the functionally related category “iron ion binding”. Likewise, genes annotated in the Mapman category “metal handling” were commonly up-regulated in response to water deficit (data not shown). In wheat roots an increase in total iron content during drought stress was detected . Besides, the level of bound iron in soybean roots was substantially higher under water stress due to sequestration by ferritin proteins . The expression of these and other metal-chelating proteins was also up-regulated in stressed maize roots . Sequestration of metal ions is beneficial as it prevents the formation of the highly toxic hydroxyl radicals via the metal-dependent Haber-Weiss reaction or the Fenton reaction . Additional studies are needed to fully elucidate the complex interaction of ROS metabolism and regulation of gene expression upon water deficit.
After 24 h of both mild and severe water deficit treatment, further enrichment of GO terms related to sugar and general carbohydrate metabolic processes (“carbohydrate metabolic process”, its child terms “polysaccharide metabolic process”, “polysaccharide catabolic process”, and “disaccharide metabolic process” as well as “glycoside metabolic process” and its child term “hydrolase activity”) was detected. Accordingly, many of the long-term water deficit responsive genes were annotated in pathways of C compound metabolisms. Accumulation of carbohydrates in response to water deficit has been reported in several species and plant organs [reviewed in . In maize primary roots, soluble carbohydrates together with proline account for osmotic adjustment which has an essential role in maintenance of root elongation at low water potentials [48, 49]. It has been reported that both increases in proline synthesis and decreases in proline oxidation occur in response to low water potentials to increase proline concentrations . Accordingly, three of four biochemical pathways conserved among water deficit responsive genes to all treatments included down-regulated proline oxidases. Besides functioning in osmotic adjustment, C compounds also participate in signaling and transcriptional and post-translational regulatory processes in metabolic and developmental programs  and the regulation of their synthesis/degradation is therefore tight and complex in response to environmental stresses.
An analysis of variance yielded differentially expressed genes in treatment main effect, time main effect, and treatment by time interactions. The highest number of genes (15,575) appeared in time main effect including more than half of all expressed genes. This suggests strong differences between root transcriptomes of seedlings differing in age as seedlings after long-term experiments were 18 h (about 20%) older and on average 1.6 cm (about 42%) longer than seedlings after short-term experiments. This is in accordance with the report of differences in protein abundance between two early stages (5 and 9 DAG) of maize root development in the inbred line B73  which was also studied in the present survey. As most of the treatment-specific genes (4,662 of 6,266) are included in time effect genes, the pairwise comparisons between treatments within each time period might be more accurate. The test for treatment by time interaction could identify genes that do not differ across water deficit conditions initially but develop differences across conditions with longer treatment periods . A comparable low number of 23 genes were assigned to this category. Hence, only a minor fraction of all water deficit responsive genes changed their expression from one time point to the other. This is in line with the large overlap observed between mild and severe deficit as well as between short and long-term treatment responsive genes in the pairwise comparisons.
A GO analysis revealed largely the same overrepresented categories among genes that display time main effects and treatment main effects. About one quarter and one half of the enriched GO terms, respectively, are related to protein biosynthesis including the biological process translation and the ribosome as cellular component. This points to a developmental function of these genes specific for time and for treatment effects as ribosomes are generally considered as housekeeping components of the cell and are involved in growth processes. Genes encoding ribosomal proteins are highly expressed in proliferating, elongating, and differentiating cells . However, it was demonstrated for Arabidopsis that expression of these genes changes in response to abiotic stress [52, 53]. Furthermore, a study in maize indicated that modifications of the translational machinery emerge in response to hypoxic stress . In the present study the same genes functionally related to ribosomes were identified for time main effects and treatment main effects. This supports the concept of these genes playing a double role in normal development/growth and in water deficit response.
In the present study, hundreds of genes were identified that are differentially expressed in response to water deficits in maize seedling primary roots. The number of responsive genes was dependent on water deficit condition and duration and increased with intensification of treatment. The products of the differentially expressed genes are known to be involved in perception and signal transduction or confer adaptation and tolerance to low water potentials. Analyses of associated GO categories and underlying biochemical pathways connect the water deficit response to ROS and carbohydrate metabolisms and signaling, revealing their complex transcriptional regulation. Further functional analyses of these genes will promote our understanding of the molecular mechanisms underlying root adaptation to water deficit and enable targeted breeding strategies. Thus, more tolerant crop varieties, that can be cultivated with stable yield in dry regions as well as on fields with uneven soil moisture, might be generated.
Plant material and water deficit stress treatment
Seeds of the maize inbred line B73 were surface sterilized and germinated in paper rolls as described in Ludwig, et al. . Four to 5 days (16 h light, 28°C; 8 h dark, 21°C) after germination seedlings with a primary root length of 2 to 4 cm were transferred to new paper rolls soaked with polyethylene glycol (PEG8000 Mr 7,300-9,000; Roth, Karlsruhe, Germany) solutions with a water potential of -0.2 MPa and -0.8 MPa for mild and severe water deficit, respectively, or distilled water for control experiments. Because manufactured PEG8000 includes a range of molar masses water potentials could not be calculated directly. Osmolyte concentrations (mol*l-1) of PEG solutions with defined mass concentrations (g*l-1) were measured with an osmometer (OSMOMAT 030-D, Gonotec GmbH, Berlin, Germany). This data was used to estimate water potentials (water potential [MPa] = concentration [mol*l-1] * gas constant [8.314 Pa*l*mol-1*K-1] * temperature [298.15 K]). Paper rolls with seedlings were incubated in aerated PEG solutions or distilled water for 6 h and 24 h (24°C, 24 h included 8 h darkness, 18°C). Before and after treatment, seedlings were photographed and primary root length was measured using WinRHIZO (http://www.regent.qc.ca/). Whole primary roots were harvested, immediately frozen in liquid nitrogen and stored at -80°C until RNA isolation. Experiments were performed in four biological replicates each consisting of 10 pooled roots.
RNA isolation and sequencing library preparation
Pooled primary roots were ground in liquid nitrogen, and RNA was extracted as previously described . RNA quality was assessed via agarose gel electrophoresis and a Bioanalyzer (Agilent RNA 6000 Nano Chip, Agilent Technologies, Santa Clara CA, USA). Only samples with an RIN [RNA integrity number;  ≥9.1 were used for downstream analyses. The cDNA libraries for Illumina sequencing were constructed in accordance with the protocol of the manufacturer (TruSeq RNA Sample Preparation, Illumina, San Diego CA, USA). For sequencing, four libraries were pooled in one lane of a flow cell. Each library per lane was indexed by one of the adapters AR001, AR008, AR010, or AR011 (Additional file 1). The indexed libraries were loaded onto a flow cell according to an incomplete block design (generated with CycDesigN, http://www.vsni.co.uk/software/cycdesign; Additional file 1). Cluster preparation and single read sequencing were performed according to the manufacturer’s instructions (HiSeq 2000, Illumina, San Diego CA, USA).
Processing and mapping of Illumina sequencing reads
Raw sequencing reads generated by the Illumina HiSeq 2000 system were initially processed and quality trimmed with SHORE (http://1001genomes.org/software/shore.html). Reads with more than 2 or 5 bases having quality scores ≤3 in the first 12 or 25 bases, respectively, were rejected . Bases with quality scores ≤5 at the 3’ end were trimmed until two succeeding bases with higher quality scores. Reads with ≥2 mismatches in adapter sequences were excluded. Only reads ≥40 bp were retained for subsequent analyses. Finally, adapter sequences were removed. Resulting reads had a length of 40 to 100 bp (≥60% of all reads were 100 bp long) and quality scores of 26 or higher at all base positions. Read mapping was performed with CLC Genomics Workbench (http://www.clcbio.com/products/clc-genomics-workbench/). All high quality reads were mapped to the maize B73 reference genome [RefGen_v2; [21, 22] allowing large gaps of up to 50 kb to span introns. At least 75% of each read had to fit with 90% similarity to the reference to be mapped. Stacked reads i.e. redundant reads sharing the same start and end coordinate, sequencing direction, and sequence were merged into one. The remaining reads were projected to the filtered gene set (FGSv2; , release 5b.60) of the B73 reference genome derived from the maize genome sequencing project (MGSP) allowing a maximum of two mismatches for reads ≤56 bp. Longer reads had to fit at least with 80% of their length thereby comprising 90% similarity. Only those reads uniquely mapping to the reference data set were subsequently used for analyses.
Statistical procedures for analyzing differential gene expression
The read counts were normalized to RPKM (reads per kilobase of exon model per million mapped reads) values  and log2 transformed to meet the assumptions of linear models. Further statistical analysis followed the empirical Bayes approach of Smyth . The mean-variance trend for log-counts was estimated and a weight assigned to each observation based on its predicted variance. The weights were then used in the linear modeling process to adjust for heteroscedasticity . To borrow strength across genes in the estimation of the residual error variance, the empirical Bayes approach of Smyth  implemented in the R-package limma was used. A linear mixed model with lane effect was applied accounting for the lanes as incomplete blocks. The lane effect is considered as random effect, thus allowing the recovery of the inter-block information. The experimental setup allowed four comparisons of control groups against different water deficit levels: mild deficit, 6 h (1), severe deficit, 6 h (2), mild deficit, 24 h (3), and severe deficit 24 h (4). After computing these contrasts, resulting p-values of each contrast were corrected for multiplicity using the FDR-approach of Benjamini and Yekutieli . Additionally, a two-way analysis of variance (ANOVA) was performed with treatment main effect, time main effect and treatment by time interaction according to the empirical Bayes approach of Smyth  implemented in the R-package limma. Computed p-values were corrected for multiplicity using the positive false discovery rate [pFDR; .
Gene Ontology (GO) and metabolic pathway analyses
GO functional categories were assigned to differentially expressed genes using the web-based agriGO software . Enriched categories were computed using singular enrichment analysis (SEA) by comparing the list of differentially expressed genes to all expressed genes as described in Du, et al. . Multiple comparison correction  was performed and FDR controlled at 5%. Through the SEACOMPARE tool analysis results were combined for cross-comparisons. Differentially expressed genes were assigned to metabolic pathways and subsequently visualized using the Mapman software [31, 65] based on the functional annotation file ZmB73_5b_FGS_cds_2011 . Additionally, genes were assigned to biochemical pathways with the web-based Plant MetGenMAP [32, 66] software according to the MaizeCyc database [version 2.1; .
Quantitative real-time PCR (qRT-PCR) analysis
To confirm gene expression levels detected by RNA-Seq, quantitative real-time PCR was performed in a Bio-Rad CFX 384TM Real-Time System (Bio-Rad, Munich, Germany) using gene-specific oligonucleotides (Additional file 9). The cDNA for qRT-PCR analyses was synthesized from 1 μg total RNA with the qScript cDNA SuperMix (Quanta Biosciences, Gaithersburg, MD, USA) using the same RNA samples as for the cDNA library construction. Each PCR reaction contained 4 μl MESA Blue qPCR™ Mastermix Plus for SYBR Assay no ROX (Eurogentec, Cologne, Germany), 1 μl cDNA sample and 100 nM gene-specific oligonucleotide primers to a final volume of 8 μl. The primer efficiency of each oligonucleotide was calculated using the following dilution series: 1, 1/2, 1/4, 1/8, 1/16, 1/32, 1/64, and 1/128. The relative expression levels of the transcripts were calculated with reference to the housekeeping gene myosin (Genbank AC: 486090G09.x1). Significant differences in gene expression levels were determined by a two-sided Student’s t-test.
We thank Jörg Hagmann (Max-Planck-Institute for Developmental Biology, Tübingen) for initial data analysis and Heiner Goldbach and Sven Berkau (University of Bonn, Bonn) for support with measuring water potentials.
This project received funding from the European Union’s Seventh Framework Programme (FP7/2007-2013) under grant agreement n° 289300.27.
- Wheeler T, von Braun J: Climate change impacts on global food security. Science. 2013, 341 (6145): 508-513. 10.1126/science.1239402.PubMedView ArticleGoogle Scholar
- Tester M, Langridge P: Breeding technologies to increase crop production in a changing world. Science. 2010, 327 (5967): 818-822. 10.1126/science.1183700.PubMedView ArticleGoogle Scholar
- Chandler VL, Brendel V: The maize genome sequencing project. Plant Physiol. 2002, 130 (4): 1594-1597. 10.1104/pp.015594.PubMed CentralPubMedView ArticleGoogle Scholar
- International Grains Council. [http://www.igc.int]
- Lobell DB, Schlenker W, Costa-Roberts J: Climate trends and global crop production since 1980. Science. 2011, 333 (6042): 616-620. 10.1126/science.1204531.PubMedView ArticleGoogle Scholar
- IPCC: Summary for Policymakers. The Physical Science Basis Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change. Edited by: Solomon S, Qin D, Manning M, Chen Z, Marquis M, Averyt KB, Tignor M, Miller HL. 2007, Cambridge, United Kingdom and New York, NY, USA: Cambridge University PressGoogle Scholar
- Boyer JS: Plant productivity and environment. Science. 1982, 218 (4571): 443-448. 10.1126/science.218.4571.443.PubMedView ArticleGoogle Scholar
- Nafziger ED, Carter PR, Graham EE: Response of corn to uneven emergence. Crop Sci. 1991, 31 (3): 811-815. 10.2135/cropsci1991.0011183X003100030053x.View ArticleGoogle Scholar
- Sharp RE, Silk WK, Hsiao TC: Growth of the maize primary root at low water potentials: I. Spatial distribution of expansive growth. Plant Physiol. 1988, 87 (1): 50-57. 10.1104/pp.87.1.50.PubMed CentralPubMedView ArticleGoogle Scholar
- Yamaguchi M, Sharp RE: Complexity and coordination of root growth at low water potentials: recent advances from transcriptomic and proteomic analyses. Plant Cell Environ. 2010, 33 (4): 590-603. 10.1111/j.1365-3040.2009.02064.x.PubMedView ArticleGoogle Scholar
- Bray EA: Plant responses to water deficit. Trends Plant Sci. 1997, 2 (2): 48-54. 10.1016/S1360-1385(97)82562-9.View ArticleGoogle Scholar
- Krasensky J, Jonak C: Drought, salt, and temperature stress-induced metabolic rearrangements and regulatory networks. J Exp Bot. 2012, 63 (4): 1593-1608. 10.1093/jxb/err460.PubMed CentralPubMedView ArticleGoogle Scholar
- Seki M, Narusaka M, Ishida J, Nanjo T, Fujita M, Oono Y, Kamiya A, Nakajima M, Enju A, Sakurai T, Satou M, Aklyama K, Taji T, Yamaguchi-Shinozaki K, Carninci P, Kawai J, Hayashizaki Y, Shinozaki K: Monitoring the expression profiles of 7000 Arabidopsis genes under drought, cold and high-salinity stresses using a full-length cDNA microarray. Plant J. 2002, 31 (3): 279-292. 10.1046/j.1365-313X.2002.01359.x.PubMedView ArticleGoogle Scholar
- Zheng J, Zhao J, Tao Y, Wang J, Liu Y, Fu J, Jin Y, Gao P, Zhang J, Bai Y, Wang G: Isolation and analysis of water stress-induced genes in maize seedlings by subtractive PCR and cDNA macroarray. Plant Mol Biol. 2004, 55 (6): 807-823. 10.1007/s11103-005-1969-9.PubMedView ArticleGoogle Scholar
- Humbert S, Subedi S, Cohn J, Zeng B, Bi Y-M, Chen X, Zhu T, McNicholas P, Rothstein S: Genome-wide expression profiling of maize in response to individual and combined water and nitrogen stresses. BMC Genomics. 2013, 14 (1): 3-10.1186/1471-2164-14-3.PubMed CentralPubMedView ArticleGoogle Scholar
- Jia J, Fu J, Zheng J, Zhou X, Huai J, Wang J, Wang M, Zhang Y, Chen X, Zhang J, Zhao J, Su Z, Lv Y, Wang G: Annotation and expression profile analysis of 2073 full-length cDNAs from stress-induced maize (Zea mays L.) seedlings. Plant J. 2006, 48 (5): 710-727. 10.1111/j.1365-313X.2006.02905.x.PubMedView ArticleGoogle Scholar
- Poroyko V, Spollen W, Hejlek L, Hernandez A, LeNoble M, Davis G, Nguyen H, Springer G, Sharp R, Bohnert H: Comparing regional transcript profiles from maize primary roots under well-watered and low water potential conditions. J Exp Bot. 2007, 58 (2): 279-289.PubMedView ArticleGoogle Scholar
- Spollen W, Tao W, Valliyodan B, Chen K, Hejlek L, Kim J-J, LeNoble M, Zhu J, Bohnert H, Henderson D, Schachtman DP, Davis GE, Springer GK, Sharp RE, Nguyen HT: Spatial distribution of transcript changes in the maize primary root elongation zone at low water potential. BMC Plant Biol. 2008, 8 (1): 32-10.1186/1471-2229-8-32.PubMed CentralPubMedView ArticleGoogle Scholar
- Yue G, Zhuang Y, Li Z, Sun L, Zhang J: Differential gene expression analysis of maize leaf at heading stage in response to water-deficit stress. Biosience Rep. 2008, 028 (3): 125-134. 10.1042/BSR20070023.View ArticleGoogle Scholar
- Metzker ML: Sequencing technologies - the next generation. Nat Rev Genet. 2010, 11 (1): 31-46. 10.1038/nrg2626.PubMedView ArticleGoogle Scholar
- MaizeSequence.org FTP site. [ftp://ftp.gramene.org/pub/gramene/maizesequence.org/]Google Scholar
- Schnable PS, Ware D, Fulton RS, Stein JC, Wei F, Pasternak S, Liang C, Zhang J, Fulton L, Graves TA, Minx P, Reily AD, Courtney L, Kruchowski SS, Tomlinson C, Strong C, Delehaunty K, Fronick C, Courtney B, Rock SM, Belter E, Du F, Kim K, Abbott RM, Cotton M, Levy A, Marchetto P, Ochoa K, Jackson SM, Gillam B: The B73 maize genome: complexity, diversity, and dynamics. Science. 2009, 326 (5956): 1112-1115. 10.1126/science.1178534.PubMedView ArticleGoogle Scholar
- Agarwal A, Koppstein D, Rozowsky J, Sboner A, Habegger L, Hillier L, Sasidharan R, Reinke V, Waterston R, Gerstein M: Comparison and calibration of transcriptome data from RNA-seq and tiling arrays. BMC Genomics. 2010, 11 (1): 383-10.1186/1471-2164-11-383.PubMed CentralPubMedView ArticleGoogle Scholar
- Marioni JC, Mason CE, Mane SM, Stephens M, Gilad Y: RNA-seq: an assessment of technical reproducibility and comparison with gene expression arrays. Genome Res. 2008, 18 (9): 1509-1517. 10.1101/gr.079558.108.PubMed CentralPubMedView ArticleGoogle Scholar
- Sharp RE, Poroyko V, Hejlek LG, Spollen WG, Springer GK, Bohnert HJ, Nguyen HT: Root growth maintenance during water deficits: physiology to functional genomics. J Exp Bot. 2004, 55 (407): 2343-2351. 10.1093/jxb/erh276.PubMedView ArticleGoogle Scholar
- Westgate ME, Boyer JS: Osmotic adjustment and the inhibition of leaf, root, stem and silk growth at low water potentials in maize. Planta. 1985, 164 (4): 540-549. 10.1007/BF00395973.PubMedView ArticleGoogle Scholar
- Zinselmeier C, Sun Y, Helentjaris T, Beatty M, Yang S, Smith H, Habben J: The use of gene expression profiling to dissect the stress sensitivity of reproductive development in maize. Field Crops Res. 2002, 75 (2–3): 111-121.View ArticleGoogle Scholar
- Bruce WB, Edmeades GO, Barker TC: Molecular and physiological approaches to maize improvement for drought tolerance. J Exp Bot. 2002, 53 (366): 13-25. 10.1093/jexbot/53.366.13.PubMedView ArticleGoogle Scholar
- NCBI Sequence Read Archive (SRA). [http://www.ncbi.nlm.nih.gov/sra]
- Nettleton D: A discussion of statistical methods for design and analysis of microarray experiments for plant scientists. Plant Cell. 2006, 18 (9): 2112-2121. 10.1105/tpc.106.041616.PubMed CentralPubMedView ArticleGoogle Scholar
- Thimm O, Bläsing O, Gibon Y, Nagel A, Meyer S, Krüger P, Selbig J, Müller LA, Rhee SY, Stitt M: Mapman: a user-driven tool to display genomics data sets onto diagrams of metabolic pathways and other biological processes. Plant J. 2004, 37 (6): 914-939. 10.1111/j.1365-313X.2004.02016.x.PubMedView ArticleGoogle Scholar
- Plant MetGenMAP. [http://bioinfo.bti.cornell.edu/cgi-bin/MetGenMAP/home.cgi]
- agriGO: GO analysis toolkit and database for the agricultural community. [http://bioinfo.cau.edu.cn/agriGO/analysis.php]
- Sinclair TR: Challenges in breeding for yield increase for drought. Trends Plant Sci. 2011, 16 (6): 289-293. 10.1016/j.tplants.2011.02.008.PubMedView ArticleGoogle Scholar
- Paschold A, Jia Y, Marcon C, Lund S, Larson NB, Yeh C-T, Ossowski S, Lanz C, Nettleton D, Schnable PS, Hochholdinger F: Complementation contributes to transcriptome complexity in maize (Zea mays L.) hybrids relative to their inbred parents. Genome Res. 2012, 22 (12): 2445-2454. 10.1101/gr.138461.112.PubMed CentralPubMedView ArticleGoogle Scholar
- Ozturk ZN, Talamé V, Deyholos M, Michalowski C, Galbraith D, Gozukirmizi N, Tuberosa R, Bohnert H: Monitoring large-scale changes in transcript abundance in drought- and salt-stressed barley. Plant Mol Biol. 2002, 48 (5–6): 551-573.View ArticleGoogle Scholar
- Kakumanu A, Ambavaram MMR, Klumas C, Krishnan A, Batlang U, Myers E, Grene R, Pereira A: Effects of drought on gene expression in maize reproductive and leaf meristem tissue revealed by RNA-seq. Plant Physiol. 2012, 160 (2): 846-867. 10.1104/pp.112.200444.PubMed CentralPubMedView ArticleGoogle Scholar
- Andjelkovic V, Thompson R: Changes in gene expression in maize kernel in response to water and salt stress. Plant Cell Rep. 2006, 25 (1): 71-79. 10.1007/s00299-005-0037-x.PubMedView ArticleGoogle Scholar
- Fan X-D, Wang J-Q, Yang N, Dong Y-Y, Liu L, Wang F-W, Wang N, Chen H, Liu W-C, Sun Y-P, Wu J-Y, Li H-Y: Gene expression profiling of soybean leaves and roots under salt, saline–alkali and drought stress by high-throughput Illumina sequencing. Gene. 2013, 512 (2): 392-402. 10.1016/j.gene.2012.09.100.PubMedView ArticleGoogle Scholar
- Zhu J, Alvarez S, Marsh EL, LeNoble ME, Cho I-J, Sivaguru M, Chen S, Nguyen HT, Wu Y, Schachtman DP, Sharp RE: Cell wall proteome in the maize primary root elongation zone. II Region-specific changes in water soluble and lightly ionically bound proteins under water deficit. Plant Physiol. 2007, 145 (4): 1533-1548. 10.1104/pp.107.107250.PubMed CentralPubMedView ArticleGoogle Scholar
- Mittler R, Vanderauwera S, Gollery M, Van Breusegem F: Reactive oxygen gene network of plants. Trends Plant Sci. 2004, 9 (10): 490-498. 10.1016/j.tplants.2004.08.009.PubMedView ArticleGoogle Scholar
- Asada K, Takahashi M: Production and scavenging of active oxygen in photosynthesis. Photoinhibition. Edited by: Kyle DJ, Osborne CB, Arntzen CJ. 1987, Amsterdam: Elsevier, 227-287.Google Scholar
- Mittler R, Vanderauwera S, Suzuki N, Miller G, Tognetti VB, Vandepoele K, Gollery M, Shulaev V, Van Breusegem F: ROS signaling: the new wave?. Trends Plant Sci. 2011, 16 (6): 300-309. 10.1016/j.tplants.2011.03.007.PubMedView ArticleGoogle Scholar
- Kar RK: Plant responses to water stress: role of reactive oxygen species. Plant Signal Behav. 2011, 6 (11): 1741-1745. 10.4161/psb.6.11.17729.PubMed CentralPubMedView ArticleGoogle Scholar
- Voothuluru P, Sharp RE: Apoplastic hydrogen peroxide in the growth zone of the maize primary root under water stress. I. Increased levels are specific to the apical region of growth maintenance. J Exp Bot. 2013, 64 (5): 1223-1233. 10.1093/jxb/ers277.PubMedView ArticleGoogle Scholar
- Price AH, Hendry GAF: Iron-catalysed oxygen radical formation and its possible contribution to drought damage in nine native grasses and three cereals. Plant Cell Environ. 1991, 14 (5): 477-484. 10.1111/j.1365-3040.1991.tb01517.x.View ArticleGoogle Scholar
- Muller B, Pantin F, Génard M, Turc O, Freixes S, Piques M, Gibon Y: Water deficits uncouple growth from photosynthesis, increase C content, and modify the relationships between C and growth in sink organs. J Exp Bot. 2011, 62 (6): 1715-1729. 10.1093/jxb/erq438.PubMedView ArticleGoogle Scholar
- Sharp RE, Hsiao TC, Silk WK: Growth of the maize primary root at Low water potentials: II. Role of growth and deposition of hexose and potassium in osmotic adjustment. Plant Physiol. 1990, 93 (4): 1337-1346. 10.1104/pp.93.4.1337.PubMed CentralPubMedView ArticleGoogle Scholar
- Voetberg GS, Sharp RE: Growth of the maize primary root at Low water potentials: III. Role of increased proline deposition in osmotic adjustment. Plant Physiol. 1991, 96 (4): 1125-1130. 10.1104/pp.96.4.1125.PubMed CentralPubMedView ArticleGoogle Scholar
- Hochholdinger F, Woll K, Guo L, Schnable PS: The accumulation of abundant soluble proteins changes early in the development of the primary roots of maize (Zea mays L.). Proteomics. 2005, 5 (18): 4885-4893. 10.1002/pmic.200402034.PubMedView ArticleGoogle Scholar
- Byrne ME: A role for the ribosome in development. Trends Plant Sci. 2009, 14 (9): 512-519. 10.1016/j.tplants.2009.06.009.PubMedView ArticleGoogle Scholar
- Sormani R, Masclaux-Daubresse C, Daniele-Vedele F, Chardon F: Transcriptional regulation of ribosome components are determined by stress according to cellular compartments in Arabidopsis thaliana. PLoS ONE. 2011, 6 (12): e28070-10.1371/journal.pone.0028070.PubMed CentralPubMedView ArticleGoogle Scholar
- Wang J, Lan P, Gao H, Zheng L, Li W, Schmidt W: Expression changes of ribosomal proteins in phosphate- and iron-deficient Arabidopsis roots predict stress-specific alterations in ribosome composition. BMC Genomics. 2013, 14 (1): 783-10.1186/1471-2164-14-783.PubMed CentralPubMedView ArticleGoogle Scholar
- Bailey-Serres J, Freeling M: Hypoxic stress-induced changes in ribosomes of maize seedling roots. Plant Physiol. 1990, 94 (3): 1237-1243. 10.1104/pp.94.3.1237.PubMed CentralPubMedView ArticleGoogle Scholar
- Ludwig Y, Zhang Y, Hochholdinger F: The maize (Zea Mays L.) AUXIN/INDOLE-3-ACETIC ACID gene family: phylogeny, synteny, and unique root-type and tissue-specific expression patterns during development. PLoS ONE. 2013, 8 (11): e78859-10.1371/journal.pone.0078859.PubMed CentralPubMedView ArticleGoogle Scholar
- Winz RA, Baldwin IT: Molecular interactions between the specialist herbivore Manduca sexta (Lepidoptera, Sphingidae) and its natural host Nicotiana attenuata. IV. Insect-induced ethylene reduces jasmonate-induced nicotine accumulation by regulating Putrescine N-Methyltransferase transcripts. Plant Physiol. 2001, 125 (4): 2189-2202. 10.1104/pp.125.4.2189.PubMed CentralPubMedView ArticleGoogle Scholar
- Schroeder A, Mueller O, Stocker S, Salowsky R, Leiber M, Gassmann M, Lightfoot S, Menzel W, Granzow M, Ragg T: The RIN: An RNA integrity number for assigning integrity values to RNA measurements. BMC Mol Biol. 2006, 7: 3-10.1186/1471-2199-7-3.PubMed CentralPubMedView ArticleGoogle Scholar
- Ewing B, Green P: Base-calling of automated sequencer traces using Phred. II. Error probabilities. Genome Res. 1998, 8 (3): 186-194.PubMedView ArticleGoogle Scholar
- Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B: Mapping and quantifying mammalian transcriptomes by RNA-seq. Nat Methods. 2008, 5 (7): 621-628. 10.1038/nmeth.1226.PubMedView ArticleGoogle Scholar
- Smyth GK: Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004, 3 (1): Article3-PubMedGoogle Scholar
- Law C: Precision Weights for Gene Expression Analysis. 2013, University of Melbourne: PhD thesisGoogle Scholar
- Benjamini Y, Yekutieli D: The control of the false discovery rate in multiple testing under dependency. Ann Stat. 2001, 29 (4): 1165-1188. 10.1214/aos/1013699998.View ArticleGoogle Scholar
- Storey JD: A direct approach to false discovery rates. J Roy Stat Soc Ser B (Stat Method). 2002, 64 (3): 479-498. 10.1111/1467-9868.00346.View ArticleGoogle Scholar
- Du Z, Zhou X, Ling Y, Zhang Z, Su Z: agriGO: a GO analysis toolkit for the agricultural community. Nucleic Acids Res. 2010, 38: W64-W70. 10.1093/nar/gkq310.PubMed CentralPubMedView ArticleGoogle Scholar
- MapMan. [http://mapman.gabipd.org/web/guest/mapman]
- Joung J-G, Corbett AM, Fellman SM, Tieman DM, Klee HJ, Giovannoni JJ, Fei Z: Plant MetGenMAP: an integrative analysis system for plant systems biology. Plant Physiol. 2009, 151 (4): 1758-1768. 10.1104/pp.109.145169.PubMed CentralPubMedView ArticleGoogle Scholar
- Metabolic pathways in maize: MaizeCyc. [http://pathway.gramene.org/gramene/maizecyc.shtml]
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.