|Home | About | Journals | Submit | Contact Us | Français|
Insight into the mechanisms for the anaerobic metabolism of aromatic compounds by the hyperthermophilic archaeon Ferroglobus placidus is expected to improve understanding of the degradation of aromatics in hot (>80°C) environments and to identify enzymes that might have biotechnological applications. Analysis of the F. placidus genome revealed genes predicted to encode enzymes homologous to those previously identified as having a role in benzoate and phenol metabolism in mesophilic bacteria. Surprisingly, F. placidus lacks genes for an ATP-independent class II benzoyl-CoA (coenzyme A) reductase (BCR) found in all strictly anaerobic bacteria, but has instead genes coding for a bzd-type ATP-consuming class I BCR, similar to those found in facultative bacteria. The lower portion of the benzoate degradation pathway appears to be more similar to that found in the phototroph Rhodopseudomonas palustris, than the pathway reported for all heterotrophic anaerobic benzoate degraders. Many of the genes predicted to be involved in benzoate metabolism were found in one of two gene clusters. Genes for phenol carboxylation proceeding through a phenylphosphate intermediate were identified in a single gene cluster. Analysis of transcript abundance with a whole-genome microarray and quantitative reverse transcriptase polymerase chain reaction demonstrated that most of the genes predicted to be involved in benzoate or phenol metabolism had higher transcript abundance during growth on those substrates vs growth on acetate. These results suggest that the general strategies for benzoate and phenol metabolism are highly conserved between microorganisms living in moderate and hot environments, and that anaerobic metabolism of aromatic compounds might be analyzed in a wide range of environments with similar molecular targets.
Ferroglobus placidus is the only hyperthermophilic microorganism and the only member of the Archaea that is known to be capable of anaerobically oxidizing aromatic compounds. F. placidus can grow by oxidizing benzoate and phenol to carbon dioxide with Fe(III) serving as the sole electron acceptor (Tor and Lovley, 2001).
It has been speculated that the anaerobic metabolism of aromatic compounds may be an important part of the carbon cycle in hot (>80°C) environments (Kashefi et al., 2000; Tor et al., 2001). Aromatic compounds have been detected in sediments associated with hydrothermal vent sites and petroleum reservoirs (Simoneit and Lonsdale, 1982; Bazylinski et al., 1989; Mangani et al., 1991; Goetz and Jannasch, 1993; Magot et al., 2000; Culotta et al., 2007). Previous studies have also shown that benzoate can be oxidized to 14CO2 in hot sediments from the same general location in Vulcano, Italy, where F. placidus was isolated (Hafenbradl et al., 1996; Tor and Lovley, 2001). Although sulfate reduction appeared to be the predominant terminal electron-accepting process in those sediments, additions of molybdate, which inhibited sulfate reduction, did not inhibit benzoate oxidation. These results suggested that benzoate oxidation in the sediments might be coupled to Fe(III) reduction, which was also a minor, but detectable, electron-accepting process in the sediments.
The metabolism of benzoate has been studied intensively in members of mesophilic Bacteria (Harwood et al., 1999; Gibson and Harwood, 2002; Boll, 2005; Carmona et al., 2009). In all anaerobic benzoate-degrading bacteria that have been studied to date, the enzyme benzoate-CoA (coenzyme A) ligase (BCL) first catalyzes the conversion of benzoate to benzoyl-CoA. The next step is the reduction of the aromatic ring to cyclohex-1,5-diene-1-carboxyl-CoA by benzoyl-CoA reductase (BCR). Two different types of BCR have been described (Figure 1). Class I BCR complexes are composed of four subunits, and accept electrons from a reduced ferredoxin protein. This reaction is ATP dependent and appears to be irreversible (Boll and Fuchs, 1995). The best-studied class I BCRs include those from the facultative denitrifying species Thauera aromatica (Boll and Fuchs, 1995), Azoarcus spp. (Lopez Barragan et al., 2004) and the phototroph Rhodopseudomonas palustris (Egland et al., 1997).
A diversity of strict anaerobic microorganisms such as Geobacter species (Wischgoll et al., 2005; Butler et al., 2007), Syntrophus aciditrophicus (Elshahed and McInerney, 2001) and Desulfococcus multivorans (Peters et al., 2004) have a class II BCR. The reaction catalyzed by this complex is reversible and ATP independent (Kung et al., 2009, 2010; Löffler et al., 2010). It has been proposed that strict anaerobes must use class II BCRs because the amount of energy available from benzoate oxidation coupled to the reduction of Fe(III), sulfate or protons is not sufficient to support the substantial energetic requirement of the ATP-dependent class I BCR reaction (Schöcke and Schink, 1999; Kung et al., 2009, 2010; Löffler et al., 2010).
Both classes of BCR complexes produce cyclohex-1,5-diene-1-carbonyl-CoA (Figure 1), which is then metabolized through a series of reactions (addition of a molecule of water to a double bond, dehydrogenation and hydrolytic ring cleavage) that generate a C7-dicarboxyl-CoA compound. Two slightly different versions of these reactions have been described (Figure 2): the BadK–BadH–BadI pathway found in the phototroph R. palustris (Harwood et al., 1999; Gibson and Harwood, 2002) and the Dch–Had–Oah pathway (Breese et al., 1998) found in all other benzoate-degrading organisms studied to date (Carmona et al., 2009). The ultimate fate of this C7-dicarboxyl-CoA compound is to be converted into acetyl-CoA by β-oxidation, either to be fully oxidized to CO2 by the tricyclic acids cycle or to be used in anabolic reactions.
The anaerobic degradation of aromatic compounds other than benzoate requires a series of reactions that convert them into the central metabolite of the pathway, benzoyl-CoA (Harwood et al., 1999; Carmona et al., 2009). In the case of phenol, two different pathways have been proposed (Figure 3). Several denitrifying species from the genera Pseudomonas, Azoarcus and Thauera, as well as the Fe(III) reducer, Geobacter metallireducens, initially phosphorylate phenol to phenylphosphate, which is then carboxylated to produce 4-hydroxybenzoate (Lack and Fuchs, 1994; Breinig et al., 2000; Rabus et al., 2005; Schleinitz et al., 2009). Clostridium hydroxybenzoicum, on the other hand, directly carboxylates phenol to 4-hydroxybenzoate (Zhang and Wiegel, 1994). Once 4-hydroxybenzoate has been formed by one of these two pathways, it is then converted to benzoyl-CoA via the sequential action of 4-hydroxybenzoate-CoA ligase (Merkel et al., 1989) and 4-hydroxybenzoyl-CoA reductase (Gibson et al., 1997).
The recent availability of the complete genome of F. placidus (http://genome.ornl.gov/microbial/fpla/) has made it feasible to evaluate the mechanisms for monoaromatic metabolism in this organism. The results presented here suggest that F. placidus metabolizes benzoate and phenol via pathways that are similar to those commonly found in Bacteria, but surprisingly does not possess the class II BCR complex, which has previously been considered to be common to strict anaerobes.
F. placidus strain AEDII12DO (DSM 10642) was obtained from the type culture collection of the Deutsche Sammlung von Mikroorganismen and Zelkuturen, Braunschweig, Germany. Strict anaerobic culturing and sampling techniques were used throughout (Miller and Wolin, 1974; Balch et al., 1979). F. placidus cells were grown with acetate (10m), benzoate (1m) or phenol (0.5m) as the electron donor, and Fe(III) citrate (56m) as the electron acceptor.
F. placidus medium was prepared as previously described (Tor and Lovley, 2001). After autoclaving, FeCl2 (1.3m), Na2SeO4 (30μgl–1), Na2WO4 (40μgl–1), APM salts (1gl–1 MgCl2, 0.23gl–1 CaCl2) (Coates et al., 1995), DL vitamins (Lovley and Phillips, 1988) and all electron donors were added to the sterilized medium from anaerobic stock solutions. Cultures were incubated under N2:CO2 (80:20) at 85°C in the dark.
Concentrations of benzoate, phenol and phenylphosphate were determined with high performance liquid chromatography (Agilent 1100 HPLC Series, Santa Clara, CA, USA), with an Altima HP C18 HL column. The eluent consisted of MeOH–H2O (60:40) and 0.1% H3PO4 and the compounds were detected at an absorbance of 280nm
Phenylphosphate, an intermediate from anaerobic phenol degradation, was also detected with thin layer chromatography. Phenol-grown cells were harvested by centrifugation (10500g; 10min) and the supernatant was filtered with a 0.2μm filter. Approximately 30μl of the supernatant was then loaded onto Whatman flexible plates (250m layer Al/Sil G/UV), 1μl at a time, and the intermediates were separated with a solvent consisting of ethanol/dichloromethane/H2O (8:1:1) as previously described (Lack and Fuchs, 1994; Schleinitz et al., 2009). The separated compounds were then visualized on an ultraviolet light box.
Operon organization of the F. placidus genome was predicted using the commercial version of the FGENESB software (V Solovyev and A Salamov, unpublished data; Softberry Inc., Mount Kisco, NY, USA; 2003–2009), as previously described for the genomes of various Geobacteraceae (Krushkal et al., 2007; Mahadevan et al., 2008; Tran et al., 2008). Sequence parameters used in Markov chain-based modeling of protein-coding genes were estimated by FGENESB via an iterative procedure using the sequence of each genome and a minimum ORF length of 100bp.
RNA for microarray and quantitative reverse transcriptase polymerase chain reaction (qRT–PCR) analyses was extracted during exponential growth of F. placidus on benzoate, acetate and phenol. F. placidus exhibits different growth rates on these three compounds. Based on previously reported data, the doubling times were 45h for acetate, 63h for benzoate and 71h for phenol (Tor et al., 2001; Tor and Lovley, 2001). In all cases, triplicate RNA samples were extracted during exponential growth.
For extraction of RNA, cultures (100ml in 156ml serum bottles) were divided into 50ml conical tubes (Falcon, BD Biosciences, Bedford, MA, USA), and cells were pelleted by centrifugation at 3000g for 15min. Pellets were then immediately frozen in liquid nitrogen and stored at –80°C.
The pellets were resuspended in 10ml HG extraction buffer preheated to 65°C. The HG extraction buffer consisted of 100m Tris–HCl, pH 8; 100m NaCl; 10m EDTA; 2.5% β-mercaptoethanol, 1% sodium dodecyl sulfate, 2% Plant RNA Isolation Aid (Ambion, Woodward, TX, USA); 5m ascorbic acid; Proteinase K (0.6mgml–1) and lysozyme (5mgml–1). The suspended cells were then dispersed into 10 2ml screw cap tubes and incubated at 65°C for 10min. After incubation, samples were placed on ice, and 2μl Superase-In (Ambion) and 0.025m CaCl2 were added. Samples were then centrifuged at 16100g for 10min, and the supernatant was transferred to new 2ml screw cap tubes. A measure of 50μl Plant RNA Isolation Aid (Ambion), 4μl linear acrylamide (5mgml–1; Ambion), 600μl of hot acidic phenol (65°C; pH 4.5) (Ambion) and 400μl of chloroform-isoamyl alcohol (24:1; Sigma, St Louis, MO, USA) were added to the supernatant. These tubes were then mixed on a Labquake rotator (Barnestead/Thermolyne, Dubuque, IA, USA) for 10min and centrifuged at 16100g for 5min. The aqueous layer was removed and transferred to new 2ml screw cap tubes, and 600μl of hot acidic (65°C; pH 4.5) phenol (Ambion) and 400μl of chloroform-isoamyl alcohol (24:1; Sigma) were added. Tubes were mixed on a rotator for 5min and centrifuged at 16100g for 10min. The aqueous layer was removed again and transferred to a new tube, and 100μl of 5 ammonium acetate (Ambion), 20μl of 5mgml–1 glycogen (Ambion) and 1ml of cold (−20°C) isopropanol (Sigma) were added.
Nucleic acids were precipitated at −30°C for 1h and pelleted by centrifugation at 16100g for 30min. The pellet was then cleaned with cold (−20°C) 70% ethanol, dried and resuspended in sterile diethylpyrocarbonate-treated water (Ambion). The resuspended pellets were combined and cleaned with the RNeasy RNA cleanup kit (Qiagen, Valencia, CA, USA) according to the manufacturer's instructions. The RNA cleanup product was then treated with DNA-free DNase (Ambion) according to the manufacturer's instructions.
High-quality RNA was extracted from these culture samples. All samples had A260/A280 ratios of 1.8–2.0, indicating that they were of high purity (Ausubel et al., 1990). In order to ensure that RNA samples were not contaminated with DNA, PCR amplification with primers targeting the 16S rRNA gene was conducted on RNA samples that had not undergone reverse transcription.
Whole-genome microarray hybridizations were carried out by Roche NimbleGen, Inc (Madison, WI, USA). The TransPlex Whole Transcriptome Amplification Kit (Sigma) was used to amplify RNA before transcriptomic analyses. RNA was obtained from three biological replicates, and triplicate technical replicates were conducted for microarray analyses. All cDNA samples were chemically labeled with Cy3 and hybridized by Nimblegen. The oligonucleotide microarrays used in this study were designed based on preliminary genome sequence data of F. placidus (accession number NC_013849) obtained from the DOE Joint Genome Institute JGI website (http://www.jgi.doe.gov).
A complete record of all oligonucleotide sequences used and raw and statistically treated data files is available in the NCBI Gene Expression Omnibus database (GEO data series numbers: GSE26421 (benzoate vs acetate), GSE26423 (phenol vs acetate) and GSE28549 (phenol vs benzoate).
Results from microarray hybridizations were analyzed with the software Array 4 Star (DNASTAR, Madison, WI, USA). P-values were determined with t-Student analysis. A gene was considered differentially expressed only if the P-value was 0.01.
Primer pairs used for qRT–PCR analysis of gene transcript abundance are provided in Supplementary Table 1. All primers were purchased from Eurofins MWG Operon (Huntsville, AL, USA) and designed according to the manufacturer's specifications (amplicon size 100–200bp). Representative products from each of these primer sets were verified by sequencing clone libraries.
Several different housekeeping genes were analyzed as potential external controls. The housekeeping genes topR, which encodes the DNA topoisomerase reverse gyrase; gyrA, which codes for the α subunit of DNA topoisomerase; and metG, which codes for methionyl-tRNA synthetase were tested. Expression of these three genes by F. placidus was stable under a variety of growth rates and conditions and none of them was differentially expressed in any of the microarray studies. Previous studies have used both metG and gyrA as controls for qRT–PCR analyses (Takle et al., 2007; Theis et al., 2007; Henssge et al., 2009). A DuraScript enhanced avian RT single-strand synthesis kit (Sigma) was used to generate cDNA as previously described (Holmes et al., 2004).
For clone library construction, PCR products were purified with the Gel Extraction Kit (Qiagen), and clone libraries were constructed with a TOPO TA cloning kit, version M (Invitrogen, Carlsbad, CA, USA) according to the manufacturer's instructions. In all, 100 plasmid inserts from each clone library were sequenced with the M13F primer at the University of Massachusetts Sequencing Facility.
Once the appropriate cDNA fragments were generated by RT–PCR, qRT–PCR amplification and detection were performed with the 7500 Real-Time PCR System (Applied Biosystems, Carlsbad, CA, USA). Optimal qRT–PCR conditions were determined using the manufacturer's guidelines. Each PCR mixture consisted of a total volume of 25μl and contained 1.5μl of the appropriate primers (stock concentrations, 15μM), and 12.5μl Power SYBR Green PCR Master Mix (Applied Biosystems). Standard curves covering eight orders of magnitude were constructed with serial dilutions of known amounts of purified cDNA quantified with a NanoDrop ND-1000 spectrophotometer at an absorbance of 260nm.
qRT–PCR analyses were conducted with 36 different genes that were differentially expressed during growth on either phenol or benzoate compared with acetate. Triplicate biological and technical replicates were done for each of these genes. The number of mRNA transcripts expressed by these genes was normalized against the average number of mRNA transcripts from the constitutively expressed housekeeping genes topR, metG and gyrA (Supplementary Tables 2 and 3). When comparisons were made between microarray and qRT–PCR studies, similar results were obtained for all three housekeeping genes (R2=0.71, 0.84 and 0.84) (Supplementary Figure 1).
In order to determine whether the F. placidus genome contained genes with homology to genes known to be involved in anaerobic benzoate and phenol metabolism in other microorganisms, comparisons were made to the characterized aromatic degradation pathway genes in T. aromatica (Breese et al., 1998; Breinig et al., 2000), Aromatoleum aromaticum EbN1 (Rabus et al., 2005), G. metallireducens (Wischgoll et al., 2005; Butler et al., 2007; Schleinitz et al., 2009) and R. palustris (Egland et al., 1997; Harwood et al., 1999). Despite the substantial phylogenetic distance between F. placidus and these mesophilic benzoate- and phenol-degrading bacteria, candidate genes for benzoate and phenol metabolism with high homology to those in the mesophilic bacteria could be identified in the F. placidus genome (Supplementary Tables 4 and 5). Many of the genes associated with benzoate metabolism are localized within two gene clusters, whereas genes associated with phenol metabolism are found in a third cluster (Figure 4). This gene organization resembles those previously described in bacteria (Carmona et al., 2009). Since F. placidus is the only hyperthermophilic archaeon known to degrade aromatic compounds, it was of interest to study these pathways in more detail. Therefore, transcriptomic studies were conducted with microarray and qRT–PCR analyses. These studies revealed that F. placidus utilizes the benzoate-CoA ligation pathway for anaerobic benzoate oxidation and the phenylphosphate/phenylcarboxylase pathway for anaerobic phenol oxidation. In addition, once phenol is converted to benzoyl-CoA, it utilizes the same benzoate-CoA ligation pathway as benzoate. Further evaluation of the pathways was based on additional analysis of gene sequences and gene expression patterns as described below.
In order to gain insight into the physiology of F. placidus cells growing with either benzoate or phenol provided as the sole carbon source, three different whole-genome DNA microarray studies were done; benzoate vs acetate, phenol vs acetate and phenol vs benzoate. In the microarray experiment comparing benzoate- and acetate-grown cells, a total of 450 genes (245 up-regulated, 205 down-regulated) exhibited at least a twofold change in expression (P-value cutoff=0.01). Many of the genes that were significantly up-regulated were involved in aromatics degradation (18 genes), energy metabolism (48 genes) or coded for hypothetical or unknown proteins (82 genes) (Supplementary Table 6). Analysis of the microarray comparing phenol- and acetate-grown cells indicated that 161 genes were differentially expressed (102 up-regulated, 59 down-regulated) under these two conditions. Among the 102 genes that were up-regulated during growth on phenol, 18 of them coded for proteins involved in aromatics degradation, 12 coded for proteins involved in energy metabolism and 42 were annotated as coding for hypothetical or unknown proteins. When phenol-grown cells were compared with benzoate-grown cells, 38 genes (10 up-regulated, 28 down-regulated) exhibited at least a fourfold difference and 60% of the up-regulated genes coded for proteins were involved in phenol degradation.
BCL catalyzes the first step of anaerobic benzoate metabolism in mesophilic bacteria and there are 12 genes in the F. placidus genome annotated as AMP-dependent synthetase/ligase genes that could potentially represent BCLs (Supplementary Table 7). Eight of the 12 are most similar to BCLs (Ferp_0091, Ferp_2180, Ferp_1044, Ferp_1567, Ferp_2239, Ferp_0083, Ferp_1484 and Ferp_0788). The other four (Ferp_1228, Ferp_1946, Ferp_0992 and Ferp_2312) are most similar to characterized phenylacetate-CoA ligases. Three of the putative BCLs had higher transcript abundance in cells grown on benzoate vs acetate: Ferp_1044, Ferp_2180 and Ferp_0788. One of the putative phenylacetate-CoA ligases, Ferp_1228, also had 2.3 times more mRNA transcripts during growth on benzoate (Supplementary Table 7).
Ferp_1044 appears to be the most likely candidate for BCL because it is located within cluster 1 (Figure 4), with other genes thought to be associated with benzoate metabolism. It has high (44%) similarity to the BCL from T. aromatica (Supplementary Table 7), and is the most highly up-regulated during growth on benzoate among the candidate BCL genes. However, in the absence of a genetic system to enable gene deletions, it cannot be definitively ruled out that one or more of the other putative BCLs might also have a role in benzoate metabolism.
Duplicate BCLs have been reported for facultative anaerobic benzoate-degrading microorganisms, Azoarcus evansii and A. aromaticum, one associated with aerobic metabolism and another with anaerobic metabolism of benzoate (Gescher et al., 2002; Rabus et al., 2005; Carmona et al., 2009). However, F. placidus is a strict anaerobe (Hafenbradl et al., 1996), and no genes for aerobic metabolism of benzoate via hydroxylation or the BOX pathway (Zaar et al., 2001; Gescher et al., 2002; Risso et al., 2009) were apparent in the genome. Multiple anaerobic BCLs have also been detected in other organisms. R. palustris has three ligases that could potentially catalyze benzoyl-CoA formation anaerobically (Egland et al., 1995). There are also multiple potential ligases in the strict anaerobe S. aciditrophicus (four potential genes), but only one is associated with a benzoate cluster (McInerney et al., 2007). In G. metallireducens, there are also four potential genes, but only Gmet_2143 (bamY) is associated with a benzoate catabolism gene cluster (Butler et al., 2007).
Surprisingly, there are no genes in F. placidus for a class II BCR, which has been proposed to catalyze the next step in benzoate metabolism in all known strict anaerobic bacteria (Löffler et al., 2010). However, there is a set of genes (Ferp_1184–1187) that encodes proteins with high homology to each of the four subunits of the class I bzd-type BCR complex (BzdNOPQ) (Supplementary Table 4) first described in Azoarcus spp. (Lopez Barragan et al., 2004). All four of these subunits are most similar to those found in A. aromaticum EbN1 (42–55% identity, 63–72% similarity). Evidence that these genes are involved in benzoate metabolism includes their localization in a cluster of other genes associated with benzoate metabolism (Figure 4) and increased expression of each of these genes during growth on benzoate vs growth on acetate (Figure 2; Supplementary Table 4). There was also an increase in the transcript abundance for these genes during growth on phenol vs growth on acetate, consistent with the expectation that once phenol is metabolized to benzoyl-CoA (Figure 3; Supplementary Table 5), further metabolism follows the same pathway as shown for benzoate metabolism (Figure 2).
BCRs use low potential ferredoxin proteins as electron donors for the initial dearomatization step (Boll and Fuchs, 1998; Harwood et al., 1999). Multiple genes coding for ferredoxin proteins were up-regulated during growth on benzoate compared with acetate in F. placidus (Supplementary Table 8), but none of these ferredoxin genes were homologous to the low potential ferredoxins found in T. aromatica (Boll and Fuchs, 1998) or R. palustris (BadB) (Boll and Fuchs, 1998; Harwood et al., 1999). However, a gene coding for a putative ferredoxin that was up-regulated fivefold according to microarray studies during growth on benzoate (Ferp_1180) is also associated with other benzoate degradation genes, including those encoding BCR, found in cluster 2 (Figure 4).
Enzymes involved in the regeneration of the reduced ferredoxin have been identified in both R. palustris and T. aromatica (Boll and Fuchs, 1998; Egland and Harwood, 1999; Dorner and Boll, 2002). In R. palustris, a putative nicotinamide adenine dinucleotide phosphate quinone oxidoreductase (BadC) is involved in ferredoxin regeneration (Egland and Harwood, 1999; Carmona et al., 2009); however, a homolog for this gene is not present in the F. placidus genome. In T. aromatica, 2-oxoglutarate:ferredoxin reductase from the citric acid cycle is responsible for regeneration of the reduced ferredoxin (Dorner and Boll, 2002). The F. placidus genome contains genes that could code for 2-oxoglutarate:ferredoxin oxidoreductase in benzoate cluster 1. Both the α and β subunits (Ferp_1033 and Ferp_1034) of this enzyme are found within this cluster and the α subunit is up-regulated 3.4-fold in benzoate-grown cells (Figure 4; Supplementary Table 9).
After the formation of cyclohex-1,5-diene-1-carboxyl-CoA by BCR, a series of reactions (addition of a molecule of water to a double bond, dehydrogenation and hydrolytic ring cleavage) generate a C7-dicarboxyl-CoA compound (Carmona et al., 2009). Two slightly different versions of these reactions have been described. One in R. palustris (Harwood et al., 1999; Gibson and Harwood, 2002) and an alternative pathway described in T. aromatica (Breese et al., 1998) and Azoarcus sp. (Lopez Barragan et al., 2004), which appears to be conserved in all previously described anaerobic benzoate-degrading heterotrophs (Carmona et al., 2009).
There is little evidence for this later pathway in F. placidus, its genome does not contain homologs for the key enzymes 6-hydroxycyclohex-1-ene-1carbonyl-CoA dehydrogenase (Had) and 6-oxocyclohex-1-ene-1-carbonyl-CoA hydratase (Oah). Instead, the sequence of water addition to a double bond followed by dehydrogenation and hydrolytic ring cleavage is more similar to that found in the autotroph R. palustris (Figure 2). This is somewhat unexpected because the BCR of F. placidus is more closely related to that of microorganisms using the T. aromatica/Azoarcus sp. pathway.
The mechanism for the formation of cyclohex-1-ene-1-carboxyl-CoA, the first unique step in the R. palustris pathway (Figure 2), has yet to be definitely determined, but it has been proposed that this results from additional activity of the BCR, with the low-potential electron donor ferredoxin (BadB) participating in multiple catalytic cycles of reduction and oxidation (Egland et al., 1997). F. placidus has several ferredoxin genes that could potentially act as electron donors that were significantly up-regulated during growth on benzoate (Supplementary Table 8), including Ferp_1180, which as noted above is within one of the clusters of benzoate-degradation genes. In R. palustris, the hydratase BadK converts cyclohex-1-ene-1-carboxyl-CoA to 2-hydroxycyclohexane-1-carboxyl-CoA. There are two putative BadK homologs in F. placidus (badK1, Ferp_1035 and badK2, Ferp_1942), but only Ferp_1035 is significantly up-regulated during growth with benzoate and it is located within one of the clusters of benzoate-degradation genes.
The F. placidus genome also contains a likely candidate for the next enzyme in the R. palustris pathway, 2-hydroxycyclohexane-1-carboyxl-CoA dehydrogenase (BadH). Ferp_1233, a NAD+-dependent secondary alcohol dehydrogenase, is predicted to fulfill this function. Ferp_1233 is not located in either of the two benzoate-specific clusters. Instead, it is associated with a cluster that contains genes potentially involved in the degradation of phenylacetate (Ferp_1228–1233), another aromatic compound that feeds into the benzoate-CoA ligation pathway in bacteria. All of the genes from this cluster were significantly up-regulated during growth on both benzoate and phenol compared with growth on acetate (Supplementary Tables 9 and 10).
Ferp_1040 is predicted to code for the ring-cleaving hydrolase, BadI. This gene is located in a cluster of benzoate-specific genes (Figure 4) and transcripts for this gene were significantly higher during growth on benzoate or phenol vs growth on acetate (Figure 2). Ferp_1040 also has homology (30% identity, 44% similarity) to Oah, an enzyme from the alternative pathway for benzoate metabolism in T. aromatica. However, the amino-acid sequence of Ferp_1040 is more similar to BadI (33% identity, 51% similarity), and when Oah, BadI and Ferp_1040 amino-acid sequences were aligned, a 49-bp deletion was apparent in BadI that was also present in Ferp_1040 (Supplementary Figure 2). These results suggest that Ferp_1040 encodes BadI. In R. palustris, the pimeloyl-CoA produced by BadI is subsequently transformed into glutaryl-CoA by a series of reactions catalyzed by enzymes encoded in the pimFABCD operon (Harrison and Harwood, 2005). Such an operon is not present in F. placidus. The large and small subunits of the pimeloyl-CoA dehydrogenase (PimCD) are best represented by Ferp_1566 and Ferp_1579, 43% and 47% similar to PimC and PimD from R. palustris, respectively (Figure 2; Supplementary Table 4). These genes are not associated with any of the aromatic-metabolism clusters. This coupled with the finding that their mRNA transcripts were not highly up-regulated in benzoate-grown cultures makes this assignment the most tenuous of all the genes predicted to be part of the benzoate degradation pathway.
The next step is hydration of the double bond of 2,3-didehydro-pimeloyl-CoA to form 3-hydroxypimeloyl-CoA. A possible candidate to catalyze this reaction is the product of the gene Ferp_1031, located in cluster 1 (Figure 4). Its expression is up-regulated during growth on benzoate and phenol and it is highly homologous to enoyl-CoA hydratases in other organisms.
Sequence analysis indicated that in F. placidus, phenol is activated by the phenylphosphate/phenylcarboxylase pathway. Genes coding for the A and B subunits of phenylphosphate synthase (ppsAB, Ferp_0085–86) and phenylphosphate carboxylase (ppcB, Ferp_0089); together with a 4-hydroxybenzoate-CoA ligase (hclA, Ferp_0091) were found in the genome located within a distinct cluster (Figure 4).
Transcript levels for Ferp_0085, Ferp_0086 and Ferp_0089, as well as other genes in the phenol-specific gene cluster, were higher during growth on phenol vs acetate and phenol vs benzoate (Figures 3 and and4),4), suggesting that this is the preferred pathway of phenol catabolism in F. placidus. Further evidence that 4-hydroxybenzoate is formed by the phosphorylation and carboxylation of phenol was provided by the fact that the intermediate phenylphosphate was detected in phenol-grown cells via thin layer chromatography and high performance liquid chromatography (Supplementary Figure 3).
The next step in phenol metabolism, synthesis of 4-hydroxybenzoyl-CoA, is most likely catalyzed by the product of Ferp_0091. This protein belongs to the BCL family, and is most similar to characterized benzoate- and 4-hydroxybenzoate-CoA ligases from other benzoate- and phenol-degrading organisms (Supplementary Table 7), particularly BclA from T. aromatica (40% identity, 57% similarity) and HbaA from R. palustris (36% identity, 58% similarity). Expression of Ferp_0091 was up-regulated during growth on phenol (Figure 3) but not during growth on benzoate, suggesting that Ferp_0091 is a phenol-specific gene that actually codes for 4-hydroxybenzoate-CoA ligase, hclA (Figure 3).
The last step in the activation of phenol, reductive dehydroxylation of 4-hydroxybenzoyl-CoA to benzoyl-CoA, is catalyzed by the 4-hydroxybenzoyl-CoA reductase (HbaBCD). This enzyme has been studied in R. palustris (Gibson et al., 1997) and has homologs in anaerobic phenol-degrading organisms such as G. metallireducens, T. aromatica and A. aromaticum EbN1, but is not present in the F. placidus genome. The enzyme performing this function in F. placidus is not known. However, it is highly likely that phenol metabolism proceeds via the formation of benzoyl-CoA as the genes involved in benzoate metabolism were highly up-regulated in cells grown on phenol vs acetate (Figure 3), and were expressed at levels comparable to those in cells grown on benzoate (Supplementary Table 11).
An alternative strategy for anaerobic phenol metabolism in which phenol is directly carboxylated to 4-hydroxybenzoate has previously been reported (Zhang and Wiegel, 1994). Two genes in the F. placidus genome that encode UbiD-like carboxylases could serve this function (Ferp_1543 and Ferp_1630). Ferp_1543 (PpcX1) was up-regulated when phenol-grown cells were compared with acetate-grown cells but not when compared with benzoate-grown cells. Ferp_1630 (PpcX2), on the other hand, was not differentially expressed under either condition (Figure 3; Supplementary Table 5). Furthermore, neither of these genes was associated with the phenol-specific gene cluster that contains the other genes involved in metabolizing the 4-hydroxybenzoate produced from this potential reaction.
The organization of genes in the phenol cluster is similar to that found in bacteria and the presence of two genes coding for putative transposases (Ferp_0074 and Ferp_0101) flanking the phenol degradation cluster (Figure 4) are consistent with recent genetic mobility. In contrast, only one of the two benzoate clusters had a transposase and only at one end of the cluster.
The clusters associated with aromatic metabolism contain additional genes with a possible role in benzoate and phenol degradation. For example, Ferp_1028, found in the first benzoate-specific gene cluster (Figure 4) codes for ketopantoate reductase, an enzyme involved in the biosynthesis of CoA, a critical cofactor in the degradation of aromatic compounds. Transcripts for this gene were significantly higher in benzoate and phenol-grown cells vs acetate-grown cells (Figure 4).
The products of Ferp_1041–1043 have some homology to BCRs (58%, 45% and 46% similar to the α, β and γ subunits of BCR from T. aromatica). These genes are also localized in the first cluster of genes with clear involvement in benzoate metabolism, and the expression of these genes is significantly higher during growth on benzoate or phenol vs growth on acetate (Figure 4). However, previously described BCRs have four subunits (Carmona et al., 2009), which suggests that Ferp_1041–1043 may code for an enzyme complex with an alternative function. Comparison to previously sequenced organisms suggests that these genes might encode a 2-hydroxyglutaryl-CoA dehydratase (Hgd) similar to that found in the fermentative bacterium, Clostridium symbiosum. Ferp_1041 is 66% similar to this enzyme's activator subunit (HgdC), and Ferp_1042 and Ferp_1043 are ca. 45% similar to this enzyme's α and β subunits (HgdAB). This enzyme is usually involved in the fermentation of amino acids and has been shown to have a small but significant homology to BCRs from T. aromatica and R. palustris, as well as to putative proteins from several Archaea (Hans et al., 1999). Given these ambiguities, further evaluation is required before an assignment of a specific role for the proteins coded by Ferp_1041–1043 in benzoate metabolism can be made.
In cluster 2, Ferp_1181–1183 were up-regulated during growth on both benzoate and phenol. These genes code for hypothetical proteins of unknown function, and Ferp_1181 appears to be unique to F. placidus (Supplementary Tables 9 and 10). In cluster 3, Ferp_0092 encodes a putative phenol carboxylase (PpcY) that is associated with anaerobic phenol degradation in other organisms such as A. aromaticum EbN1, T. aromatica, G. metallireducens and G. daltonii (Schleinitz et al., 2009). Its closest relative is Gmet_2105 (58% identity/78% similarity). This gene was first identified in T. aromatica and is referred to as ORF8. Although its expression is up-regulated during growth on phenol vs acetate (fourfold) and phenol vs benzoate (2.80-fold), its role in phenol degradation is unclear (Breinig et al., 2000). The roles of Ferp_0090 and Ferp_0094, also up-regulated during growth on phenol, are unknown as well.
There are at least two other clusters of genes whose expression was up-regulated during growth on aromatic compounds; Ferp_1227–1233 and Ferp_1683–1689. Ferp_1227–1233 includes some genes coding for enzymes associated with aromatic metabolism; Ferp_1228 encodes a putative phenylacetate-CoA ligase, Ferp_1233 encodes BadH and Ferp_1230–1232 code for redox proteins. The genes in the Ferp_1683-1689 cluster code primarily for redox proteins. The role of the majority of genes in this cluster in benzoate and phenol degradation is uncertain. However, the fact that they are up-regulated during growth on both aromatic compounds suggests that they belong to a shared pathway, possibly the lower steps of degradation that lead to the formation of acetyl-CoA.
These results suggest that although F. placidus is a hyperthermophilic archaeon, it metabolizes benzoate and phenol via pathways previously described in mesophilic bacteria. However, several unifying features of benzoate metabolism that have been observed in the mesophilic bacteria that have been studied to date, do not apply to F. placidus.
Most surprising is the finding that F. placidus lacks genes for the ATP-independent class II BCRs that have previously been found in all strict anaerobes. It has been postulated that strict anaerobes must use ATP-independent class II BCRs because the consumption of two ATPs with class I BCRs would preclude a sufficient net gain of energy to support growth under strict anaerobic conditions (Kung et al., 2009, 2010; Löffler et al., 2010). However, the fact that F. placidus is able to metabolize benzoate anoxically with an ATP-consuming class I BCR indicates that this rule does not apply to all anaerobic organisms and might be restricted to mesophilic bacteria. The energy available for anaerobic oxidation of benzoate coupled to Fe(III) oxide reduction at 25°C (−92kJ/electron equivalent) is comparable to that available at 85°C (−99kJ/electron equivalent) (Tor and Lovley, 2001), indicating that this alone is not sufficient to account for differences in the metabolic pathways utilized by F. placidus and a mesophilic bacterium such as G. metallireducens.
In addition to providing further insights into the metabolism of aromatic compounds in F. placidus, this study has identified key genes whose expression can be indicative of ongoing metabolism of aromatic compounds at high temperatures. The physiological characteristics of hyperthermophilic microorganisms are often used to infer aspects of evolution of microbial physiology as well as for understanding reactions that might have been prevalent on early Earth. Thus, further study of the diverse microbes capable of anaerobically oxidizing aromatic compounds at high temperatures and their enzymes seems warranted.
We thank Julia Krushkal and Yanhua Qu for their contribution with operon analysis, Muktak Aklujkar for helpful discussions on this manuscript and Joy Ward for her excellent technical support. This research was supported by the Office of Science (BER), U.S. Department of Energy, Cooperative Agreement No. DE-FC02-02ER63446.
Supplementary Information accompanies the paper on The ISME Journal website (http://www.nature.com/ismej)