Previous Article | Next Article ![]()
Applied and Environmental Microbiology, April 2007, p. 2491-2497, Vol. 73, No. 8
0099-2240/07/$08.00+0 doi:10.1128/AEM.02589-06
Copyright © 2007, American Society for Microbiology. All Rights Reserved.


Department of Biological Sciences, University of Warwick, Coventry CV4 7AL, United Kingdom
Received 7 November 2006/ Accepted 13 February 2007
|
|
|---|
|
|
|---|
Spectroscopic studies with species of Sulfolobus and Metallosphaera have shown a differential expression of respiratory complexes dependent on the growth substrate. The current view of the aerobic respiratory chains of the most-studied species for which genomes are available (Sulfolobus acidocaldarius, Sulfolobus solfataricus, and Sulfolobus tokodaii) is that a SoxABCD complex (a heme aa3-CuB oxidase) and a SoxM (heme bb3-CuB) supercomplex serve as the terminal aerobic oxidases (13, 14, 19). Whereas SoxABCD is a quinol oxidase, SoxM in conjunction with the SoxH subunit, which has a CuA center binding motif, is thought to receive electrons from SoxE, sulfocyanin (13). The additional Rieske iron-sulfur proteins, SoxF and SoxL (22), as well as the CbsAB/SoxL2N complex (9), are assumed to be intermediate cytochrome bc1-like quinol oxidases. A further type of heme, aa3 quinol oxidase, was found in Acidianus ambivalens (20), where it apparently represents the entire respiratory system of the organism. The genomes of the sequenced Sulfolobus species possess open reading frames (ORFs) whose products have high similarities to the A. ambivalens subunit I protein of this oxidase, DoxB.
An absorption peak at 573 nm was found as a novel feature of cytochrome spectra of thermoacidophilic archaea (Sulfolobus metallicus, Metallosphaera sedula, and Acidianus brierleyi) grown on ferrous iron or pyrite (3, 6). The expression of soxN, soxL2, cbsA, soxM, and soxB in M. sedula grown on different substrates indicated cbsA to be the only gene out of these five that was highly expressed in pyrite-grown cells (10), and it was suggested that its gene product could be related to the 573-nm cytochrome absorption peak. However, since no iron-expressed genes for terminal oxidases were identified and expression of cbsA during growth on sulfur and yeast extract was also relatively high in comparison to that during growth on pyrite, it can be concluded that the genetic basis of iron oxidation in such organisms remains largely unknown.
The sequences, transcriptional organization, and energy substrate-dependent expression of some of the genes upregulated during growth of S. metallicus on ferrous iron are described in this paper, and the organism's sulfur oxygenase-reductase gene is identified. These genes were found using micro-representational-difference analysis (mRDA), a subtractive hybridization approach (5), starting with cDNAs from S. metallicus grown with the different substrates, ferrous iron and sulfur.
|
|
|---|
RNA extraction and DNase treatment.
Total RNA was extracted using the TRIzol reagent (Invitrogen, Carlsbad, CA), following the manufacturer's instructions. Contaminating DNA was removed with RNase-free DNase I treatment (New England Biolabs, Ipswich, MA), and the sample was considered to be DNA free when 0.5 to 1 µg of RNA used as a template in a PCR did not yield any products after 35 cycles.
mRDA.
Experiments were performed with tester cDNA from ferrous iron-grown cells and driver cDNA from sulfur-grown cells and with the reverse combination of these cDNAs. Double-stranded cDNAs were synthesized using the Universal RiboClone cDNA synthesis system (Promega, Madison, WI) with the times for first- and second-strand syntheses extended to 2 h and 4 h, respectively. Syntheses of amplicons and two rounds of subtractive hybridizations were carried out as described by Becker et al. (5). The difference products obtained after the second round of subtractive hybridization (DP2s) were shotgun cloned using the TOPO-TA cloning kit (Invitrogen). Clones were screened by colony PCR using vector primers M13f/r, and a number of inserts of the appropriate sizes were selected for sequencing.
DNA extraction and inverse PCR.
A series of treatments with neutral-pH sodium dodecyl sulfate (SDS), proteinase K, high salt, and hexadecyltrimethylammonium bromide (2) was applied to S. metallicus to extract DNA. RNA was removed by treatment with pancreatic RNase to allow an accurate spectrophotometric determination of the DNA concentration. The DNA was digested with a number of restriction enzymes and recircularized at a concentration of 5 ng µl1. Inverse PCR (17) was conducted using these religated DNA samples as templates.
RT-PCR and real-time RT-PCR (quantitative RT-PCR [qRT-PCR]).
First-strand cDNA synthesis for reverse transcription-PCR (RT-PCR) was performed using the ImProm-II reverse transcription system (Promega) using random hexamers and 5 mM MgCl2. For RT-PCR, a cDNA amount corresponding to 50 ng RNA was used in a 25-µl reaction volume with 30 cycles using 2.5 mM MgCl2, 0.2 mM of each deoxynucleoside triphosphate, 0.2 µM of each primer (Table 1), and 0.02 U µl1 of Platinum Taq DNA polymerase (Invitrogen). The cycles comprised 0.5 min at 94°C, annealing for 1 min at approximately 2°C below the lower primer melting temperature, and an extension of 1 min per 1,000-bp product length at 72°C. The program was preceded by a 2-min initial denaturation at 94°C and concluded with 5 min at 72°C. Primer combinations Orf1f and FoxDr2, FoxDf2 and FoxCr2, FoxCf2 and FoxBr, FoxBf and FoxGr2, FoxGf2 and FoxHr2, FoxHf2 and FoxIr, FoxIf and FoxJr, FoxJf and FoxEr2, FoxEf and FoxFr2, FoxEf and FoxAr2, and FoxAf2 and Orf2r were used to evaluate the transcriptional organization of the described ORFs.
|
View this table: [in a new window] |
TABLE 1. Oligonucleotides used in this study
|
2D blue native/SDS-PAGE analysis of S. metallicus membrane extracts.
Protein extracts were prepared from S. metallicus by suspension of the cells in ACA750 buffer (24), cell lysis by sonication, and removal of the cell debris by centrifugation at 15,000 x g for 10 min. Membranes were pelleted by centrifugation for 30 min at 100,000 x g, washed once in ACA750 buffer, and solubilized for 1 h at 4°C in ACA750 buffer containing 1% (wt/vol) n-dodecyl-ß-D-maltoside. Solubilized membrane proteins were separated on a 5 to 15% (wt/vol) acrylamide blue native gel prior to second-dimension (2D) SDS-polyacrylamide gel electrophoresis (PAGE) (12% [wt/vol] acrylamide) following described protocols (24). Heme staining was performed using o-dianisidine (23).
Mass spectroscopy and bioinformatic analyses.
Sequences obtained from mRDA fragments were analyzed by using BLASTN and BLASTX searches (1) against GenBank (http://www.ncbi.nlm.nih.gov/BLAST/), and percent identity and similarity values of amino acid sequences from full ORFs were obtained using the Fasta3 algorithm (18) of the EBI database (http://www.ebi.ac.uk/fasta33/#). Theoretical molecular weights and isoelectric points were obtained using the ExPASy Compute pI/Mw tool (http://us.expasy.org/tools/pi_tool.html), and prediction of transmembrane helices used the TMHMM Server, v. 2.0 (http://www.cbs.dtu.dk/services/TMHMM-2.0/).
Selected Coomassie blue-stained protein spots were excised from 2D gels. Tryptic-digest fragments were analyzed by liquid chromatography electrospray ionization tandem mass spectrometry with a Q-ToF instrument fitted with an in-line capillary LC system (Waters MS Technologies). Peptide mass fingerprints were analyzed using the Mascot v2.1 (MatrixScience) search engine aimed at the NCBInr database and the in-house database (consisting of UniProt sequences and the S. metallicus genome fragment sequenced in this work).
Nucleotide sequence accession numbers.
The sequences of the S. metallicus DSM 6482T sulfur oxygenase-reductase gene, the genomic DNA region containing the fox gene cluster, and the partial aspartate oxidase gene have been deposited in the GenBank database with accession numbers EF040586, DQ813663, and EF040585, respectively.
|
|
|---|
![]() View larger version (58K): [in a new window] |
FIG. 1. Agarose gel electrophoresis of the cDNA amplicon (Amp) and the first (DP1) and second (DP2) difference products from the mRDA of ferrous iron- versus sulfur-grown S. metallicus. Tester cDNA was from (a) ferrous iron-grown cells and (b) sulfur-grown cells. M, 100-bp ladder.
|
With cDNA from ferrous iron-grown cells as the tester sample and from sulfur-grown cells as the driver sample (Fig. 1a), 8 different sequences were found in 25 clones of the DP2 fragments. All deduced amino acid sequences showed highest percent identities to hypothetical proteins of S. tokodaii, and with one exception, these were to five hypothetical proteins, ST2587, ST2595, and ST2598 to ST2600, whose encoding genes are clustered in the S. tokodaii genome.
S. metallicus genes similar to genes ST2584 to ST2603 of S. tokodaii.
A 15-kb region of S. metallicus DNA was sequenced using a combination of inverse PCR and conventional PCR to bridge the gaps between sequences obtained from the mRDA DP2 products that were derived from the cDNA of ferrous iron-grown cells. Most of the identified ORFs were found to be very similar to ORFs in S. tokodaii (Table 2). ORFs ST2586, ST2596, and ST2602 within the ST2584-to-ST2603 ORF cluster in S. tokodaii had no counterparts in the S. metallicus sequence (Fig. 2). There are significant differences in the orders of the similar ORFs in the two organisms, and two of them in S. metallicus appear to be gene fusions, comprising the equivalent of ORFs ST2592 and ST2593 and ORFs ST2598 and ST2599 (Fig. 2). The genes in S. metallicus are here referred to as fox genes, for ferrous iron oxidation, in recognition of a presumed link to growth on ferrous iron. The amino acid sequences deduced from the S. metallicus ORFs indicate all proteins except FoxH would have at least one potential transmembrane region (Table 2). With the exception of FoxJ, with 35% and 37% identity to hypothetical proteins ST1860 of S. tokodaii and SSO2676 of S. solfataricus, respectively, only relatively low percent identities of the fox gene products to other than the S. tokodaii ORFs were detected. Nevertheless, these might have indicated possible functional assignments for FoxA, FoxB, and FoxC, as already noted by Kappler et al. with reference to their counterparts in S. tokodaii (10). FoxC has about 25% identity to cytochrome b558/566 subunits A (CbsA) from S. tokodaii (ST1664), S. acidocaldarius (Saci_1858), and S. solfataricus (SSO2801). FoxA and FoxB identities to cytochrome c oxidase subunits I and II were, for FoxA, 24% to DoxB of Acidianus ambivalens (P94117) and 26% to DoxB of S. solfataricus (SSO0044) and for FoxB, 30% to PoxB of Pyrobaculum oguniense (BAC56131) and 32% to SoxH of S. acidocaldarius (Q59825). Additionally, FoxG has a polyferredoxin domain (COG0348) and FoxH contains an N-terminal Fe-S binding domain (COG1146), as well as a C-terminal CBS domain (COG2905), which is frequently found in proteins involved in signal transduction processes.
|
View this table: [in a new window] |
TABLE 2. Characteristics of deduced S. metallicus proteins and their sequence identities and similarities to corresponding predicted proteins of S. tokodaii
|
![]() View larger version (18K): [in a new window] |
FIG. 2. A comparison of ORFs in genome fragments of S. metallicus (ORFs of >300 bp) or S. tokodaii. S. tokodaii ORFs and corresponding fox genes are connected by diagonal lines, and the S. tokodaii gene numbers are given below the S. metallicus gene names.
|
Expression of fox and sor genes under different growth conditions.
The fox genes that were relatively highly expressed in ferrous iron-grown versus sulfur-grown S. metallicus were also highly expressed in pyrite-grown cells, where ferrous iron was available from the mineral (Fig. 3). foxA, foxB, foxC, and foxE were most highly expressed, and foxA most of all. foxG and foxH most noticeably had higher expression levels with pyrite than when ferrous iron was added to the medium, although their final expression levels were quite low. orf2 did not appear to be differentially regulated by ferrous iron and sulfur, and orf1 had higher expression levels with sulfur as the substrate than with ferrous iron or pyrite, which is consistent with the assumption that these are not involved in the same functional task as the other genes. The highest expression level of the sor gene was observed during growth on sulfur (Fig. 3), with a level comparable to that of the fox genes that were highly expressed during growth on ferrous iron. The expression level was lower with pyrite than with sulfur but still higher than with ferrous iron.
![]() View larger version (31K): [in a new window] |
FIG. 3. Relative expression levels of ORFs of the sequenced S. metallicus genome fragment and of the sulfur oxygenase-reductase gene (sor) assessed by qRT-PCR and normalized against 16S rRNA gene expression in individual cultures. Mean values and standard deviations are from analyses of triplicate cultures grown autotrophically on ferrous iron, pyrite, or sulfur.
|
|
View larger version (8K): [in a new window] |
FIG. 4. Investigation of the transcriptional organization of S. metallicus ORFs using RT-PCR with primers situated in adjacent ORFs. Products produced from cDNAs of a ferrous iron-grown culture (Fe) or a pyrite-grown culture (P) or from genomic DNA (G) are shown.
|
Identification of FoxA in membranes of Sulfolobus metallicus.
2D gel analyses of S. metallicus membrane proteins were performed to investigate differential protein expression with different substrates. Blue native PAGE of membrane proteins from pyrite-grown cells showed the presence of a complex with an apparent molecular weight of 270 kDa which was not seen in the membrane proteins from sulfur-grown cells. This complex separated into at least three subunits with apparent molecular sizes around 60, 50, and 27 kDa in SDS-PAGE (Fig. 5). The spot containing protein(s) of apparent relative molecular masses of about 50 kDa (Fig. 5) was analyzed by tandem mass spectrometry of tryptic digest fragments and yielded four different peptides (8 to 19 residues in length) with sequences identical to fragments of the deduced sequence of FoxA. Other individual peptides were identified in the sample, but no proteins in databases were indicated by more than one of these other peptides. A major protein spot appearing in the same relative position on a comparable gel of membrane proteins from ferrous iron-grown S. tokodaii was also analyzed (not shown). Five different peptides (9 to 21 residues) with sequences identical to those of fragments of the FoxA equivalent (ST2595) in S. tokodaii were found in database searches. One other protein, a hypothetical protein from the S. tokodaii genome, was indicated by two peptide sequences from this sample. The migration of these putative FoxA proteins in SDS-PAGE was not predicted from their indicated size of approximately 67 kDa (Table 2). The protein spot with an apparent relative mass of 60 kDa (circled in Fig. 5) showed heme-staining activity after SDS-PAGE.
![]() View larger version (53K): [in a new window] |
FIG. 5. Coomassie blue-stained 2D blue native/SDS-PAGE gels of membrane preparations of Sulfolobus metallicus cells grown on pyrite (a) and sulfur (b). The putative FoxA is indicated as one of three spots (circled) that represent a possible protein complex present in pyrite-grown cells but not sulfur-grown cells. The positions of molecular mass markers (in kDa) are shown on the borders of the gels.
|
|
|
|---|
The RT-PCR experiments investigating the transcriptional organization of the foxD-foxA locus are not entirely conclusive. foxA appears to form a monocistronic mRNA. Only foxB and foxC form a definite transcriptional unit, which may also include foxD. Additional bicistronic mRNAs may be formed by foxG and foxH, foxI and foxJ, and foxE and foxF. The possible transcriptional links between the ORFs would coincide with these genes having the same order in S. metallicus and S. tokodaii. However, the weak RT-PCR products observed here may also have been formed due to an incomplete termination of transcription. Since the expression level of foxC was approximately double that of foxB, it is possible that foxC also has an additional, independent promoter.
A sulfur oxygenase reductase was expected in S. metallicus based on its similar growth characteristics and phylogenetic relationship to the Acidianus species with which the enzyme has been studied. The high level of transcription of sor in sulfur-grown S. metallicus in contrast to the low level in iron-grown cells indicated a growth substrate-dependent transcriptional regulation, an observation also made by Kletzin (11) when RNA from aerobically and anaerobically grown cells of Acidianus ambivalens was compared.
The subtractive hybridization technique used in this study does not lead to an exhaustive sampling of all differentially regulated ORFs. The probability of targeting a differentially regulated gene depends on the formation of restriction fragments which have a suitable size for PCR amplification. Therefore, large ORFs and a group of adjacent ORFs which are all differentially regulated will have a higher likelihood of detection in mRDA than short and functionally isolated ORFs. For example, the ORFs foxE, foxC, and foxB were not detected directly by subtractive hybridization. Similarly, only the sor gene of the genes involved in the sulfur oxidation pathway was identified here. In addition, any constitutively expressed genes involved in iron or sulfur oxidation would not have been revealed by the subtractive hybridization method.
The 16S rRNA gene is often used for normalization of transcriptional data obtained by real-time PCR, but the data here should be regarded as semiquantitative and considered with caution, because it is not known if the 16S rRNA reference gene expression was at the same level during growth with the three substrates. However, the proteomic data showing differential FoxA expression supports the transcriptional trends described here. The gene for an antioxidant protein of the alkyl hydroperoxide reductase family which was shown to be a major cell protein induced in ferrous iron-oxidizing S. metallicus (7) was also found to be upregulated in iron- and pyrite-grown cells but not in sulfur-grown cells, using the same normalization procedure (data not shown).
A possible link between the CbsA protein of the cbsBA-soxL2N gene cluster in M. sedula and the previously noted cytochrome with an absorption peak at 573 nm, which was induced during ferrous iron oxidation, has been considered (10). However, the observed gene expression patterns suggest a significant involvement of the fox genes in iron oxidation, and one of their products may be a more likely candidate for this cytochrome. foxA and foxB homologues have been found in the type strain of M. sedula and in another iron-oxidizing Metallosphaera species, strain J1 (Bathe and Norris, unpublished data). However, additional genetic and protein biochemical work, including the investigation of the redox properties of the fox gene products, is required for elucidation of the potential roles of the individual genes and their products in ferrous iron-oxidizing, thermoacidophilic archaea.
The sequences and expression patterns of the genes described here can be used in development of functional gene probes for investigation of the abundance, diversity, and activity of iron- and sulfur-oxidizing archaea in natural environments and commercial biomining operations. Such probes should also be useful in the assessment of the relative activities of iron and sulfur oxidation in the degradation of different mineral sulfides where, for example, sulfur oxidation is considered to be more critical in efficient dissolution of the relatively refractory chalcopyrite than in dissolution of pyrite. The intermediate transcription level of sor in pyrite-grown cells together with the overall high levels of the foxD-foxA genes might indicate that the dominant metabolic activity on this substrate is iron oxidation but with a certain extent of simultaneous sulfur oxidation. However, the threshold substrate concentrations necessary for induction of the genes, as well as their time-dependent expression during growth on mineral sulfides, are unknown.
This work was supported by EU FP6 BioMinE project NMP2-CT-2005-500329 and by BHPBilliton PLC.
Published ahead of print on 23 February 2007. ![]()
Present address: Institut für Ingenieurbiologie und Biotechnologie des Abwassers, Universität Karlsruhe (TH), Am Fasanengarten, 76128 Karlsruhe, Germany. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»