Genome-wide analysis of the Populus Hsp90 gene family reveals differential expression patterns, localization, and heat stress responses
© Zhang et al.; licensee BioMed Central Ltd. 2013
Received: 25 January 2013
Accepted: 30 July 2013
Published: 5 August 2013
Members of the heat shock protein 90 (Hsp90) class of proteins are evolutionarily conserved molecular chaperones. They are involved in protein folding, assembly, stabilization, activation, and degradation in many normal cellular processes and under stress conditions. Unlike many other well-characterized molecular chaperones, Hsp90s play key roles in signal transduction, cell-cycle control, genomic silencing, and protein trafficking. However, no systematic analysis of genome organization, gene structure, and expression compendium has been performed in the Populus model tree genus to date.
We performed a comprehensive analysis of the Populus Hsp90 gene family and identified 10 Populus Hsp90 genes, which were phylogenetically clustered into two major groups. Gene structure and motif composition are relatively conserved in each group. In Populus trichocarpa, we identified three paralogous pairs, among which the PtHsp90-5a/PtHsp90-5b paralogous pair might be created by duplication of a genome segment. Subcellular localization analysis shows that PtHsp90 members are localized in different subcellular compartments. PtHsp90-3 is localized both in the nucleus and in the cytoplasm, PtHsp90-5a and PtHsp90-5b are in chloroplasts, and PtHsp90-7 is in the endoplasmic reticulum (ER). Furthermore, microarray and semi-quantitative real-time RT-PCR analyses show that a number of Populus Hsp90 genes are differentially expressed upon exposure to various stresses.
The gene structure and motif composition of PtHsp90s are highly conserved among group members, suggesting that members of the same group may also have conserved functions. Microarray and RT-PCR analyses show that most PtHsp90s were induced by various stresses, including heat stress. Collectively, these observations lay the foundation for future efforts to unravel the biological roles of PtHsp90 genes.
KeywordsExpression analysis Gene family Gene structure Hsp90 Phylogenetic analysis Populus
Plants are exposed to various environmental stresses. Primary stresses such as high light intensity, heat shock, drought, chilling, salinity, and chemical pollutants act simultaneously on plants, causing cell injury and producing secondary stresses such as osmotic and oxidative stresses . Plants cannot avoid exposure to these factors, but adapted morphologically and physiologically by some mechanisms. Biosynthesis of many proteins called “stress proteins” is induced to protect cells from these harmful stimuli .
Heat shock proteins (Hsps) are responsible for protein folding, assembly, translocation, and degradation in many normal cellular processes. They stabilize proteins and membranes, and can assist in protein refolding under stress conditions. They also play a crucial role in protecting plants from stresses by reestablishing normal protein conformations and thus cellular homeostasis . Plant Hsps are classified into five families according to their molecular size: Hsp100, Hsp90, Hsp70, Hsp60, and small Hsps (sHsps). They have been well characterized in a few model plants such as the tomato, Arabidopsis, and rice [3, 4].
Hsp90s are a class of chaperone proteins that are highly conserved in prokaryotes and all eukaryotes. They are the major species of molecular chaperones and require ATP for their functions . Although Hsp90s are expressed in most organisms, their expression increases in response to stresses. Distinct from many other well-characterized molecular chaperones, Hsp90s display considerable specificity for their client proteins. Most of their known substrates are signal-transduction proteins such as steroid hormone receptors and signaling kinases . Although the major function of Hsp90s is to assist protein folding, they play key roles in signal transduction, cell-cycle control, protein degradation, genomic silencing, and protein trafficking [5, 6]. Expression of Hsp90 in Arabidopsis is developmentally regulated and is responsive to heat, cold, salinity, heavy metals, phytohormones, and light and dark transitions [4, 7]. Tobacco NbHsp90-1 and Arabidopsis AtHsp90-2 confer pathogen resistance by reacting to resistance proteins (R proteins), which are signal receptors from the pathogen [8, 9]. In addition, Hsp90s interact with the 26S proteasome and play a key role in its ATP-dependent assembly and maintenance in budding yeast . To fulfill their cellular roles, Hsp90s cooperate with other chaperones to form a multiprotein chaperone complex . Moreover, Hsp90s also act as buffers to phenotypic changes and are portrayed as “capacitors for evolution” .
Hsp90s are encoded by multiple genes. They consist of conserved N-terminal and C-terminal domains that are joined by a charged linker region that varies in length. Genes encoding cytosol-, ER-, and plastid-localized Hsp90 proteins have been characterized in several plant species . In the Arabidopsis genome, seven Hsp90 family members have been identified. Sequence analyses of Arabidopsis Hsp90 family genes have revealed two major subfamilies. AtHsp90-1–4 proteins containing the C-terminal pentapeptide MEEVD form the cytoplasmic subfamily; AtHsp90-5–7 form the other subfamily. AtHsp90-5 and AtHsp90-7 are localized in chloroplasts  and the endoplasmic reticulum (ER) , respectively. AtHsp90-6 is localized in mitochondria . Overexpression of cytosolic AtHsp90-2, chloroplast-localized AtHsp90-5, and ER-localized AtHsp90-7 reduces tolerance to salt and drought stresses, but improves tolerance to high concentrations of Ca2+. The induction of ABA-responsive genes is delayed by overexpression of cytosolic AtHsp90-2, but is hardly affected by overexpression of AtHsp90-5 and AtHsp90-7 under conditions of salt and drought stress, which implies that different cellular compartment-localized Hsp90s in Arabidopsis might contribute to responses to abiotic stresses by different functional mechanisms, probably through ABA- or Ca2+-dependent pathways .
The Populus genus comprises woody plants that are important to humans and animals. Completion of the P. trichocarpa genome sequence in 2006 rendered it a model species for research on trees , providing an opportunity to analyze and further understand Hsp90s. To determine the structure-function relationship of Hsp90s in the Populus genus, we performed detailed systematic analyses of genome organization, gene structure, and expression compendium. We report the comprehensive genomic identification and phylogenetic analysis of all 10 members of the Hsp90 gene family in the Populus genus, as well as their expression profiles in different tissues and their responses under heat stress. Our results provide a framework for further functional investigations of these genes.
Results and discussion
Identification of the Hsp90 gene family in P. trichocarpa and other plant species
Hsp90 genes families in Arabidopsis , Populus , and rice
Mol. Wt. (kDa), Length (aa), pI
Gene length, ORF, introns
81.18, 705, 4.95
2606, 2118, 3
N: 6, C: 5, Ch: 2
80.06, 699, 4.95
2511, 2187, 2
C: 11, P: 2
80.05, 699, 4.95
2502, 2100, 2
C: 9, P: 2, N: 1, M: 1
80.14, 699, 4.96
2466, 2100, 2
C: 9, P: 4
88.66, 780, 4.93
3927, 2343, 18
Ch: 12, N: 1
90.57, 799, 5.26
4360, 2400, 19
M: 9.5, Ch_M: 6.5, Ch: 2.5
94.2, 823, 4.94
3950, 2472, 14
E.R.: 7, V: 4, N: 1
Chr04: 6143799 - 6146941(-)
80.72, 703, 5
2566, 2112, 3
C: 7, N: 3, P: 2, M: 1
Chr17: 15369800 - 15373207(+)
80.75, 703, 4.98
2696, 2112, 3
C: 8, N: 2, P: 2, M: 1
Chr01: 49936276 - 49939430(+)
80.03, 699, 4.95
2767, 2100, 2
C: 6, P: 2, E.R.: 2, Ch: 1, N: 1, M: 1
Chr16: 167702 - 171110(+)
79.95, 698, 4.94
2906, 2097, 3
N: 7, C: 5, Ch: 2
Chr01: 29295547 - 29298729(+)
80.86, 706, 4.88
2734, 2121, 2
C: 7, N: 2, P: 2, Ch: 1, M: 1
Chr06: 219117 - 222439(+)
80.01, 699, 4.91
2734, 2100, 2
N: 4, C: 4, P: 3, Ch: 1, V: 1
Chr08: 7190800 - 7197203(+)
90.16, 791, 4.92
5884, 2376, 17
Ch: 10, N: 2, P: 1
Chr10: 15004039 - 15010595(-)
90.15, 793, 4.96
5723, 2382, 17
Chr14: 13043761 - 13049852(-)
82.95, 723, 5.53
4850, 2172, 13
M: 7.5, Ch_M: 7, Ch: 5.5
Chr05: 24712685 - 24717968(-)
94.05, 823, 4.85
4926, 2472, 14
E.R.: 5, V: 5, P: 2, Ch: 1
Chr4: 483241 - 486065 (+)
80.25, 703, 5.04
2353, 2112, 2
C: 11, Ch: 1, P: 1
Chr8: 24719086 - 24723553 (-)
80.19, 699, 4.99
4002, 2100, 2
C: 8, P: 2, Ch: 1, N: 1, M: 1
Chr9: 18514572 - 18518316 (-)
80.2, 699, 4.97
3232, 2100, 2
N: 7, C: 4, Ch: 2
Chr9: 18535746 - 18541109 (-)
94.2, 830, 5.15
5364, 2493, 3
N: 8, C: 3, Ch: 1, P: 1
Chr8: 24124838 - 24129488 (-)
85.85, 761, 4.97
4285, 2286, 19
Ch: 12.5, Ch_M: 7
Chr9: 18150618 - 18155512 (-)
89.22, 791, 5.03
4359, 2376, 18
Ch: 12.5, Ch_M: 7
Chr12: 19921576 - 19927766 (-)
91.48, 811, 5.21
5680, 2436, 19
Ch: 12, M: 1
Chr6: 30444411 - 30450497 (-)
93.04, 812, 4.89
5669, 2439, 14
E.R.: 5, V: 5, P: 2, Ch: 1
To investigate the evolutionary relationships of Hsp90 proteins from different plants, we identified Hsp90 genes from seven other plant species, including the moss Physcomitrella patens, the monocotyledonous angiosperms Oryza sativa, Sorghum bicolor, and Brachypodium distachyon, the dicotyledonous angiosperms Arabidopsis thaliana, Vitis vinifera, and Medicago truncatula. All angiosperm genomes, as well as the moss genome, contain Hsp90 genes. The number of Hsp90s identified is seven in A. thaliana, ten in P. trichocarpa, five in V. vinifera, five in M. truncatula, eight in O. sativa, seven in S. bicolor, eight in B. distachyon, and ten in P. patens. Additional file 3 provides a complete list of all Hsp90 genes identified in the present study.
Phylogenetic analyses of the Hsp90 gene family
Numbers of Hsp90s within each plant species
It is more accurate to reflect an evolutionary relationship by using conserved domain sequences . Therefore we also constructed the phylogenetic tree with the conserved Hsp90 motif sequences from Populus, Arabidopsis, and rice using the maximum likelihood method with 1000 bootstrap replicates (Additional file 6). The resulted phylogenetic tree is consistent with the one generated based on the full length protein sequences.
Gene structure and conserved motifs of Hsp90 genes in Populus, Arabidopsis, and rice
To further investigate the structural diversity of Hsp90 genes in Populus, Arabidopsis, and rice, we first constructed a separate phylogenetic tree using the full-length Hsp90 protein sequences from these three species. The Hsp90 proteins are classified into two groups as described above (Figure 2A). Then we analyzed the exon/intron organization in the coding sequence of each Hsp90 gene (Figure 2B). In general, the positions of some spliceosomal introns are conserved in orthologous genes. In many cases, conservation of exon/intron organization or gene structure in paralogous genes is high and sufficient to reveal the evolutionary relationship between introns . In the present study, Hsp90 gene family members within the same group shared similar gene structures in terms of intron number or exon length (Figure 2B). Hsp90 group I comprises the cytosolic Hsp90s whose members have two or three introns, while group II comprises organelle-type Hsp90s, which have 13–19 introns (Figure 2 and Table 1). The gene structure difference between group I and group II Hsp90 might associate with their functions in different biological processes in subcellular compartments. We also investigated intron phases with respect to codons. The intron phases are remarkably well conserved among group members, while the intron arrangements and intron phases are strikingly distinct between groups (Figure 2). This may lend support to the results of phylogenetic and genome duplication analyses. We further examined the exon/intron organization of paralogous pairs of Hsp90 genes to explore traceable intron gain or loss within these genes. Three paralogous pairs in Populus (PtHsp90-1a/1b, PtHsp90-4a/4b, and PtHsp90-5a/5b) show conserved exon/intron structures in terms of intron number or gene length, while OsHsp90-5a shows a single intron gain event during the structural evolution of the OsHsp90-5a/5b paralogous pair. Interestingly, OsHsp90-4 has an additional C-terminal exon compared to other members of the group I Hsp90.
Next, we predicted the major domains of these proteins in all three species using Pfam and CDD . All of the proteins contain a HATPase_c superfamily domain and a Hsp90 family domain (Additional file 1). Although the tools we used are suitable for defining the presence or absence of recognizable domains, they are unable to recognize smaller individual motifs and more divergent patterns. Thus, the program MEME was used to further study the diversification of these proteins . Twenty distinct motifs were identified (Figure 2C). Details of the 20 motifs are presented in Additional file 5. Most of the closely related members have common motif composition, suggesting possible functional similarity among these Hsp90 proteins (Figure 2C). Motif 2 and 6 (corresponding to the HATPase_c superfamily domain at the N-terminus) are found in all Hsp90 proteins from the species we examined. It has reported that both ATP binding and hydrolysis are required for Hsp90 function in vivo[23, 24]. Noticeably, motif 20 (representing the LEA_6 subdomain) is only found in OsHsp90-4. This additional LEA_6 subdomain might explain the specific ability of OsHsp90-4 to acclimatize to various stresses.
Chromosomal location and gene duplication of Hsp90 genes in Populus, Arabidopsis, and rice
Chromosomal mapping of the gene loci shows that the 10 PtHsp90 genes are distributed unevenly among nine chromosomes (Additional file 7). Two PtHsp90 genes are localized on chromosome I, and one is localized on each of chromosome IV, V, VI, VIII, X, XIV, XVI, and XVII. Gene duplication events are thought to occur frequently in organismal evolution [25, 26]. Previous studies report that the Populus genome has experienced at least two genome-wide duplication events (eurosid and salicoid), followed by a series of chromosomal reorganizations involving reciprocal tandem/terminal fusions and translocations . To investigate the possible relationship between Hsp90 genes and segmental chromosome duplication, we also compared the locations of Hsp90 genes in duplicated chromosomal blocks that were previously identified in Populus, Arabidopsis, and rice [17, 27, 28]. Their distributions are shown in Additional file 7 (Populus), Additional file 8 (Arabidopsis), and Additional file 9 (rice). The results suggest that segmental duplication and transposition events are not the major factors that led to the expansion of the Populus Hsp90 gene family. It may be that dynamic changes occurred following segmental duplication and led to the loss of many of the duplicated Hsp90 genes.
A search for duplicated genes using the Plant Genome Duplication Database (PGDD; http://chibba.agtec.uga.edu/duplication/) revealed the existence of three gene pairs (PtHsp90-1a/PtHsp90-1b, PtHsp90-3/PtHsp90-4b, and PtHsp90-5a/PtHsp90-5b) in P. trichocarpa (Additional file 10A) and two pairs (OsHsp90-2/OsHsp90-3 and OsHsp90-5a/OsHsp90-5b) in O. sativa (Additional file 10B). Interestingly, PtHsp90-4a was not assigned as a duplicated gene with PtHsp90-3 and PtHsp90-4b, indicating that PtHsp90-4a had experienced intensive recombination events after the recent duplication with PtHsp90-4b, which led to the great divergence in its adjacent regions. Of the three Hsp90 pairs in Populus that we examined, only one pair, PtHsp90-5a/PtHsp90-5b, remained in a conserved position in segmental duplicated blocks (Additional file 7), suggesting that only this paralogous pair survived during the evolutionary process after chromosome duplication event.
Subcellular localization of Populus Hsp90 proteins
Differential expression patterns of Hsp90 genes in Populus
Differential stress responses of Hsp90 genes in Populus
To verify the expression profiles of PtHsp90 genes in response to heat stress, qRT-PCR analysis was performed for four selected PtHsp90 genes under heat stress (Additional file 12C-F). Notably, PtHsp90-3, PtHsp90-4a and PtHsp90-5a are induced 3 h after 37°C pretreatment and significantly induced 3 h after 45°C treatment in leaves (Additional file 12C-E). The expression of PtHsp90-5b is also induced by heat stress, but the induction is not that dramatic compared with that of the other PtHsp90 genes in both 37°C pretreatment (2-fold) and 45°C treatment (3-fold) (Additional file 12F). In addition, we found that the paralogous pair PtHsp90-5a/PtHsp90-5b shared the same expression profile in different tissues but were different under wounding and heat stresses (Figures 5 and 6). We then analyzed the promotors (2000 bp upstream of the start codon) of PtHsp90-5a and PtHsp90-5b using PlantCARE . The sequence of the promotors share a low sequence identity (43.3%) and two heat shock elements (HSE) exist in the promotor of PtHsp90-5a while none in PtHsp-5b (data not shown), which may contribute to the different expression pattern of the two genes. These results suggest different response mechanisms of PtHsp90 members may exist under heat stress, and provide significant insights into their functions.
We performed a comprehensive analysis of the Populus Hsp90 gene family covering phylogeny, chromosomal location, gene structure, subcellular localization, expression profiling, and heat stress responses. A total of 10 full-length Hsp90 genes were identified in the Populus genome, all of which are clustered into two distinct groups. Exon/intron structure and motif compositions are found to be relatively conserved in each subgroup. The Populus genome contains three paralogous Hsp90 gene pairs, but only PtHsp90-5a/PtHsp90-5b is located in conserved positions in duplicated blocks, suggesting that it may be derived from a segmental duplication event during evolution. Furthermore, subcellular localization analysis revealed that PtHsp90 members are localized in different organelles. In addition, comparative expression profile analysis of Populus Hsp90s revealed that Hsp90s may play various conserved roles in different biological processes in plants. Although the functions of PtHsp90s remain largely unknown and many experiments are needed to determine their precise functions, our phylogenetic and expression analyses of the Populus Hsp90 gene family establishes a solid foundation for future comprehensive functional analyses of PtHsp90s.
Database searching and sequence retrieval
To identify potential members of the Populus Hsp90 gene family, we performed multiple database searches. Published Arabidopsis Hsp90 gene sequences  were retrieved and used as queries in BLAST searches against the Poplar Genome Database (http://www.phytozome.net/poplar.php, release 3.0). BLAST searches were also performed against the poplar genomes at the National Center for Biotechnology Information (NCBI, http://www.ncbi.nlm.nih.gov) and Phytozome (http://www.phytozome.net). Rice Hsp90 gene sequences were downloaded from the Rice Genome Annotation Project Database (http://rice.plantbiology.msu.edu/, release 7). Sequences of M. truncatula, S. bicolor, B. distachyon, V. vinifera, and P. patens were downloaded from Phytozome (http://www.phytozome.net). Local BLAST searches were performed using Arabidopsis Hsp90 protein sequences as queries to identify Hsp90 sequences in these plant species. All of the sequences were manually analyzed to confirm the presence of HATPase and Hsp90 domains using InterProScan (http://www.ebi.ac.uk/Tools/pfa/iprscan/).
WoLF PSORT (http://wolfpsort.org) was used to predict protein subcellular localization. The pI and molecular weight were estimated using the Compute pI/Mw tool from ExPASy (http://web.expasy.org/compute_pi).
Multiple sequence alignment of the full-length protein sequences was performed using ClustalX2 (version 2.1) . A maximum likelihood (ML) phylogenetic tree was constructed using PhyML (v3.0) with the JTT amino acid substitution model, 1000 bootstrap replicates, estimated proportions of invariable sites, four rate categories, estimated gamma distribution parameters, and an optimized starting BIONJ tree [42, 43].
Chromosomal location and gene structure of the Hsp90 genes
The chromosomal locations of the Hsp90 genes were determined using the Populus genome browser (http://www.phytozome.net/poplar). Information on intron/exon structure was collected from the genome annotations of P. trichocarpa from NCBI and Phytozome.
Gene structure analysis
The exon and intron structures of individual Hsp90 genes were illustrated using Gene Structure Display Server (GSDS, http://gsds.cbi.pku.edu.cn/)  by aligning the cDNA sequences with the corresponding genomic DNA sequences from Phytozome.
Conserved motif analysis
Functional motifs or domains of PtHsp90 protein sequences were analyzed using PROSITE and the Conserved Domain database. MEME (http://meme.sdsc.edu)  was used to identify motifs in candidate sequences. MEME was run locally with the following parameters: number of repetitions = any, maximum number of motifs = 20, and optimum motif width = 30 to 70 residues.
Transient expression and imaging
Transient expression in Nicotiana benthamiana lower leaf epidermal cells was performed as described by Zheng et al.  with slight modifications. Plants were cultivated under short-day conditions (8 h light/16 h dark). When the agrobacterium culture reached the stationary growth phase at 28°C with agitation, cells were pelleted and resuspended in infiltration buffer (100 μM acetosyringone in 10 mM MgCl2).
Publicly available microarray data analyses
For abiotic and hormonal treatments, Affymetrix microarray data available in the NCBI GEO database under the series accession numbers GSE26199 (heat stress), GSE17230 (drought stress) and GSE16786 were analyzed [39, 46, 47]. GSE16786 is composed of the following five subsets: GSE14893 (nitrogen limitation, genotype 1979), GSE14515 (nitrogen limitation, genotype 3200), GSE16783 (1 week after leaf wounding), GSE16785 (90 h after leaf wounding), and GSE16773 (methyl jasmonate-elicited suspension cell cultures). The Affymetrix CEL files representing different abiotic and hormonal treatments were downloaded from the GEO database and preprocessed using GeneSpring GX (V11.5) software (Agilent Technologies). The data were normalized using the GCRMA algorithm and then log transformed. The averages were calculated. After normalization and log transformation of data for all of the Populus genes presented on the chip, the log signal intensity values for Populus probe IDs corresponding to the Hsp90 gene models (v1.1) were extracted as a subset for further analyses. Expression was shown as fold change in experimental treatment samples relative to control samples. Tab-delimited files for the average log signal intensity values were imported into Genesis (v1.75) to generate the heat maps .
Probe sets corresponding to PtHsp90 genes were identified using the online Probe Match tool POParray (http://aspendb.uga.edu/poparray). For probe sets matching several Populus Hsp90 gene models, only these exhibiting consistently high hybridization signals across multiple samples were considered.
Plant material and growth conditions
Plant materials were collected from clonally propagated 1-year-old hybrid poplar (P. alba × P. glandulosa) clones (84K) grown in a growth chamber under long-day conditions (16 h light/8 h dark) at 23–25°C. Poplar saplings were subjected to heat treatment. Briefly, chamber was heated to 37°C for 3 h (pretreatment), returned to 23°C for 2 h, heated to 45°C for 3 h (treatment), and then allowed to recover for 2 h. Two biological replicates were performed. Leaves from three different plants were harvested at seven selected time points during heat stress treatment, frozen immediately in liquid nitrogen, and stored at −80°C for further analysis.
RNA isolation and semi-quantitative RT-PCR
Total RNA was extracted using the RNeasy Plant Mini Kit (Qiagen) with on-column treatment with RNase-free DNase I (Qiagen) to remove any contamination of genomic DNA according to the manufacturer’s instructions. RNA integrity was verified by 2% agar gel electrophoresis. First-strand cDNA synthesis was carried out with approximately 1 μg RNA using the SuperScript III reverse transcription kit (Invitrogen) and random primers according to the manufacturer’s procedure. Primers with melting temperatures of 58–60°C, lengths of 20–27 bp, and amplicon lengths of 160–260 bp were designed using Primer3 software (http://frodo.wi.mit.edu/primer3/input.htm). All primer sequences are listed in Additional file 13.
Real-time PCR was conducted on 7500 Real Time PCR System (Applied Biosystems, CA, USA) using SYBR Premix Ex Taq™ Kit (TaKaRa, Tokyo, Japan). Reactions were prepared in a total volume of 20 μl containing: 10 μl of 2×SYBR Premix, 2 μl of cDNA template, 0.4 μl of each specific primer to a final concentration of 200 nM. The reactions were performed as the following conditions: initial denaturation step of 95°C for 30 s followed by two-step thermal cycling profile of denaturation at 95°C for 10 s, and combined primer annealing/extension at 60°C for 34 s for 40 cycles. Negative PCR control without templates was performed for each primer pair. To verify the specificity of each primer pair, a melting curve analysis was performed ranging from 60°C to 95°C with temperature increasing steps of 0.06°C/s (5 acquisitions per °C) at the end of each run. The final threshold cycle (Ct) values were the mean of eight values including two biological replicates for each treatment and four technical replicates. The PtActin gene was used as an internal control.
This work was supported by the National Basic Research Program of China [2012CB114500] and the National Natural Science Foundation of China  to ML, and the National High Technology Research and Development Program of China [2013AA102702] to JC.
- Wang W, Vinocur B, Shoseyov O, Altman A: Role of plant heat-shock proteins and molecular chaperones in the abiotic stress response. Trends Plant Sci. 2004, 9 (5): 244-252. 10.1016/j.tplants.2004.03.006.View ArticlePubMedGoogle Scholar
- Gupta SC, Sharma A, Mishra M, Mishra RK, Chowdhuri DK: Heat shock proteins in toxicology: how close and how far?. Life Sci. 2010, 86 (11–12): 377-384.View ArticlePubMedGoogle Scholar
- Hu W, Hu G, Han B: Genome-wide survey and expression profiling of heat shock proteins and heat shock factors revealed overlapped and stress specific response under abiotic stresses in rice. Plant Sci. 2009, 176 (4): 583-590. 10.1016/j.plantsci.2009.01.016.View ArticlePubMedGoogle Scholar
- Krishna P, Gloor G: The Hsp90 family of proteins in Arabidopsis thaliana. Cell Stress Chaperone. 2001, 6 (3): 238-246. 10.1379/1466-1268(2001)006<0238:THFOPI>2.0.CO;2.View ArticleGoogle Scholar
- Young JC, Moarefi I, Hartl FU: Hsp90: a specialized but essential protein-folding tool. J Cell Biol. 2001, 154 (2): 267-273. 10.1083/jcb.200104079.PubMed CentralView ArticlePubMedGoogle Scholar
- Richter K, Buchner J: Hsp90: chaperoning signal transduction. J Cell Physiol. 2001, 188 (3): 281-290. 10.1002/jcp.1131.View ArticlePubMedGoogle Scholar
- Milioni D, Hatzopoulos P: Genomic organization of hsp90 gene family in Arabidopsis. Plant Mol Biol. 1997, 35 (6): 955-961. 10.1023/A:1005874521528.View ArticlePubMedGoogle Scholar
- Liu Y, Burch-Smith T, Schiff M, Feng S, Dinesh-Kumar SP: Molecular chaperone Hsp90 associates with resistance protein N and its signaling proteins SGT1 and Rar1 to modulate an innate immune response in plants. J Biol Chem. 2004, 279 (3): 2101-2108. 10.1074/jbc.M310029200.View ArticlePubMedGoogle Scholar
- Hubert DA, Tornero P, Belkhadir Y, Krishna P, Takahashi A, Shirasu K, Dangl JL: Cytosolic HSP90 associates with and modulates the Arabidopsis RPM1 disease resistance protein. EMBO J. 2003, 22 (21): 5679-5689. 10.1093/emboj/cdg547.PubMed CentralView ArticlePubMedGoogle Scholar
- Imai J, Maruya M, Yashiroda H, Yahara I, Tanaka K: The molecular chaperone Hsp90 plays a role in the assembly and maintenance of the 26S proteasome. EMBO J. 2003, 22 (14): 3557-3567. 10.1093/emboj/cdg349.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang Z, Quick MK, Kanelakis KC, Gijzen M, Krishna P: Characterization of a plant homolog of hop, a cochaperone of hsp90. Plant Physiol. 2003, 131 (2): 525-535. 10.1104/pp.011940.PubMed CentralView ArticlePubMedGoogle Scholar
- Rutherford SL, Lindquist S: Hsp90 as a capacitor for morphological evolution. Nature. 1998, 396 (6709): 336-342. 10.1038/24550.View ArticlePubMedGoogle Scholar
- Cao D, Froehlich JE, Zhang H, Cheng CL: The chlorate-resistant and photomorphogenesis-defective mutant cr88 encodes a chloroplast-targeted HSP90. Plant J. 2003, 33 (1): 107-118. 10.1046/j.1365-313X.2003.016011.x.View ArticlePubMedGoogle Scholar
- Ishiguro S, Watanabe Y, Ito N, Nonaka H, Takeda N, Sakai T, Kanaya H, Okada K: SHEPHERD is the Arabidopsis GRP94 responsible for the formation of functional CLAVATA proteins. EMBO J. 2002, 21 (5): 898-908. 10.1093/emboj/21.5.898.PubMed CentralView ArticlePubMedGoogle Scholar
- Prassinos C, Haralampidis K, Milioni D, Samakovli D, Krambis K, Hatzopoulos P: Complexity of Hsp90 in organelle targeting. Plant Mol Biol. 2008, 67 (4): 323-334. 10.1007/s11103-008-9322-8.View ArticlePubMedGoogle Scholar
- Song H, Zhao R, Fan P, Wang X, Chen X, Li Y: Overexpression of AtHsp90.2, AtHsp90.5 and AtHsp90.7 in Arabidopsis thaliana enhances plant sensitivity to salt and drought stresses. Planta. 2009, 229 (4): 955-964. 10.1007/s00425-008-0886-y.View ArticlePubMedGoogle Scholar
- Tuskan GA, Difazio S, Jansson S, Bohlmann J, Grigoriev I, Hellsten U, Putnam N, Ralph S, Rombauts S, Salamov A: The genome of black cottonwood, Populus trichocarpa (Torr. & Gray). Science. 2006, 313 (5793): 1596-1604. 10.1126/science.1128691.View ArticlePubMedGoogle Scholar
- Marchler-Bauer A, Lu S, Anderson JB, Chitsaz F, Derbyshire MK, DeWeese-Scott C, Fong JH, Geer LY, Geer RC, Gonzales NR: CDD: a conserved domain database for the functional annotation of proteins. Nucleic Acids Res. 2011, 39: D225-229. 10.1093/nar/gkq1189.PubMed CentralView ArticlePubMedGoogle Scholar
- Marchler-Bauer A, Anderson JB, Chitsaz F, Derbyshire MK, DeWeese-Scott C, Fong JH, Geer LY, Geer RC, Gonzales NR, Gwadz M: CDD: specific functional annotation with the conserved domain database. Nucleic Acids Res. 2009, 37: D205-210. 10.1093/nar/gkn845.PubMed CentralView ArticlePubMedGoogle Scholar
- Dutta R, Inouye M: GHKL, An emergent ATPase/kinase superfamily. Trends Biochem Sci. 2000, 25 (1): 24-28. 10.1016/S0968-0004(99)01503-0.View ArticlePubMedGoogle Scholar
- Hardison RC: A brief history of hemoglobins: plant, animal, protist, and bacteria. Proc Natl Acad Sci U S A. 1996, 93 (12): 5675-5679. 10.1073/pnas.93.12.5675.PubMed CentralView ArticlePubMedGoogle Scholar
- Bailey TL, Williams N, Misleh C, Li WW: MEME: discovering and analyzing DNA and protein sequence motifs. Nucleic Acids Res. 2006, 34: W369-373. 10.1093/nar/gkl198.PubMed CentralView ArticlePubMedGoogle Scholar
- Obermann WM, Sondermann H, Russo AA, Pavletich NP, Hartl FU: In vivo function of Hsp90 is dependent on ATP binding and ATP hydrolysis. J Cell Biol. 1998, 143 (4): 901-910. 10.1083/jcb.143.4.901.PubMed CentralView ArticlePubMedGoogle Scholar
- Panaretou B, Prodromou C, Roe SM, O'Brien R, Ladbury JE, Piper PW, Pearl LH: ATP binding and hydrolysis are essential to the function of the Hsp90 molecular chaperone in vivo. EMBO J. 1998, 17 (16): 4829-4836. 10.1093/emboj/17.16.4829.PubMed CentralView ArticlePubMedGoogle Scholar
- Kent WJ, Baertsch R, Hinrichs A, Miller W, Haussler D: Evolution's cauldron: duplication, deletion, and rearrangement in the mouse and human genomes. Proc Natl Acad Sci U S A. 2003, 100 (20): 11484-11489. 10.1073/pnas.1932072100.PubMed CentralView ArticlePubMedGoogle Scholar
- Mehan MR, Freimer NB, Ophoff RA: A genome-wide survey of segmental duplications that mediate common human genetic variation of chromosomal architecture. Hum Genomics. 2004, 1 (5): 335-344. 10.1186/1479-7364-1-5-335.PubMed CentralView ArticlePubMedGoogle Scholar
- Blanc G, Hokamp K, Wolfe KH: A recent polyploidy superimposed on older large-scale duplications in the Arabidopsis genome. Genome Res. 2003, 13 (2): 137-144. 10.1101/gr.751803.PubMed CentralView ArticlePubMedGoogle Scholar
- Guyot R, Keller B: Ancestral genome duplication in rice. Genome. 2004, 47 (3): 610-614. 10.1139/g04-016.View ArticlePubMedGoogle Scholar
- Yoshida T, Ohama N, Nakajima J, Kidokoro S, Mizoi J, Nakashima K, Maruyama K, Kim J-M, Seki M, Todaka D: Arabidopsis HsfA1 transcription factors function as the main positive regulators in heat shock-responsive gene expression. Mol Genet Genomics. 2011, 286 (5–6): 321-332.View ArticlePubMedGoogle Scholar
- Haseloff J, Siemering KR, Prasher DC, Hodge S: Removal of a cryptic intron and subcellular localization of green fluorescent protein are required to mark transgenic Arabidopsis plants brightly. Proc Natl Acad Sci USA. 1997, 94 (6): 2122-2127. 10.1073/pnas.94.6.2122.PubMed CentralView ArticlePubMedGoogle Scholar
- Lin Y, Cheng CL: A chlorate-resistant mutant defective in the regulation of nitrate reductase gene expression in Arabidopsis defines a new HY locus. Plant Cell. 1997, 9 (1): 21-35.PubMed CentralView ArticlePubMedGoogle Scholar
- Cao D, Lin Y, Cheng CL: Genetic interactions between the chlorate-resistant mutant cr88 and the photomorphogenic mutants cop1 and hy5. Plant Cell. 2000, 12 (2): 199-210.PubMed CentralPubMedGoogle Scholar
- Felts SJ, Owen BA, Nguyen P, Trepel J, Donner DB, Toft DO: The hsp90-related protein TRAP1 is a mitochondrial protein with distinct functional properties. J Biol Chem. 2000, 275 (5): 3305-3312. 10.1074/jbc.275.5.3305.View ArticlePubMedGoogle Scholar
- Gesualdi NM, Chirico G, Pirozzi G, Costantino E, Landriscina M, Esposito F: Tumor necrosis factor-associated protein 1 (TRAP-1) protects cells from oxidative stress and apoptosis. Stress. 2007, 10 (4): 342-350. 10.1080/10253890701314863.View ArticleGoogle Scholar
- Takabatake R, Ando Y, Seo S, Katou S, Tsuda S, Ohashi Y, Mitsuhara I: MAP kinases function downstream of HSP90 and upstream of mitochondria in TMV resistance gene N-mediated hypersensitive cell death. Plant Cell Physiol. 2007, 48 (3): 498-510. 10.1093/pcp/pcm021.View ArticlePubMedGoogle Scholar
- Fletcher JC: Shoot and floral meristem maintenance in Arabidopsis. Annu Rev Plant Biol. 2002, 53 (1): 45-66. 10.1146/annurev.arplant.53.092701.143332.View ArticlePubMedGoogle Scholar
- Sangster TA, Queitsch C: The HSP90 chaperone complex, an emerging force in plant development and phenotypic plasticity. Curr Opin Plant Biol. 2005, 8 (1): 86-92. 10.1016/j.pbi.2004.11.012.View ArticlePubMedGoogle Scholar
- Barrett T, Edgar R: Gene expression omnibus: microarray data storage, submission, retrieval, and analysis. Methods Enzymol. 2006, 411: 352-369.PubMed CentralView ArticlePubMedGoogle Scholar
- Weston DJ, Karve AA, Gunter LE, Jawdy SS, Yang X, Allen SM, Wullschleger SD: Comparative physiology and transcriptional networks underlying the heat shock response in Populus trichocarpa, Arabidopsis thaliana and Glycine max. Plant Cell Environ. 2011, 34 (9): 1488-1506. 10.1111/j.1365-3040.2011.02347.x.View ArticlePubMedGoogle Scholar
- Lescot M, Déhais P, Thijs G, Marchal K, Moreau Y, Van de Peer Y, Rouzé P, Rombauts S: PlantCARE, a database of plant cis-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences. Nucleic Acids Res. 2002, 30 (1): 325-327. 10.1093/nar/30.1.325.PubMed CentralView ArticlePubMedGoogle Scholar
- Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R: Clustal W and Clustal X version 2.0. Bioinformatics. 2007, 23 (21): 2947-2948. 10.1093/bioinformatics/btm404.View ArticlePubMedGoogle Scholar
- Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003, 52 (5): 696-704. 10.1080/10635150390235520.View ArticlePubMedGoogle Scholar
- Jones DT, Taylor WR, Thornton JM: The rapid generation of mutation data matrices from protein sequences. Comput Appl Biosci. 1992, 8 (3): 275-282.PubMedGoogle Scholar
- Guo AY, Zhu QH, Chen X, Luo JC: GSDS: a gene structure display server. Yi chuan. 2007, 29 (8): 1023-1026. 10.1360/yc-007-1023.View ArticlePubMedGoogle Scholar
- Zheng H, Camacho L, Wee E, Batoko H, Legen J, Leaver CJ, Malho R, Hussey PJ, Moore I: A Rab-E GTPase mutant acts downstream of the Rab-D subclass in biosynthetic membrane traffic to the plasma membrane in tobacco leaf epidermis. Plant Cell. 2005, 17 (7): 2020-2036. 10.1105/tpc.105.031112.PubMed CentralView ArticlePubMedGoogle Scholar
- Cohen D, Bogeat-Triboulot MB, Tisserant E, Balzergue S, Martin-Magniette ML, Lelandais G, Ningre N, Renou JP, Tamby JP, Le Thiec D: Comparative transcriptomics of drought responses in Populus: a meta-analysis of genome-wide expression profiling in mature leaves and root apices across two genotypes. BMC Genomics. 2010, 11 (1): 630-10.1186/1471-2164-11-630.PubMed CentralView ArticlePubMedGoogle Scholar
- Yuan Y, Chung JD, Fu X, Johnson VE, Ranjan P, Booth SL, Harding SA, Tsai CJ: Alternative splicing and gene duplication differentially shaped the regulation of isochorismate synthase in Populus and Arabidopsis. Proc Natl Acad Sci U S A. 2009, 106 (51): 22020-22025. 10.1073/pnas.0906869106.PubMed CentralView ArticlePubMedGoogle Scholar
- Sturn A, Quackenbush J, Trajanoski Z: Genesis: cluster analysis of microarray data. Bioinformatics. 2002, 18 (1): 207-208. 10.1093/bioinformatics/18.1.207.View ArticlePubMedGoogle 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.