Transcriptome analysis of heat stress response in switchgrass (Panicum virgatumL.)
© Li et al.; licensee BioMed Central Ltd. 2013
Received: 7 July 2013
Accepted: 3 October 2013
Published: 6 October 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.
KeywordsSwitchgrass Biofuel Microarray Heat stress Transcriptome
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
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.
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/).
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.
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 common heat responsive genes in switchgrass, rice, maize and wheat
PviUT ID annotation
CS 1 sh
CS 1 h
CS 24 sh
CS 24 h
TAM 1 sh
TAM 1 h
TAM 24 sh
TAM 24 h
AP13CTG24376 universal stress protein
KanlCTG39278 DnaK family protein
AP13ISTG73695 glycosyl hydrolases family 17
KanlSGLT49533 activator of 90 kDa HSP ATPase homolog
AP13CTG04492 chaperone protein clpB 1
KanlCTG05390 heat shock 22 kDa protein, mitochondrial precursor
AP13CTG00593 T-complex protein
KanlCTG37554 Putative transposon
AP13CTG06779 T-complex protein
AP13CTG59854 expressed protein
AP13ISTG34202 activator of 90 kDa HSP ATPase homolog
AP13CTG25439 peptidyl-prolyl isomerase
AlamCTG07708 co-chaperone GrpE protein
AP13CTG01830 T-complex protein
AP13CTG30798 HIT zinc finger domain protein
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.
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.
- Hitchcock AS, Chase A: Manual of the Grass of the United State. Washington, D. C: U S Government Printing Office; 1950.Google Scholar
- 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.View ArticleGoogle Scholar
- Porter CL: An analysis of variation between upland and lowland switchgrass panicum virgatum L in Central Oklahoma. Ecology. 1966, 47: 980-992.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Gunter LE, Tuskan GA, Wullschleger SD: Diversity among populations of switchgrass based on RAPD markers. Crop Sci. 1996, 36: 1017-1022.View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Narasimhamoorthy B, Saha MC, Swaller T, Bouton JH: Genetic diversity in switchgrass collections assessed by EST-SSR markers. Bioenergy Res. 2008, 1: 136-146.View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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-PubMedPubMed CentralView ArticleGoogle Scholar
- 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-PubMedPubMed CentralView ArticleGoogle Scholar
- 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-PubMedPubMed CentralView ArticleGoogle Scholar
- 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-PubMedPubMed CentralView ArticleGoogle Scholar
- Young HA, Lanzatella CL, Sarath G, Tobias CM: Chloroplast genome variation in upland and lowland switchgrass. PloS One. 2011, 6 (8): e23980-PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Xi Y, Ge Y, Wang ZY: Genetic transformation of switchgrass. Methods Mol Biol. 2009, 581: 53-59.PubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Kandel TP, Wu Y, Kakani VG: Growth and yield responses of switchgrass ecotypes to temperature. Am J Plant Sci. 2013, 4: 1173-1180.View ArticleGoogle Scholar
- Peck SC, Teisberg TJ: CETA: a model for carbon emissions trajectory assessment. Energy J. 1992, 13: 55-77.View ArticleGoogle Scholar
- Jagadish SV, Craufurd PQ, Wheeler TR: High temperature stress and spikelet fertility in rice (Oryza sativa L.). J Exp Bot. 2007, 58: 1627-1635.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Semenov MA, Halford NG: Identifying target traits and molecular mechanisms for wheat breeding under a changing climate. J Exp Bot. 2009, 60: 2791-2804.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Casler MD, Vogel KP, Taliaferro CM, Wynia RL: Latitudinal adaptation of switchgrass populations. Crop Sci. 2004, 44: 293-303.Google Scholar
- Casler MD, Boe AR: Cultivar X environment interactions in switchgrass. Crop Sci. 2003, 43: 2226-2233.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.xGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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-432PubMedPubMed CentralView ArticleGoogle Scholar
- Larkindale J, Vierling E: Core genome responses involved in acclimation to high temperature. Plant Physiol. 2008, 146: 748-761.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.0049652PubMedPubMed CentralView ArticleGoogle Scholar
- Du Z, Zhou X, Ling Y, Zhang Z, Su Z: agriGO: a GO analysis toolkit for the agricultural community. Nucleic Acids Res. 2010, 38 (Web Server issue): W64-70.PubMedPubMed CentralView ArticleGoogle Scholar
- McClung CR, Davis SJ: Ambient thermometers in plants: from physiological outputs towards mechanisms of thermal sensing. Curr Biol. 2010, 20: R1086-1092.PubMedView ArticleGoogle Scholar
- Ruelland E, Zachowski A: How plants sense temperature. Environ Exp Bot. 2010, 69: 225-232.View ArticleGoogle Scholar
- Mittler R, Finka A, Goloubinoff P: How do plants feel the heat?. Trends Biochem Sci. 2012, 37: 118-125.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- Mahalingam R, Fedoroff N: Stress response, cell death and signalling: the many faces of reactive oxygen species. Physiol Plant. 2003, 119: 56-68.View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Morimoto RI, Santoro MG: Stress-inducible responses and heat shock proteins: new pharmacologic targets for cytoprotection. Nat Biotechnol. 1998, 16: 833-838.PubMedView ArticleGoogle Scholar
- Apel K, Hirt H: Reactive oxygen species: metabolism, oxidative stress, and signal transduction. Ann Rev Plant Biol. 2004, 55: 373-399.View ArticleGoogle Scholar
- Tipping AJ, McPherson MJ: Cloning and molecular analysis of the pea seedling copper amine oxidase. J Biol Chem. 1995, 270: 16939-16946.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Gholizadeh A: Heterologous expression of stress-responsive DUF538 domain containing protein and its morpho-biochemical consequences. Protein J. 2011, 30: 351-358.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralGoogle Scholar
- 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.View ArticleGoogle Scholar
- Hill JE, Hemmingsen SM: Arabidopsis thaliana type I and II chaperonins. Cell Stress Chaperones. 2001, 6: 190-200.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Jia XY, He LH, Jing RL, Li RZ: Calreticulin: conserved protein and diverse functions in plants. Physiol Plantarum. 2009, 136: 127-138.View ArticleGoogle Scholar
- Supek F, Bosnjak M, Skunca N, Smuc T: REVIGO summarizes and visualizes long lists of gene ontology terms. PloS One. 2011, 6: e21800-PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedGoogle Scholar
- Waters ER: The evolution, function, structure, and expression of the plant sHSPs. J Exp Bot. 2013, 64: 391-403.PubMedView ArticleGoogle Scholar
- Lund AA, Blum PH, Bhattramakki D, Elthon TE: Heat-stress response of maize mitochondria. Plant Physiol. 1998, 116: 1097-1110.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- Padidam M, Reddy VS, Beachy RN, Fauquet CM: Molecular characterization of a plant mitochondrial chaperone GrpE. Plant Mol Biol. 1999, 39: 871-881.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.