Previous Article | Next Article ![]()
Journal of Bacteriology, January 2006, p. 202-210, Vol. 188, No. 1
0021-9193/06/$08.00+0 doi:10.1128/JB.188.1.202-210.2006
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
Laboratory of Applied Microbiology, Marine Biotechnology Institute, Kamaishi, Iwate 026-0001,1 National Institute of Bioscience and Human-Technology, Agency of Industrial Science and Technology, Tsukuba, Ibaraki 305-8566,2 Department of Environmental Systems Engineering, Nagaoka University of Technology, Nagaoka, Niigata 940-2188, Japan3
Received 31 August 2005/ Accepted 17 October 2005
|
|
|---|
|
|
|---|
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 mol1, 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 mol1). 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.
![]() View larger version (32K): [in a new window] |
FIG. 1. Proposed methylmalonyl-CoA pathway (modified from reference 37). PCT, propionate CoA transferase; POT, propionyl-CoA:oxaloacetate transcarboxylase; MCM, methylmalonyl-CoA mutase; SCS, succinyl-CoA synthetase; SDH, succinate dehydrogenase; FHT, fumarate hydratase (fumarase); MDH, malate dehydrogenase; PDH, pyruvate dehydrogenase.
|
|
|
|---|
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 1e100, 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 ml1 DNase I, 0.25 mg ml1 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.
|
|
|---|
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).
|
View this table: [in a new window] |
TABLE 1. Identification of ORFs relevant to the methylmalonyl-CoA pathway in P. thermopropionicum
|
![]() View larger version (35K): [in a new window] |
FIG. 2. Reconstructed MMC pathway in P. thermopropionicum based on the genomic data. See the explanations for abbreviations in the legend for Fig. 1 and footnote b of Table 1. When there is more than one ORF, the numbers of ORFs found in the genome are presented in parentheses. Fed(red) and Fed(ox), reduced and oxidized ferredoxin, respectively; FR, fumarate reductase; EA, electron acceptor.
|
|
View larger version (6K): [in a new window] |
FIG. 3. Genetic organization of the mmc cluster. For ORF names, refer to Table 1. The transposase gene (black arrow), regulatory gene (white arrow), and mmc structural genes (gray arrows) are indicated. T, terminator.
|
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.
![]() View larger version (36K): [in a new window] |
FIG. 4. 2D-GE patterns of soluble proteins in P. thermopropionicum cells recovered from coculture with M. thermautotrophicus. The growth substrates were propionate (A) and butanol (B). The circled protein spots were excised and identified (Table 2).
|
|
View this table: [in a new window] |
TABLE 2. Characterization of major protein spots on the 2D gels
|
|
|
|---|
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.
This work was supported by New Energy and Industrial Technology Development Organization (NEDO).
|
|
|---|
54 (
N) transcription factor. J. Bacteriol. 182:4129-4136.
30 (
H). J. Bacteriol. 170:1054-1062.
H: functional analysis and comparative genomics. J. Bacteriol. 179:7135-7155.
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»