Reconstruction and Regulation of the Central Catabolic Pathway in the Thermophilic Propionate-Oxidizing Syntroph Pelotomaculum thermopropion
http://www.100md.com
细菌学杂志 2006年第1期
Laboratory of Applied Microbiology, Marine Biotechnology Institute, Kamaishi, Iwate 026-0001,National Institute of Bioscience and Human-Technology, Agency of Industrial Science and Technology, Tsukuba, Ibaraki 305-8566,Department of Environmental Systems Engineering, Nagaoka University of Technology, Nagaoka, Niigata 940-2188, Japan
ABSTRACT
Obligate anaerobic bacteria fermenting volatile fatty acids in syntrophic association with methanogenic archaea share the intermediate bottleneck step in organic-matter decomposition. These organisms (called syntrophs) are biologically significant in terms of their growth at the thermodynamic limit and are considered to be the ideal model to address bioenergetic concepts. We conducted genomic and proteomic analyses of the thermophilic propionate-oxidizing syntroph Pelotomaculum thermopropionicum to obtain the genetic basis for its central catabolic pathway. Draft sequencing and subsequent targeted gap closing identified all genes necessary for reconstructing its propionate-oxidizing pathway (i.e., methylmalonyl coenzyme A pathway). Characteristics of this pathway include the following. (i) The initial two steps are linked to later steps via transferases. (ii) Each of the last three steps can be catalyzed by two different types of enzymes. It was also revealed that many genes for the propionate-oxidizing pathway, except for those for propionate coenzyme A transferase and succinate dehydrogenase, were present in an operon-like cluster and accompanied by multiple promoter sequences and a putative gene for a transcriptional regulator. Proteomic analysis showed that enzymes in this pathway were up-regulated when grown on propionate; of these enzymes, regulation of fumarase was the most stringent. We discuss this tendency of expression regulation based on the genetic organization of the open reading frame cluster. Results suggest that fumarase is the central metabolic switch controlling the metabolic flow and energy conservation in this syntroph.
INTRODUCTION
Methane fermentation is a complex microbiological process and involves metabolic interactions among a variety of microorganisms. These microorganisms are classified into at least four trophic groups, namely, primary fermenting bacteria, secondary fermenting bacteria, hydrogenotrophic methanogenic archaea, and acetoclastic methanogenic archaea (39). Among these microorganisms, secondary fermenting bacteria oxidize volatile fatty acids and alcohols in syntrophic association with methanogenic archaea and are thus called syntrophic bacteria (or syntrophs). Despite plenty of genetic and genomic information being available for primary fermenting bacteria (18, 30) and methanogenic archaea (11, 13, 41), such information is currently scarce for syntrophic bacteria.
A significant feature of syntrophic bacteria is their growth at the thermodynamic limit. For instance, the Gibbs free energy (G°) change in syntrophic propionate oxidation is approximately –25 kJ mol–1, which is less than the energy needed for synthesizing one ATP molecule (44). Syntrophs and methanogens share this energy for their growth. More surprisingly, Jackson and McInerney (22) have shown that substrate oxidation by syntrophs proceeds at values close to the thermodynamic equilibrium (G° 0 kJ mol–1). Therefore, one can deduce that these bacteria should have extremely efficient catabolic systems (22), and their catabolic pathway is of wide biological interest.
Propionate-oxidizing syntrophs have been isolated from mesophilic (4, 16, 27, 47) and thermophilic (21, 34) methanogenic ecosystems. Two propionate-oxidizing pathways have been proposed for mesophilic syntrophic bacteria, i.e., the methylmalonyl-coenzyme A (CoA) pathway (see Fig. 1) and a pathway via a six-carbon intermediate metabolite. The methylmalonyl-CoA (MMC) pathway has first been proposed for methanogenic freshwater sediment (38) and later confirmed for syntrophic propionate-oxidizing bacteria on the basis of the results of 13C-nuclear magnetic resonance analyses of intermediate metabolites (33) and enzyme activity measurements (19). The second pathway has been proposed for Smithella propionica (9), which produces acetate and butyrate via a six-carbon intermediate metabolite. In addition, a recent study has detected activities of several enzymes in the MMC pathway in a thermophilic syntrophic propionate-oxidizing bacterium (34). However, there has been no information regarding the genetics of syntrophic propionate oxidation.
The present study was conducted to gain the genetic basis for the central catabolic pathway in a thermophilic syntrophic propionate-oxidizing bacterium, Pelotomaculum thermopropionicum strain SI (20, 21). Recent advances in genomics coupled to database-associated bioinformatics have allowed scientists to identify microbial metabolic pathways on the basis of genome sequences (18, 26, 40). Accordingly, we performed genomic and proteomic analyses of P. thermopropionicum, which revealed novel insights into how this organism controls the central metabolic flow and energy production.
MATERIALS AND METHODS
Strains and growth conditions. P. thermopropionicum strain SI (DSM 13744) was grown in cultures alone and in cocultures with Methanothermobacter thermautotrophicus H (DSM 1053) using a culture medium described elsewhere (21). The medium was supplemented with 0.1% Bacto yeast extract (Difco) and an appropriate fermentation substrate at 20 mM. Cultivation was conducted at 55°C under an atmosphere of N2 plus CO2 (80/20 [vol/vol]) without shaking.
DNA isolation and genome-size analysis. P. thermopropionicum was grown alone in 1 liter of medium supplemented with fumarate as the sole substrate. Cells were harvested at the late exponential growth phase by centrifugation, and the total DNA was extracted by the method of Marmur (29). To check for the presence of plasmids, the purified DNA was electrophoresed in agarose gels by the standard procedure (37). To determine the genome size, pulsed-field gel electrophoresis was conducted using a contour-clamped homogeneous electric field (CHEF DRIII; Bio-Rad Laboratories) according to the manufacturer's instructions.
Library construction, shotgun sequencing, and assembly. The procedures for draft genome sequencing, namely, shotgun and fosmid library construction, sequencing, and contig assembly, were conducted by Dragon Genomics. Two types of shotgun library were constructed, i.e., a plasmid library containing 2- to 3-kb inserts and a fosmid library containing 30- to 50-kb inserts. For constructing a plasmid library, the extracted DNA was mechanically sheared, blunted, and ligated into HincII-digested pUC118 (Takara). For constructing a fosmid library, the extracted DNA was fractured by the phenol-chloroform treatment, fragments of approximately 33 to 48 kb were recovered by electrophoresis, and they were ligated into pCC1FOS (Epicenter) using a copy control fosmid library production kit (Epicenter). Clones were selected from the plasmid (15,360 clones) and fosmid (960 clones) libraries and subjected to sequencing from both ends by using a DYEnamic ET dye terminator kit (Amersham Bioscience) and a MegaBASE 4000 sequencer (Amersham Bioscience). Raw sequencing data (quality values in the Phred analysis [CodonCode] of more than 15) were assembled to generate "contigs" by using the Paracel genome assembler with the CAP4 algorithm (Paracel). Physically linked contigs were tentatively connected to generate "scaffolds." A scaffold was named CWW_CXX, where CWW and CXX were identification (ID) numbers of the first and last contigs. Gaps were closed where needed by direct sequencing of an PCR-recovered intercontig fragment.
Gene prediction and annotation. The assembled sequences were first analyzed by the GRIMMER 2.10 program (10) trained with the genome sequence of Clostridium acetobutylicum ATCC 824 (30). Open reading frames (ORFs) detected were subjected to the BLAST search (49) against the Swiss-Prot (3) and COG (Clusters of Orthologous Groups of proteins) (42) databases to select the ORFs that satisfy the following criteria: sizes larger than 50 amino acids, E values to best-match sequences smaller than 1e–100, best-match sequences being not hypothetical or putative proteins, and the initial and last amino acids being identical to those of the best-match proteins. The GRIMMER program was next trained using the obtained ORFs and again used for analyzing the assembled sequences. After this procedure was repeated three times in total, selected ORFs were subjected to homology search against the COG database using the BLASTp program (1) and the nonredundant GenBank database using the BLASTx program (1). Results of the homology search were linked to the scaffold information and stored in Annotation Viewer (Dragon Genomics) and FileMaker pro version 3.0 (Claris). We also constructed an P. thermopropionicum ORF database available for the BLASTp search. The naming scheme for an ORF was CWW_CXX_YY_ZZ (ORF ID number), where CWW_CXX was the ID number of a scaffold, while YY and ZZ were the base numbers of the first and last nucleotides of the ORF.
Functions of ORFs were also analyzed by using the InterPro program (2) linked to the PROSITE, PRINTS, Pfam, SMART, TIGRFAMs, PIR SuperFamily, and ProDom databases to search for conserved domains and motifs and to validate prediction of gene function. The secondary structure and membrane topology were analyzed using the nnPredict (25) and TMpred (17) programs, respectively. Transcription regulation sequences were primarily analyzed using the DBTBS database constructed for analyzing the Bacillus subtilis genome (28), and the results were manually checked. Terminator sequences were analyzed using the FindTerm program (SoftBerry). Molecular weights and isoelectric points (pI) of proteins were predicted using the Compute pI/Mw tool (14). Alignment of amino acid sequences was conducted by using the ClustalX program (45) for visually inspecting conservation of motifs and functionally important amino acids. For phylogenetic analyses, a neighbor-joining tree (36) was constructed by using the njplot software in the ClustalX program and plotted by using TREEVIEW (32).
Proteomics. P. thermopropionicum cells were grown in cocultures (500 ml) with M. thermautotrophicus as described above, harvested by centrifugation at 11,000 x g and 4°C for 10 min, and washed two times using 5 mM sodium phosphate buffer (pH 7.0). Each cell pellet was suspended in 25 ml of an Percoll solution (Amersham Bioscience) supplemented with 0.15 M NaCl and subjected to centrifugation at 10,000 x g and 4°C for 30 min. After the centrifugation, two distinct bands were seen; microscopic observation showed that the upper band was mostly comprised of cells of P. thermopropionicum. We carefully collected this fraction using a pipette, and Percoll particles were removed by washing with the phosphate buffer. The harvested cells were disrupted by sonication (Sonifier 250; Branson), treated with an RNase-DNase solution (1 mg ml–1 DNase I, 0.25 mg ml–1 RNase A, 0.5 M Tris-HCl [pH 7.0], and 50 mM MgCl2), and lysed in a lysis solution {5 M urea, 2 M thiourea, 4% 3-[(3-cholamidopropyl)dimethylammonio]-1-propane-sulfonate, 80 mM dithiothreitol, and 4% IPG buffer (Amersham Bioscience)}. The resultant lysate was centrifuged at 30,000 x g and 4°C for 20 min, and the supernatant was collected. The protein concentration in this supernatant solution was determined using a two-dimensional (2D) Quant kit (Amersham Bioscience).
Two-dimensional gel electrophoresis (2D-GE) was performed using an investigator 2D electrophoresis system (Genomic Solution). Immobiline DryStrip pH 4 to 7 (Amersham Bioscience) was used for the first-dimensional isoelectric focusing, to which approximately 50 μg of protein was loaded. The second-dimensional separation was performed in 12.5% polyacrylamide gels (PAGs) containing 0.1% sodium dodecyl sulfate (SDS). SDS-PAG electrophoresis molecular weight standard low range (Bio-Rad) was used as a size marker. The gels were stained using Sypro Ruby (Molecular Probes) and visualized using the FM BIOII Multi-view system (Hitachi). We performed 2D-GE more than three times for each culture condition and confirmed the reproducibility. Protein spots were quantified using an PDQuest system (Bio-Rad).
For the N-terminal sequencing (NTS) identification, protein spots on an SDS-PAG were transferred to a ProBlott membrane (Applied Biosystems) using an electro membrane blotter (Nippon Eido). Spots were excised from the membrane, and an N-terminal sequence was determined using a 494cLC protein sequencer (Applied Biosystems). For the peptide mass fingerprinting (PMF) identification, a gel was stained using a silver stain kit (Wako). Protein spots were excised and subjected to in-gel trypsin digestion by the method of Katayama et al. (24). Time-of-flight mass spectrometry was performed using a Voyager DE-PRO time-flight mass spectrometer (Applied Biosystems) and -cyano-4-hydroxycinnamic acid (Sigma) as a matrix. Peptide mass peaks were calibrated using the Data Explorer DM software (Applied Biosystems). Mascot search (Matrix Science) was performed using the P. thermopropionicum draft sequence database for identifying protein spots.
Nucleotide sequence accession numbers. The draft genome sequence of P. thermopropionicum has been deposited in the DDBJ, EMBL, and NCBI databases under accession numbers BAAC01000001 to BAAC01000195. The accession numbers for nucleotide sequences of Pct3, MmcD, and OdcAB are AB221127, AB221128, and AB221129, respectively.
RESULTS
Draft genome sequencing. Pulsed-field gel electrophoresis determined the genome size of P. thermopropionicum to be approximately 3 Mb (data not shown). The genome was single and circular, and no plasmid was found. We performed shotgun sequencing of its genome after constructing plasmid (approximately 2-kb inserts) and fosmid (approximately 42-kb inserts) libraries. A total of 28,453 reads (approximately 17 Mb in total, x5.5 coverage) generated 105 scaffolds and an assembled sequence of 2,810,264 bp, which was estimated to cover up to 94% of the total genome. Its G+C content was in good agreement with an estimate in a hybridization test (52.8% (21)). In the draft sequence, we found 2,988 of ORFs, of which approximately 70% could be grouped into the COG (42) categories (data not shown). On the basis of the current draft data, we constructed ORF databases available for keyword and BLAST searches (these databases will be provided on request).
Reconstruction of the central catabolic pathway. P. thermopropionicum is capable of anaerobic growth on propionate, propanol, butanol, ethanol, and lactate in cocultures with a hydrogenotrophic methanogen, such as Methanothermobacter thermautotrophicus H, while it can also grow on fumarate and pyruvate in culture alone (21). It is known that this bacterium produces acetate and propionate (at a molar ratio of 3:1) from pyruvate (21). Since the available substrates and fermentation products are either constituents of the MMC pathway or compounds directly linked to this pathway, we assumed that P. thermopropionicum likely utilized the MMC pathway (Fig. 1) for metabolizing propionate.
No genetic information had been available for the MMC pathway. However, since genes for enzymes functionally identical to those in this pathway have been identified in other organisms, such as genes for enzymes in the citrate cycle in aerobic bacteria and those in the propionate fermentation pathway in Propionibacterium acnes (5), we were able to search for putative genes for enzymes in the MMC pathway. We picked ORFs of interest by keyword and BLAST search against the P. thermopropionicum ORF databases, for which we selected search terms and sequences according to information obtained from the KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway database (23). The selected ORFs were further analyzed for the conservation of signature motifs and important amino acids and, in some cases, similarities in the overall structures. ORFs in the draft genome sequence for Pct3, MmcD, and OdcA (Table 1) were found over two contigs, whose sequences were later completed and independently deposited to the databases (information available from authors upon request).
Table 1 summarizes results of the ORF search, showing that ORFs for all enzymes in the MMC pathway (Fig. 1) were identified. Of the enzymes in the MMC pathway that are different from previously known enzymes, methylmalonyl-CoA mutase (MCM) and fumarate hydratase (FHT) (fumarase) seem to be composed of two different subunits (information available from authors upon request). The two-subunit FHT enzyme was expressed in Escherichia coli, and its fumarate-transforming activity was detected at 55°C (our unpublished result). We also found an ORF for methylmalonyl-CoA epimerase (MCE) [which catalyzes isomerization between (R)-methylmalonyl-CoA and (S)-methylmalonyl-CoA], which should be necessary prior to step III (Fig. 1). Pyruvate:ferredoxin oxidoreductase (POR), rather than pyruvate dehydrogenase (PDH), was found in the P. thermopropionicum genome (Table 1), which is known to utilize low-potential ferredoxin as an electron acceptor. It was also found that each of the last three steps in this pathway, oxaloacetate decarboxylation (step VIII), acetyl-CoA formation (step IX), and acetyl-CoA hydrolysis (step X), can be catalyzed by two different enzymes (Table 1). According to the genomic information, we reconstructed the MMC pathway in P. thermopropionicum (Fig. 2).
An interesting feature found in the sequence analysis was that many ORFs for the enzymes in the MMC pathway, except for those for propionate CoA transferase (PCT) and SDH, formed a tight cluster on the genome (mmcBCDEFGHIJKLM in Fig. 3), which we named the mmc cluster. We also found two ORFs (tps and mmcA) upstream of this cluster (Table 1 and Fig. 3). MmcA is homologous (37% identical in total amino acids) to RocR, an NtrC/NifA family transcriptional regulator in B. subtilis (7). Analyses of two long intergenic regions (mmcA-mmcB [274 bp] and mmcC-mmcD [452 bp]) performed using the DBTBS database (28) found that there existed a consensus sequence for a sigma L-dependent promoter (8) in the mmcA-mmcB intergenic region, while consensus sequences for sigma H-dependent promoters (12) and sigma A-dependent promoters (48) were present in the mmcC-mmcD region (data not shown).
Peripheral pathways. We also found ORFs that showed substantial homologies to genes for enzymes catalyzing peripheral pathways of the MMC pathways (Table 1 and Fig. 2), such as lactate dehydrogenase (LDH), 2-oxoglutarate:ferredoxin oxidoreductase (KOR), and phosphoenolpyruvate synthetase (PPS). LDH is considered to be used for growth of this bacterium on lactate, while KOR and PPS are important for linking the MMC pathway to biosyntheses of cellular components.
P. thermopropionicum is capable of growth on several short-chain alcohols, such as ethanol, propanol, and butanol. The end product from butanol is butyrate, while propionate produced from propanol is further transformed to acetate. Our growth test showed that acetate and propionate were produced from ethanol (data not shown). In the genome of P. thermopropionicum, multiple ORFs for alcohol dehydrogenase (ADH) and aldehyde:ferredoxin oxidoreductase (AOR) and one ORF for alcohol/aldehyde dehydrogenase (AAD) were identified (Table 1 and Fig. 2). A combination of ADH and AOR produces acetate from ethanol, while AAD directly produces acetyl-CoA. Substrates for these ADHs and AORs need to be experimentally identified.
We found ORFs for four hydrogenases (three Fe-only hydrogenases and one NiFeSe hydrogenase [Table 1]) in the draft genome; these ORFs may be important for scavenging reducing equivalents produced in the MMC pathway and alcohol metabolism as molecular hydrogen. As electron carriers, six ORFs for ferredoxins were found (Table 1), while there was no ORF for a soluble c-type cytochrome. A complete set of ORFs for bacterial F-type ATPase was also identified (Table 1).
Proteomic analyses. In order to detect proteins encoded by ORFs in the mmc cluster, we carried out 2D-GE proteomic analysis of soluble proteins in cells grown on propionate or butanol (Fig. 4). In this analysis, the amounts of proteins loaded on a gel were optimized for evaluating the levels of expression of major soluble proteins. The butanol culture was chosen as the control, because its catabolism is not linked to the MMC pathway. We used the draft genome data to identify major protein spots by PMF and/or NTS. As summarized in Table 2, we were able to detect subunits of propionyl-CoA:oxaloacetase transcarboxylase (POT) (McmH), MCM (MmcE), succinyl-CoA synthetase (SCS) (MmcD1), malate dehydrogenase (MDH) (MmcK), FHT (MmcB), and succinate dehydrogenase (SDH) (Sdh1C) as major spots in the gel. MCE and POR were undetected, probably because their sizes were out of the range of the gel. However, we assume that they were also expressed, since the genes of these enzymes were present in the same transcriptional unit (Fig. 3). In addition, we could not detect the predicted PCT (Pct3 in Table 1), whereas a subunit of an alternative putative CoA transferase (spot 17) was detected. This transferase is homologous (31% identical in amino acids) to GctA, a subunit of glutaconate CoA-transferase from Acidaminococcus fermentans (CAA57199). Although spot 17 was not specifically expressed in the propionate culture, we also need to consider a possibility that this protein functions as PCT.
We next determined the relative intensity of each protein spot (the intensity of a spot was divided by the sum of intensities of all spots in a gel) and compared them between the propionate and butanol cultures (Table 2). MmcK and MmcB each appeared as two spots (spots 15 and 16 and spots 18 and 19, respectively), whose average ratios in the relative intensity were 2.3 and 4.8, respectively. This analysis showed that the enzymes in the MMC pathway were generally up-regulated in the propionate culture and that the up-regulation of MmcB was the most stringent. These tendencies of changes in the spot intensity were also observed when the propionate culture was compared with the propanol or ethanol culture (data not shown).
DISCUSSION
In the present study, we carried out draft sequencing of the genome of P. thermopropionicum to obtain the genetic basis for its central metabolic pathway. In combination with the results of the proteomic analysis, we conclude that this thermophilic syntroph utilizes the MMC pathway for propionate oxidation. The draft sequencing and subsequent targeted gap closing provided the complete sequences for genes encoding the enzymes necessary for reconstructing the MMC pathway, which will assist in subsequent studies for biochemical characterization of the enzymes and molecular analyses of transcriptional regulation mechanisms. Below, we discuss some aspects of the central catabolic pathway in P. thermopropionicum on the basis of current genomic information.
A similar pathway has been identified as a propionate-fermenting pathway in the genome of Propionibacterium acnes (5), while we found differences between the MMC pathway in P. thermopropionicum and the propionate-fermenting pathway in P. acnes. First, P. acnes utilizes a malic enzyme catalyzing the conversion between malate and pyruvate (malate + NAD+ = pyruvate + CO2 + NADH) for steps VII and VIII. Second, PCT is used for step I in P. thermopropionicum, while the same reaction is catalyzed by propionyl-CoA synthetase in P. acnes.
The first two steps in the MMC pathways in P. thermopropionicum (steps I and II in Fig. 1) are catalyzed by transferases, i.e., PCT and POT. These enzymes require counterpart substrates for their activities (CoA derivatives, e.g., acetyl-CoA, for PCT and oxaloacetate for POT), while these counterparts are also the downstream intermediate metabolites in the MMC pathway. We think that this is one of the reasons why growth of P. thermopropionicum was so slow. As reported previously (20), growth of this bacterium on propionate in coculture with a hydrogenotrophic methanogen occurred after a long lag period of 20 days. It is likely that this lag period is necessary for P. thermopropionicum to accumulate these counterpart substrates, and propionate oxidation is initiated only after the concentrations of these substrates reach certain levels. Another possibility for the long lag period would be that the cells need much time for the expression of the catabolic enzymes under the energy-limiting conditions.
Proteomic analyses (Fig. 4 and Table 2) showed that a subunit of fumarase (MmcB) was up-regulated in the propionate culture compared to the butanol culture. Other proteins identified to be encoded in the mmc cluster were also up-regulated in the propionate culture, although up-regulation of these proteins was weaker than that of MmcB. Northern blot analysis confirmed this tendency of gene expression regulation (our unpublished results). We think that this tendency reflected the genetic organization of the mmc cluster. As shown in Fig. 3, the sequence analysis suggested that transcription of mmcBC is solely regulated by the sigma L-dependent promoter, while that of mmcDEFGHIJKLM is also under the control of the sigma A and H-dependent promoters. Sigma A is known to be the major sigma factor in B. subtilis (48), which corresponds to sigma 70 in gram-negative bacteria (15). Sigma A has been considered to control expression of approximately 250 operons in B. subtilis, many of which encode housekeeping proteins (the DBTBS database). Sigma H is a nonessential sigma factor involved in the expression of vegetative and early stationary-phase genes in B. subtilis (12). We found 19 ORFs for sigma factors in the draft genome of P. thermopropionicum; these ORFs include 1 ORF for sigma A, 2 ORFs for sigma H, and 1ORF for sigma L (data not shown). These findings suggest that mmcDEFGHIJKLM are transcribed at certain levels under a wide range of growth conditions. In contrast, sigma L is a member of the sigma 54 family of bacterial sigma factors, which are characterized by their cooperative action with enhancer-binding proteins whose function requires nucleotide hydrolysis (6). Many sigma 54-controlled genes are specifically transcribed in the presence of enhancer molecules (6). In the case of the mmc cluster, it is likely that MmcA functions as an enhancer-binding protein, although we have not yet identified what molecules (or conditions) actually enhanced the expression of mmcBC. Since MmcA includes the PAS domain (43), we deduce that expression of MmcB is enhanced by a global physiological status, e.g., the membrane potential (43), rather than a specific substrate of the catabolic pathway. Up-regulation of the other MMC proteins may also be governed by the sigma L-dependent promoter, although its regulation may be weakened by the other promoters.
As discussed above, the genome analyses suggested that the particular genetic organization of P. thermopropionicum has evolved to accomplish the hierarchical transcriptional regulation of its central metabolic pathway. We think that the stringent regulation of fumarase is physiologically expedient in terms of multiple roles of fumarate in that system. First, since fumarate is an intermediate metabolite of the propionate-oxidizing pathway, fumarase should be up-regulated together with other MMC enzymes, when the bacterium metabolizes propionate. In addition, fumarate is a growth substrate in the culture of P. thermopropionicum alone, where P. thermopropionicum primarily oxidizes it to acetate, while this organism also utilizes it as an electron acceptor to form succinate (namely, fumarate disproportionation [46]) (Fig. 2). As shown in Fig. 4 and Table 2, the expression of fumarase was conditional, while that of fumarate reductase (i.e., Sdh1C, spot 20 in Fig. 4) did not show much change. It is beneficial for P. thermopropionicum to mainly utilize fumarate as the electron acceptor when the MMC pathway does not efficiently work, because fumarate respiration can produce the membrane proton gradient. Proteomics (Table 2 and Fig. 4) showed that P. thermopropionicum expressed a significant amount of ATPase, suggesting that this organism can utilize reducing equivalents produced by fermentation for fumarate respiration and oxidative phosphorylation. It is worth noting that two putative SDHs are present in this organism (Table 1); their roles will be addressed in future biochemical studies. All together, we suggest that fumarase acts as the central metabolic switch controlling the metabolic flow and energy conservation in this syntrophic bacterium.
In conclusion, the genomic analyses of P. thermopropionicum offer a glimpse into the sophisticated catabolic system of syntrophic bacteria. Owing to its unique genetic organization, the mmc cluster may also be interesting for scientists who investigate the evolutionary aspects of operon formation (31, 35). We are currently conducting genetic and biochemical experiments for more-quantitative evaluations of the contribution of each promoter under different growth conditions.
ACKNOWLEDGMENTS
We thank Yasuo Igarashi (University of Tokyo) for helpful advice and continuous encouragement. Acknowledgment is also made to Youhei Urushida for N-terminal amino acid sequencing and Reiko Hirano for technical assistance.
This work was supported by New Energy and Industrial Technology Development Organization (NEDO).
REFERENCES
Altschul, S., T. Madden, A. Schaffer, J. Zhang, Z. Zhang, W. Miller, and D. J. Lipman. 1997. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 25:3389-3402.
Apweiler, R., T. K. Attwood, A. Bairoch, A. Bateman, E. Birney, M. Biswas, P. Bucher, L. Cerutti, F. Corpet, M. D. Croning, R. Durbin, L. Falquet, W. Fleischmann, J. Gouzy, H. Hermjakob, N. Hulo, I. Jonassen, D. Kahn, A. Kanapin, Y. Karavidopoulou, R. Lopez, B. Marx, N. J. Mulder, T. M. Oinn, M. Pagni, F. Servant, C. J. Sigrist, and E. M. Zdobnov. 2001. The InterPro database, an integrated documentation resource for protein families, domains and functional sites. Nucleic Acids Res. 29:37-40.
Boeckmann, B., A. Bairoch, R. Apweiler, M. Blatter, A. Estreicher, E. Gasteiger, M. Martin, K. Michoud, C. O'Donovan, I. Phan, S. Pilbout, and M. Schneider. 2003. The SWISS-PROT protein knowledge base and its supplement TrEMBL in 2003. Nucleic Acids Res. 31:365-370.
Boone, D. R., and M. P. Bryant. 1980. Propionate-degrading bacterium, Syntrophobacter wolinii sp. nov., gen. nov., from methanogenic ecosystems. Appl. Environ. Microbiol. 40:626-632.
Brüggemann, H., A. Henne, F. Hoster, H. Liesegang, A. Wiezer, A. Strittmatter, S. Hujer, P. Durre, and G. Gottschalk. 2004. The complete genome sequence of Propionibacterium acnes, a commensal of human skin. Science 305:671-673.
Buck, M., M. T. Gallegos, D. J. Studholme, Y. Guo, and J. D. Gralla. 2000. The bacterial enhancer-dependent 54 (N) transcription factor. J. Bacteriol. 182:4129-4136.
Calogero, S., R. Gardan, P. Glaser, J. Schweizer, G. Rapoport, and M. Debarbouille. 1994. RocR, a novel regulatory protein controlling arginine utilization in Bacillus subtilis, belongs to the NtrC/NifA family of transcriptional activators. J. Bacteriol. 176:1234-1241.
Debarbouille, M., I. Martin-Verstraete, F. Kunst, and G. Rapoport. 1991. The Bacillus subtilis sigL gene encodes an equivalent of sigma 54 from Gram-negative bacteria. Proc. Natl. Acad. Sci. USA 88:9092-9096.
De Bok, F. A., A. J. Stams, C. Dijkema, and D. R. Boone. 2001. Pathway of propionate oxidation by a syntrophic culture of Smithella propionica and Methanospirillum hungatei. Appl. Environ. Microbiol. 67:1800-1804.
Delcher, A. L., D. Harmon, S. Kasif, O. White, and S. L. Salzberg. 1999. Improved microbial gene identification with GLIMMER. Nucleic Acids Res. 27:4636-4641.
Deppenmeier, U., A. Johann, T. Hartsch, R. Merkl, R. A. Schmitz, R. Martinez-Arias, A. Henne, A. Wiezer, S. Baumer, C. Jacobi, H. Bruggemann, T. Lienard, A. Christmann, M. Bomeke, S. Steckel, A. Bhattacharyya, A. Lykidis, R. Overbeek, H. P. Klenk, R. P. Gunsalus, H. J. Fritz, and G. Gottschalk. 2002. The genome of Methanosarcina mazei: evidence for lateral gene transfer between Bacteria and Archaea. J. Mol. Microbiol. Biotechnol. 4:453-461.
Dubnau, E., J. Weir, G. Nair, L. Carter III, C. Moran, Jr., and I. Smith. 1988. Bacillus sporulation gene spo0H codes for 30 (H). J. Bacteriol. 170:1054-1062.
Galagan, J. E., C. Nusbaum, A. Roy, M. G. Endrizzi, P. Macdonald, W. FitzHugh, S. Calvo, R. Engels, S. Smirnov, D. Atnoor, A. Brown, N. Allen, J. Naylor, N. Stange-Thomann, K. DeArellano, R. Johnson, L. Linton, P. McEwan, K. McKernan, J. Talamas, A. Tirrell, W. Ye, A. Zimmer, R. D. Barber, I. Cann, D. E. Graham, D. A. Grahame, A. M. Guss, R. Hedderich, C. Ingram-Smith, H. C. Kuettner, J. A. Krzycki, J. A. Leigh, W. Li, J. Liu, B. Mukhopadhyay, J. N. Reeve, K. Smith, T. A. Springer, L. A. Umayam, O. White, R. H. White, E. Conway de Macario, J. G. Ferry, K. F. Jarrell, H. Jing, A. J. Macario, I. Paulsen, M. Pritchett, K. R. Sowers, R. V. Swanson, S. H. Zinder, E. Lander, W. W. Metcalf, and B. Birren. 2002. The genome of M. acetivorans reveals extensive metabolic and physiological diversity. Genome Res. 12:532-542.
Gasteiger, E., C. Hoogland, A. Gattiker, S. Duvaud, M. R. Wilkins, R. D. Appel, and A. Bairoch. 2005. Protein identification and analysis tools on the ExPASy server, p. 571-607. In J. M. Walker (ed.), The proteomics protocols handbook. Humana Press, Totowa, N.J.
Gitt, M. A., L. F. Wang, and R. H. Doi. 1985. A strong sequence homology exists between the major RNA polymerase sigma factors of Bacillus subtilis and Escherichia coli. J. Biol. Chem. 260:7178-7185.
Harmsen, H. J., B. L. van Kuijk, C. M. Plugge, A. D. Akkermans, W. M. de Vos, and A. J. Stams. 1998. Syntrophobacter fumaroxidans sp. nov., a syntrophic propionate-degrading sulfate-reducing bacterium. Int. J. Syst. Bacteriol. 48:1383-1387.
Hofmann, K., and W. Stoffel. 1993. TMPred—a database of membrane spanning protein segments. Biol. Chem. Hoppe-Seyler 347:166.
Hong, S. H., J. S. Kim, S. Y. Lee, Y. H. In, S. S. Choi, J. K. Rih, C. H. Kim, H. Jeong, C. G. Hur, and J. J. Kim. 2004. The genome sequence of the capnophilic rumen bacterium Mannheimia succiniciproducens. Nat. Biotechnol. 22:1275-1281.
Houwen, F. P., J. Plokker, A. J. M. Stams, and A. J. B. Zehnder. 1990. Enzymatic evidence for involvement of the methylmalonyl-CoA pathway in propionate oxidation by Syntrophobacter wolinii. Arch. Microbiol. 155:52-55.
Imachi, H., Y. Sekiguchi, Y. Kamagata, A. Ohashi, and H. Harada. 2000. Cultivation and in situ detection of a thermophilic bacterium capable of oxidizing propionate in syntrophic association with hydrogenotrophic methanogens in a thermophilic methanogenic granular sludge. Appl. Environ. Microbiol. 66:3608-3615.
Imachi, H., Y. Sekiguchi, Y. Kamagata, A. Ohashi, and H. Harada. 2002. Pelotomaculum thermopropionicum gen. nov., sp. nov., an anaerobic, thermophilic, syntrophic propionate-oxidizing bacterium. Int. J. Syst. Evol. Microbiol. 52:1729-1735.
Jackson, B. E., and M. J. McInerney. 2002. Anaerobic microbial metabolism can proceed close to thermodynamic limits. Nature 415:454-456.
Kanehisa, M., and S. Goto. 2000. KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 28:27-30.
Katayama, H., T. Nagasu, and Y. Oda. 2001. Improvement of in-gel digestion protocol for peptide mass fingerprinting by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry. Rapid Commun. Mass Spectrom. 15:1416-1421.
Kneller, D. G., F. E. Cohen, and R. Langridge. 1990. Improvements in protein secondary structure prediction by an enhanced neural network. J. Mol. Biol. 214:171-182.
Larimer, F. W., P. Chain, L. Hauser, J. Lamerdin, S. Malfatti, L. Do, M. L. Land, D. A. Pelletier, J. T. Beatty, A. S. Lang, F. R. Tabita, J. L. Gibson, T. E. Hanson, C. Bobst, J. L. T. Torres, C. Peres, F. H. Harrison, J. Gibson, and C. S. Harwood. 2004. The genome sequence of the metabolically versatile photosynthetic bacterium Rhodopseudomonas palustris. Nat. Biotechnol. 22:55-61.
Liu, Y., D. L. Balkwill, H. C. Aldrich, G. R. Drake, and D. R. Boone. 1999. Characterization of the anaerobic propionate-degrading syntrophs Smithella propionica gen. nov., sp. nov. and Syntrophobacter wolinii. Int. J. Syst. Bacteriol. 49:545-556.
Makita, Y., M. Nakao, N. Ogasawara, and K. Nakai. 2004. DBTBS: database of transcriptional regulation in Bacillus subtilis and its contribution to comparative genomics. Nucleic Acids Res. 32:D75-D77.
Marmur, J. 1961. A procedure for the isolation of deoxyribonucleic acid from microorganisms. J. Mol. Biol. 3:208-218.
Nlling, J., G. Breton, M. V. Omelchenko, K. S. Makarova, Q. Zeng, R. Gibson, H. M. Lee, J. Dubois, D. Qiu, J. Hitti, Y. I. Wolf, R. L. Tatusov, F. Sabathe, L. Doucette-Stamm, P. Soucaille, M. J. Daly, G. N. Bennett, E. V. Koonin, and D. R. Smith. 2001. Genome sequence and comparative analysis of the solvent-producing bacterium Clostridium acetobutylicum. J. Bacteriol. 183:4823-4838.
Omerchenko, M. V., K. S. Makarova, Y. I. Wolf, I. B. Rogozin, and E. V. Koonin. 2003. Evolution of mosaic operons by horizontal gene transfer and gene displacement in situ. Genome Biol. 4:R55.
Page, R. D. M. 1996. TREEVIEW: an application to display phylogenetic trees on personal computers. Comput. Appl. Biosci. 12:357-358.
Plugge, C. M., C. Dijkema, and A. J. M. Stams. 1993. Acetyl-CoA cleavage pathway in a syntrophic propionate oxidizing bacterium grown on fumarate in the absence of methanogens. FEMS Microbiol. Lett. 110:71-76.
Plugge, C. M., M. Balk, and A. J. M. Stams. 2002. Desulfotomaculum thermobenzoicum subsp. thermosyntrophicum subsp. nov., a thermophilic, syntrophic, propionate-oxidizing, spore-forming bacterium. Int. J. Syst. Evol. Microbiol. 52:391-399.
Price, M. N., K. H. Huang, A. P. Arkin, and E. J. Alm. 2005. Operon formation is driven by co-regulation and not by horizontal gene transfer. Genome Res. 15:809-819.
Saitou, N., and M. Nei. 1987. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol. Biol. Evol. 4:406-425.
Sambrook, J., E. F. Fritsch, and T. Maniatis. 1989. Molecular cloning: a laboratory manual, 2nd ed. Cold Spring Harbor Laboratory, Cold Spring Harbor, N.Y.
Schink, B. 1985. Mechanism and kinetics of succinate and propionate degradation in anoxic freshwater sediments and sewage sludge. J. Gen. Microbiol. 131:643-650.
Schink, B. 1997. Energetics of syntrophic cooperation in methanogenic degradation. Microbiol. Mol. Biol. Rev. 61:262-280.
Siebers, B., B. Tjaden, K. Michalke, C. Drr, H. Ahmed, M. Zaparty, P. Gordon, C. W. Sensen, A. Zibat, H. P. Klenk, S. C. Schuster, and R. Hensel. 2004. Reconstruction of the central carbohydrate metabolism of Thermoproteus tenax by use of genomic and biochemical data. J. Bacteriol. 186:2179-2194.
Smith, D. R., L. A. Doucette-Stamm, C. Deloughery, H. Lee, J. Dubois, T. Aldredge, R. Bashirzadeh, D. Blakely, R. Cook, K. Gilbert, D. Harrison, L. Hoang, P. Keagle, W. Lumm, B. Pothier, D. Qiu, R. Spadafora, R. Vicaire, Y. Wang, J. Wierzbowski, R. Gibson, N. Jiwani, A. Caruso, D. Bush, and J. N. Reeve. 1997. Complete genome sequence of Methanobacterium thermoautotrophicum H: functional analysis and comparative genomics. J. Bacteriol. 179:7135-7155.
Tatusov, R. L., N. D. Fedorova, J. D. Jackson, A. R. Jacobs, B. Kiryutin, E. V. Koonin, D. M. Krylov, R. Mazumder, S. L. Mekhedov, A. N. Nikolskaya, B. S. Rao, S. Smirnov, A. V. Sverdlov, S. Vasudevan, Y. I. Wolf, J. J. Yin, and D. A. Natale. 2003. The COG database: an updated version includes eukaryotes. BMC Bioinformatics 4:41.
Taylor, B. L., and I. B. Zhulin. 1999. PAS domains: internal sensors of oxygen, redox potential, and light. Microbiol. Mol. Biol. Rev. 63:479-506.
Thauer, R. K., K. Jungermann, and K. Decker. 1977. Energy conservation in chemotrophic anaerobic bacteria. Bacteriol. Rev. 41:100-180.
Thompson, J. D., T. J. Gibson, F. Plewniak, F. Jeanmougin, and D. G. Higgins. 1997. The ClustalX windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 24:4876-4882.
Van Kuijk, B. L., E. Schlosser, and A. J. Stams. 1998. Investigation of the fumarate metabolism of the syntrophic propionate-oxidizing bacterium strain MPOB. Arch. Microbiol. 169:346-352.
Wallrabenstein, C., E. Hauschild, and B. Schink. 1995. Syntrophobacter pfennigii sp. nov., new syntrophically propionate-oxidizing anaerobe growing in pure culture with propionate and sulfate. Arch. Microbiol. 164:346-352.
Wang, L. F., and R. H. Doi. 1986. Nucleotide sequence and organization of Bacillus subtilis RNA polymerase major sigma (sigma 43) operon. Nucleic Acids Res. 14:4293-4307.
Zhang, J., and T. L. Madden. 1997. PowerBLAST: a new network BLAST application for interactive or automated sequence analysis and annotation. Genome Res. 7:649-656.(Tomoyuki Kosaka, Taku Uch)
ABSTRACT
Obligate anaerobic bacteria fermenting volatile fatty acids in syntrophic association with methanogenic archaea share the intermediate bottleneck step in organic-matter decomposition. These organisms (called syntrophs) are biologically significant in terms of their growth at the thermodynamic limit and are considered to be the ideal model to address bioenergetic concepts. We conducted genomic and proteomic analyses of the thermophilic propionate-oxidizing syntroph Pelotomaculum thermopropionicum to obtain the genetic basis for its central catabolic pathway. Draft sequencing and subsequent targeted gap closing identified all genes necessary for reconstructing its propionate-oxidizing pathway (i.e., methylmalonyl coenzyme A pathway). Characteristics of this pathway include the following. (i) The initial two steps are linked to later steps via transferases. (ii) Each of the last three steps can be catalyzed by two different types of enzymes. It was also revealed that many genes for the propionate-oxidizing pathway, except for those for propionate coenzyme A transferase and succinate dehydrogenase, were present in an operon-like cluster and accompanied by multiple promoter sequences and a putative gene for a transcriptional regulator. Proteomic analysis showed that enzymes in this pathway were up-regulated when grown on propionate; of these enzymes, regulation of fumarase was the most stringent. We discuss this tendency of expression regulation based on the genetic organization of the open reading frame cluster. Results suggest that fumarase is the central metabolic switch controlling the metabolic flow and energy conservation in this syntroph.
INTRODUCTION
Methane fermentation is a complex microbiological process and involves metabolic interactions among a variety of microorganisms. These microorganisms are classified into at least four trophic groups, namely, primary fermenting bacteria, secondary fermenting bacteria, hydrogenotrophic methanogenic archaea, and acetoclastic methanogenic archaea (39). Among these microorganisms, secondary fermenting bacteria oxidize volatile fatty acids and alcohols in syntrophic association with methanogenic archaea and are thus called syntrophic bacteria (or syntrophs). Despite plenty of genetic and genomic information being available for primary fermenting bacteria (18, 30) and methanogenic archaea (11, 13, 41), such information is currently scarce for syntrophic bacteria.
A significant feature of syntrophic bacteria is their growth at the thermodynamic limit. For instance, the Gibbs free energy (G°) change in syntrophic propionate oxidation is approximately –25 kJ mol–1, which is less than the energy needed for synthesizing one ATP molecule (44). Syntrophs and methanogens share this energy for their growth. More surprisingly, Jackson and McInerney (22) have shown that substrate oxidation by syntrophs proceeds at values close to the thermodynamic equilibrium (G° 0 kJ mol–1). Therefore, one can deduce that these bacteria should have extremely efficient catabolic systems (22), and their catabolic pathway is of wide biological interest.
Propionate-oxidizing syntrophs have been isolated from mesophilic (4, 16, 27, 47) and thermophilic (21, 34) methanogenic ecosystems. Two propionate-oxidizing pathways have been proposed for mesophilic syntrophic bacteria, i.e., the methylmalonyl-coenzyme A (CoA) pathway (see Fig. 1) and a pathway via a six-carbon intermediate metabolite. The methylmalonyl-CoA (MMC) pathway has first been proposed for methanogenic freshwater sediment (38) and later confirmed for syntrophic propionate-oxidizing bacteria on the basis of the results of 13C-nuclear magnetic resonance analyses of intermediate metabolites (33) and enzyme activity measurements (19). The second pathway has been proposed for Smithella propionica (9), which produces acetate and butyrate via a six-carbon intermediate metabolite. In addition, a recent study has detected activities of several enzymes in the MMC pathway in a thermophilic syntrophic propionate-oxidizing bacterium (34). However, there has been no information regarding the genetics of syntrophic propionate oxidation.
The present study was conducted to gain the genetic basis for the central catabolic pathway in a thermophilic syntrophic propionate-oxidizing bacterium, Pelotomaculum thermopropionicum strain SI (20, 21). Recent advances in genomics coupled to database-associated bioinformatics have allowed scientists to identify microbial metabolic pathways on the basis of genome sequences (18, 26, 40). Accordingly, we performed genomic and proteomic analyses of P. thermopropionicum, which revealed novel insights into how this organism controls the central metabolic flow and energy production.
MATERIALS AND METHODS
Strains and growth conditions. P. thermopropionicum strain SI (DSM 13744) was grown in cultures alone and in cocultures with Methanothermobacter thermautotrophicus H (DSM 1053) using a culture medium described elsewhere (21). The medium was supplemented with 0.1% Bacto yeast extract (Difco) and an appropriate fermentation substrate at 20 mM. Cultivation was conducted at 55°C under an atmosphere of N2 plus CO2 (80/20 [vol/vol]) without shaking.
DNA isolation and genome-size analysis. P. thermopropionicum was grown alone in 1 liter of medium supplemented with fumarate as the sole substrate. Cells were harvested at the late exponential growth phase by centrifugation, and the total DNA was extracted by the method of Marmur (29). To check for the presence of plasmids, the purified DNA was electrophoresed in agarose gels by the standard procedure (37). To determine the genome size, pulsed-field gel electrophoresis was conducted using a contour-clamped homogeneous electric field (CHEF DRIII; Bio-Rad Laboratories) according to the manufacturer's instructions.
Library construction, shotgun sequencing, and assembly. The procedures for draft genome sequencing, namely, shotgun and fosmid library construction, sequencing, and contig assembly, were conducted by Dragon Genomics. Two types of shotgun library were constructed, i.e., a plasmid library containing 2- to 3-kb inserts and a fosmid library containing 30- to 50-kb inserts. For constructing a plasmid library, the extracted DNA was mechanically sheared, blunted, and ligated into HincII-digested pUC118 (Takara). For constructing a fosmid library, the extracted DNA was fractured by the phenol-chloroform treatment, fragments of approximately 33 to 48 kb were recovered by electrophoresis, and they were ligated into pCC1FOS (Epicenter) using a copy control fosmid library production kit (Epicenter). Clones were selected from the plasmid (15,360 clones) and fosmid (960 clones) libraries and subjected to sequencing from both ends by using a DYEnamic ET dye terminator kit (Amersham Bioscience) and a MegaBASE 4000 sequencer (Amersham Bioscience). Raw sequencing data (quality values in the Phred analysis [CodonCode] of more than 15) were assembled to generate "contigs" by using the Paracel genome assembler with the CAP4 algorithm (Paracel). Physically linked contigs were tentatively connected to generate "scaffolds." A scaffold was named CWW_CXX, where CWW and CXX were identification (ID) numbers of the first and last contigs. Gaps were closed where needed by direct sequencing of an PCR-recovered intercontig fragment.
Gene prediction and annotation. The assembled sequences were first analyzed by the GRIMMER 2.10 program (10) trained with the genome sequence of Clostridium acetobutylicum ATCC 824 (30). Open reading frames (ORFs) detected were subjected to the BLAST search (49) against the Swiss-Prot (3) and COG (Clusters of Orthologous Groups of proteins) (42) databases to select the ORFs that satisfy the following criteria: sizes larger than 50 amino acids, E values to best-match sequences smaller than 1e–100, best-match sequences being not hypothetical or putative proteins, and the initial and last amino acids being identical to those of the best-match proteins. The GRIMMER program was next trained using the obtained ORFs and again used for analyzing the assembled sequences. After this procedure was repeated three times in total, selected ORFs were subjected to homology search against the COG database using the BLASTp program (1) and the nonredundant GenBank database using the BLASTx program (1). Results of the homology search were linked to the scaffold information and stored in Annotation Viewer (Dragon Genomics) and FileMaker pro version 3.0 (Claris). We also constructed an P. thermopropionicum ORF database available for the BLASTp search. The naming scheme for an ORF was CWW_CXX_YY_ZZ (ORF ID number), where CWW_CXX was the ID number of a scaffold, while YY and ZZ were the base numbers of the first and last nucleotides of the ORF.
Functions of ORFs were also analyzed by using the InterPro program (2) linked to the PROSITE, PRINTS, Pfam, SMART, TIGRFAMs, PIR SuperFamily, and ProDom databases to search for conserved domains and motifs and to validate prediction of gene function. The secondary structure and membrane topology were analyzed using the nnPredict (25) and TMpred (17) programs, respectively. Transcription regulation sequences were primarily analyzed using the DBTBS database constructed for analyzing the Bacillus subtilis genome (28), and the results were manually checked. Terminator sequences were analyzed using the FindTerm program (SoftBerry). Molecular weights and isoelectric points (pI) of proteins were predicted using the Compute pI/Mw tool (14). Alignment of amino acid sequences was conducted by using the ClustalX program (45) for visually inspecting conservation of motifs and functionally important amino acids. For phylogenetic analyses, a neighbor-joining tree (36) was constructed by using the njplot software in the ClustalX program and plotted by using TREEVIEW (32).
Proteomics. P. thermopropionicum cells were grown in cocultures (500 ml) with M. thermautotrophicus as described above, harvested by centrifugation at 11,000 x g and 4°C for 10 min, and washed two times using 5 mM sodium phosphate buffer (pH 7.0). Each cell pellet was suspended in 25 ml of an Percoll solution (Amersham Bioscience) supplemented with 0.15 M NaCl and subjected to centrifugation at 10,000 x g and 4°C for 30 min. After the centrifugation, two distinct bands were seen; microscopic observation showed that the upper band was mostly comprised of cells of P. thermopropionicum. We carefully collected this fraction using a pipette, and Percoll particles were removed by washing with the phosphate buffer. The harvested cells were disrupted by sonication (Sonifier 250; Branson), treated with an RNase-DNase solution (1 mg ml–1 DNase I, 0.25 mg ml–1 RNase A, 0.5 M Tris-HCl [pH 7.0], and 50 mM MgCl2), and lysed in a lysis solution {5 M urea, 2 M thiourea, 4% 3-[(3-cholamidopropyl)dimethylammonio]-1-propane-sulfonate, 80 mM dithiothreitol, and 4% IPG buffer (Amersham Bioscience)}. The resultant lysate was centrifuged at 30,000 x g and 4°C for 20 min, and the supernatant was collected. The protein concentration in this supernatant solution was determined using a two-dimensional (2D) Quant kit (Amersham Bioscience).
Two-dimensional gel electrophoresis (2D-GE) was performed using an investigator 2D electrophoresis system (Genomic Solution). Immobiline DryStrip pH 4 to 7 (Amersham Bioscience) was used for the first-dimensional isoelectric focusing, to which approximately 50 μg of protein was loaded. The second-dimensional separation was performed in 12.5% polyacrylamide gels (PAGs) containing 0.1% sodium dodecyl sulfate (SDS). SDS-PAG electrophoresis molecular weight standard low range (Bio-Rad) was used as a size marker. The gels were stained using Sypro Ruby (Molecular Probes) and visualized using the FM BIOII Multi-view system (Hitachi). We performed 2D-GE more than three times for each culture condition and confirmed the reproducibility. Protein spots were quantified using an PDQuest system (Bio-Rad).
For the N-terminal sequencing (NTS) identification, protein spots on an SDS-PAG were transferred to a ProBlott membrane (Applied Biosystems) using an electro membrane blotter (Nippon Eido). Spots were excised from the membrane, and an N-terminal sequence was determined using a 494cLC protein sequencer (Applied Biosystems). For the peptide mass fingerprinting (PMF) identification, a gel was stained using a silver stain kit (Wako). Protein spots were excised and subjected to in-gel trypsin digestion by the method of Katayama et al. (24). Time-of-flight mass spectrometry was performed using a Voyager DE-PRO time-flight mass spectrometer (Applied Biosystems) and -cyano-4-hydroxycinnamic acid (Sigma) as a matrix. Peptide mass peaks were calibrated using the Data Explorer DM software (Applied Biosystems). Mascot search (Matrix Science) was performed using the P. thermopropionicum draft sequence database for identifying protein spots.
Nucleotide sequence accession numbers. The draft genome sequence of P. thermopropionicum has been deposited in the DDBJ, EMBL, and NCBI databases under accession numbers BAAC01000001 to BAAC01000195. The accession numbers for nucleotide sequences of Pct3, MmcD, and OdcAB are AB221127, AB221128, and AB221129, respectively.
RESULTS
Draft genome sequencing. Pulsed-field gel electrophoresis determined the genome size of P. thermopropionicum to be approximately 3 Mb (data not shown). The genome was single and circular, and no plasmid was found. We performed shotgun sequencing of its genome after constructing plasmid (approximately 2-kb inserts) and fosmid (approximately 42-kb inserts) libraries. A total of 28,453 reads (approximately 17 Mb in total, x5.5 coverage) generated 105 scaffolds and an assembled sequence of 2,810,264 bp, which was estimated to cover up to 94% of the total genome. Its G+C content was in good agreement with an estimate in a hybridization test (52.8% (21)). In the draft sequence, we found 2,988 of ORFs, of which approximately 70% could be grouped into the COG (42) categories (data not shown). On the basis of the current draft data, we constructed ORF databases available for keyword and BLAST searches (these databases will be provided on request).
Reconstruction of the central catabolic pathway. P. thermopropionicum is capable of anaerobic growth on propionate, propanol, butanol, ethanol, and lactate in cocultures with a hydrogenotrophic methanogen, such as Methanothermobacter thermautotrophicus H, while it can also grow on fumarate and pyruvate in culture alone (21). It is known that this bacterium produces acetate and propionate (at a molar ratio of 3:1) from pyruvate (21). Since the available substrates and fermentation products are either constituents of the MMC pathway or compounds directly linked to this pathway, we assumed that P. thermopropionicum likely utilized the MMC pathway (Fig. 1) for metabolizing propionate.
No genetic information had been available for the MMC pathway. However, since genes for enzymes functionally identical to those in this pathway have been identified in other organisms, such as genes for enzymes in the citrate cycle in aerobic bacteria and those in the propionate fermentation pathway in Propionibacterium acnes (5), we were able to search for putative genes for enzymes in the MMC pathway. We picked ORFs of interest by keyword and BLAST search against the P. thermopropionicum ORF databases, for which we selected search terms and sequences according to information obtained from the KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway database (23). The selected ORFs were further analyzed for the conservation of signature motifs and important amino acids and, in some cases, similarities in the overall structures. ORFs in the draft genome sequence for Pct3, MmcD, and OdcA (Table 1) were found over two contigs, whose sequences were later completed and independently deposited to the databases (information available from authors upon request).
Table 1 summarizes results of the ORF search, showing that ORFs for all enzymes in the MMC pathway (Fig. 1) were identified. Of the enzymes in the MMC pathway that are different from previously known enzymes, methylmalonyl-CoA mutase (MCM) and fumarate hydratase (FHT) (fumarase) seem to be composed of two different subunits (information available from authors upon request). The two-subunit FHT enzyme was expressed in Escherichia coli, and its fumarate-transforming activity was detected at 55°C (our unpublished result). We also found an ORF for methylmalonyl-CoA epimerase (MCE) [which catalyzes isomerization between (R)-methylmalonyl-CoA and (S)-methylmalonyl-CoA], which should be necessary prior to step III (Fig. 1). Pyruvate:ferredoxin oxidoreductase (POR), rather than pyruvate dehydrogenase (PDH), was found in the P. thermopropionicum genome (Table 1), which is known to utilize low-potential ferredoxin as an electron acceptor. It was also found that each of the last three steps in this pathway, oxaloacetate decarboxylation (step VIII), acetyl-CoA formation (step IX), and acetyl-CoA hydrolysis (step X), can be catalyzed by two different enzymes (Table 1). According to the genomic information, we reconstructed the MMC pathway in P. thermopropionicum (Fig. 2).
An interesting feature found in the sequence analysis was that many ORFs for the enzymes in the MMC pathway, except for those for propionate CoA transferase (PCT) and SDH, formed a tight cluster on the genome (mmcBCDEFGHIJKLM in Fig. 3), which we named the mmc cluster. We also found two ORFs (tps and mmcA) upstream of this cluster (Table 1 and Fig. 3). MmcA is homologous (37% identical in total amino acids) to RocR, an NtrC/NifA family transcriptional regulator in B. subtilis (7). Analyses of two long intergenic regions (mmcA-mmcB [274 bp] and mmcC-mmcD [452 bp]) performed using the DBTBS database (28) found that there existed a consensus sequence for a sigma L-dependent promoter (8) in the mmcA-mmcB intergenic region, while consensus sequences for sigma H-dependent promoters (12) and sigma A-dependent promoters (48) were present in the mmcC-mmcD region (data not shown).
Peripheral pathways. We also found ORFs that showed substantial homologies to genes for enzymes catalyzing peripheral pathways of the MMC pathways (Table 1 and Fig. 2), such as lactate dehydrogenase (LDH), 2-oxoglutarate:ferredoxin oxidoreductase (KOR), and phosphoenolpyruvate synthetase (PPS). LDH is considered to be used for growth of this bacterium on lactate, while KOR and PPS are important for linking the MMC pathway to biosyntheses of cellular components.
P. thermopropionicum is capable of growth on several short-chain alcohols, such as ethanol, propanol, and butanol. The end product from butanol is butyrate, while propionate produced from propanol is further transformed to acetate. Our growth test showed that acetate and propionate were produced from ethanol (data not shown). In the genome of P. thermopropionicum, multiple ORFs for alcohol dehydrogenase (ADH) and aldehyde:ferredoxin oxidoreductase (AOR) and one ORF for alcohol/aldehyde dehydrogenase (AAD) were identified (Table 1 and Fig. 2). A combination of ADH and AOR produces acetate from ethanol, while AAD directly produces acetyl-CoA. Substrates for these ADHs and AORs need to be experimentally identified.
We found ORFs for four hydrogenases (three Fe-only hydrogenases and one NiFeSe hydrogenase [Table 1]) in the draft genome; these ORFs may be important for scavenging reducing equivalents produced in the MMC pathway and alcohol metabolism as molecular hydrogen. As electron carriers, six ORFs for ferredoxins were found (Table 1), while there was no ORF for a soluble c-type cytochrome. A complete set of ORFs for bacterial F-type ATPase was also identified (Table 1).
Proteomic analyses. In order to detect proteins encoded by ORFs in the mmc cluster, we carried out 2D-GE proteomic analysis of soluble proteins in cells grown on propionate or butanol (Fig. 4). In this analysis, the amounts of proteins loaded on a gel were optimized for evaluating the levels of expression of major soluble proteins. The butanol culture was chosen as the control, because its catabolism is not linked to the MMC pathway. We used the draft genome data to identify major protein spots by PMF and/or NTS. As summarized in Table 2, we were able to detect subunits of propionyl-CoA:oxaloacetase transcarboxylase (POT) (McmH), MCM (MmcE), succinyl-CoA synthetase (SCS) (MmcD1), malate dehydrogenase (MDH) (MmcK), FHT (MmcB), and succinate dehydrogenase (SDH) (Sdh1C) as major spots in the gel. MCE and POR were undetected, probably because their sizes were out of the range of the gel. However, we assume that they were also expressed, since the genes of these enzymes were present in the same transcriptional unit (Fig. 3). In addition, we could not detect the predicted PCT (Pct3 in Table 1), whereas a subunit of an alternative putative CoA transferase (spot 17) was detected. This transferase is homologous (31% identical in amino acids) to GctA, a subunit of glutaconate CoA-transferase from Acidaminococcus fermentans (CAA57199). Although spot 17 was not specifically expressed in the propionate culture, we also need to consider a possibility that this protein functions as PCT.
We next determined the relative intensity of each protein spot (the intensity of a spot was divided by the sum of intensities of all spots in a gel) and compared them between the propionate and butanol cultures (Table 2). MmcK and MmcB each appeared as two spots (spots 15 and 16 and spots 18 and 19, respectively), whose average ratios in the relative intensity were 2.3 and 4.8, respectively. This analysis showed that the enzymes in the MMC pathway were generally up-regulated in the propionate culture and that the up-regulation of MmcB was the most stringent. These tendencies of changes in the spot intensity were also observed when the propionate culture was compared with the propanol or ethanol culture (data not shown).
DISCUSSION
In the present study, we carried out draft sequencing of the genome of P. thermopropionicum to obtain the genetic basis for its central metabolic pathway. In combination with the results of the proteomic analysis, we conclude that this thermophilic syntroph utilizes the MMC pathway for propionate oxidation. The draft sequencing and subsequent targeted gap closing provided the complete sequences for genes encoding the enzymes necessary for reconstructing the MMC pathway, which will assist in subsequent studies for biochemical characterization of the enzymes and molecular analyses of transcriptional regulation mechanisms. Below, we discuss some aspects of the central catabolic pathway in P. thermopropionicum on the basis of current genomic information.
A similar pathway has been identified as a propionate-fermenting pathway in the genome of Propionibacterium acnes (5), while we found differences between the MMC pathway in P. thermopropionicum and the propionate-fermenting pathway in P. acnes. First, P. acnes utilizes a malic enzyme catalyzing the conversion between malate and pyruvate (malate + NAD+ = pyruvate + CO2 + NADH) for steps VII and VIII. Second, PCT is used for step I in P. thermopropionicum, while the same reaction is catalyzed by propionyl-CoA synthetase in P. acnes.
The first two steps in the MMC pathways in P. thermopropionicum (steps I and II in Fig. 1) are catalyzed by transferases, i.e., PCT and POT. These enzymes require counterpart substrates for their activities (CoA derivatives, e.g., acetyl-CoA, for PCT and oxaloacetate for POT), while these counterparts are also the downstream intermediate metabolites in the MMC pathway. We think that this is one of the reasons why growth of P. thermopropionicum was so slow. As reported previously (20), growth of this bacterium on propionate in coculture with a hydrogenotrophic methanogen occurred after a long lag period of 20 days. It is likely that this lag period is necessary for P. thermopropionicum to accumulate these counterpart substrates, and propionate oxidation is initiated only after the concentrations of these substrates reach certain levels. Another possibility for the long lag period would be that the cells need much time for the expression of the catabolic enzymes under the energy-limiting conditions.
Proteomic analyses (Fig. 4 and Table 2) showed that a subunit of fumarase (MmcB) was up-regulated in the propionate culture compared to the butanol culture. Other proteins identified to be encoded in the mmc cluster were also up-regulated in the propionate culture, although up-regulation of these proteins was weaker than that of MmcB. Northern blot analysis confirmed this tendency of gene expression regulation (our unpublished results). We think that this tendency reflected the genetic organization of the mmc cluster. As shown in Fig. 3, the sequence analysis suggested that transcription of mmcBC is solely regulated by the sigma L-dependent promoter, while that of mmcDEFGHIJKLM is also under the control of the sigma A and H-dependent promoters. Sigma A is known to be the major sigma factor in B. subtilis (48), which corresponds to sigma 70 in gram-negative bacteria (15). Sigma A has been considered to control expression of approximately 250 operons in B. subtilis, many of which encode housekeeping proteins (the DBTBS database). Sigma H is a nonessential sigma factor involved in the expression of vegetative and early stationary-phase genes in B. subtilis (12). We found 19 ORFs for sigma factors in the draft genome of P. thermopropionicum; these ORFs include 1 ORF for sigma A, 2 ORFs for sigma H, and 1ORF for sigma L (data not shown). These findings suggest that mmcDEFGHIJKLM are transcribed at certain levels under a wide range of growth conditions. In contrast, sigma L is a member of the sigma 54 family of bacterial sigma factors, which are characterized by their cooperative action with enhancer-binding proteins whose function requires nucleotide hydrolysis (6). Many sigma 54-controlled genes are specifically transcribed in the presence of enhancer molecules (6). In the case of the mmc cluster, it is likely that MmcA functions as an enhancer-binding protein, although we have not yet identified what molecules (or conditions) actually enhanced the expression of mmcBC. Since MmcA includes the PAS domain (43), we deduce that expression of MmcB is enhanced by a global physiological status, e.g., the membrane potential (43), rather than a specific substrate of the catabolic pathway. Up-regulation of the other MMC proteins may also be governed by the sigma L-dependent promoter, although its regulation may be weakened by the other promoters.
As discussed above, the genome analyses suggested that the particular genetic organization of P. thermopropionicum has evolved to accomplish the hierarchical transcriptional regulation of its central metabolic pathway. We think that the stringent regulation of fumarase is physiologically expedient in terms of multiple roles of fumarate in that system. First, since fumarate is an intermediate metabolite of the propionate-oxidizing pathway, fumarase should be up-regulated together with other MMC enzymes, when the bacterium metabolizes propionate. In addition, fumarate is a growth substrate in the culture of P. thermopropionicum alone, where P. thermopropionicum primarily oxidizes it to acetate, while this organism also utilizes it as an electron acceptor to form succinate (namely, fumarate disproportionation [46]) (Fig. 2). As shown in Fig. 4 and Table 2, the expression of fumarase was conditional, while that of fumarate reductase (i.e., Sdh1C, spot 20 in Fig. 4) did not show much change. It is beneficial for P. thermopropionicum to mainly utilize fumarate as the electron acceptor when the MMC pathway does not efficiently work, because fumarate respiration can produce the membrane proton gradient. Proteomics (Table 2 and Fig. 4) showed that P. thermopropionicum expressed a significant amount of ATPase, suggesting that this organism can utilize reducing equivalents produced by fermentation for fumarate respiration and oxidative phosphorylation. It is worth noting that two putative SDHs are present in this organism (Table 1); their roles will be addressed in future biochemical studies. All together, we suggest that fumarase acts as the central metabolic switch controlling the metabolic flow and energy conservation in this syntrophic bacterium.
In conclusion, the genomic analyses of P. thermopropionicum offer a glimpse into the sophisticated catabolic system of syntrophic bacteria. Owing to its unique genetic organization, the mmc cluster may also be interesting for scientists who investigate the evolutionary aspects of operon formation (31, 35). We are currently conducting genetic and biochemical experiments for more-quantitative evaluations of the contribution of each promoter under different growth conditions.
ACKNOWLEDGMENTS
We thank Yasuo Igarashi (University of Tokyo) for helpful advice and continuous encouragement. Acknowledgment is also made to Youhei Urushida for N-terminal amino acid sequencing and Reiko Hirano for technical assistance.
This work was supported by New Energy and Industrial Technology Development Organization (NEDO).
REFERENCES
Altschul, S., T. Madden, A. Schaffer, J. Zhang, Z. Zhang, W. Miller, and D. J. Lipman. 1997. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 25:3389-3402.
Apweiler, R., T. K. Attwood, A. Bairoch, A. Bateman, E. Birney, M. Biswas, P. Bucher, L. Cerutti, F. Corpet, M. D. Croning, R. Durbin, L. Falquet, W. Fleischmann, J. Gouzy, H. Hermjakob, N. Hulo, I. Jonassen, D. Kahn, A. Kanapin, Y. Karavidopoulou, R. Lopez, B. Marx, N. J. Mulder, T. M. Oinn, M. Pagni, F. Servant, C. J. Sigrist, and E. M. Zdobnov. 2001. The InterPro database, an integrated documentation resource for protein families, domains and functional sites. Nucleic Acids Res. 29:37-40.
Boeckmann, B., A. Bairoch, R. Apweiler, M. Blatter, A. Estreicher, E. Gasteiger, M. Martin, K. Michoud, C. O'Donovan, I. Phan, S. Pilbout, and M. Schneider. 2003. The SWISS-PROT protein knowledge base and its supplement TrEMBL in 2003. Nucleic Acids Res. 31:365-370.
Boone, D. R., and M. P. Bryant. 1980. Propionate-degrading bacterium, Syntrophobacter wolinii sp. nov., gen. nov., from methanogenic ecosystems. Appl. Environ. Microbiol. 40:626-632.
Brüggemann, H., A. Henne, F. Hoster, H. Liesegang, A. Wiezer, A. Strittmatter, S. Hujer, P. Durre, and G. Gottschalk. 2004. The complete genome sequence of Propionibacterium acnes, a commensal of human skin. Science 305:671-673.
Buck, M., M. T. Gallegos, D. J. Studholme, Y. Guo, and J. D. Gralla. 2000. The bacterial enhancer-dependent 54 (N) transcription factor. J. Bacteriol. 182:4129-4136.
Calogero, S., R. Gardan, P. Glaser, J. Schweizer, G. Rapoport, and M. Debarbouille. 1994. RocR, a novel regulatory protein controlling arginine utilization in Bacillus subtilis, belongs to the NtrC/NifA family of transcriptional activators. J. Bacteriol. 176:1234-1241.
Debarbouille, M., I. Martin-Verstraete, F. Kunst, and G. Rapoport. 1991. The Bacillus subtilis sigL gene encodes an equivalent of sigma 54 from Gram-negative bacteria. Proc. Natl. Acad. Sci. USA 88:9092-9096.
De Bok, F. A., A. J. Stams, C. Dijkema, and D. R. Boone. 2001. Pathway of propionate oxidation by a syntrophic culture of Smithella propionica and Methanospirillum hungatei. Appl. Environ. Microbiol. 67:1800-1804.
Delcher, A. L., D. Harmon, S. Kasif, O. White, and S. L. Salzberg. 1999. Improved microbial gene identification with GLIMMER. Nucleic Acids Res. 27:4636-4641.
Deppenmeier, U., A. Johann, T. Hartsch, R. Merkl, R. A. Schmitz, R. Martinez-Arias, A. Henne, A. Wiezer, S. Baumer, C. Jacobi, H. Bruggemann, T. Lienard, A. Christmann, M. Bomeke, S. Steckel, A. Bhattacharyya, A. Lykidis, R. Overbeek, H. P. Klenk, R. P. Gunsalus, H. J. Fritz, and G. Gottschalk. 2002. The genome of Methanosarcina mazei: evidence for lateral gene transfer between Bacteria and Archaea. J. Mol. Microbiol. Biotechnol. 4:453-461.
Dubnau, E., J. Weir, G. Nair, L. Carter III, C. Moran, Jr., and I. Smith. 1988. Bacillus sporulation gene spo0H codes for 30 (H). J. Bacteriol. 170:1054-1062.
Galagan, J. E., C. Nusbaum, A. Roy, M. G. Endrizzi, P. Macdonald, W. FitzHugh, S. Calvo, R. Engels, S. Smirnov, D. Atnoor, A. Brown, N. Allen, J. Naylor, N. Stange-Thomann, K. DeArellano, R. Johnson, L. Linton, P. McEwan, K. McKernan, J. Talamas, A. Tirrell, W. Ye, A. Zimmer, R. D. Barber, I. Cann, D. E. Graham, D. A. Grahame, A. M. Guss, R. Hedderich, C. Ingram-Smith, H. C. Kuettner, J. A. Krzycki, J. A. Leigh, W. Li, J. Liu, B. Mukhopadhyay, J. N. Reeve, K. Smith, T. A. Springer, L. A. Umayam, O. White, R. H. White, E. Conway de Macario, J. G. Ferry, K. F. Jarrell, H. Jing, A. J. Macario, I. Paulsen, M. Pritchett, K. R. Sowers, R. V. Swanson, S. H. Zinder, E. Lander, W. W. Metcalf, and B. Birren. 2002. The genome of M. acetivorans reveals extensive metabolic and physiological diversity. Genome Res. 12:532-542.
Gasteiger, E., C. Hoogland, A. Gattiker, S. Duvaud, M. R. Wilkins, R. D. Appel, and A. Bairoch. 2005. Protein identification and analysis tools on the ExPASy server, p. 571-607. In J. M. Walker (ed.), The proteomics protocols handbook. Humana Press, Totowa, N.J.
Gitt, M. A., L. F. Wang, and R. H. Doi. 1985. A strong sequence homology exists between the major RNA polymerase sigma factors of Bacillus subtilis and Escherichia coli. J. Biol. Chem. 260:7178-7185.
Harmsen, H. J., B. L. van Kuijk, C. M. Plugge, A. D. Akkermans, W. M. de Vos, and A. J. Stams. 1998. Syntrophobacter fumaroxidans sp. nov., a syntrophic propionate-degrading sulfate-reducing bacterium. Int. J. Syst. Bacteriol. 48:1383-1387.
Hofmann, K., and W. Stoffel. 1993. TMPred—a database of membrane spanning protein segments. Biol. Chem. Hoppe-Seyler 347:166.
Hong, S. H., J. S. Kim, S. Y. Lee, Y. H. In, S. S. Choi, J. K. Rih, C. H. Kim, H. Jeong, C. G. Hur, and J. J. Kim. 2004. The genome sequence of the capnophilic rumen bacterium Mannheimia succiniciproducens. Nat. Biotechnol. 22:1275-1281.
Houwen, F. P., J. Plokker, A. J. M. Stams, and A. J. B. Zehnder. 1990. Enzymatic evidence for involvement of the methylmalonyl-CoA pathway in propionate oxidation by Syntrophobacter wolinii. Arch. Microbiol. 155:52-55.
Imachi, H., Y. Sekiguchi, Y. Kamagata, A. Ohashi, and H. Harada. 2000. Cultivation and in situ detection of a thermophilic bacterium capable of oxidizing propionate in syntrophic association with hydrogenotrophic methanogens in a thermophilic methanogenic granular sludge. Appl. Environ. Microbiol. 66:3608-3615.
Imachi, H., Y. Sekiguchi, Y. Kamagata, A. Ohashi, and H. Harada. 2002. Pelotomaculum thermopropionicum gen. nov., sp. nov., an anaerobic, thermophilic, syntrophic propionate-oxidizing bacterium. Int. J. Syst. Evol. Microbiol. 52:1729-1735.
Jackson, B. E., and M. J. McInerney. 2002. Anaerobic microbial metabolism can proceed close to thermodynamic limits. Nature 415:454-456.
Kanehisa, M., and S. Goto. 2000. KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 28:27-30.
Katayama, H., T. Nagasu, and Y. Oda. 2001. Improvement of in-gel digestion protocol for peptide mass fingerprinting by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry. Rapid Commun. Mass Spectrom. 15:1416-1421.
Kneller, D. G., F. E. Cohen, and R. Langridge. 1990. Improvements in protein secondary structure prediction by an enhanced neural network. J. Mol. Biol. 214:171-182.
Larimer, F. W., P. Chain, L. Hauser, J. Lamerdin, S. Malfatti, L. Do, M. L. Land, D. A. Pelletier, J. T. Beatty, A. S. Lang, F. R. Tabita, J. L. Gibson, T. E. Hanson, C. Bobst, J. L. T. Torres, C. Peres, F. H. Harrison, J. Gibson, and C. S. Harwood. 2004. The genome sequence of the metabolically versatile photosynthetic bacterium Rhodopseudomonas palustris. Nat. Biotechnol. 22:55-61.
Liu, Y., D. L. Balkwill, H. C. Aldrich, G. R. Drake, and D. R. Boone. 1999. Characterization of the anaerobic propionate-degrading syntrophs Smithella propionica gen. nov., sp. nov. and Syntrophobacter wolinii. Int. J. Syst. Bacteriol. 49:545-556.
Makita, Y., M. Nakao, N. Ogasawara, and K. Nakai. 2004. DBTBS: database of transcriptional regulation in Bacillus subtilis and its contribution to comparative genomics. Nucleic Acids Res. 32:D75-D77.
Marmur, J. 1961. A procedure for the isolation of deoxyribonucleic acid from microorganisms. J. Mol. Biol. 3:208-218.
Nlling, J., G. Breton, M. V. Omelchenko, K. S. Makarova, Q. Zeng, R. Gibson, H. M. Lee, J. Dubois, D. Qiu, J. Hitti, Y. I. Wolf, R. L. Tatusov, F. Sabathe, L. Doucette-Stamm, P. Soucaille, M. J. Daly, G. N. Bennett, E. V. Koonin, and D. R. Smith. 2001. Genome sequence and comparative analysis of the solvent-producing bacterium Clostridium acetobutylicum. J. Bacteriol. 183:4823-4838.
Omerchenko, M. V., K. S. Makarova, Y. I. Wolf, I. B. Rogozin, and E. V. Koonin. 2003. Evolution of mosaic operons by horizontal gene transfer and gene displacement in situ. Genome Biol. 4:R55.
Page, R. D. M. 1996. TREEVIEW: an application to display phylogenetic trees on personal computers. Comput. Appl. Biosci. 12:357-358.
Plugge, C. M., C. Dijkema, and A. J. M. Stams. 1993. Acetyl-CoA cleavage pathway in a syntrophic propionate oxidizing bacterium grown on fumarate in the absence of methanogens. FEMS Microbiol. Lett. 110:71-76.
Plugge, C. M., M. Balk, and A. J. M. Stams. 2002. Desulfotomaculum thermobenzoicum subsp. thermosyntrophicum subsp. nov., a thermophilic, syntrophic, propionate-oxidizing, spore-forming bacterium. Int. J. Syst. Evol. Microbiol. 52:391-399.
Price, M. N., K. H. Huang, A. P. Arkin, and E. J. Alm. 2005. Operon formation is driven by co-regulation and not by horizontal gene transfer. Genome Res. 15:809-819.
Saitou, N., and M. Nei. 1987. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol. Biol. Evol. 4:406-425.
Sambrook, J., E. F. Fritsch, and T. Maniatis. 1989. Molecular cloning: a laboratory manual, 2nd ed. Cold Spring Harbor Laboratory, Cold Spring Harbor, N.Y.
Schink, B. 1985. Mechanism and kinetics of succinate and propionate degradation in anoxic freshwater sediments and sewage sludge. J. Gen. Microbiol. 131:643-650.
Schink, B. 1997. Energetics of syntrophic cooperation in methanogenic degradation. Microbiol. Mol. Biol. Rev. 61:262-280.
Siebers, B., B. Tjaden, K. Michalke, C. Drr, H. Ahmed, M. Zaparty, P. Gordon, C. W. Sensen, A. Zibat, H. P. Klenk, S. C. Schuster, and R. Hensel. 2004. Reconstruction of the central carbohydrate metabolism of Thermoproteus tenax by use of genomic and biochemical data. J. Bacteriol. 186:2179-2194.
Smith, D. R., L. A. Doucette-Stamm, C. Deloughery, H. Lee, J. Dubois, T. Aldredge, R. Bashirzadeh, D. Blakely, R. Cook, K. Gilbert, D. Harrison, L. Hoang, P. Keagle, W. Lumm, B. Pothier, D. Qiu, R. Spadafora, R. Vicaire, Y. Wang, J. Wierzbowski, R. Gibson, N. Jiwani, A. Caruso, D. Bush, and J. N. Reeve. 1997. Complete genome sequence of Methanobacterium thermoautotrophicum H: functional analysis and comparative genomics. J. Bacteriol. 179:7135-7155.
Tatusov, R. L., N. D. Fedorova, J. D. Jackson, A. R. Jacobs, B. Kiryutin, E. V. Koonin, D. M. Krylov, R. Mazumder, S. L. Mekhedov, A. N. Nikolskaya, B. S. Rao, S. Smirnov, A. V. Sverdlov, S. Vasudevan, Y. I. Wolf, J. J. Yin, and D. A. Natale. 2003. The COG database: an updated version includes eukaryotes. BMC Bioinformatics 4:41.
Taylor, B. L., and I. B. Zhulin. 1999. PAS domains: internal sensors of oxygen, redox potential, and light. Microbiol. Mol. Biol. Rev. 63:479-506.
Thauer, R. K., K. Jungermann, and K. Decker. 1977. Energy conservation in chemotrophic anaerobic bacteria. Bacteriol. Rev. 41:100-180.
Thompson, J. D., T. J. Gibson, F. Plewniak, F. Jeanmougin, and D. G. Higgins. 1997. The ClustalX windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 24:4876-4882.
Van Kuijk, B. L., E. Schlosser, and A. J. Stams. 1998. Investigation of the fumarate metabolism of the syntrophic propionate-oxidizing bacterium strain MPOB. Arch. Microbiol. 169:346-352.
Wallrabenstein, C., E. Hauschild, and B. Schink. 1995. Syntrophobacter pfennigii sp. nov., new syntrophically propionate-oxidizing anaerobe growing in pure culture with propionate and sulfate. Arch. Microbiol. 164:346-352.
Wang, L. F., and R. H. Doi. 1986. Nucleotide sequence and organization of Bacillus subtilis RNA polymerase major sigma (sigma 43) operon. Nucleic Acids Res. 14:4293-4307.
Zhang, J., and T. L. Madden. 1997. PowerBLAST: a new network BLAST application for interactive or automated sequence analysis and annotation. Genome Res. 7:649-656.(Tomoyuki Kosaka, Taku Uch)