Pankaj Barah
Naresh D. Jayavelu
John Mundy
Atle M. Bones- 1Cell Molecular Biology and Genomics Group, Department of Biology, Norwegian University of Science and Technology, Trondheim, Norway
- 2Department of Chemical Engineering, Norwegian University of Science and Technology, Trondheim, Norway
- 3Department of Biology, University of Copenhagen, Copenhagen, Denmark
In the scenario of global warming and climate change, heat stress is a serious threat to crop production worldwide. Being sessile, plants cannot escape from heat. Plants have developed various adaptive mechanisms to survive heat stress. Several studies have focused on diversity of heat tolerance levels in divergent Arabidopsis thaliana (A. thaliana) ecotypes, but comprehensive genome scale understanding of heat stress response in plants is still lacking. Here we report the genome scale transcript responses to heat stress of 10 A. thaliana ecotypes (Col, Ler, C24, Cvi, Kas1, An1, Sha, Kyo2, Eri, and Kond) originated from different geographical locations. During the experiment, A. thaliana plants were subjected to heat stress (38°C) and transcript responses were monitored using Arabidopsis NimbleGen ATH6 microarrays. The responses of A. thaliana ecotypes exhibited considerable variation in the transcript abundance levels. In total, 3644 transcripts were significantly heat regulated (p < 0.01) in the 10 ecotypes, including 244 transcription factors and 203 transposable elements. By employing a systems genetics approach- Network Component Analysis (NCA), we have constructed an in silico transcript regulatory network model for 35 heat responsive transcription factors during cellular responses to heat stress in A. thaliana. The computed activities of the 35 transcription factors showed ecotype specific responses to the heat treatment.
Introduction
Climate change is increasingly viewed as a current and future cause of hunger and poverty (Lobell et al., 2011; Wheeler and von Braun, 2013). In the scenario of global climatic change, different environmental stresses are severe threats to agricultural production worldwide (Brown and Funk, 2008; Ahuja et al., 2010). Among all stress conditions, elevated temperature is seen as the most serious threat to crop production (Wheeler et al., 2000; Ciais et al., 2005; Semenov and Shewry, 2011). Recurrent heat stress also affects disease resistance in plants by suppressing plant immunity, as plant heat stress and defense responses share important mediators such as calcium ions and heat shock proteins (HSPs) (Lee et al., 2012). Climate data suggest that heat waves became more common during the twentieth century (Stott et al., 2004). Recently, Bita et al. reviewed the effects of high temperature stress on physiology, biochemistry, and gene regulation pathways in plants leading to catastrophic loss of crop productivity (Bita and Gerats, 2013). Transient or continuous high temperatures cause a range of morphological, physiological, and biochemical changes in plants which affect growth and development and may lead to a drastic reduction in economic yield (Richter et al., 2010). Plants are highly sensitive to temperature and can differentiate minute variations of as little as 1°C (Mittler et al., 2012). Upon exposure to heat stress, seed germination, and photosynthetic efficiency decline (Endo et al., 2009). Considering the predicted severity of changing climatic situation, dissecting the molecular basis of heat stress responses in plants, and identifying key components of the heat stress sensing and signal transduction pathways, are becoming major concern of present time (Bita and Gerats, 2013; Qu et al., 2013). Such knowledge could be used toward developing plants and crops with enhanced tolerance to heat stress (Zhang et al., 2000; Mittler and Blumwald, 2010).
Environmental stress is a key factor driving the genome regulation, evolutionary history, and geographical distribution of organisms including plants (Alonso-Blanco et al., 2009). Intraspecific natural variation or within-species phenotypic variation caused by spontaneous, favorable mutations contribute to the local adaptations of plants (Mitchell-Olds and Schmitt, 2006). Such natural variation in crop plants has been exploited by human society for the selection of developmental traits and physiological features beneficial for agriculture (Weigel and Nordborg, 2005; Doebley et al., 2006). Additionally, studying natural variation in wild species can tell us about the molecular basis of phenotypic differences related to plant adaptation to diverse natural environments (Borevitz and Nordborg, 2003). There have been very few studies conducted till date focusing on the diversity of heat tolerance in phenotypically divergent ecotypes (Alonso-Blanco and Koornneef, 2000; Larkindale et al., 2005; Al-Quraan et al., 2012). Thus, the molecular basis of the natural variation during heat stress response in plants at genome scale is not fully understood yet (Yeh et al., 2012).
Transcriptomics, proteomics and metabolomics approaches have been frequently used to identify heat stress-responsive genes, proteins, and metabolites in plants (Kaplan et al., 2004; Jagadish et al., 2010; Pecinka et al., 2010; Weston et al., 2011; Zou et al., 2011b; Rocco et al., 2013). Transcript profiling is a major tool to identify genes exhibiting transcriptional regulation in response to changing environmental conditions. For such studies in plants, A. thaliana remains a model system (Somerville and Koornneef, 2002). Variation in experimental conditions and protocols makes it difficult to extract and compare information from data sets produced by individual laboratories (Moreau et al., 2003). To overcome such problems, 10 ecotypes of A. thaliana were subjected to 5 individual stress treatments and 6 combinations of these stress treatments under the same experimental set up and profiling protocols (Rasmussen et al., 2013). We have considered all the heat experiments conducted on 10 ecotypes from this published dataset (GEO accessionGSE41935) to explore genome-scale transcriptomic response signatures of A. thaliana during heat stress treatment. Being highly dynamic in nature, any biological system changes in response to environmental and genetic perturbations. Differential dynamic network mapping facilitates the exploration of previously unknown interactions (Ideker and Krogan, 2012). While the A. thaliana genome has ~1922 TFs (Guo et al., 2005), experimentally confirmed regulatory relations are available for less than 100 TFs, as per information from the AGRIS database version updated in September, 2012 (Davuluri et al., 2003). Tirosh et al. (Tirosh and Barkai, 2011) have explained how regulatory relationships can also be deduced from patterns of evolutionary divergence in molecular properties such as gene expression (Keurentjes et al., 2007). To compensate the lack of information on transcription factor activity at the genome-scale, computational algorithms have been developed to identify regulatory modules and their condition-specific regulators from gene expression data (Alter et al., 2000; Segal et al., 2003; Herrgard et al., 2004). Network Component Analysis (NCA) is such an approach, which has been successfully implemented in species including A. thaliana to determine both the activities and regulatory influences for a set of transcription factors on target genes (Liao et al., 2003; Kao et al., 2004; Wang et al., 2011). Using the NCA method, we have predicted ecotype specific regulatory relationships which generated new information toward understanding the natural variation in heat response pattern among different ecotypes of the model plant A. thaliana.
Results
Different Transcriptome Signatures of 10 Arabidopsis ecotypes responding to heat stress
To cover a wide array of phenotypic variations, 10 natural accessions of A. thaliana representing their originally reported habitats from 16 to 56.5° north latitudes were selected during the ERA-PG Multi-stress project. These accessions or ecotypes were- Cvi (Cape Verde Islands), Kas-1 (Kashmir, India), Kyo-2 (Kyoto, Japan), Sha (Shakdara, Tadjikistan), Col-0(Columbia, USA), Kond (Kondara, Tadjikistan), C24 (Coimbra, Portugal), Ler (Landsberg, Poland), An-1 (Antwerpe, Belgium), Eri (Erigsboda, Sweden) (details in Table 1). We chose a cut-off p ≤ 0.01 to define a gene as differentially stress regulated. Using the results from the 10 ecotypes, we examined the differences in transcript abundences that occurred during early hours of heat treatment (38°C). The results (Table 1 and Figure 1) indicated that the A. thaliana ecotypes have visibly different transcriptome level signatures in response to heat stress. Variable numbers of transcripts were up or down regulated among the ecotypes (Table 1). Kas-1 (797) and Cvi (776) exhibited higher numbers of differentially regulated transcripts while Col-0 (143) had comparatively few differentially regulated transcripts. A unified list of 3644 differentially regulated transcripts (p < 0.01) was generated from the 10 ecotypes (Table S1A.) Surprisingly, 3114 (85%) transcripts were differentially regulated in only one of the 10 ecotypes. Figure 2 displays fold change values (treatment vs. control) calculated from normalized expression index for the top 1000 most significant genes from the 10 ecotypes. Global observation of the heat map indicates differentially regulated transcriptome signatures in response to heat treatment in the 10 ecotypes. The significant list of differentially regulated transcripts includes most of the previously documented heat regulated genes including Hsps (heat shock proteins) and Hsfs (heat shock transcription factors) (Swindell et al., 2007).
Table 1. Summary of the ecotypes and their gene expression pattern during heat stress.
Figure 1. Numbers of differentially regulated transcripts in each of the 10 ecotypes at significance level p ≤ 0.01. Ecotypes are on the x axis and numbers of differentially regulated transcripts on the y axis. Blue bar represents total number of differentially regulated transcripts, red bar the number of positively regulated (up) transcripts and green bar represents number of negatively regulated (down) transcripts.
Figure 2. Fold change values (treatment vs. control) calculated from normalized expression index for top 1000 significant genes from all the 10 ecotypes. Hierarchical clustering (HCL) was performed with Pearson correlation using average linkage method and 10,000 bootstrapping for the top 1000 heat regulated transcripts based on fold-change ratios compared to their respective controls.
Ecotype Specific Heat Regulated Transcript Lists Contain Many Transcription Factors (TFs) and Transposable Elements (TEs)
The unified list of 3644 differentially regulated transcripts during the heat stress contained 244 TFs (annotated in Table S1B). Only AT5G57660 (CONSTANS-like 5 zinc finger family protein) was significantly (p ≤ 0.01) upregulated in all of the 10 ecotypes. Two other TFs, AT4G25480 (Dehydration response element B1A) and AT5G24470 (Arabidopsis pseudo-response regulator 5), were significantly upregulated in 9 ecotypes. MBF1C/AT3G24500 (multiprotein bridging factor 1C) was significantly down-regulated in 8 ecotypes (Suzuki et al., 2008). Among others, 70 TFs were significantly regulated in 2 ecotypes and 62, TFs were significantly regulated only in one of the ecotypes. The differentially regulated TFs included members of prominent TF families such as ABF3,ADOF, AFO, AGL, NAC, AP1, AP2, Prr5, ARF, bZIP, HSF, IDD, MYB, BLJ, DNAJ, JAZ, MYB, PHD finger, WRKY, C2H2 zinc finger etc. Such differential regulation of diverse TF families was obvious from the fact that apart from heat shock protein induction, other pathways involving ethylene, salicylic acid (SA), and trehalose were shown to play crucial roles in plant thermotolerance (Larkindale and Knight, 2002; Larkindale et al., 2005).
The Nimblgen12-plex Arabidopsis microarray chip included 3822 transposable element (TE) probes. Of them, 203 TEs were differentially regulated during heat stress (Table S1C). Except for 5, TEs, the rest were differentially regulated in single ecotypes. The distribution of the differentially regulated TEs in ten ecotypes were: Col-0 (10), Ler (24), Cvi (24), Eri (18), Kas2 (27), Kond (23), Kyo2 (30), C24 (11), Sha (27), and An1 (14).
Gene Set Enrichment Analysis (GSEA) Indicates Activation of Diverse Processes
To investigate functionally over-represented gene ontology categories, BinGO software was used on the list of 3644 differentially regulated transcripts from the 10 ecotypes. No annotations were retrieved for 60 genes which were eliminated from the final analysis. In total, 82 statistically significant gene ontology categories were detected, including many parent categories such as response to stimulus, stress, biotic stimulus, abiotic stimulus etc. (Table S2). Apart from these global terms, genes showing significant variation in mRNA levels in A. thaliana during heat stress were mainly belong to categories like response to heat, temperature stimulus, water deprivation, light stimulus, wounding, osmotic stress, oxidative stress, salt stress, and protein folding etc. The rest of the differentially regulated genes covered various functions, such as transcription, translation, signaling, metabolism, and general stress response. These results indicated that, during exposure to heat stress, plants extensively reprogrammed gene expression, to limit damage caused by high temperatures.
Hsp Genes Exhibit Differential Expression Patterns in Arabidopsis Ecotypes during Heat Stress
A list of total 145 Hsps was generated containing the term “heat shock protein” as per annotations available from TAIR10 database (Table S3). Among them, 31 Hsps were significantly (p = 0.01) differentially regulated in at least one of the 10 ecotypes. Most of them were encoded for DNAJ heat shock N-terminal domain-containing proteins. Other upregulated members were HSP70, HSP21, HSP17, HSP18 etc. None of the 31 significant HSPs were expressed in a similar pattern across all 10 ecotypes, which indicated differentially regulated activity profiles of them across A. thaliana ecotypes during heat stress responses (Figure 3).
Figure 3. Heat map of log2 fold change values of the 31 significantly regulated (p ≤ 0.01) HSP genes in 10 A. thaliana ecotypes. Genes and ecotypes were clustered using Pearson's correlation coefficient with average linkage method. The P-values and log2 fold-change values associated with all of the 145 Hsps are provided in Table S3.
Re-Construction of a Transcriptional Regulatory Network During the Heat Stress Response in A. thaliana
By looking at the differential expression levels of a large number of TFs during the heat experiments, we wanted to explore the pattern of regulatory interactions between the TFs and their corresponding target genes (TGs) in the 10 A. thaliana ecotypes using a benchmarked algorithm, NCA (Liao et al., 2003; Wang et al., 2011; Barah et al., 2013). Simple correlation between the expression profile of a transcription factor and its targets is not obvious, and simple clustering based methods have not been very successful in deciphering them (Qian et al., 2003). The key assumption during predictions of interactions between TFs and their target genes using gene expression data is that high dimensional mRNA expression profiles contain hidden regulatory signals which can be decomposed to low-dimensional regulatory signals driven through an interacting network (Holter et al., 2000; Carrera et al., 2009). The lower dimensional regulatory signals can be represented as a bipartite networked system of the transcription factors and the target genes in which the gene expression levels are transformed into weighted functions of the intracellular states corresponding to the activity of the transcription factors.
The NCA algorithm requires two inputs to calculate the hidden regulatory activity profiles: a series of gene expression profiles and a pre-defined regulatory network. A list of 1922, TFs in A. thaliana genome were collected from the Database of Arabidopsis Transcription Factors (DATF) (Guo et al., 2005), The Arabidopsis Gene Regulatory Information Server (AGRIS) (Yilmaz et al., 2011), and the Plant Transcription Factor Database (PlantTFDB) (Riano-Pachon et al., 2007). A list of 59 previously known heat regulated transcription factors was generated from the Gene Ontology database (Ashburner et al., 2000) under the annotation category “response to heat” or containing the term “heat shock factor.” The list of differentially regulated TFs in our transcriptome data contains 35 out the curated list of 59 heat responsive factors. A bipartite co-regulatory network (Alvarez and Woolf, 2011) was constructed from the gene expression values based on correlation-coefficient threshold ≥ 0.8 between the 35 heat regulated TFs and 1294, TGs (Table S4). The resulting network contained 1947 connections. Of them, 687 connections were activations (positive) and 1260 were repressions (negative). Few of the TFs in the network are highly connected (hubs), which supports the scale-free behavior of the predicted TF-TG network (Albert, 2005). The number of connections for each of the TFs is listed in Table 2. This co-regulatory network model was further used as an input to the NCA algorithm to predict the activities of the TFs based on differential expression profiles (log2 fold change values) of their linked TGs (Figure 4). Noticeable variation was observed in the activity profiles of the 35 TFs among the 10 ecotypes.
Table 2. Number of predicted regulatory connections for each of the TFs.
Figure 4. Predicted the activities of the 35 TFs among the 10 ecotypes. The NCA algorithm predicts the activities of the TFs based on differential expression profiles (log2 fold change values) of their linked TGs. The predicted activity profiles of the 35 heat regulated TFs shows variations in the 10 A. thaliana ecotypes.
The predicted activity profiles of the 35 heat regulated TFs clearly show the ecotype specific activities in the 10 A. thaliana ecotypes. For example, transcription factor AT5G02810 (PRR7) was highly active in the Kond ecotype. We identified both multi responsive (active in more than one ecotype) and unique responsive transcription factors (active only in one specific ecotype). The detailed results are provided in Table 3. The majority of the ecotype specific transcription factors were active in Cvi ecotype in response to heat treatment. Multi responsive transcription factors are mostly active in Kond, An-1 and Sha. The transcription factor AT1G74950 (TIFY10B) is highly responsive in Kond, An-1 and Sha.
Table 3. Ecotype specific transcriptional activity profiles of the 35 heat responsive TFs.
Discussion
Here we undertook an experiment to analyze the natural variation in genome-scale heat stress response in 10 A. thaliana ecotypes at a single time point (3 h) of gene expression measurement. The analysis indicated that the 10 A. thaliana ecotypes had significantly different transcriptome level signatures in response to heat stress. It raises question about global acceptability of results generated from previously conducted plant stress experiments based only on Col-0 and Ler as model ecotypes.
Among the differentially heat regulated transcripts, 85% showed ecotype specific expression patterns. Heat shock proteins or molecular chaperones were the most prominent group within the list of differentially regulated list of transcripts. Apart from common, heat stress related functional categories, GSEA of the differentially regulated transcripts uncovered many functional categories related to other stress response. Profound transcriptional reprogramming during heat stress involves extensive regulation of transcription that affects a large part of the whole transcriptome (Zeller et al., 2009; Zou et al., 2011a).
The differential expression of 243 TFs among the 10 ecotypes indicates a complex level of transcriptional regulation during the exposure of plants to heat stress. Due to the lack of experimentally validated transcriptional regulatory information in A. thaliana, an in-silico transcript regulatory network model during cellular responses to heat stress in A. thaliana was constructed using the homogeneous gene expression data. The predicted activities of the heat regulated TFs showed significant variations among 10 ecotypes (Figure 4). The observed differential regulatory activities among the heat regulated TFs might contribute to high temperature acclimation of the specific ecotypes. Swindell et al. (2007) reported that multiple stress treatments interact with HSF and HSP response pathways to varying extents, suggesting that there is a basis of cross-tolerance in plant species through a shock response network. Expression of HSPs confers heat stress tolerance in plants that leads to improved photosynthesis, assimilate partitioning, water and nutrient use efficiency, and membrane stability (Wahid et al., 2007). The function of HSPs in enhancing stress tolerance may vary among genotypes and also depends on the nature of the stress imposed upon the cell. Such quantitative variation in the gene expression among the Hsp genes in the 10 ecotypes is clearly visible from Figure 3. Heat stress leads to direct denaturation of cellular proteins. Earlier, some in vitro data indicated that HSPs acted as molecular chaperones to prevent thermal aggregation of proteins by binding non-native intermediates which could then be refolded in an ATP-dependent fashion by other chaperones (Lee and Vierling, 2000). Therefore, the molecular chaperone activity of the HSPs may contribute to high temperature tolerance via prevention of protein misfolding and removal of non-native aggregations.
The 203 differentially regulated transposable elements (TEs) among the 10 ecotypes may play an important role in genome adapting to local climatic temperatures (Fedoroff, 2012). In a recent review, (Lisch, 2013) summarize the impact of stress activated retrotransposons on genome evolution in plants. Natural populations can show diverse responses when exposed to adverse environmental conditions because of their genetic variation as well as because of their epigenetic variations. Only a few studies have reported that stress responses in plants affect epigenetic regulation and require specific epigenetic regulators (Chinnusamy and Zhu, 2009). For example, UV, cold, and heat stress result in the reactivation of silent transgenes and endogenous transposable elements, although without reductions in DNA methylation and repressive histone marks (Grativol et al., 2012; Popova et al., 2013). Pecinka et al. (2010) showed that several repetitive elements of A. thaliana are under epigenetic regulation by transcriptional gene silencing at ambient temperatures and become activated by prolonged heat stress. A change in the epigenetic state of TEs by heat stress might also contribute to regulatory activities for adjacent genes. Recently, Wang et al. (2013) demonstrated that both TE sequence polymorphisms and the presence of linked TEs are positively correlated with intraspecific variation in gene expression. Some of the differentially regulated TEs in our heat experiments may therefore, be potentially interesting targets to explore diversity of heat stress responses among different ecotypes. Further targeted experiments in this direction can explore the molecular details of any potential role of these TEs on genomic adaptation of the ecotypes to their local environment.
Materials and Methods
Microarray Data
We have considered all the heat stress microarray experiments conducted on 10 ecotypes during the ERA-PG Multi-stress project (Rasmussen et al., 2013) to explore genome-scale transcriptome response signatures of A. thaliana during heat stress (microarray data available at GEO with the accession GSE41935). All experiments of ERA-PG Multistress project were performed in environmentally controlled rooms at the plant growth facilities at RISØ DTU National Laboratory for Sustainable Energy (Roskilde, Denmark). A pilot study using wild type Col and Ler plants was set up to find the appropriate conditions at sub-lethal doses (Rasmussen et al., 2013). These initial observations indicated that an optimal time before the onset of a phenotypic response (e.g., wilting, dehydration) while avoiding tissue damage was 3 h. 10 A. thaliana wild ecotypes (Table 1) were grown in soil under long day photoperiod and 24°C in a greenhouse setting for one generation to amplify homogeneous seed for all different genotypes. The seeds were then sown into trays and grown in a Conviron growth chamber (Winnipeg, Manitoba, Canada) under a 12h/12h photoperiod, 24°C and standard A. thaliana growth conditions. 3 week-old plants were then placed for 3 h in the environmentally controlled growth rooms that were preset to heat stress conditions (38°C). Triplicated (biological) trays with the wild type controls were subject to the heat stress. After the stress treatments, leaf samples were collected and promptly frozen in liquid nitrogen for subsequent microarray experiments.
Statistical Analysis of the Data
Resulting data from the microarray experiments was pre-processed using the RMA (Irizarry et al., 2003) implementation in the oligo package (Carvalho and Irizarry, 2010) in R programming platform (R Core Team, 2012). Gene annotation was acquired from TAIR10 (Lamesch et al., 2012) using the BioMart data mining tool (Guberman et al., 2011). Differentially expressed genes between control and treated plants were identified using t-test (p < 0.01). Genotype specific responses to stress were identified by the interaction effect from a Two-Way ANOVA (Kerr et al., 2000; Cui and Churchill, 2003) of the genotype and treatment effect (p < 0.01). The union of stress responsive genes were further used for network-based analysis. Heat maps were plotted using TM4 microarray software suite (Saeed et al., 2006).
Gene Set Enrichment Analysis (GSEA)
The Biological Networks Gene Ontology (BiNGO) tool (Maere et al., 2005), an open-source Java tool, was used to determine Gene Ontology (GO) terms (Ashburner et al., 2000) that were significantly overrepresented in our differentially regulated gene lists (p-values were Bonferroni corrected).
Network Component Analysis and Network Reconstruction
Network component analysis is a computational method for reconstructing the hidden regulatory signals (TFAs-Transcription Factor Activities) from gene expression data with known connectivity in terms of matrix decomposition (Liao et al., 2003; Galbraith et al., 2006). The algorithm for NCA analysis is implemented in MATLAB by Liao and colleagues and is online for download (www.ee.ucla.edu/~riccardo/NCA/nca.html). With NCA as a reconstruction method, we predicted significant TFs and connectivity strength on target genes and TFAs of TFs.
Author Contributions
Atle M. Bones and John Mundy conceived the Multi-Stress project. Pankaj Barah developed the concept of the current manuscript, performed bioinformatics analyses and drafted the manuscript. Naresh D. Jayavelu performed the NCA analysis. John Mundy led the ERA-NET PG Multi-Stress project, his laboratory generated all sample RNA/cDNAs. Atle M. Bones coordinated the overall development of the manuscript. All authors contributed toward improvement of the manuscript and have read and approved the manuscript.
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.
Acknowledgments
We sincerely acknowledge the funding from ERA-NET Plant Genomics (www.erapg.org/) for Multi-Stress project. The arrays experiments were run at the DTU Multi-Assay Core (DMAC) thanks to an infrastructure grant from the Technical University of Denmark. Pankaj Barah was supported by the Biotechnology and Functional Genomics (FUGE) program of the Norwegian Research Council through grants NFR 184147/S10, 185173/V40, 184146/S10 and Naresh D. Jayavelu was supported by the Norwegian Research Council, through grant number 7017430.
Supplementary Material
The Supplementary Material for this article can be found online at: www.frontiersin.org/Crop_Science_and_Horticulture/10.3389/fpls.2013.00532/abstract
Table S1 | List of differentially regulated transcripts during heat stress. (A) List of all the transcripts whose expressions were altered, during heat stress in the 10 ecotypes. Both p-values and log2 fold change values have been included along with short annotations. (B) List of 244 transcription factors (TFs) whose expressions were altered, during heat stress in the 10 ecotypes. Both p-values and log2 fold change values have been included. (C) List of 203 transposable elements (TEs) whose expressions were altered, during heat stress in the 10 ecotypes. Both p-values and log2 fold change values have been included.
Table S2 | Results of Gene Set Enrichment Analysis (GSEA) using BinGO software.
Table S3 | List of Heat Shock Protein coding genes (Hsps). List of all the 145 Hsps with corresponding p-values and log2 fold change values during the heat stress treatment in 10 ecotypes.
Table S4 | TF-TG regulatory bipartite connections predicted using NCA algorithm. TF-TG regulatory bipartite connections predicted using NCA algorithm based on their expression profiles, using Pearson correlation coefficient threshold (PCC) ≥0.80. Annotations and functional roles (TF/TG) are also provided.
References
Ahuja, I., de Vos, R. C., Bones, A. M., and Hall, R. D. (2010). Plant molecular stress responses face climate change. Trends Plant Sci. 15, 664–674. doi: 10.1016/j.tplants.2010.08.002
Albert, R. (2005). Scale-free networks in cell biology. J. Cell. Sci. 118, 4947–4957. doi: 10.1242/jcs.02714
Alonso-Blanco, C., Aarts, M. G. M., Bentsink, L., Keurentjes, J. J. B., Reymond, M., Vreugdenhil, D., et al. (2009). What has natural variation taught us about plant development, physiology, and adaptation? Plant Cell 21, 1877–1896. doi: 10.1105/tpc.109.068114
Alonso-Blanco, C., and Koornneef, M. (2000). Naturally occurring variation in Arabidopsis: an underexploited resource for plant genetics. Trends Plant Sci. 5, 22–29. doi: 10.1016/S1360-1385(99)01510-1
Al-Quraan, N. A., Locy, R. D., and Singh, N. K. (2012). Heat and cold stresses phenotypes of Arabidopsis thaliana calmodulin mutants: regulation of gamma-aminobutyric acid shunt pathway under temperature stress. Int. J. Plant Biol. 3:e2. doi: 10.4081/pb.2012.e2
Alter, O., Brown, P. O., and Botstein, D. (2000). Singular value decomposition for genome-wide expression data processing and modeling. Proc. Natl. Acad Sci. U.S.A. 97, 10101–10106. doi: 10.1073/pnas.97.18.10101
Alvarez, A., and Woolf, P. J. (2011). Partially observed bipartite network analysis to identify predictive connections in transcriptional regulatory networks. BMC Syst. Biol. 5:86. doi: 10.1186/1752-0509-5-86
Ashburner, M., Ball, C. A., Blake, J. A., Botstein, D., Butler, H., Cherry, J. M., et al. (2000). Gene ontology: tool for the unification of biology. The gene ontology consortium. Nat. Genet. 25, 25–29. doi: 10.1038/75556
Barah, P., Jayavelu, N. D., Rasmussen, S., Nielsen, H. B., Mundy, J., and Bones, A. M. (2013). Genome-scale cold stress response regulatory networks in ten Arabidopsis thaliana ecotypes. BMC Genomics 14:722. doi: 10.1186/1471-2164-14-722
Bita, C. E., and Gerats, T. (2013). Plant tolerance to high temperature in a changing environment: scientific fundamentals and production of heat stress-tolerant crops. Front. Plant Sci. 4:273. doi: 10.3389/fpls.2013.00273
Borevitz, J. O., and Nordborg, M. (2003). The impact of genomics on the study of natural variation in Arabidopsis. Plant Physiol. 132, 718–725. doi: 10.1104/pp.103.023549
Brown, M. E., and Funk, C. C. (2008). Climate. Food security under climate change. Science 319, 580–581. doi: 10.1126/science.1154102
Carrera, J., Rodrigo, G., Jaramillo, A., and Elena, S. F. (2009). Reverse-engineering the Arabidopsis thaliana transcriptional network under changing environmental conditions. Genome Biol. 10, R96. doi: 10.1186/gb-2009-10-9-r96
Carvalho, B. S., and Irizarry, R. A. (2010). A framework for oligonucleotide microarray preprocessing. Bioinformatics 26, 2363–2367. doi: 10.1093/bioinformatics/btq431
Chinnusamy, V., and Zhu, J. K. (2009). Epigenetic regulation of stress responses in plants. Curr. Opin. Plant Biol. 12, 133–139. doi: 10.1016/j.pbi.2008.12.006
Ciais, P., Reichstein, M., Viovy, N., Granier, A., Ogee, J., Allard, V., et al. (2005). Europe-wide reduction in primary productivity caused by the heat and drought in 2003. Nature 437, 529–533. doi: 10.1038/nature03972
Cui, X., and Churchill, G. A. (2003). Statistical tests for differential expression in cDNA microarray experiments. Genome Biol. 4, 210. doi: 10.1186/gb-2003-4-4-210
Davuluri, R. V., Sun, H., Palaniswamy, S. K., Matthews, N., Molina, C., Kurtz, M., et al. (2003). AGRIS: arabidopsis gene regulatory information server, an information resource of Arabidopsis cis-regulatory elements and transcription factors. BMC Bioinform. 4:25. doi: 10.1186/1471-2105-4-25
Doebley, J. F., Gaut, B. S., and Smith, B. D. (2006). The molecular genetics of crop domestication. Cell 127, 1309–1321. doi: 10.1016/j.cell.2006.12.006
Endo, M., Tsuchiya, T., Hamada, K., Kawamura, S., Yano, K., Ohshima, M., et al. (2009). High temperatures cause male sterility in rice plants with transcriptional alterations during pollen development. Plant Cell Physiol. 50, 1911–1922. doi: 10.1093/pcp/pcp135
Fedoroff, N. V. (2012). Presidential address. Transposable elements, epigenetics, and genome evolution. Science 338, 758–767. doi: 10.1126/science.338.6108.758
Galbraith, S. J., Tran, L. M., and Liao, J. C. (2006). Transcriptome network component analysis with limited microarray data. Bioinformatics 22, 1886–1894. doi: 10.1093/bioinformatics/btl279
Grativol, C., Hemerly, A. S., and Ferreira, P. C. G. (2012). Genetic and epigenetic regulation of stress responses in natural plant populations. Biochim. Biophys. Acta 1819, 176–185. doi: 10.1016/j.bbagrm.2011.08.010
Guberman, J. M., Ai, J., Arnaiz, O., Baran, J., Blake, A., Baldock, R., et al. (2011). BioMart Central Portal: an open database network for the biological community. Database (Oxford) 2011, bar041. doi: 10.1093/database/bar041
Guo, A., He, K., Liu, D., Bai, S., Gu, X., Wei, L., et al. (2005). DATF: a database of Arabidopsis transcription factors. Bioinformatics 21, 2568–2569. doi: 10.1093/bioinformatics/bti334
Herrgard, M. J., Covert, M. W., and Palsson, B. O. (2004). Reconstruction of microbial transcriptional regulatory networks. Curr. Opin. Biotechnol. 15, 70–77. doi: 10.1016/j.copbio.2003.11.002
Holter, N. S., Mitra, M., Maritan, A., Cieplak, M., Banavar, J. R., and Fedoroff, N. V. (2000). Fundamental patterns underlying gene expression profiles: simplicity from complexity. Proc. Natl. Acad. Sci. U.S.A. 97, 8409–8414. doi: 10.1073/pnas.150242097
Ideker, T., and Krogan, N. J. (2012). Differential network biology. Mol. Syst. Biol. 8, 565. doi: 10.1038/msb.2011.99
Irizarry, R. A., Bolstad, B. M., Collin, F., Cope, L. M., Hobbs, B., and Speed, T. P. (2003). Summaries of Affymetrix GeneChip probe level data. Nucleic Acids Res. 31, e15. doi: 10.1093/nar/gng015
Jagadish, S. V. K., Muthurajan, R., Oane, R., Wheeler, T. R., Heuer, S., Bennett, J., et al. (2010). Physiological and proteomic approaches to address heat tolerance during anthesis in rice (Oryza sativa L.). J. Exp. Bot. 61, 143–156. doi: 10.1093/jxb/erp289
Kao, K. C., Yang, Y. L., Boscolo, R., Sabatti, C., Roychowdhury, V., and Liao, J. C. (2004). Transcriptome-based determination of multiple transcription regulator activities in Escherichia coli by using network component analysis. Proc. Natl. Acad. Sci. U.S.A. 101, 641–646. doi: 10.1073/pnas.0305287101
Kaplan, F., Kopka, J., Haskell, D. W., Zhao, W., Schiller, K. C., Gatzke, N., et al. (2004). Exploring the temperature-stress metabolome of Arabidopsis. Plant Physiol. 136, 4159–4168. doi: 10.1104/pp.104.052142
Kerr, M. K., Martin, M., and Churchill, G. A. (2000). Analysis of variance for gene expression microarray data. J. Comput. Biol. 7, 819–837. doi: 10.1089/10665270050514954
Keurentjes, J. J. B., Fu, J. Y., Terpstra, I. R., Garcia, J. M., van den Ackerveken, G., Snoek, L. B., et al. (2007). Regulatory network construction in Arabidopsis by using genome-wide gene expression quantitative trait loci. Proc. Natl. Acad. Sci. U.S.A. 104, 1708–1713. doi: 10.1073/pnas.0610429104
Lamesch, P., Berardini, T. Z., Li, D., Swarbreck, D., Wilks, C., Sasidharan, R., et al. (2012). The Arabidopsis Information Resource (TAIR): improved gene annotation and new tools. Nucleic Acids Res. 40, D1202–D1210. doi: 10.1093/nar/gkr1090
Larkindale, J., Hall, J. D., Knight, M. R., and Vierling, E. (2005). Heat stress phenotypes of Arabidopsis mutants implicate multiple signaling pathways in the acquisition of thermotolerance. Plant Physiol. 138, 882–897. doi: 10.1104/pp.105.062257
Larkindale, J., and Knight, M. R. (2002). Protection against heat stress-induced oxidative damage in Arabidopsis involves calcium, abscisic acid, ethylene, and salicylic acid. Plant Physiol. 128, 682–695. doi: 10.1104/pp.010320
Lee, G. J., and Vierling, E. (2000). A small heat shock protein cooperates with heat shock protein 70 systems to reactivate a heat-denatured protein. Plant Physiol. 122, 189–197. doi: 10.1104/pp.122.1.189
Lee, J. H., Yun, H. S., and Kwon, C. (2012). Molecular communications between plant heat shock responses and disease resistance. Mol. Cells 34, 109–116. doi: 10.1007/s10059-012-0121-3
Liao, J. C., Boscolo, R., Yang, Y. L., Tran, L. M., Sabatti, C., and Roychowdhury, V. P. (2003). Network component analysis: reconstruction of regulatory signals in biological systems. Proc. Natl. Acad. Sci. U.S.A. 100, 15522–15527. doi: 10.1073/pnas.2136632100
Lobell, D. B., Schlenker, W., and Costa-Roberts, J. (2011). Climate trends and global crop production since 1980. Science 333, 616–620. doi: 10.1126/science.1204531
Lisch, D. (2013). How important are transposons for plant evolution? Nat. Rev. Genet. 14, 49–61. doi: 10.1038/nrg3374
Maere, S., Heymans, K., and Kuiper, M. (2005). BiNGO: a Cytoscape plugin to assess overrepresentation of gene ontology categories in biological networks. Bioinformatics 21, 3448–3449. doi: 10.1093/bioinformatics/bti551
Mitchell-Olds, T., and Schmitt, J. (2006). Genetic mechanisms and evolutionary significance of natural variation in Arabidopsis. Nature 441, 947–952. doi: 10.1038/nature04878
Mittler, R., and Blumwald, E. (2010). Genetic engineering for modern agriculture: challenges and perspectives. Annu. Rev. Plant Biol. 61, 443–462. doi: 10.1146/annurev-arplant-042809-112116
Mittler, R., Finka, A., and Goloubinoff, P. (2012). How do plants feel the heat? Trends Biochem. Sci. 37, 118–125. doi: 10.1016/j.tibs.2011.11.007
Moreau, Y., Aerts, S., de Moor, B., de Strooper, B., and Dabrowski, M. (2003). Comparison and meta-analysis of microarray data: from the bench to the computer desk. Trends Genet. 19, 570–577. doi: 10.1016/j.tig.2003.08.006
Pecinka, A., Dinh, H. Q., Baubec, T., Rosa, M., Lettner, N., and Scheid, O. M. (2010). Epigenetic regulation of repetitive elements is attenuated by prolonged heat stress in arabidopsis. Plant Cell 22, 3118–3129. doi: 10.1105/tpc.110.078493
Popova, O. V., Dinh, H. Q., Aufsatz, W., and Jonak, C. (2013). The RdDM pathway is required for basal heat tolerance in arabidopsis. Mol Plant 6, 396–410. doi: 10.1093/mp/sst023
Qian, J., Lin, J., Luscombe, N. M., Yu, H., and Gerstein, M. (2003). Prediction of regulatory networks: genome-wide identification of transcription factor targets from gene expression data. Bioinformatics 19, 1917–1926. doi: 10.1093/bioinformatics/btg347
Qu, A. L., Ding, Y. F., Jiang, Q., and Zhu, C. (2013). Molecular mechanisms of the plant heat stress response. Biochem. Biophys. Res. Commun. 432, 203–207. doi: 10.1016/j.bbrc.2013.01.104
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
R Core Team (2012). R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing.
Riano-Pachon, D. M., Ruzicic, S., Dreyer, I., and Mueller-Roeber, B. (2007). PlnTFDB: an integrative plant transcription factor database. BMC Bioinformatics 8:42. doi: 10.1186/1471-2105-8-42
Richter, K., Haslbeck, M., and Buchner, J. (2010). The heat shock response: life on the verge of death. Mol. Cell 40, 253–266. doi: 10.1016/j.molcel.2010.10.006
Rocco, M., Arena, S., Renzone, G., Scippa, G. S., Lomaglio, T., Verrillo, F., et al. (2013). Proteomic analysis of temperature stress-responsive proteins in Arabidopsis thaliana rosette leaves. Mol. Biosyst. 9, 1257–1267. doi: 10.1039/c3mb70137a
Saeed, A. I., Hagabati, N. K., Braisted, J. C., Liang, W., Sharov, V., Howe, E. A., et al. (2006). TM4 microarray software suite. Methods Enzymol. 411, 134–193. doi: 10.1016/S0076-6879(06)11009-5
Segal, E., Shapira, M., Regev, A., Pe'er, D., Botstein, D., Koller, D., et al. (2003). Module networks: identifying regulatory modules and their condition-specific regulators from gene expression data. Nat. Genet. 34, 166–176. doi: 10.1038/ng1165
Semenov, M. A., and Shewry, P. R. (2011). Modeling predicts that heat stress, not drought, will increase vulnerability of wheat in Europe. Sci. Rep. 1:66. doi: 10.1038/srep00066
Somerville, C., and Koornneef, M. (2002). A fortunate choice: the history of Arabidopsis as a model plant. Nat. Rev. Genet. 3, 883–889. doi: 10.1038/nrg927
Stott, P. A., Stone, D. A., and Allen, M. R. (2004). Human contribution to the European heatwave of 2003. Nature 432, 610–614. doi: 10.1038/nature03089
Suzuki, N., Bajad, S., Shuman, J., Shulaev, V., and Mittler, R. (2008). The transcriptional co-activator MBF1c is a key regulator of thermotolerance in Arabidopsis thaliana. J. Biol. Chem. 283, 9269–9275. doi: 10.1074/jbc.M709187200
Swindell, W. R., Huebner, M., and Weber, A. P. (2007). Transcriptional profiling of Arabidopsis heat shock proteins and transcription factors reveals extensive overlap between heat and non-heat stress response pathways. BMC Genomics 8:125. doi: 10.1186/1471-2164-8-125
Tirosh, I., and Barkai, N. (2011). Inferring regulatory mechanisms from patterns of evolutionary divergence. Mol. Syst. Biol. 7, 530. doi: 10.1038/msb.2011.60
Wahid, A., Gelani, S., Ashraf, M., and Foolad, M. R. (2007). Heat tolerance in plants: an overview. Environ. Exp. Bot. 61, 199–223. doi: 10.1016/j.envexpbot.2007.05.011
Wang, J. G., Qiu, X. J., Li, Y. H., Deng, Y. P., and Shi, T. L. (2011). A transcriptional dynamic network during Arabidopsis thaliana pollen development. BMC Syst. Biol. 5. doi: 10.1186/1752-0509-5-S3-S8
Wang, X., Weigel, D., and Smith, L. M. (2013). Transposon variants and their effects on gene expression in Arabidopsis. PLoS Genet. 9, 1–13. doi: 10.1371/journal.pgen.1003255
Weigel, D., and Nordborg, M. (2005). Natural variation in arabidopsis. How do we find the causal genes? Plant Physiol. 138, 567–568. doi: 10.1104/pp.104.900157
Weston, D. J., Karve, A. A., Gunter, L. E., Jawdy, S. S., Yang, X. H., Allen, S. M., et al. (2011). Comparative physiology and transcriptional networks underlying the heat shock response in Populus trichocarpa, Arabidopsis thaliana and Glycine max. Plant Cell Environ. 34, 1488–1506. doi: 10.1111/j.1365-3040.2011.02347.x
Wheeler, T., and von Braun, J. (2013). Climate change impacts on global food security. Science 341, 508–513. doi: 10.1126/science.1239402
Wheeler, T. R., Craufurd, P. Q., Ellis, R. H., Porter, J. R., and Prasad, P. V. V. (2000). Temperature variability and the yield of annual crops. Agric. Ecosyst. Environ. 82, 159–167. doi: 10.1016/S0167-8809(00)00224-3
Yeh, C. H., Kaplinsky, N. J., Hu, C., and Charng, Y. Y. (2012). Some like it hot, some like it warm: phenotyping to explore thermotolerance diversity. Plant Sci. 195, 10–23. doi: 10.1016/j.plantsci.2012.06.004
Yilmaz, A., Mejia-Guerra, M. K., Kurz, K., Liang, X., Welch, L., and Grotewold, E. (2011). AGRIS: the Arabidopsis Gene Regulatory Information Server, an update. Nucleic Acids Res. 39, D1118–D1122. doi: 10.1093/nar/gkq1120
Zeller, G., Henz, S. R., Widmer, C. K., Sachsenberg, T., Ratsch, G., Weigel, D., et al. (2009). Stress-induced changes in the Arabidopsis thaliana transcriptome analyzed using whole-genome tiling arrays. Plant J. 58, 1068–1082. doi: 10.1111/j.1365-313X.2009.03835.x
Zhang, J. X., Klueva, N. Y., Wang, Z., Wu, R., Ho, T. H. D., and Nguyen, H. T. (2000). Genetic engineering for abiotic stress resistance in crop plants. In Vitro Cell. Dev. Biol. Plant 36, 108–114. doi: 10.1007/s11627-000-0022-6
Zou, C., Sun, K., Mackaluso, J. D., Seddon, A. E., Jin, R., Thomashow, M. F., et al. (2011a). Cis-regulatory code of stress-responsive transcription in Arabidopsis thaliana. Proc. Natl. Acad. Sci. U.S.A. 108, 14992–14997. doi: 10.1073/pnas.1103202108
Keywords: heat stress, natural variation, microarray transcriptional profiling, regulatory networks, systems biology, Arabidopsis thaliana, ecotypes
Citation: Barah P, Jayavelu ND, Mundy J and Bones AM (2013) Genome scale transcriptional response diversity among ten ecotypes of Arabidopsis thaliana during heat stress. Front. Plant Sci. 4:532. doi: 10.3389/fpls.2013.00532
Received: 17 October 2013; Accepted: 10 December 2013;
Published online: 26 December 2013.
Edited by:
Craita Bita, Radboud University Nijmegen, NetherlandsReviewed by:
Soren K. Rasmussen, University of Copenhagen, DenmarkChenglin Chai, University of Missouri-Columbia, USA
Copyright © 2013 Barah, Jayavelu, Mundy and Bones. 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: Atle M. Bones, Cell Molecular Biology and Genomics Group, Department of Biology, Norwegian University of Science and Technology, Hoegskoleringen 5, N-7491 Trondheim, Norway e-mail: atle.bones@bio.ntnu.no