- Research article
- Open Access
Transcriptome analysis of heat stress response in switchgrass (Panicum virgatumL.)
BMC Plant Biology volume 13, Article number: 153 (2013)
Global warming predictions indicate that temperatures will increase by another 2-6°C by the end of this century. High temperature is a major abiotic stress limiting plant growth and productivity in many areas of the world. Switchgrass (Panicum virgatum L.) is a model herbaceous bioenergy crop, due to its rapid growth rate, reliable biomass yield, minimal requirements of water and nutrients, adaptability to grow on marginal lands and widespread distribution throughout North America. The effect of high temperature on switchgrass physiology, cell wall composition and biomass yields has been reported. However, there is void in the knowledge of the molecular responses to heat stress in switchgrass.
We conducted long-term heat stress treatment (38°/30°C, day/night, for 50 days) in the switchgrass cultivar Alamo. A significant decrease in the plant height and total biomass was evident in the heat stressed plants compared to controls. Total RNA from control and heat stress samples were used for transcriptome analysis with switchgrass Affymetrix genechips. Following normalization and pre-processing, 5365 probesets were identified as differentially expressed using a 2-fold cutoff. Of these, 2233 probesets (2000 switchgrass unigenes) were up-regulated, and 3132 probesets (2809 unigenes) were down-regulated. Differential expression of 42 randomly selected genes from this list was validated using RT-PCR. Rice orthologs were retrieved for 78.7% of the heat stress responsive switchgrass probesets. Gene ontology (GOs) enrichment analysis using AgriGO program showed that genes related to ATPase regulator, chaperone binding, and protein folding was significantly up-regulated. GOs associated with protein modification, transcription, phosphorus and nitrogen metabolic processes, were significantly down-regulated by heat stress.
Plausible connections were identified between the identified GOs, physiological responses and heat response phenotype observed in switchgrass plants. Comparative transcriptome analysis in response to heat stress among four monocots – switchgrass, rice, wheat and maize identified 16 common genes, most of which were associated with protein refolding processes. These core genes will be valuable biomarkers for identifying heat sensitive plant germplasm since they are responsive to both short duration as well as chronic heat stress treatments, and are also expressed in different plant growth stages and tissue types.
Switchgrass (Panicum virgatum L.) is an economically important, warm-season, and widely adapted C4 perennial grass . Due to its rapid growth rate, reliable biomass yield across locations, minimal requirements of water and nutrients, adaptability to growth on marginal lands and widespread distribution throughout North America, switchgrass has been selected in 1992 by the U. S. Department of Energy (USDOE) as a model herbaceous bioenergy crop for the development of renewable feed stock resource to produce transportation fuel . Based on morphology and habitat preference, switchgrass has been classified into two groups: lowland ecotype and upland ecotype . Lowland ecotypes are mostly tetraploid (2n = 4× = 36), and generally adapted to wet areas with milder winter temperatures, while upland ecotypes are mainly octaploid (2n = 8× = 72) or hexaploid (2n = 6× = 54), and thrive well in drier and colder areas. Lowland plants are usually taller and have longer and wider leaf blades, fewer tillers per plant, larger stem diameter and later in heading and flowering compared with upland plants [2, 4, 5].
Switchgrass germplasm collection and breeding for increasing biomass and conversion from cellulosic feedstock to ethanol have become a high priority [4, 6, 7]. In recent years, molecular markers have been extensively used to examine variation in switchgrass germplasm [5, 8–13]. Switchgrass genetic linkage map has been established [14, 15]. Three bacterial artificial chromosome (BAC) libraries have been generated [16, 17]. Switchgrass expressed sequence tag (EST) database is being populated using traditional and RNA-Seq based technologies [18–21]. A switchgrass Affymetrix gene chip has been made available to the community through the DOE Bioenergy Science Center (BESC) co-operative efforts . The procedure for switchgrass transformation has been established [22–24]. However, compared with rice, a well-studied monocot, studies on switchgrsass responses to biotic and abiotic stresses are limited [25–27].
Global warming predictions indicate that temperatures will increase another 2- 6°C by the end of this century . High temperature can retard plant growth, development and yield, therefore, agriculture will be seriously affected by global warming in the future [29–31]. Based on the climate change observation, Behrman et al.  forecast that there will be substantial variation in switchgrass productivity within regions inside United States and over time. The southern United States, the main region of switchgrass production at present, is predicted to have the largest decrease in biomass in 2080–2090, due to the increased temperature and decreased precipitation . Several studies reported the effect of high temperature on switchgrass focusing on physiology, composition and yields [27, 33–36]. In the present study, we conducted transcriptome analysis using Affymetrix gene chips to elucidate the transcriptional changes in response to heat stress in Switchgrass Alamo, a lowland cultivar, extensively grown as biofuel feedstock. Comparative transcriptome analysis with other monocots identified a core set of 16 common heat stress responsive genes. The identified genes will provide rational candidates for germplasm screening to enhance switchgrass heat tolerance.
Results and discussion
Phenotypic responses to heat stress in switchgrass
Switchgrass plants grown under the optimal 28°/20°C (day/night) condition produced more foliage and were nearly twice the height of the plants that were under extreme heat stress at 38°/30°C (day/night). Plants that were subjected to a moderate heat stress treatment of 33°/25°C (day/night) showed only a slight reduction in the plant height (Figure 1). In a recent report on heat stress response in four switchgrass cultivars including Alamo, it was reported that stem elongation rate, leaf elongation rate and plant height were greatly impaired, ultimately lowering the growth and biomass of switchgrass cultivars at 38°/30°C (day/night) . This reduction of total biomass by nearly 50% in response to elevated temperatures is a concern in the wake of the predicted increases in the global temperatures and particularly in the south-western USA, an ideal location for switchgrass production for biofuels. In fact, during the summer of 2011–2012 there were more than 100 days with more than 38°C temperature in the State of Oklahoma (http://www.mesonet.org).
Transcriptome analysis in response to heat stress in switchgrass
A switchgrass affymetrix array containing 122,868 probesets corresponding to 110,208 unigene transcripts has been recently developed . We used these gene chips to examine changes in switchgrass transcriptome following the heat stress treatment. The correlation between the two biological replicate experiments was 0.97. Using the arbitrary 2-fold cutoff about 4.2% of the probesets on the array was identified as heat-stress responsive. In several other plants about 5% of the transcriptome has been reported to be heat responsive [37–39], indicating that switchgrass heat stress transcriptional response is quantitively comparable to other plants.
Among the 5164 differentially expressed probesets in response to heat in switchgrass, 2076 probesets corresponding to 2002 switchgrass unigenes were induced, and 3088 probesets corresponding to 2809 switchgrass unigenes were repressed. In the heat transcriptome studies in maize, wheat, and rice, the number of induced genes was 3–6 times more than the number of repressed genes [38, 40, 41]. We surmise that the reason for this difference may be due to the differences in the heat stress treatment regimes. In the long-term heat stress imposed in our studies the observed transcript changes reflect the acclimative response. Whereas in the short duration heat stress in the other studies mentioned above, the transcriptional changes reflect the more active defense response.
To validate the transcriptional pattern identified by microarray analysis, expression of 42 differentially expressed genes was analyzed using RT-PCR. The results showed that majority of the tested genes (38 genes) followed the expression pattern observed in the gene chip experiments (Figure 2, Additional file 1), validating the array results. In the few cases where there was a discrepancy it was observed that the hybridization intensities associated with those genes were weak for either the control or the heat stressed sample.
Gene ontology analysis of switchgrass heat responsive transcripts
To investigate the biological significance of the genes regulated by heat stress in switchgrass it is important to have the gene ontology (GO) descriptions i.e., detailed annotations of gene function, biological process it is involved, and cellular location of the gene product. Since switchgrass genes have not been well annotated yet, rice orthologs of the differently regulated transcripts were identified. The best rice transcripts matching the switchgrass probesets (E-value > 1e10-5 and with at least 100 HSPs) were retrieved from the switchgrass genomics database maintained by the Samuel Roberts Noble Foundation (http://switchgrassgenomics.noble.org). Rice orthologs that showed 98-100% homology with the switchgrass probesets was identified for 4062 (78.7%) heat stress responsive switchgrass probesets. Among these, 1158 unique rice orthologs, representing 1478 switchgrass probesets were up-regulated, and 1857 unique rice orthologs, representing 2584 switchgrass probesets were down-regulated by heat stress. These unique rice orthologs were subjected to singular enrichment analysis (SEA) in agriGO to identify enriched GOs . SEA is designed to identify enriched GO terms in a list of microarray probe sets or gene identifiers. Finding enriched GO terms corresponds to finding enriched biological facts, and term enrichment level is judged by comparing query list to a background population from which the query list is derived. In this study the background query list comprised of 54,971 annotated rice genes from the MSU 6.1 version (http://rice.plantbiology.msu.edu/).
Among the heat stress induced transcripts, GOs associated with ATPase regulator, protein folding, chaperone binding and catalytic activity were significant (Figure 3 and Additional file 2). Heat stress affects the stability of various proteins, membranes, RNA species and cytoskeleton structures. In order to counter the imbalance, survive and continue grow at a higher temperature, plants have to reprogram their transcriptome, proteome, metabolome and lipidome by altering the transcripts, proteins, metabolites and lipids [43–45]. It is well known Heat shock proteins (HSPs) and other chaperones are induced by various stresses. They play an important role in protein-protein interactions such as folding, assisting in proper protein conformation, stabilizing partially unfolded proteins and prevention of unwanted protein aggregation. The induction of HSPs expression is one of the common heat responsive mechanisms in all organisms [46–50]. Thus the enrichment of the GO associated with ATPase regulator, protein folding and chaperone binding in switchgrass is not surprising. There were 50 heat shock proteins and other chaperones that were strongly induced by heat stress in switchgrass, about 2/3 of these induced genes were also observed in rice panicle (Additional file 3). Since a large number of proteins will be misfolded under heat stress, there is a need for these chaperones to assist in the repair and/or salvage process. This is an energy demanding activity and apart from ATP requires the assistance of nucleotide exchange factor proteins such as GrpE . Upregulation of the ATP synthase subunit along with three different GrpE genes in switchgrass is consistent with the findings in other plant species [38, 40, 41]. A switchgrass heat shock protein (ortholog of Os04g01740) was induced more than 66-fold by heat stress, and this pattern was further validated by RT-PCR (Figure 2). This gene was also induced 102 times in rice panicle under heat stress . Other common heat stress response proteins included DnaK proteins (Os01g62290, Os03g11910 and Os03g16920), DnaJ protein (Os05g48810), and heat shock protein (Os02g52150, Os05g44340 and Os01g04370) (Additional file 3).
Under the GO category of catalytic activity, isomerase, transferase and oxidoreductase activities were enriched in response to heat stress in switchgrass (Additional file 2). The GO for oxido-reductase activity formed the largest group with 67 genes. Many of the genes in this group were associated with oxidative stress that is caused due to excess accumulation of ROS and/or due to insufficient antioxidant defenses [45, 52–55]. A superoxide dismutase (SOD) (Os08g44770), two different amine oxidases and an alcohol oxidase were identified and could be involved in the generation of ROS during heat stress in switchgrass . Six different peroxidases were induced, among them were cytosolic, stromal and peroxisomal ascorbate peroxidase and glutathione peroxidase. Interestingly, lack of chloroplastic ascorbate peroxidase resulted in enhanced heat stress tolerance in Arabidopsis plants . This suggests that chloroplastic ROS, may play different role in switchgrass heat stress response. On the same lines, induction of thylakoid APX may be an important factor leading to heat sensitivity in Alamo plants. Cytosolic ascorbate peroxidases are important for protecting the organelles especially chloroplast ROS-scavenging system in Arabidopsis and lack of this enzyme can lead to increased protein oxidation . The induction of two different peptide methionine sulfoxide reductase, a well-known enzyme vital for repairing oxidatively damaged proteins especially photosynthetic antennae , implies a severe oxidative stress due to high temperature in switchgrass.
Eleven different cytochrome P450 genes were identified in response to heat stress in switchgrass. Upregulation of two genes for CYP71A1 suggest the increased synthesis of indole alkaloids such as secologanin . Induction of three different peroxiredoxins and a thioredoxin in switchgrass is again indicative of the oxidative stress induced by high temperature. Peroxiredoxins are important for reducing hydrogen peroxide and alkyl-hydroperoxides  and has been reported in response to heat stress in other plants . In Arabidopsis, a thioredoxin was shown to form low to oligomeric protein structures and also high molecular weight complexes in response to heat stress . Furthermore, the low weight structures showed disulfide reductase activity while the higher complexes exhibited chaperone function . Interestingly, several targets of cytosolic thioredoxin , such as ascorbate peroxidase, malate dehydrogenase, glyceroldehyde-3-phosphate dehydrogenase, alcohol dehydrogenase were up regulated in response to heat stress in switchgrass. Further characterization of the switchgrass thioredoxin will provide a better understanding for the role of this protein in the heat stress response. A switchgrass DUF538 domain containing protein (ortholog of Os01g11240) expressed at very low level under normal condition was strongly induced by heat stress (56.7-fold). The function of DUF538 protein is unknown, however it has been proposed as putative candidate for the common stress related proteins in the plant system . Exogenously applying maltose-binding fusion protein (MBP-DUF538) on the leaves of tobacco can elevate activities of redox enzymes including catalase, peroxidase, polyphenol oxidase and phenyalanine ammonia lyase . We speculate that the switchgrass DUF538 protein is involved in ROS detoxification during heat stress.
The GO for the unfolded protein response consisted of 15 genes and six of them were annotated as T-complex proteins. The first member of this complex was identified in the hyperthermophilic archeon, Pyrodictium occultum, in response to heat stress and was shown to have ATPase activity . Detailed structural analysis revealed that this enzyme is a complex of eight subunits and the fact that it is heat-inducible led to the name thermosome . Careful sequence analysis of these switchgrass T-complex genes (data not shown) indicated that they were part of the group II chaperonin or the CCT complex .
Just as proteins become misfolded during heat stress, it is conceivable that the secondary structures of mRNAs can be disrupted in response to heat stress. We identified five different DEAD box RNA helicases that were up regulated in response to heat stress in switchgrass. Several studies on DEAD box RNA helicases in response to cold and salinity stress has been previously reported [70–74]. A rice DEAD box RNA helicase, OSABP, was strongly repressed in response to heat and cold stress . A detailed analysis of the switchgrass DEAD box RNA helicases induced in response to heat merits further attention as these proteins may function as RNA chaperones.
Heat stress induces inward calcium flux, the increased calcium ion (Ca2+) level can in turn regulate multiple signaling pathway in plants . The inward flux of calcium can activate several calcium-dependent protein kinases (CDPKs), which can, in turn, activate multiple mitogen-activated protein kinases (MAPKs) . A calcium-transporting ATPase (Os04g51610), a calcium-dependent protein kinase CPK1 adapter protein (Os06g50146), calcineurin B (Os01g39770), and four calmodulin dependent protein kinases were induced in switchgrass by heat stress. In Arabidopsis, calmodulin AtCaM3 is required for heat stress signaling and is involved in the activation of WRKY and HSF transcription factors [78–81]. Similarly, the WRKY transcription factor (Os03g55164) and four HSFs associated with MAPK signaling were induced by heat stress in switchgrass. Surprisingly, a calreticulin precursor (Os07g14270) is the highest induced gene (166.4-fold) and was validated by RT-PCR (Figure 2 and Additional file 1). Calreticulin (CRT) is a key Ca2+-binding protein mainly resident in the endoplasmic reticulum in plants. CRT plays important roles in a variety of cellular processes including Ca2+ signaling, protein folding and as a key alleviator of endoplasmic reticulum stress [82, 83]. CRT mRNA and protein are upregulated in response to cold stress, salt stress and exogenous phytohormones  and this study revealed that CRT is also induced in the heat stress response.
In proportion with the larger number of repressed genes in response to heat stress in switchgrass, more enriched GO categories were identified (Additional file 4 and Additional file 5). In order to reduce the number of GO terms, enriched GO categories with false discovery rates (FDR) < 0.05 from AgriGO analysis were submitted to the REVIGO (REduce and Visualize GO) program . Using the Uniprot database as background and the default semantic similarity measure (Simrel), this analysis clearly showed that biological processes associated with metabolism, cellular homeostasis, cell death, regulation of transcription and transporters were significantly over-represented among the genes repressed by heat stress in switchgrass (Figure 4).
All the genes involved in glycolysis were strongly repressed by heat stress in switchgrass. The strong repression of the genes involved in glucose metabolism strongly supports the significant reduction reported in the cellulose and hemicellulose levels in response to heat stress in switchgrass cultivars . The down regulation of metabolism in general, nitrogen and phosphorus metabolism in particular will have pronounced impact on vegetative growth (Figure 1). This is also supported by the observation that the dry shoot biomass in switchgrass cultivars was reduced by nearly 66% in response to heat stress . These studies clearly demonstrate that under higher temperatures, genes associated with carbon fixation is down regulated, leading to reduced biomass. The reduction in the growth and development is probably brought into effect by the strong repression of a gamut of transcription factor families including 19 WRKYs, 13 NAMs (No Apical Meristem), 11 Myb TFs, nine AP2 domain containing TFs, five AUX/IAA type TFs (Additional file 6). Based on the predictions of higher temperatures in areas where switchgrass will be used as the main bioenergy crop our transcriptome studies indicate significant repression of carbon fixation processes that ultimately will have a negative impact on cellulosic bioethanol production.
Comparative transcriptome analysis of heat stress response in monocots
In order to identify the unique responses and commonalities in response to heat stress among monocots we identified three other transcriptome studies in rice , wheat  and maize . We examined transcriptome profiles from these studies using publicly available datasets. The switchgrass heat response transcriptome had the largest number of unique genes (2272). In contrast to the short-term heat stress treatments (few hours), long-term heat treatments in switchgrass (for up to 50 days) may be providing a totally different snap shot of the transcriptome. Among unique differentially expressed genes, 1462 were down regulated and 886 genes were up. Among the enriched GOs of repressed switchgrass-specific genes, redox homeostasis, regulation of transcription, transport and ubiquitination was particularly interesting since these were up regulated in the short-term heat treatment studies in rice . In response to 24-h heat treatment in wheat, (considered as long-term heat stress), two different NADPH oxidases were induced while these genes were repressed in a short-term one hour treatment . In contrast to short-term stress, we speculate that the long-term heat regime lead to excess oxidative stress that causes a redox imbalance. Similarly, short-term heat response invokes a sleuth of transcription factors in wheat and maize [38, 40], while the long-term heat stress in switchgrass leads to repression of transcription factors that are associated with growth and development. This suggests that key processes involved in signaling and/or adaptation based on transcriptome are very different for short and long term heat stress treatments. However, it is important to realize that many of these genes exists as gene families and specific gene family members may be specific for certain stress and/or developmental stage as has been shown in maize and wheat [38, 40]. GO enrichment analysis using the up regulated gene identifiers uniquely responsive to heat stress in switchgrass did not reveal any new categories at the level 3 terms i.e. two levels more specific than the top-level, molecular function or biological process.
Sixteen genes were identified in all the four monocots (Figure 5 & Table 1) despite the fact that the heat stress treatment conditions and the plant growth stages used in these studies were extremely different. Among them, only two genes were down-regulated (Table 1). Interestingly, the expression patterns of these 16 genes are very similar in the four different plant species (Table 1). The significant GO category of the common up-regulated genes is related with protein folding and unfolding, a common biochemical response to short-term and long-term heat stress.
Members of 'Activator of Hsp90 ATPase’ bind to the molecular chaperone HSP82 and stimulate its ATPase activity . Small heat shock protein, sHsps, are small stress induced proteins with monomeric masses between 12–43 kDa, are believed to be ATP-independent chaperones that prevent aggregation and are important in refolding in combination with other Hsps . GroEL_like type I chaperonin are involved in productive folding of proteins and in plants are called as cpn60 . With the aid of co-chaperonin GroES, GroEL encapsulates non-native substrate proteins inside the cavity of the GroEL-ES complex and promotes folding by using energy derived from ATP hydrolysis. Chaperonin 10 kDa subunit (cpn10 or GroES) cooperates with chaperonin 60 (cpn60 or GroEL), an ATPase, to assist the folding and assembly of proteins and is found in cytosol, as well as in the matrix of mitochondria and chloroplasts. It forms heptameric rings with a dome-like structure, forming a lid to the large cavity of the tetradecameric cpn60 cylinder and thereby tightly regulating release and binding of proteins to the cpn60 surface. The 60 kDa chaperonin alpha subunit is a part of the T-complex proteins important in the unfolded protein response of cytoskeletal proteins-actin and tubulins and also other proteins . Aha1 is one of several co-chaperones, which regulate the dimeric chaperone Hsp90. Hsp90, Aha1, and other accessory proteins interact in a chaperone cycle driven by ATP binding and hydrolysis. Aha1 promotes dimerization of the N-terminal domains of Hsp90, and stimulates its low intrinsic ATPase activity . Aha1 may regulate the dwell time of Hsp90 with client proteins. Aha1 may act as either a negative or positive regulator of chaperone-dependent activation, depending on the client protein . GrpE is the adenine nucleotide exchange factor of DnaK (Hsp70)-type ATPases and is important for thermo-tolerance to chronic heat stress in plants . The GrpE dimer binds to the ATPase domain of Hsp70 catalyzing the dissociation of ADP, which enables rebinding of ATP, one step in the Hsp70 reaction cycle in protein folding and can direct incompetent “client” proteins towards degradation . ATP-dependent Clp protease, ClpB, has been described as HSP101 in plants and is indispensable for basal thermotolerance and negatively impacts root growth  though in the absence of stress, this protein is dispensable for normal growth and development.
The universal stress protein (Usp) is a small cytoplasmic protein whose expression is enhanced when the cell is exposed to stress agents. Usp enhances the rate of cell survival during prolonged exposure to such conditions, and may provide a general “stress endurance” activity . FKBP-type peptidyl-prolyl cis-trans isomerase with a Tetratricopeptide repeat domain is involved in chaperone, cell-cycle, transcription, and protein transport complexes. We speculate that in response to heat stress it most likely behaves as a chaperone. Identification of a gene similar to beta-1,3-glucanase suggests that similar to the biotic stresses, modification of cell walls is an important component of the heat stress response. L-Asparaginase type 2-like enzymes are important for nitrogen remobilization and seed production . Asparaginase has been shown to be important for low temperature response in soybeans . Asparaginase as a common component in the heat stress response in plants may be important for remobilizing the scanty nitrogen reserves for ensuring seed development.
The present study identified significantly altered transcripts in switchgrass under chronic heat stress. Using a comparative transcriptome analysis we identified 16 common genes in the heat stress response in plants. Since these genes are expressed both under short and long-term heat stress treatments, they are ideal biomarkers for screening germplasm for thermo-tolerance. This will be valuable for developing new plant types that can adapt and thrive well under high temperatures that are predicted for the future.
Plant growth and treatment
Switchgrass (Panicum virgatum L.) Alamo cultivar seeds were sown in pots (0.2 meter diameter × 0.45 m tall) filled with pure, fine sand soil and grown in growth chambers (Conviron Ltd., Winnipeg, Canada) at 28°/20°C (day/night) with a photoperiod of 14 h/10 h (day/night). After germination, switchgrass seedlings were thinned to four plants per pot. Thirty days after sowing, switchgrass seedlings were subjected to heat stress treatment of 38°/30°C (day/night) for 50 days. Leaf samples from plants growing under control condition (28°/20°C) and heat stress (38°/30°C) were harvested and snap frozen in liquid nitrogen for RNA isolation at the end of treatment.
RNA isolation and switchgrass affymetrix genechip hybridizations
Total RNA from switchgrass leaves was isolated using RNEasy Plant Isolation kits (Qiagen Valencia, CA, USA). The quality of the RNA was tested using the BioAnalyzer (Agilent, Palo Alto, CA, USA). About 10 μg RNA was used for probe labeling according to procedure of GeneChip 3’ IVT Express Kit (Affymetrix, Santa Clara, CA). Labeled probe was hybridized to Affymetrix switchgrass cDNA chip containing more than 120,000 probe sets. Hybridizations were done in duplicates using RNA prepared from two biological replicate samples.
Microarray data analysis
Microarray data normalization was conducted by robust multi-array average (RMA) provided with Expression Console . Probesets showing a normalized transcript level of more than 28 were selected for further analysis. The reason for this cut-off is that values obtained from probe sets with 'absent’ values, as determined by the Affymetrix software, consistently exhibited log 2 normal distribution with an upper boundary below 5 (log228 = 4.8). Furthermore, the 20 negative controls probe sets detected no transcripts when hybridized to switchgrass RNA using this threshold. Probesets showing more than two-fold change (heat stress/control) in expression were considered as differentially regulated genes.
Analysis of enriched gene ontologies
The best rice transcripts matching the switchgrass probesets (E-value > 1×10-5 and >100 HSPs) were retrieved from the switchgrass genomics database maintained by the Samuel Roberts Noble Foundation (http://switchgrassgenomics.noble.org). Rice orthologs of switchgrass genes were input for Gene ontology analysis using agriGO (http://bioinfo.cau.edu.cn/agriGO/) and REVIGO (http://revigo.irb.hr/)  software.
Gene expression validation
One microgram total RNA was used for cDNA synthesis with superscript reverse transcriptase II (Invitrogen). A 1:10 dilution of cDNA was used for PCR amplification. Forty-two differentially expressed genes were selected for RT- PCR analysis. Amplification of actin gene was used as the reference for normalization. The gene list, primer sequence, and corresponding rice orthologs are listed in Additional file 1.
A detailed analysis of the rice heat responsive transcriptome was published recently . In the rice study 3364 heat responsive probesets, represented 3213 rice genes, of which 2451 were identified as unique loci. These 2451 unique differentially expressed genes were used for comparative analysis with the genes identified in switchgrass analysis. Nearly 6560 differentially expressed probesets were reported in response to heat stress treatment in wheat . Rice orthologs of the wheat probesets were retrieved using the 'microarray platform translator’ tool in PLEXdb . Using this tool, 5288 rice probesets were identified while for 1272 probesets there were no orthologs using the default BLAST parameters. The rice affymetrix probe identifiers were then fed into the model genome interrogator tool in the PLEXdb to retrieve the corresponding rice locus identifiers. This analysis yielded 1831 unique rice loci. A heat transcriptome analysis in maize using the long oligo array reported 1081 differentially expressed genes . Using the PLEXdb microarray platform translator tool 887 corresponding rice orthologs were retrieved and this corresponded to 633 unique rice loci. A four-way comparison of the heat responsive transcriptomes of rice, wheat, maize and switchgrass was undertaken.
Availability of supporting data
The microarray data sets supporting the results of this article are available in the ArrayExpress repository under the accession number E-MTAB-1897.
Hitchcock AS, Chase A: Manual of the Grass of the United State. Washington, D. C: U S Government Printing Office; 1950.
Sanderson MA, Adler PR, Boateng AA, Casler MD, Sarath G: Switchgrass as a biofuels feedstock in the USA. Can J Plant Sci. 2006, 86: 1315-1325.
Porter CL: An analysis of variation between upland and lowland switchgrass panicum virgatum L in Central Oklahoma. Ecology. 1966, 47: 980-992.
Casler MD, Tobias CM, Kaeppler SM, Buell CR, Wang ZY, Cao P, Schmutz J, Ronald P: The switchgrass genome: tools and strategies. Plant Genome. 2011, 4: 273-282.
Gunter LE, Tuskan GA, Wullschleger SD: Diversity among populations of switchgrass based on RAPD markers. Crop Sci. 1996, 36: 1017-1022.
Marra M, Keene T, Skousen J, Griggs T: Switchgrass yield on reclaimed surface mines for bioenergy production. J Environ Qual. 2013, 42 (3): 696-703.
Brudecki G, Cybulska I, Rosentrater K: Optimization of clean fractionation process applied to switchgrass to produce pulp for enzymatic hydrolysis. Bioresour Technol. 2013, 131: 101-112.
Cortese LM, Honig J, Miller C, Bonos SA: Genetic diversity of twelve switchgrass populations using molecular and morphological markers. Bioenergy Res. 2010, 3: 262-271.
Huang SX, Su XJ, Haselkorn R, Gornicki P: Evolution of switchgrass (Panicum virgatum L.) based on sequences of the nuclear gene encoding plastid acetyl-CoA carboxylase. Plant Sci. 2003, 164: 43-49.
Hultquist SJ, Vogel KP, Lee DJ, Arumuganathan K, Kaeppler S: Chloroplast DNA and nuclear DNA content variations among cultivars of switchgrass, Panicum virgatum L. Crop Sci. 1996, 36: 1049-1052.
Narasimhamoorthy B, Saha MC, Swaller T, Bouton JH: Genetic diversity in switchgrass collections assessed by EST-SSR markers. Bioenergy Res. 2008, 1: 136-146.
Todd J, Wu YQ, Wang Z, Samuels T: Genetic diversity in tetraploid switchgrass revealed by AFLP marker polymorphisms. Genet Mol Res. 2011, 10: 2976-2986.
Zalapa JE, Price DL, Kaeppler SM, Tobias CM, Okada M, Casler MD: Hierarchical classification of switchgrass genotypes using SSR and chloroplast sequences: ecotypes, ploidies, gene pools, and cultivars. Theor Appl Genet. 2011, 122: 805-817.
Okada M, Lanzatella C, Saha MC, Bouton J, Wu R, Tobias CM: Complete switchgrass genetic maps reveal subgenome collinearity, preferential pairing and multilocus interactions. Genetics. 2010, 185: 745-760.
Liu L, Wu Y, Wang Y, Samuels T: A high-density simple sequence repeat-based genetic linkage map of switchgrass. G3. 2012, 2: 357-370.
Saski CA, Li Z, Feltus FA, Luo H: New genomic resources for switchgrass: a BAC library and comparative analysis of homoeologous genomic regions harboring bioenergy traits. BMC Genomics. 2011, 12: 369-
Sharma MK, Sharma R, Cao P, Jenkins J, Bartley LE, Qualls M, Grimwood J, Schmutz J, Rokhsar D, Ronald PC: A genome-wide survey of switchgrass genome structure and organization. PloS One. 2012, 7 (4): e33892-
Ersoz ES, Wright MH, Pangilinan JL, Sheehan MJ, Tobias C, Casler MD, Buckler ES, Costich DE: SNP Discovery with EST and NextGen sequencing in switchgrass (Panicum virgatum L.). PloS One. 2012, 7 (9): e44112-
Wang YX, Zeng X, Iyer NJ, Bryant DW, Mockler TC, Mahalingam R: Exploring the switchgrass transcriptome using second-generation sequencing technology. PloS One. 2012, 7 (3): e34225-
Young HA, Lanzatella CL, Sarath G, Tobias CM: Chloroplast genome variation in upland and lowland switchgrass. PloS One. 2011, 6 (8): e23980-
Zhang JY, Lee YC, Torres-Jerez I, Wang M, Yin Y, Chou WC, He J, Shen H, Srivastava AC, Pennacchio C: Development of an integrated transcript sequence database and a gene expression atlas for gene discovery and analysis in switchgrass (Panicum virgatum L.). Plant J. 2013, 74: 160-173.
Fu C, Sunkar R, Zhou C, Shen H, Zhang JY, Matts J, Wolf J, Mann DG, Stewart CN, Tang Y: Overexpression of miR156 in switchgrass (Panicum virgatum L.) results in various morphological alterations and leads to improved biomass production. Plant Biotechnol J. 2012, 10: 443-452.
Mann DG, Lafayette PR, Abercrombie LL, King ZR, Mazarei M, Halter MC, Poovaiah CR, Baxter H, Shen H, Dixon RA: Gateway-compatible vectors for high-throughput gene functional analysis in switchgrass (Panicum virgatum L.) and other monocot species. Plant Biotechnol J. 2012, 10: 226-236.
Xi Y, Ge Y, Wang ZY: Genetic transformation of switchgrass. Methods Mol Biol. 2009, 581: 53-59.
Kakani VG, Reddy KR: Temperature response of C4 species big bluestem (Andropogon gerardii) is modified by growing carbon dioxide concentration. Environ Exp Bot. 2007, 61: 281-290.
Kakani VG, Surabhi GK, Reddy KR: Photosynthesis and fluorescence responses of C4 plant an- dropogon gerardii acclimated to temperature and carbon dioxide. Photosynthetica. 2009, 46: 420-430.
Kandel TP, Wu Y, Kakani VG: Growth and yield responses of switchgrass ecotypes to temperature. Am J Plant Sci. 2013, 4: 1173-1180.
Peck SC, Teisberg TJ: CETA: a model for carbon emissions trajectory assessment. Energy J. 1992, 13: 55-77.
Jagadish SV, Craufurd PQ, Wheeler TR: High temperature stress and spikelet fertility in rice (Oryza sativa L.). J Exp Bot. 2007, 58: 1627-1635.
Qu A-L, Ding Y-F, Jiang Q, Zhu C: Molecular mechanisms of the plant heat stress response. Biochem Biophys Res Commun. 2013, 432: 203-207.
Semenov MA, Halford NG: Identifying target traits and molecular mechanisms for wheat breeding under a changing climate. J Exp Bot. 2009, 60: 2791-2804.
Behrman KD, Kiniry JR, Winchell M, Juenger TE, Keitt TH: Spatial forecasting of switchgrass productivity under current and future climate change scenarios. Ecol Appl. 2013, 23: 73-85.
Casler MD, Vogel KP, Taliaferro CM, Wynia RL: Latitudinal adaptation of switchgrass populations. Crop Sci. 2004, 44: 293-303.
Casler MD, Boe AR: Cultivar X environment interactions in switchgrass. Crop Sci. 2003, 43: 2226-2233.
Balasko JA, Smith D: Influence of temperature and nitrogen fertilization on the growth and composition of switchgrass (Panicum virgatum L.) and timothy (Phleum pratense L.) at anthesis. Agron J. 1971, 63: 853-857.
Hartman JC, Nipper JB: Physiological and growth responses of switchgrass (Panicum virgatum L.) in native stands under passive Air temperature manipulation. GCB Bioenergy. 2012, doi:10.1111/j.1757-1707.2012.01204.x
Finka A, Mattoo RU, Goloubinoff P: Meta-analysis of heat- and chemically upregulated chaperone genes in plant and human cells. Cell Stress Chaperones. 2011, 16: 15-31.
Qin DD, Wu HY, Peng HR, Yao YY, Ni ZF, Li ZX, Zhou CL, Sun QX: Heat stress-responsive transcriptome analysis in heat susceptible and tolerant wheat (Triticum aestivum L.) by using wheat genome array. BMC Genomics. 2008, 9: 432-doi:10.1186/1471-2164-9-432
Larkindale J, Vierling E: Core genome responses involved in acclimation to high temperature. Plant Physiol. 2008, 146: 748-761.
Fernandes J, Morrow DJ, Casati P, Walbot V: Distinctive transcriptome responses to adverse environmental conditions in Zea mays L. Plant Biotechnol J. 2008, 6: 782-798.
Zhang XW, Li JP, Liu AL, Zou J, Zhou XY, Xiang JH, Rerksiri W, Peng Y, Xiong XY, Chen XB: Expression profile in rice panicle: insights into heat response mechanism at reproductive stage. PloS One. 2012, 7 (11): e49652-doi:10.1371/journal.pone.0049652
Du Z, Zhou X, Ling Y, Zhang Z, Su Z: agriGO: a GO analysis toolkit for the agricultural community. Nucleic Acids Res. 2010, 38 (Web Server issue): W64-70.
McClung CR, Davis SJ: Ambient thermometers in plants: from physiological outputs towards mechanisms of thermal sensing. Curr Biol. 2010, 20: R1086-1092.
Ruelland E, Zachowski A: How plants sense temperature. Environ Exp Bot. 2010, 69: 225-232.
Mittler R, Finka A, Goloubinoff P: How do plants feel the heat?. Trends Biochem Sci. 2012, 37: 118-125.
Queitsch C, Hong SW, Vierling E, Lindquist S: Heat shock protein 101 plays a crucial role in thermotolerance in Arabidopsis. Plant Cell. 2000, 12: 479-492.
Su PH, Li HM: Arabidopsis stromal 70-kD heat shock proteins are essential for plant development and important for thermotolerance of germinating seeds. Plant Physiol. 2008, 146: 1231-1241.
Yamada K, Fukao Y, Hayashi M, Fukazawa M, Suzuki I, Nishimura M: Cytosolic HSP90 regulates the heat shock response that is responsible for heat acclimation in Arabidopsis thaliana. J Biol Chem. 2007, 282: 37794-37804.
Dafny-Yelin M, Tzfira T, Vainstein A, Adam Z: Non-redundant functions of sHSP-CIs in acquired thermotolerance and their role in early seed development in Arabidopsis. Plant Mol Biol. 2008, 67: 363-373.
Kotak S, Larkindale J, Lee U, Von Koskull-Doring P, Vierling E, Scharf KD: Complexity of the heat stress response in plants. Curr Opin Plant Biol. 2007, 10: 310-316.
Hu C, Lin SY, Chi WT, Charng YY: Recent gene duplication and subfunctionalization produced a mitochondrial GrpE, the nucleotide exchange factor of the Hsp70 complex, specialized in thermotolerance to chronic heat stress in Arabidopsis. Plant Physiol. 2012, 158: 747-758.
Mahalingam R, Fedoroff N: Stress response, cell death and signalling: the many faces of reactive oxygen species. Physiol Plant. 2003, 119: 56-68.
Volkov RA, Panchuk II, Mullineaux PM, Schoffl F: Heat stress-induced H2O2 is required for effective expression of heat shock genes in Arabidopsis. Plant Mol Biol. 2006, 61: 733-746.
Morimoto RI, Santoro MG: Stress-inducible responses and heat shock proteins: new pharmacologic targets for cytoprotection. Nat Biotechnol. 1998, 16: 833-838.
Apel K, Hirt H: Reactive oxygen species: metabolism, oxidative stress, and signal transduction. Ann Rev Plant Biol. 2004, 55: 373-399.
Tipping AJ, McPherson MJ: Cloning and molecular analysis of the pea seedling copper amine oxidase. J Biol Chem. 1995, 270: 16939-16946.
Miller G, Suzuki N, Rizhsky L, Hegie A, Koussevitzky S, Mittler R: Double mutants deficient in cytosolic and thylakoid ascorbate peroxidase reveal a complex mode of interaction between reactive oxygen species, plant development, and response to abiotic stresses. Plant Physiol. 2007, 144: 1777-1785.
Davletova S, Rizhsky L, Liang HJ, Zhong SQ, Oliver DJ, Coutu J, Shulaev V, Schlauch K, Mittler R: Cytosolic ascorbate peroxidase 1 is a central component of the reactive oxygen gene network of Arabidopsis. Plant Cell. 2005, 17: 268-281.
Laugier E, Tarrago L, Dos Santos CV, Eymery F, Havaux M, Rey P: Arabidopsis thaliana plastidial methionine sulfoxide reductases B, MSRBs, account for most leaf peptide MSR activity and are essential for growth under environmental constraints through a role in the preservation of photosystem antennae. Plant J. 2010, 61: 271-282.
Irmler S, Schroder G, St-Pierre B, Crouch NP, Hotze M, Schmidt J, Strack D, Matern U, Schroder J: Indole alkaloid biosynthesis in Catharanthus roseus: new enzyme activities and identification of cytochrome P450CYP72A1 as secologanin synthase. Plant J. 2000, 24: 797-804.
Dietz KJ, Jacob S, Oelze ML, Laxa M, Tognetti V, De Miranda SM, Baier M, Finkemeier I: The function of peroxiredoxins in plant organelle redox metabolism. J Exp Bot. 2006, 57: 1697-1709.
Durand TC, Sergeant K, Carpin S, Label P, Morabito D, Hausman JF, Renaut J: Screening for changes in leaf and cambial proteome of Populus tremula x P. alba under different heat constraints. J Plant Physiol. 2012, 169: 1698-1718.
Park UH, Han HS, Um E, An XH, Kim EJ, Um SJ: Redox regulation of transcriptional activity of retinoic acid receptor by thioredoxin glutathione reductase (TGR). Biochem Biophys Res Commun. 2009, 390: 241-246.
Yamazaki D, Motohashi K, Kasama T, Hara Y, Hisabori T: Target proteins of the cytosolic thioredoxins in Arabidopsis thaliana. Plant Cell Physiol. 2004, 45: 18-27.
Gholizadeh A, Kohnehrouz BB: Identification of DUF538 cDNA clone from Celosia cristata expressed sequences of nonstressed and stressed leaves. Russ J Plant Physl+. 2010, 57: 247-252.
Gholizadeh A: Heterologous expression of stress-responsive DUF538 domain containing protein and its morpho-biochemical consequences. Protein J. 2011, 30: 351-358.
Phipps BM, Hoffmann A, Stetter KO, Baumeister W: A novel atpase complex selectively accumulated upon heat-shock is a major cellular-component of thermophilic archaebacteria. Embo J. 1991, 10: 1711-1722.
Phipps BM, Typke D, Hegerl R, Volker S, Hoffmann A, Stetter KO, Baumeister W: Structure of a molecular chaperone from a thermophilic archaebacterium. Nature. 1993, 361: 475-477.
Hill JE, Hemmingsen SM: Arabidopsis thaliana type I and II chaperonins. Cell Stress Chaperones. 2001, 6: 190-200.
Gong ZZ, Dong CH, Lee H, Zhu JH, Xiong LM, Gong DM, Stevenson B, Zhu JK: A DEAD box RNA helicase is essential for mRNA export and important for development and stress responses in Arabidopsis. Plant Cell. 2005, 17: 256-267.
Gong ZZ, Lee H, Xiong LM, Jagendorf A, Stevenson B, Zhu JK: RNA helicase-like protein as an early regulator of transcription factors for plant chilling and freezing tolerance. PNAS. 2002, 99: 11507-11512.
Chung E, Cho CW, Yun BH, Choi HK, So HA, Lee SW, Lee JH: Molecular cloning and characterization of the soybean DEAD-box RNA helicase gene induced by low temperature and high salinity stress. Gene. 2009, 443: 91-99.
Seki M, Narusaka M, Abe H, Kasuga M, Yamaguchi-Shinozaki K, Carninci P, Hayashizaki Y, Shinozaki K: Monitoring the expression pattern of 1300 Arabidopsis genes under drought and cold stresses by using a full-length cDNA microarray. Plant Cell. 2001, 13: 61-72.
Nakamura T, Muramoto Y, Yokota S, Ueda A, Takabe T: Structural and transcriptional characterization of a salt-responsive gene encoding putative ATP-dependent RNA helicase in barley. Plant Sci. 2004, 167: 63-70.
Macovei A, Vaid N, Tula S, Tuteja N: A new DEAD-box helicase ATP-binding protein (OsABP) from rice is responsive to abiotic stress. Plant Signal Behav. 2012, 7: 1138-1143.
Reddy AS, Ali GS, Celesnik H, Day IS: Coping with stresses: roles of calcium- and calcium/calmodulin-regulated gene expression. Plant Cell. 2011, 23: 2010-2032.
Sangwan V, Orvar BL, Beyerly J, Hirt H, Dhindsa RS: Opposite changes in membrane fluidity mimic cold and heat stress activation of distinct plant MAP kinase pathways. Plant J. 2002, 31: 629-638.
Zhang W, Zhou RG, Gao YJ, Zheng SZ, Xu P, Zhang SQ, Sun DY: Molecular and genetic evidence for the key role of AtCaM3 in heat-shock signal transduction in Arabidopsis. Plant Physiol. 2009, 149: 1773-1784.
Liu HT, Gao F, Li GL, Han JL, Liu DL, Sun DY, Zhou RG: The calmodulin-binding protein kinase 3 is part of heat-shock signal transduction in Arabidopsis thaliana. Plant J. 2008, 55: 760-773.
Li S, Zhou X, Chen L, Huang W, Yu D: Functional characterization of Arabidopsis thaliana WRKY39 in heat stress. Mol Cells. 2010, 29: 475-483.
Liu HC, Liao HT, Charng YY: The role of class A1 heat shock factors (HSFA1s) in response to heat and other stresses in Arabidopsis. Plant Cell Environ. 2011, 34: 738-751.
Christensen A, Svensson K, Persson S, Jung J, Michalak M, Widell S, Sommarin M: Functional characterization of Arabidopsis calreticulin1a: A key alleviator of endoplasmic reticulum stress. Plant Cell Physiol. 2008, 49: 912-924.
Jia XY, He LH, Jing RL, Li RZ: Calreticulin: conserved protein and diverse functions in plants. Physiol Plantarum. 2009, 136: 127-138.
Supek F, Bosnjak M, Skunca N, Smuc T: REVIGO summarizes and visualizes long lists of gene ontology terms. PloS One. 2011, 6: e21800-
Nadeau K, Das A, Walsh CT: Hsp90 chaperonins possess ATPase activity and bind heat shock transcription factors and peptidyl prolyl isomerases. J Biol Chem. 1993, 268: 1479-1487.
Waters ER: The evolution, function, structure, and expression of the plant sHSPs. J Exp Bot. 2013, 64: 391-403.
Lund AA, Blum PH, Bhattramakki D, Elthon TE: Heat-stress response of maize mitochondria. Plant Physiol. 1998, 116: 1097-1110.
Retzlaff M, Hagn F, Mitschke L, Hessling M, Gugel F, Kessler H, Richter K, Buchner J: Asymmetric activation of the hsp90 dimer by its cochaperone aha1. Mol Cell. 2010, 37: 344-354.
Ran FL, Gadura N, Michels CA: Hsp90 Cochaperone Aha1 is a negative regulator of the saccharomyces MAL activator and acts early in the chaperone activation pathway. J Biol Chem. 2010, 285: 13850-13862.
Padidam M, Reddy VS, Beachy RN, Fauquet CM: Molecular characterization of a plant mitochondrial chaperone GrpE. Plant Mol Biol. 1999, 39: 871-881.
Nieto-Sotelo J, Martinez LM, Ponce G, Cassab GI, Alagon A, Meeley RB, Ribaut JM, Yang RY: Maize HSP101 plays important roles in both induced and basal thermotolerance and primary root growth. Plant Cell. 2002, 14: 1621-1633.
Nystrom T, Neidhardt FC: Expression and role of the universal stress protein, UspA, of Escherichia coli during growth arrest. Mol Microbiol. 1994, 11: 537-544.
Credali A, Garcia-Calderon M, Dam S, Perry J, Diaz-Quintana A, Parniske M, Wang TL, Stougaard J, Vega JM, Marquez AJ: The K + -dependent asparaginase, NSE1, is crucial for plant growth and seed production in lotus japonicus. Plant Cell Physiol. 2013, 54: 107-118.
Cho CW, Lee HJ, Chung E, Kim KM, Kim JEHJI, Chung J, Ma YZ, Fukui K, Lee DW, Kim DH, et al: Molecular characterization of the soybean L-asparaginase gene induced by low temperature stress. Mol Cells. 2007, 23: 280-286.
Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U, Speed TP: Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics. 2003, 4: 249-264.
Dash S, Van Hemert J, Hong L, Wise RP, Dickerson JA: PLEXdb: gene expression resources for plants and plant pathogens. Nucleic Acids Res. 2012, 40 (Database issue): D1194-1201.
This research is supported by the National Science Foundation under Grant No. EPS-0814361. The authors thank the Oklahoma State University Recombinant DNA core facility for technical assistance in this project.
RM received grant support. RM, YW and YL designed the experiment. VGK grew switchgrass plants, performed the heat stress treatments and collected the tissue samples. YT performed the affymetrix genechip hybridizations, conducted the preliminary data analysis and assisted in submitting the microarray data to ArrayExpress. YL and YW validated gene expression. YL and RM analyzed the array results and prepared the manuscript. All authors reviewed and approved the final manuscript.
Electronic supplementary material
About this article
Cite this article
Li, Y., Wang, Y., Tang, Y. et al. Transcriptome analysis of heat stress response in switchgrass (Panicum virgatumL.). BMC Plant Biol 13, 153 (2013) doi:10.1186/1471-2229-13-153
- Heat stress