- freely available
Int. J. Mol. Sci. 2012, 13(7), 7963-7979; doi:10.3390/ijms13077963
Published: 27 June 2012
Abstract: Heat shock protein 90 (Hsp90) is a molecular chaperone highly conserved across the species from prokaryotes to eukaryotes. Hsp90 is essential for cell viability under all growth conditions and is proposed to act as a hub of the signaling network and protein homeostasis of the eukaryotic cells. By interacting with various client proteins, Hsp90 is involved in diverse physiological processes such as signal transduction, cell mobility, heat shock response and osmotic stress response. In this research, we cloned the dshsp90 gene encoding a polypeptide composed of 696 amino acids from the halotolerant unicellular green algae Dunaliella salina. Sequence alignment indicated that DsHsp90 belonged to the cytosolic Hsp90A family. Further biophysical and biochemical studies of the recombinant protein revealed that DsHsp90 possessed ATPase activity and existed as a dimer with similar percentages of secondary structures to those well-studied Hsp90As. Analysis of the nucleotide sequence of the cloned genomic DNA fragment indicated that dshsp90 contained 21 exons interrupted by 20 introns, which is much more complicated than the other plant hsp90 genes. The promoter region of dshsp90 contained putative cis-acting stress responsive elements and binding sites of transcriptional factors that respond to heat shock and salt stress. Further experimental research confirmed that dshsp90 was upregulated quickly by heat and salt shock in the D. salina cells. These findings suggested that dshsp90 might serve as a component of the early response system of the D. salina cells against environmental stresses.
The unicellular green algae Dunaliella is the most halotolerant eukaryote known, and can survive in media of a wide range of chemical compositions and salt concentrations ranging from 0.05 M to saturation. The genus Dunaliella, a member of the Polyblepharidaceae family, was first discovered by Michel Felix Dunal in 1838 and named after him in 1905 . During the past over 100 years, several well-known Dunaliella species, including D. salina, have been used as model organisms to investigate the molecular mechanism of saline adaption [1,2]. The morphology of Dunaliella is similar to another model organism Chlamydomonas reinhardtii except that the former one contains no rigid polysaccharide cell wall to facilitate the rapid alternation of cell volume upon osmotic pressure. In addition to the exceptional halotolerance, some Dunaliella species also accumulate high concentrations of carotenoids and glycerol, which makes them widely used as the cell factory for the commercial production of β-carotene and glycerol [3,4]. The accumulation of carotenoids may be important for Dunaliella to fight against excess blue light or UV irradiation , while glycerol is rapidly converted from starch to maintain the homeostasis when Dunaliella cells are under high-salt stress .
The morphological and metabolic changes of Dunaliella cells during osmotic stress have been well-established . In brief, when subject to hyperosmotic or hypoosmotic shock, an alteration occurs in the Dunaliella cell shape, volume and intracellular ion strength within minutes, followed by a balance of the osmolarity via adjustment of intracellular glycerol concentration achieved within 2–3 hours, and finally the cells adapt to the new situation by synthesizing salt-induced proteins. Undoubtedly, the regulation of intracellular glycerol concentration is the most important mechanism of Dunaliella halotolerance [6,7], and the key enzymes involved in glycerol synthesis including glycerol 3-phosphate dehydrogenase1 (GPD1) have been cloned and studied [2,8]. However, the high-osmolarity glycerol (HOG) pathway, which is well-characterized in the hyperosmotic response of Saccharomyces cerevisiae , has not been identified in Dunaliella. Moreover, the mechanism of Dunaliella halotolerance is much more complex since the cells have to adjust their metabolism state to adapt to the new environment. Recently, several proteomic or genomic studies have revealed that many genes are regulated when Dunaliella cells are subjected to salinity stress [10–13]. Among them, the proteins involved in protein turnover are found to be up-regulated . This discovery is not surprising since the conversion of the Dunaliella cells’ life style requires the synthesis of new proteins, the degradation of unnecessary proteins and the stabilization of proteins against stresses by protein quality control systems. Particularly, two heat shock proteins (Hsps) with high homology to Hsp70 and Hsp90 have been identified by proteomic analysis .
Hsp90, a molecular chaperone highly conserved across the species from prokaryotes to eukaryotes, is abundantly expressed in higher eukaryotes and is essential for cell viability under all growth conditions. Hsp90 exists as a dimer and each subunit contains three domains: the N-terminal ATPase domain (ND) and the C-terminal dimerization domain (CD) connected by the M-domain (MD) . Hsp90 is proposed to be a hub of the signaling network and protein homeostasis of eukaryotic cells [15–17]. Compared with the other general Hsps, Hsp90 specifically interact with a subset of proteins, and more than 200 substrate (also called client) proteins have been identified. These client proteins are involved in diverse physiological/pathological processes such as signal transduction, cell mobility, tumorigenesis, steroid signaling, innate immunity, heat shock response and telomere maintenance [15–18]. The client proteins, Hsp90, Hsp70 and the regulatory proteins assemble into a large multichaperone machine to regulate the client protein function and turnover [14,19]. In Saccharomyces cerevisiae, Hsp90 is found to be required for the high osmotic stress response and regulates a novel signaling pathway via the co-chaperone Cdc37p parallel to the HOG pathway [20,21]. Considering that Hsp70- and Hsp90-like proteins were upregulated in Dunaliella under halophilic conditions  and the importance of Hsp90 in the stress response of eukaryotic cells, it is possible that Dunaliella Hsp90 (DsHsp90) is involved in the halotolerance of Dunaliella. In this research, we cloned the dshsp90 gene, monitored its salt-induced expression and characterized the biochemical and biophysical properties of the recombinant DsHsp90. The results herein suggested that DsHsp90 belongs to the Hsp90A family and is involved in the adaption of Dunaliella against extreme environmental conditions.
2. Results and Discussion
2.1. Cloning and Gene Structure of Dshsp90 from D. Salina
The construction of the cDNA library and the generation of the ESTs were as described previously , and a 566 bp fragment was identified to be homologous to Chlamydomonas reinhardtii Hsp90A (CrHsp90A). The full length gene was obtained by 5′ and 3′ RACEs to amplify the flanking regions using the primers listed in Table S1, which resulted in a 2395 bp cDNA including an ATG start codon, a TAA stop codon, a 41 bp 5′ untranslated region (UTR) and a 263 bp 3′ UTR. The cDNA open reading frame encoded a protein composed of 696 amino acids (Accession No. JQ735968) with a predicted pI of 4.96 and an estimated molecular weight of 79619.2 Da. Hsp90 is a family of proteins with the size ranging from 80 to 110 kDa, and can be divided into five families: cytoplasmic Hsp90 (Hsp90A), endoplasmic reticulum Hsp90 (Hsp90B), chloroplast Hsp90 (Hsp90C), mitochondria Hsp90 (TRAP) and bacteria Hsp90 (HtpG) . The phylogenetic tree (Figure 1) constructed according to the sequences of various Hsp90s clearly indicated that DsHsp90 belongs to the cytoplasmic Hsp90, and has a close evolutionary relationship to the Hsp90As from algae.
Sequence alignment indicated that DsHsp90 shared the highest homology with CrHsp90A (83.3%) and Hsp90 from Volvox carteri f. nagariensis (81.8%). The plant Hsp90As are highly conserved and the similarity of DsHsp90 to the other plant Hsp90s was more than 75%. As shown in Figure 2, the three structural/functional domains (ND, MD and CD) as well as the charged linker region were well-conserved in DsHsp90. Consistent with the phylogenetic tree analysis in Figure 1, DsHsp90 belonging to the cytosolic Hsp90 family was also indicated by the absence of any predictable organelle-specific N-terminal transit peptide sequence. Moreover, the five conserved signature sequences of cytosolic Hsp90As were well-conserved in DsHsp90. Motifs II-V were almost identical to those identified previously , while the Phe conserved in the other plant Hsp90As (NKEIFLRELISNASDALDKIR) was replaced by another aromatic residue Trp in motif I of DsHsp90 (NKEIWLRELISNASDALDKIR). Consistent with previous sequence analysis of Hsp90As, the charged linker and the C-terminus were the most variable regions in plant Hsp90s, which provided Hsp90 specificity in binding client proteins. Notably, the highly conserved motif MEEVD at the C-terminus characteristic of the Hsp90A family, which is crucial to the binding of co-chaperones containing tetratricopeptide repeats to cytosolic Hsp90s , was replaced by MEDVD in DsHsp90. Nonetheless, the similar properties of the acidic amino acid residues Asp and Glu implied that the motif MEDVD in DsHsp90 also maintained the conserved function of the C-terminus in Hsp90As.
Analysis of the nucleotide sequence of the cloned genomic DNA fragment (Accession No. JQ735969) indicated that dshsp90 was 8.79 kb (Figure S1A), which might be the longest genomic sequence of Hsp90s (Figure 3). A comparison of the genomic sequence with the cDNA sequence of dshsp90 revealed that the ORF of dshsp90 contained 21 exons interrupted by 20 introns with an average size of 313 bp, which is longer than the mean length of higher plant introns (249 bp)  but smaller than that of C. reinhardtii (373 bp). The smallest intron was 108 bp (intron 19), while the largest one was 707 bp (intron 15). The 21 exons spread over the 8.79 kb region and ranged in the size from 35 bp to 241 bp with an average size of 110 bp (Table 1), which is much smaller than the mean length of 183 bp in higher plants . Six of the splice junctions were type I, nine were type II, and five were type 0. There were 7 splice types found in dshsp90, and only 2 splice junctions followed the typical canonical consensus GT-AG dinucleotide sequence (Table 1). The boundary structures of the 18 other introns were distributed among 8 TG-GG type, 3 TA-GG type, 2 GC-GT type, 2 TC-GG type, 1 GT-GA type, 1 AA-GT type and 1 GC-GG type. The types of splice site of dshsp90 were much more complicated than the other algae hsp90s. For example, the 7 introns of hsp90 from C. reinhardtii contain 3 TG-GG type, 2 GT-AG type, 1 GT-GG type and 1 GT-GC type, while the 9 introns of hsp90 from V. carteri have 3 TG-GG, 2 GT-AG, 2 TA-GG, 1 AG-GT and 1 AG-TA borders. Alternate splicing of exons can occur between splice junctions of the same type. To examine whether alternate splicing occurred in dshsp90, Northern blot analysis was performed and only a single band could be obtained, implying there was a single specie of mRNA (Figure S2). Moreover, Western blot analysis also indicated that only one form of protein was expressed in the D. salina cells under various conditions (data not shown).
Alignment of the exon sequence of the dshsp90 gene with the other plant hsp90 genes showed that the dshsp90 gene has the most complicated gene structure. Although D. salina and C. reinhardtii have a close evolutionary relationship and DsHsp90 and CrHsp90A shared a high sequence homology of 83.3%, the gene structures between these two algae were quite different (Figure 3). The gene of hsp90 from C. reinhardtii has much fewer introns (only 7). The other plant hsp90 genes also contained less than 10 introns. For example, there are 9 introns for hsp90 from Volvox carteri f. nagariensis, 4 for Coprinopsis cinerea okayama and Schizophyllum commune, 3 for Micromonas and Selaginella moellendorffii, 2 for Oryza sativa, Arabidopsis thaliana and Ostreococcus lucimarinus, and only 1 for Physcomitrella patenshave. Although the debate regarding the origin of introns in eukaryotic cells has lasted for a long time between the introns-early and introns-late concepts, a hypothesis combining both concepts has been proposed to explain the mechanism of intron evolution in eukaryotes, which involves intron gain, intron loss and intron slide [25–27]. Moreover, the genomes of single-cell eukaryotes are usually intron-poor when compared to the intron-rich multicellular organisms . However, the genomes of algae seem to be intron-rich although the mean intron length of the multicellular V. carteri is longer than the single-cellular C. reinhardtii [29,30]. Compared with hsp90s from C. reinhardtii and V. carteri, the hsp90 gene from D. salina was unusually complicated, and all of the three well-structured domains of DsHsp90 were interrupted by several introns. It is unclear yet how the intron number and size differ in the genomes of D. salina and C. reinhardtii, and whether the large number of introns also exists in the other genes from D. salina. Nonetheless, the extremely large number of the introns of dshsp90 suggested that it might be a highly expressed gene containing more introns with regulatory roles since the highly expressed genes in plants are usually less compact and contain more and long introns compared to the low expression genes .
To identify the potential regulatory factors involved in dshsp90 gene expression, sequence analysis was performed of the 5′ upstream region of the dshsp90 gene. The regulatory elements were analyzed by the online software plantCARE , and the results are shown in Figure S3. The promoter region of dshsp90 beginning at −781 bp contained several canonical cis-acting elements such as CAAT box and TATA box to ensure the transcription start. Putative cis-elements responding to light, gibberellin, MeJA, abscisic acid (ABA) and cold stimulation were also identified. The existence of response elements to ABA and MYB binding site and the lack of the other stress-responsive elements such as DRE/CRT suggested that the expression of dshsp90 might be induced by stresses via the ABA-dependent pathway . The putative transcription factor (TF) binding site was also analyzed by TFSEARCHER , which indicated that the promoter of dshsp90 had potential heat shock transcription factor (HSF) binding sites. It is worth noting that three Alfin1 putative binding sites were identified with the motif of GTGGNG or GNGGTG (Figure 4). Alfin1 has been shown to be a novel TF that binds to the promoter of salt-inducible genes . Thus the in silico analysis suggested that dshsp90 might be involved in the salt response of the D. salina cells.
2.2. Biophysical and Biochemical Characterization of Recombinant DsHsp90
The analysis of the amino acid sequence indicated that DsHsp90 belongs to the Hsp90A family. To confirm this proposal, the recombinant DsHsp90 was produced in E. coli, and was purified to a homogenous state with only one single peak in the SEC elution profile (Figure 5A). From the SEC elution volume, the apparent molecular weight (MW) of DsHsp90 was calculated to be 191 kDa, which was slightly larger than the actual MW value of the dimeric DsHsp90 (159 kDa). This is consistent with the fact that the shape of the Hsp90A molecules deviated from the sphere shape , which will affect its position in the SEC profile. Nonetheless, it is clear that DsHsp90 mainly existed as a dimer in solution.
The structural features of DsHsp90 were evaluated by CD spectroscopy, which revealed that the recombinant DsHsp90 was well-folded containing both α-helical and β-sheet structures (Figure 5B). The percentages of the secondary structures were calculated by CDpro , and the three algorithms used for prediction gave rather consistent results. The predicted values are quite consistent with the crystal structures of Hsp90As. For example, the structure of yeast Hsp90 is composed of 32.3% α-helix and 310-helix, 16.7% β-sheet and 50.9% other secondary structures . Since Hsp90s are ATPases, we also examined whether recombinant DsHsp90 possessed ATPase activity using the malachite green reagent method . The data shown in Figure 5C clearly indicated that DsHsp90 could exhibit ATPase activity as revealed by the significant increase in the absorbance at 630 nm when compared to the controls that either lacked the substrate ATP or the enzyme DsHsp90. The three-dimensional structure of DsHsp90 was modeled using SWISS-MODEL using the yeast Hsp90 crystal structure  as the template. The modeled structure of DsHsp90 (Figure 5D) shared similar structural features to the Hsp90A family.
2.3. Relative Expression of the Dshsp90 Gene
The in silico analysis above indicated that dshsp90 might be a heat- or salt-inducible gene. To verify this proposal, real-time quantitative PCR was used to analyze the relative expressions of dshsp90 under heat and salt stress. As shown in Figure 6A, the transcriptional level of DsHsp90 increased immediately after heat shock, and reached its maximum after 30 min heat shock treatment. The maximum level of the dshsp90 transcript was about 7-fold of that without heat treatment. This observation confirmed that dshsp90 responded early upon heat shock, which might facilitate to stabilize its client proteins.
The D. salina cells have been shown to grow optimally in ~2 M NaCl , and a salt shock was performed by transferring the living conditions of the D. salina cells from 2 M to 4 M NaCl. Similar to the heat shock results, the amount of the dshsp90 transcripts increased immediately upon the salt stress and reached maximum after 30 min treatment (Figure 6B). The expression of dshsp90 was maintained at a high level for the duration of our experimental time of 24 h. Meanwhile, the expression level of dsgpd was not altered significantly during the first 6 h cultivation in 4 M NaCl, increased about three-fold at 12 h, and dropped to the normal level at 24 h. One possible explanation of the minor increase of the expression level of dsgpd during the first 6 h is that the osmotic-stress induced gene expression usually occurs 12–24 h after the salt shock , while the other possible reason is that the Dunaliella cells have already accumulated considerable amounts of glycerol to fight against the osmotic pressure under 2 M NaCl . Nonetheless, the observations in Figure 6 suggested that dshsp90 was an early response gene to heat or salt stress. The high dshsp90 expression level under stress conditions also implied that DsHsp90 was essential for the Dunaliella cells to survive under extreme conditions.
3. Experimental Section
Tris, sodium dodecylsulfate (SDS), isopropyl-1-thio-β-d-galactopyranoside (IPTG), ATP, ammonium molybdate, malachite green oxalate, hydrochloride acid and tween-20 were Sigma products. All other chemicals were local products of analytical grade or higher. All restriction enzymes used in this study were purchased from TaKaRa (Shiga, Japan). The SMART™ RACE cDNA Amplification Kit was from Clontech (Palo Alto, CA, USA). The PMD18-T vector, LA Tag and RNAiso Plus kit was from TaKaRa (Shiga, Japan).
3.2. Algae and Growth Conditions
The green algae D. salina strain was bought from the Institute of Hydrobiology, the Chinese Academy of Sciences. The D. Salina cells were grown at 25 °C, illuminated at 4500 lux with a long-day condition (18 h light and 6 h dark) in medium comprising 2 M NaCl, 50 mM NaHCO3, 2.5 mM NaNO3, 5 mM MgSO4, 0.2 mM KH2 PO4, 6 μM EDTA, 2 μM FeCl3, 0.2 mM CaCl2, 7 μM MnCl2, 1 μM ZnSO4, 1 μM Co(NO3)2 and 1 μM CuSO4. The pH was adjusted to 7.5 using 2 M HCl. The Dualiella salina cells used in this study were collected at their logarithmic phase.
3.3. Gene Cloning and Sequence Analysis
The genome of D. salina was extracted by cetyltrimethylammonium bromide (CTAB) using the standard protocol , and the RNA was isolated by RNAiso Plus. The full length cDNA was obtained using the SMARTer™ RACE cDNA Amplification Kit. The whole cDNA was divided into two segments with an overlap, and two pairs of primers were designed accordingly to amplify these two segments by step-down PCR method (Table S1). PCR was run using the following conditions: 1cycle of denaturation at 95 °C/4 min followed by 30 three-round cycles of amplification (95 °C/30 sec, 59 °C, 56 °C, 53 °C/30 sec (5/5/20 cycles), 72 °C/3 min), finally extension at 72 °C/10 min. The PCR products were constructed to the pMD19-T vector, and were confirmed by DNA sequencing (Invitrogen, Shanghai, China). The boundaries of the introns and exons were identified by the Vector NTI 11 software. The promoter sequences were cloned by adaptor PCR using the genome walking method [41,42], and the GSP primers were designed according to the 5′-UTR region. The sequence was analyzed using plantCARE [32,43].
3.4. Material Treatments and Total RNA Isolation
50 mL algae cells in the logarithmic phase were collected for further processing. The heat shock treatment was performed by transferring the cells from 25 °C to 37 °C with the temperature controlled by water bath. The samples used for further time-course analysis were gathered after 15 min, 30 min, 60 min, 90 min and 120 min heat shock treatment. The salt shock was carried out by adding 5.844 g NaCl to 50 mL algae cell solutions, and the final concentration of NaCl was 4 M. The collected cells were treated by liquid nitrogen before total RNA isolation. RNA was isolated by RNAiso Plus, and contaminating DNA was degraded by treating each sample with RNase-free DNase. The total RNA was quantified by optical density and the quality was evaluated by gel electrophoresis by the intactness of the 28S and 18S rRNA subunit. One μg total RNA was reverse transcribed using ReverTra Ace (Toyobo, Japan) and random primer according to the manufacturers’ recommendations. The obtained cDNA was diluted 1:20 with H2O prior to the quantitative PCR analysis.
3.5. Quantitative PCR with SYBR Green
The relative gene expressions of DsHsp90 and DsGPD (Accession No.AY845323.1) in D. salina were analyzed by real-time quantitative reverse-transcription PCR. The relative abundance of the ds18S gene (Accession No. EF195157.1) was also determined and used as the internal control. For each transcript, a standard curve was constructed using the purified PCR product generated for each specific primer pair. The 20 μL reaction solutions contained 4 μL diluted cDNA, 10 pmol primer, 4 μL H2O and 10 μL SYBR Green real time PCR master mix (TOYOBO, Japan). The real time qPCR was run on iQ5(Bio-Rad) with 1 cycle of denaturation at 95 °C/10 min, followed by 40 two-segment cycles of amplification (95 °C/30 sec, 55 °C (18S and GPD) or 57 °C (Hsp90)/15 sec, 72 °C/15 sec). The fluorescence was measured during PCR and one three-segment cycle of product melting (95 °C/1 min, 55 °C/30 sec, 95 °C/30 sec). The baseline adjustment method of the Bio-RadiQ5 software was used to determine the Ct in each reaction. A melting curve was constructed for each primer pair to verify the presence of one gene-specific peak and the absence of primer dimer. All samples were repeated three times.
3.6. Protein Expression and Purification
After digestion of the PCR product and of the plasmid pET28a with NdeI and BamHI, the dshsp90 gene was ligated into the expression vector pET28a (Novagen). The six-His Tag sequence of pET28a vector was fused to the N-terminus of the open reading frame for further purification. The recombinant plasmids were then transformed into E. coli BL21(DE3) (Novagen). The recombinant strain was inoculated into 5 mL of Luria–Bertani medium supplemented with 50 μg/mL kanamycin and grown overnight at 200 r/min and 37 °C. The cultures were diluted (1:100) in the same medium and grown at 37 °C to reach an OD value of 0.6–0.8, and then the recombinant protein expression was induced by the addition of 1 mM IPTG. After 4 h induction at 37 °C, the cultures were harvested and sonicated. The His-tagged recombinant proteins in the soluble fractions were first purified by the metal-chelated affinity chromatography using a 1 mL Pharmacia Ni-NTA His-Bind column, followed by a further step of purification by Mono Q 5/50 GL anion-exchange chromatography. The final products were collected after purification by size exclusion chromatography (SEC) using a Superdex G200 column equipped on an ÄKTA purification system. The buffer used for gel filtration contained 40 mM HEPES and 25 mM KCl, pH 7.5. The purified proteins were found to be homogeneous by SDS-polyacrylamide gel electrophoresis (SDS-PAGE) and SEC analysis, and the purity of the final products was above 98%.
3.7. ATPase Activity Assay
The ATPase activity of DsHsp90 was determined by the liberation of inorganic phosphate using a modified procedure  of the malachite green reagent method . In brief, the working reagent to terminate the reaction was freshly prepared by mixing 130 mM Malachite Green, 7.5% (w/v) ammonium molybdate and 11% (v/v) Tween-20 with a volume ratio of 10:2.5:0.2. The reaction solution contained the enzyme in 40 mM HEPES and 100 mM KCl, pH 7.5, with the addition of 5 mM ATP. The reaction was terminated after 30 min reaction by adding 400 μL hydrochloride acid, 500 μL water and 250 μL working reagent in 100 μL reaction solutions. After 10 min equilibrium, the absorbance at 630 nm was recorded.
The samples used for spectroscopic experiments were prepared by dissolving the protein in 50 mM Tris-HCl buffer containing 25 mM KCl, pH 7.5, with a final protein concentration of 0.2 mg/mL. Far-UV circular dichroism (CD) spectra were recorded on a Jasco 715 spectrophotometer with a 1 mm pathlength cell. The resultant spectrum was baseline-subtracted versus the buffer spectrum. The percentages of the secondary structures were determined from the CD spectrum with a wavelength range of 200–250 nm using the CONTINLL, SELCON3 and CDSSTR algorithms within the CDPro analytical software . The results were presented as the average ± standard error calculated from the values of the three methods.
3.9. Homology Modeling
The automated protein structure homology modeling was performed using the online software SWISS-MODEL in the project mode . The crystal structure of Hsp90 from Saccharomyces cerevisiae (PDB ID: 2CG9)  was used as the template structure. The model structure of DsHsp90 was manipulated and rendered in PyMol (Available online: http://www.pymol.org/).
In this research, we cloned the dshsp90 gene encoding a polypeptide composed of 696 amino acids from the halotolerant unicellular green algae D. salina. Sequence alignment indicated that DsHsp90 shares high homology to the cytosolic Hsp90 family with conserved motifs characteristic of Hsp90A. Further biophysical and biochemical studies of the recombinant protein revealed that DsHsp90 possessed ATPase activity and existed as a dimer with similar percentages of secondary structures to those well-studied Hsp90As. Thus the gene was named dshsp90 according to its structure and function relationship to the Hsp90A family. Analysis of the nucleotide sequence of the cloned genomic DNA fragment indicated that dshsp90 contained 21 exons interrupted by 20 introns, which is much more complicated than the other plant hsp90 genes. Such a large number of introns of dshsp90 suggested that it might be a highly expressed gene containing more introns with regulatory roles. The promoter region of dshsp90 contained putative cis-acting stress-responsive elements and binding sites of transcriptional factors that respond to heat shock and salt stress. Further experimental study confirmed that the dshsp90 transcripts were induced early by heat shock and salt shock, suggesting that dshsp90 was involved in the early response of the D. salina cells against environmental stresses. The characterization of dshsp90 provides a starting point for further investigation of the role of dshsp90 in the haloadaption mechanism of D. salina.
The authors thank Yang Zhou at Sichuan University for his help with the bioinformatics analysis. This investigation was supported by grant 30870498 from the National Natural Science Foundation of China, grant 20070610177 from the Ministry of Education, China, grant 08ZQ026-080 from Sichuan Province and funds from the State Key Laboratory of Biomembrane and Membrane Biotechnology.
- Oren, A. A hundred years of Dunaliella research: 1905–2005. Saline Syst 2005, 1. [Google Scholar] [CrossRef]
- Chen, H.; Jiang, J.G. Osmotic responses of Dunaliella to the changes of salinity. J. Cell Physiol 2009, 219, 251–258. [Google Scholar]
- Hosseini Tafreshi, A.; Shariati, M. Dunaliella biotechnology: Methods and applications. J. Appl. Microbiol 2009, 107, 14–35. [Google Scholar]
- Lamers, P.P.; Janssen, M.; de Vos, R.C.; Bino, R.J.; Wijffels, R.H. Exploring and exploiting carotenoid accumulation in Dunaliella salina for cell-factory applications. Trends Biotechnol 2008, 26, 631–638. [Google Scholar]
- Ben-Amotz, A.; Shaish, A.; Avron, M. Mode of action of the massively accumulated b-carotene of Dunaliella bardawil in protecting the alga against damage by excess irradiation. Plant Physiol 1989, 91, 1040–1043. [Google Scholar]
- Ben-Amotz, A.; Avron, M. The role of glycerol in the osmotic regulation of the halophilic alga Dunaliella parva. Plant Physiol 1973, 51, 875–878. [Google Scholar]
- Borowitzka, L.J.; Brown, A.D. The salt relations of marine and halophilic species of the unicellular green alga, Dunaliella. The role of glycerol as a compatible solute. Arch. Mikrobiol 1974, 96, 37–52. [Google Scholar]
- He, Q.; Qiao, D.; Bai, L.; Zhang, Q.; Yang, W.; Li, Q.; Cao, Y. Cloning and characterization of a plastidic glycerol 3-phosphate dehydrogenase cDNA from Dunaliella salina. J. Plant Physiol 2007, 164, 214–220. [Google Scholar]
- O’Rourke, S.M.; Herskowitz, I.; O’Shea, E.K. Yeast go the whole HOG for the hyperosmotic response. Trends Genet 2002, 18, 405–412. [Google Scholar]
- Jia, Y.; Xue, L.; Li, J.; Liu, H. Isolation and proteomic analysis of the halotolerant alga Dunaliella salina flagella using shotgun strategy. Mol. Biol. Rep 2010, 37, 711–716. [Google Scholar]
- Katz, A.; Waridel, P.; Shevchenko, A.; Pick, U. Salt-induced changes in the plasma membrane proteome of the halotolerant alga Dunaliella salina as revealed by blue native gel electrophoresis and nano-LC-MS/MS analysis. Mol. Cell. Proteomics 2007, 6, 1459–1472. [Google Scholar]
- Liska, A.J.; Shevchenko, A.; Pick, U.; Katz, A. Enhanced photosynthesis and redox energy production contribute to salinity tolerance in Dunaliella as revealed by homology-based proteomics. Plant Physiol 2004, 136, 2806–2817. [Google Scholar]
- Zhao, R.; Cao, Y.; Xu, H.; Lv, L.; Qiao, D.; Cao, Y. Analysis of expressed sequence tags from the green alga Dunaliella salina (chlorophyta). J. Phycol 2011, 47, 1454–1460. [Google Scholar]
- Pearl, L.H.; Prodromou, C. Structure and mechanism of the Hsp90 molecular chaperone machinery. Annu. Rev. Biochem 2006, 75, 271–294. [Google Scholar]
- Rutherford, S.; Knapp, J.R.; Csermely, P. Hsp90 and developmental networks. Adv. Exp. Med. Biol 2007, 594, 190–197. [Google Scholar]
- Zhao, R.; Houry, W.A. Molecular interaction network of the Hsp90 chaperone system. Adv. Exp. Med. Biol 2007, 594, 27–36. [Google Scholar]
- Taipale, M.; Jarosz, D.F.; Lindquist, S. HSP90 at the hub of protein homeostasis: Emerging mechanistic insights. Nat. Rev. Mol. Cell Biol 2010, 11, 515–528. [Google Scholar]
- DeZwaan, D.C.; Freeman, B.C. HSP90 manages the ends. Trends Biochem. Sci 2010, 35, 384–391. [Google Scholar]
- Wandinger, S.K.; Richter, K.; Buchner, J. The Hsp90 chaperone machinery. J. Biol. Chem 2008, 283, 18473–18477. [Google Scholar]
- Yang, X.X.; Maurer, K.C.; Molanus, M.; Mager, W.H.; Siderius, M.; van der Vies, S.M. The molecular chaperone Hsp90 is required for high osmotic stress response in Saccharomyces cerevisiae. FEMS Yeast Res 2006, 6, 195–204. [Google Scholar]
- Yang, X.X.; Hawle, P.; Bebelman, J.P.; Meenhuis, A.; Siderius, M.; van der Vies, S.M. Cdc37p is involved in osmoadaptation and controls high osmolarity-induced cross-talk via the MAP kinase Kss1p. FEMS Yeast Res 2007, 7, 796–807. [Google Scholar]
- Li, G.; Liu, M.; Jiang, Y.; Qiao, D.R.; Cao, Y. Construction and functional gene screening of cDNA library of Dunaliella salina. J. Trop. Subtrop. Bot 2004, 12, 74–78. [Google Scholar]
- Gupta, R.S. Phylogenetic analysis of the 90 kD heat shock family of protein sequences and an examination of the relationship among animals, plants, and fungi species. Mol. Biol. Evol 1995, 12, 1063–1073. [Google Scholar]
- Hawkins, J.D. A survey on intron and exon lengths. Nucleic Acids Res 1988, 16, 9893–9908. [Google Scholar]
- Gilbert, W.; Glynias, M. On the ancient nature of introns. Gene 1993, 135, 137–144. [Google Scholar]
- Rogozin, I.B.; Sverdlov, A.V.; Babenko, V.N.; Koonin, E.V. Analysis of evolution of exon-intron structure of eukaryotic genes. Brief. Bioinform 2005, 6, 118–134. [Google Scholar]
- Roy, S.W.; Gilbert, W. Rates of intron loss and gain: Implications for early eukaryotic evolution. Proc. Natl. Acad. Sci. USA 2005, 102, 5773–5778. [Google Scholar]
- Mourier, T.; Jeffares, D.C. Eukaryotic intron loss. Science 2003, 300. [Google Scholar] [CrossRef]
- Merchant, S.S.; Prochnik, S.E.; Vallon, O.; Harris, E.H.; Karpowicz, S.J.; Witman, G.B.; Terry, A.; Salamov, A.; Fritz-Laylin, L.K.; Marechal-Drouard, L.; et al. The Chlamydomonas genome reveals the evolution of key animal and plant functions. Science 2007, 318, 245–250. [Google Scholar]
- Prochnik, S.E.; Umen, J.; Nedelcu, A.M.; Hallmann, A.; Miller, S.M.; Nishii, I.; Ferris, P.; Kuo, A.; Mitros, T.; Fritz-Laylin, L.K.; et al. Genomic analysis of organismal complexity in the multicellular green alga Volvox carteri. Science 2010, 329, 223–226. [Google Scholar]
- Ren, X.Y.; Vorst, O.; Fiers, M.W.; Stiekema, W.J.; Nap, J.P. In plants, highly expressed genes are the least compact. Trends Genet 2006, 22, 528–532. [Google Scholar]
- Lescot, M.; Dehais, P.; Thijs, G.; Marchal, K.; Moreau, Y.; van de Peer, Y.; Rouze, 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, 325–327. [Google Scholar]
- Yamaguchi-Shinozaki, K.; Shinozaki, K. Organization of cis-acting regulatory elements in osmotic- and cold-stress-responsive promoters. Trends Plant Sci 2005, 10, 88–94. [Google Scholar]
- TFSEARCH: Searching Transcription Factor Binding Sites, version 1.3; Computational Biology Research Center, AIST: Japan, 2002.
- Bastola, D.R.; Pethe, V.V.; Winicov, I. Alfin1, a novel zinc-finger protein in alfalfa roots that binds to promoter elements in the salt-inducible MsPRP2 gene. Plant Mol. Biol 1998, 38, 1123–1135. [Google Scholar]
- Ali, M.M.; Roe, S.M.; Vaughan, C.K.; Meyer, P.; Panaretou, B.; Piper, P.W.; Prodromou, C.; Pearl, L.H. Crystal structure of an Hsp90-nucleotide-p23/Sba1 closed chaperone complex. Nature 2006, 440, 1013–1017. [Google Scholar]
- Sreerama, N.; Woody, R.W. Estimation of protein secondary structure from circular dichroism spectra: Comparison of CONTIN, SELCON, and CDSSTR methods with an expanded reference set. Anal. Biochem 2000, 287, 252–260. [Google Scholar]
- Itaya, K.; Ui, M. A new micromethod for the colorimetric determination of inorganic phosphate. Clin. Chim. Acta 1966, 14, 361–366. [Google Scholar]
- Brock, T.D. Salinity and the ecology of Dunaliella from great salt lake. J. Gen. Microbiol 1975, 89, 285–292. [Google Scholar]
- Stewart, C.N., Jr; Via, L.E. A rapid CTAB DNA isolation technique useful for RAPD fingerprinting and other PCR applications. Biotechniques 1993, 14, 748–750. [Google Scholar]
- Siebert, P.D.; Chenchik, A.; Kellogg, D.E.; Lukyanov, K.A.; Lukyanov, S.A. An improved PCR method for walking in uncloned genomic DNA. Nucleic Acids Res 1995, 23, 1087–1088. [Google Scholar]
- Noguchi, A.; Takekawa, N.; Einarsdottir, T.; Koura, M.; Noguchi, Y.; Takano, K.; Yamamoto, Y.; Matsuda, J.; Suzuki, O. Chromosomal mapping and zygosity check of transgenes based on flanking genome sequences determined by genomic walking. Exp. Anim 2004, 53, 103–111. [Google Scholar]
- plantCARE, available online: http://bioinformatics.psb.ugent.be/webtools/plantcare/html/ accessed on 27 June 2012.
- Fisher, D.K.; Higgins, T.J. A sensitive, high-volume, colorimetric assay for protein phosphatases. Pharm. Res 1994, 11, 759–763. [Google Scholar]
- Arnold, K.; Bordoli, L.; Kopp, J.; Schwede, T. The SWISS-MODEL workspace: A web-based environment for protein structure homology modelling. Bioinformatics 2006, 22, 195–201. [Google Scholar]
|Table 1. Nucleotide sequence at the splice junction sequences and sizes of exons and introns of dshsp90.|
|Exon No||Exon Size(bp)||Intron Size(bp)||Boundary Sequence||Junction Type|
© 2012 by the authors; licensee Molecular Diversity Preservation International, Basel, Switzerland. This article is an open-access article distributed under the terms and conditions of the Creative Commons Attribution license (http://creativecommons.org/licenses/by/3.0/).