Skip to main content

Genome analysis of Desulfotomaculum kuznetsovii strain 17T reveals a physiological similarity with Pelotomaculum thermopropionicum strain SIT


Desulfotomaculum kuznetsovii is a moderately thermophilic member of the polyphyletic spore-forming genus Desulfotomaculum in the family Peptococcaceae. This species is of interest because it originates from deep subsurface thermal mineral water at a depth of about 3,000 m. D. kuznetsovii is a rather versatile bacterium as it can grow with a large variety of organic substrates, including short-chain and long-chain fatty acids, which are degraded completely to carbon dioxide coupled to the reduction of sulfate. It can grow methylotrophically with methanol and sulfate and autotrophically with H2 + CO2 and sulfate. For growth it does not require any vitamins. Here, we describe the features of D. kuznetsovii together with the genome sequence and annotation. The chromosome has 3,601,386 bp organized in one contig. A total of 3,567 candidate protein-encoding genes and 58 RNA genes were identified. Genes of the acetyl-CoA pathway, possibly involved in heterotrophic growth with acetate and methanol, and in CO2 fixation during autotrophic growth are present. Genomic comparison revealed that D. kuznetsovii shows a high similarity with Pelotomaculum thermopropionicum. Genes involved in propionate metabolism of these two strains show a strong similarity. However, main differences are found in genes involved in the electron acceptor metabolism.


Desulfotomaculum kuznetsovii strain 17T (VKM B-1805; DSM 6115) is a moderately thermophilic sulfate-reducing bacterium isolated from deep subsurface thermal mineral water [1]. It grows with a wide range of substrates, including organic acids, such as long-chain fatty acids, short-chain fatty acids (butyrate, propionate, acetate), lactate, pyruvate, fumarate and succinate as well as ethanol and methanol. These substrates are degraded to CO2 coupled to sulfate reduction. The strain is also able to grow autotrophically with H2/CO2 and sulfate and to ferment pyruvate and fumarate. For growth, D. kuznetsovii has no vitamin requirement.

Desulfotomaculum is a genus of Gram-positive, spore-forming anaerobes that is phylogenetically and physiologically very diverse. The genus is poorly studied physiologically, while its members are known to play an important role in the carbon and sulfur cycle in a variety of often adverse environments. The genus is divided phylogenetically into different sub-groups [2,3]. To get a thorough understanding of the evolutionary relationship of the different Desulfotomaculum sub-groups and the physiology of the individual species, it is important to have genome sequence information. Here, we present a summary of the features of D. kuznetsovii strain 17T, together with the description of the complete genomic sequencing and annotation. Moreover, we describe a physiological and genomic comparison of D. kuznetsovii strain 17T and Pelotomaculum thermopropionicum strain SIT, because phylogenetically P. thermopropionicum is the closest related organism with validly published name that has a completely sequenced genome. However, the two strains have different physiological traits. For example, P. thermopropionicum is not able to grow by sulfate reduction, but is able to grow in syntrophy with methanogens. D. kuznetsovii lacks this ability. By comparing the genomes of the two bacteria we were able to identify the main similarities and differences.

Classification and features

D. kuznetsovii is a member of the phylum Firmicutes. Phylogenetic analysis of the 16S rRNA genes of D. kuznetsovii shows that it clusters in Desulfotomaculum cluster 1. This cluster not only contains Desulfotomaculum species, but also members of the genera Sporotomaculum, Cryptanaerobacter and Pelotomaculum. D. kuznetsovii is part of sub-group 1c together with D. solfataricum, D. luciae, D. thermosubterraneum, D. salinum, D. australicum, and D. thermocisternum, while Pelotomaculum species belong to sub-group 1h (Figure 1) [2].

Figure 1.
figure 1

Neighbor joining tree based on 16S rRNA sequences showing the phylogenetic affiliation of Desulfotomaculum and related species divided in the subgroups of Desulfotomaculum cluster 1. D. kuznetsovii is printed in bold type. The sequences of different Thermotogales were used as outgroup, but were pruned from the tree. Closed circles represent bootstrap values between 75 and 100%. The scale bar represents 10% sequences difference.

D. kuznetsovii cells are rod-shaped (1.0–1.4 x 3.5–5 εm) with rounded ends and peritrichous flagella [Figure 2]. Spores of D. kuznetsovii are spherical (1.3 εm in diameter) and centrally located causing swelling of the cells. D. kuznetsovii grows between 50 and 85°C, but the optimal growth temperature is 60–65°C. The substrates D. kuznetsovii can grow with are completely oxidized to CO2. Suitable electron acceptors are sulfate, thiosulfate and sulfite. D. kuznetsovii is also able to grow by fermentation of pyruvate and fumarate. A summary of the classification and general features of D. kuznetsovii is presented in Table 1 [1].

Figure 2.
figure 2

Scanning electron microscopic photograph of D. kuznetsovii.

Table 1. Classification and general features of D. kuznetsovii DSM 6115 according to the MIGS recommendations [4].

Genome sequencing and annotation

Genome project history

D. kuznetsovii was selected for sequencing in the DOE Joint Genome Institute Community Sequencing Program 2009, proposal 300132_795700 ‘Exploring the genetic and physiological diversity of Desulfotomaculum species’, because of its phylogenetic position in one of the Desulfotomaculum sub-groups, its important role in bioremediation, and its ability to use propionate, acetate and methanol for growth. The genome project is listed in the Genome OnLine Database (GOLD) [20] as project Gc01781, and the complete genome sequence was deposited in Genbank. Sequencing, finishing and annotation of the D. kuznetsovii genome were performed by the DOE Joint Genome Institute (JGI). A summary of the project information is shown in Table 2.

Table 2. Genome sequencing project information

Growth conditions and DNA isolation

D. kuznetsovii was grown anaerobically at 60oC in bicarbonate buffered medium with propionate and sulfate as substrates [1]. DNA of cell pellets was isolated using the standard DOE-JGI CTAB method recommended by the DOE Joint Genome Institute (JGI, Walnut Creek, CA, USA). In short, cells were resuspended in TE (10 mM tris; 1 mM EDTA, pH 8.0). Subsequently, cells were lysed using lysozyme and proteinase K, and DNA was extracted and purified using CTAB and phenol:chloroform:isoamylalcohol extractions. After precipitation in 2-propanol and washing in 70% ethanol, the DNA was resuspended in TE containing RNase. Following a quality and quantity check using agarose gel electrophoresis in the presence of ethidium bromide, and spectrophotometric measurement using a NanoDrop ND-1000 spectrophotometer (NanoDrop® Technologies, Wilmington, DE, USA).

Genome sequencing and assembly

The genome was sequenced using a combination of Illumina and 454 sequencing platforms. All general aspects of library construction and sequencing can be found at the JGI website [21]. Pyrosequencing reads were assembled using the Newbler assembler (Roche).

The initial Newbler assembly consisting of 81 contigs in five scaffolds was converted into a phrap [22] assembly by making fake reads from the consensus, to collect the read pairs in the 454 paired end library. Illumina GAii sequencing data (570.2 Mb) was assembled with Velvet [23] and the consensus sequences were shredded into 1.5 kb overlapped fake reads and assembled together with the 454 data. The 454 draft assembly was based on 134.6 Mb 454 draft data and all of the 454 paired end data. Newbler parameters are -consed -a 50 -l 350 -g -m -ml 20. The Phred/Phrap/Consed software package [22] was used for sequence assembly and quality assessment in the subsequent finishing process. After the shotgun stage, reads were assembled with parallel phrap (High Performance Software, LLC). Possible mis-assemblies were corrected with gapResolution [21], Dupfinisher [24], or sequencing cloned bridging PCR fragments with subcloning. Gaps between contigs were closed by editing in Consed, by PCR and by Bubble PCR primer walks (J.-F. Chang, unpublished). A total of 400 additional reactions and one shatter library were necessary to close gaps and to raise the quality of the finished sequence. Illumina reads were also used to correct potential base errors and increase consensus quality using a software Polisher developed at JGI [25]. The error rate of the completed genome sequence is less than 1 in 100,000. Together, the combination of the Illumina and 454 sequencing platforms provided 188.8 × coverage of the genome. The final assembly contained 323,815 pyrosequence and 15,594,144 Illumina reads.

Genome annotation

Genes were identified using Prodigal [26] as part of the Oak Ridge National Laboratory genome annotation pipeline, followed by a round of manual curation using the JGI GenePRIMP pipeline [27]. The predicted CDSs were translated and used to search the National Center for Biotechnology Information (NCBI) nonredundant database, UniProt, TIGR-Fam, Pfam, PRIAM, KEGG, COG, and InterPro databases. Additional gene prediction analysis and functional annotation was performed within the Integrated Microbial Genomes - Expert Review (IMG-ER) platform [28].

Genome properties and genome comparison with other strains

The genome of D. kuznetsovii consists of a circular chromosome of 3,601,386 bp with 54.88% GC content (Table 3 and Figure 3). Pseudogenes comprise 4.66% of the genes identified. Of the 3,625 genes predicted, 3,567 are protein-coding genes of which 2,560 are assigned to COG functional categories. The distribution of these genes into COG functional categories is presented in Table 4.

Figure 3.
figure 3

Graphical map of the chromosome of D. kuznetsovii. From outside to the center: Genes on the forward strand (colored by COG categories), Genes on the reverse strand (colored by COG categories), RNA genes (tRNAs green, rRNAs red, other RNAs black), GC content, GC skew.

Table 3. Genome statistics
Table 4. Number of genes associated with the general COG functional categories

The genome of D. kuznetsovii has 58 RNA genes of which, three are 16S rRNA genes. This is one more than the previously described rrnA and rrnB [29]. These two rRNA genes contained two large inserts. One at the variable 5′terminal region and one at the variable 3′terminal region. The main differences between the two rRNA genes were found in these inserts. These inserts were hypothesized to be involved in the operation of ribosomes at high temperatures. However, more research is needed to assess the function of these inserts. All three rRNA genes of D. kuznetsovii have a size of approximately 1,700 nucleotides. This suggests that the third rRNA gene might also contain inserts. Alignment of the 16S rRNA genes confirmed the presence of inserts in all three 16S rRNA genes (data not shown).

BLAST analysis [30,31] of the genes of D. kuznetsovii against genes in the KEGG Sequence Similarity DataBase revealed similarity with other Desulfotomaculum strains (Table 5), D. acetoxidans, D. carboxydivorans, “D. reducens” and D. ruminis, but interestingly also with non-Desulfotomaculum strains. D. kuznetsovii contains 873 genes with high similarity to genes of Pelotomaculum thermopropionicum, which is more than to any of the sequenced Desulfotomaculum species. Moreover, we identified the conserved proteins of D. kuznetsovii across three related fully sequenced species (Table 6). The bidirectional best blast hits showed that despite the smaller genome of P. thermopropionicum it contained more homologous predicted proteins with D. kuznetsovii (1,406) compared to D. acetoxidans (1,309) and “D. reducens” (1330). This suggests a strong physiological similarity between D. kuznetsovii and P. thermopropionicum.

Table 5. Taxonomic distribution of the top KEGG hits of D. kuznetsovii genes based on BLAST against KEGG database.
Table 6. Proteins of D. kuznetsovii conserved across three related species with fully sequenced genomes.

Insights into the genome

Involvement of the acetyl-coA pathway in growth with acetate and methanol

D. kuznetsovii oxidizes acetate completely to CO2. The pathway of acetate degradation has not been studied yet, but sulfate reducers may employ the tricarboxylic acid (TCA) cycle or the acetyl-CoA pathway for acetate degradation, as exemplified by Desulfobacter postgatei and Desulfobacca acetoxidans, respectively [32]. Most genes predicted to code for enzymes of the TCA cycle are present in the genome of D. kuznetsovii, but genes with similarity to those coding for an ATP-dependent citrate synthase and isocitrate dehydrogenase are missing. This suggests that the TCA cycle is not complete and that the TCA cycle enzymes have mainly an anabolic function or a function in other catabolic pathways, such as the propionate degradation pathway. Genes with similarity to those coding for enzymes involved in the acetyl-CoA pathway are all present in the genome of D. kuznetsovii (Figure 4), which suggests its involvement in acetate oxidation. However, there are no genes similar to those that code for acetate kinase and phosphate acetyltransferase present in the genome. The reaction from acetate to acetyl-CoA is likely performed by acetyl-CoA synthetase (Desku_1241).

Figure 4.
figure 4

Pathway of acetate oxidation to CO2 by D. kuznetsovii. Enzymes in this figure are in bold italic and their locus tags are included. Genes with the locus tags Desku_1488 and Desku_1490 putatively code for the small subunit and the large subunit of the iron-sulfur protein, respectively. This protein is involved in transferring the methyl from acetyl-CoA to tetrahydrofolate. Abbreviations: A-CoA S, acetyl-CoA synthetase; AcsA, carbon-monoxide dehydrogenase; AcsB, acetyl-CoA synthase; CFeSP, iron-sulfur protein; CH3, methyl; THF, tetrahydrofolate; MeTr, methyltransferase.

D. acetoxidans is an acetate-oxidizing Desulfotomaculum species, positioned in sub-group 1e (Figure 1), that also uses the acetyl-CoA pathway for acetate oxidation to CO2 [33]. The genes involved in acetate oxidation in D. acetoxidans are similar to those in D. kuznetsovii, but there are some exceptions. The genome of D. acetoxidans does not contain a gene that putatively codes for acetyl-CoA synthetase, similar to D. kuznetsovii, but contains genes that putatively code for an acetate kinase and a phosphate acetyltransferase [34]. Additionally, putative carbon-monoxide dehydrogenase complex coding genes involved in the acetyl-CoA pathway show differences between the two Desulfotomaculum species. D. kuznetsovii lacks a ferredoxin coding gene that is located between cooC (Desku_1493) and acsE (Desku_1487), which in contrast is present in the genome of D. acetoxidans (Dtox_1273). Moreover, three genes similar to heterodisulfide reductase encoding genes (Desku_1486-1484) are located upstream of acsE in D. kuznetsovii, which is not the case in the genome of D. acetoxidans.

Methanol metabolism

Growth of D. kuznetsovii with methanol and sulfate was studied [35]. In that study the activity of methyltransferase, an enzyme that is involved in methanol metabolism in methanogens and acetogens [36,37], could not be assessed, while low activities of an alcohol dehydrogenase could be measured. An alcohol dehydrogenase with a molecular mass of 42 kDa was partially purified and showed activity with methanol [35]. The genome of D. kuznetsovii contains several alcohol dehydrogenase genes (Desku_0165, 0619, 0624, 0628, 2955, 3082) that each code for an enzyme with a size of approximately 42 kDa. In the genome, genes with similarity to those coding for a methanol methyltransferase mtaA (Desku_0050, 0055, 0060), mtaB (Desku_0051) and mtaC (Desku_0048, 0049, 0052, 0056) were also found, suggesting a methanol metabolism as described in Moorella thermoacetica [36]. Further studies are needed to obtain information about the diversity of the methanol-degradation pathways in D. kuznetsovii.

Comparison of D. kuznetsovii and P. thermopropionicum genomes

Genomic comparison revealed that a large number of D. kuznetsovii genes show similarity to genes of Pelotomaculum thermopropionicum, a syntrophic propionate-oxidizing thermophile (Table 5 and 6). Interestingly, among them are genes that putatively code for enzymes involved in propionate metabolism (Table 7). Moreover, the genetic organization of the methylmalonyl-CoA (mmc) cluster in the genome of both bacteria is similar (Figure 5). However, D. kuznetsovii lacks tps, mmcA and mmcM in the mmc cluster. mmcA codes for a response regulator and mmcM for pyruvate ferredoxin oxidoreductase.

Figure 5.
figure 5

Gene organization of the mmc cluster in D. kuznetsovii and P. thermopropionicum. Names of the genes can be found in table 6, except for tps, which is a transposase gene.

Table 7. Genes in D. kuznetsovii that are annotated as enzymes involved in propionate metabolism.

Based on 16S rRNA gene sequences, D. kuznetsovii and P. thermopropionicum group in cluster group c and h of the Desulfotomaculum cluster 1, respectively (see Figure 1). P. thermopropionicum is known for its ability to grow with propionate and ethanol in syntrophic association with methanogens. It is not able to grow by sulfate respiration, despite the presence of sulfate reduction genes in the genome [38]. In contrast, D. kuznetsovii is able to grow with propionate (Figure 6) and ethanol with sulfate. However, in the absence of sulfate, it cannot grow in syntrophic association with methanogens. Therefore, differences are expected in genes coding for hydrogenases, formate dehydrogenases, and those involved in sulfate reduction.

Figure 6.
figure 6

Propionate degradation pathway in D. kuznetsovii based on genomic data. Enzymes are depicted in bold italic. Next to these enzymes are the possible encoding genes, and their locus tags. GCT, Glutaconate CoA-transferase; MCD, Methylmalonyl-CoA decarboxylase; PCC, Propionyl-CoA carboxylase; MCE, Methylmalonyl-CoA epimerase; MCM, Methylmalonyl-CoA mutase; SCS, Succinyl-CoA synthetase; SDH, Succinate dehydrogenase; FHT, Fumarase; MDH, Malate dehydrogenase; OAD, Oxaloacetate decarboxylase; PFL, Pyruvate formate lyase.

Sulfate reduction genes: Figure 7 depicts the sulfate reduction pathway of the two strains. In the genome of D. kuznetsovii two genes (Desku_2103; Desku_3527) are annotated as phosphoadenosine phosphosulfate reductase encoding genes whose corresponding proteins might be involved in assimilatory sulfate metabolism. The P. thermopropionicum genome lacks these genes [39]. Instead, the P. thermopropionicum genome contains an adenylylsulfate kinase gene (PTH_0238). In the dissimilatory sulfate reduction pathway, the two strains both have genes that code for enzymes to reduce sulfate to H2S. However, P. thermopropionicum is missing the gene that codes for an adenylylsulfate reductase beta subunit, which is present in the D. kuznetsovii genome (Desku_1073). Moreover, the gene labeled as a dissimilatory sulfite reductase (dsr) alpha and beta subunit in the P. thermopropionicum genome (PTH_0242) is not similar to dsrA or dsrB from D. kuznetsovii or any other Desulfotomaculum strain. However, it has high similarity to the dsrC gene from D. kuznetsovii, indicating that it is not a dsrA or dsrB gene but a dsrC gene (data not shown). Therefore, the inability of P. thermopropionicum to grow by sulfate respiration is most likely caused by the absence of an adenylylsulfate reductase beta subunit encoding gene and the dsrAB genes.

Figure 7.
figure 7

Sulfate reduction pathway of D. kuznetsovii and P. thermopropionicum. Depicted in green are genes that code for sulfate reduction enzymes that are present in the genome. Dashed arrows indicate the presence of a subunit encoding gene, but not the presence of all genes required for the enzyme. Dashed dotted arrows are used when no genes were found for the reaction. Abbreviations: APS, adenylylsulfate; DP, diphosphate; PAPS, 3′-Phosphoadenylyl-sulfate (PAPS); redA, reduced acceptor; oxA, oxidized acceptor.

Hydrogenase and formate dehydrogenase genes: Schut and Adams (2009) [40] showed that the trimeric [FeFe]-hydrogenase from Thermotoga maritima oxidizes NADH and ferredoxin simultaneously to produce H2. Similar bifurcating / confurcating [FeFe]-hydrogenases and formate dehydrogenases are present in Syntrophobacter fumaroxidans and P. thermopropionicum [41]. Both generate NADH and ferredoxin during propionate degradation via the methylmalonyl-CoA pathway and might use confurcating hydrogenases and formate dehydrogenases to drive the unfavorable re-oxidation of NADH (E0’=-320mV) by the exergonic re-oxidation of ferredoxin (E0’=-398mV) to produce hydrogen (E0’= -414mV) or formate (E0’= -432mV) that are subsequently transferred to hydrogen and formate scavenging methanogens. Additionally, up-regulation of genes encoding hydrogenases and formate dehydrogenases in P. thermopropionicum was shown during syntrophic growth [42]. The P. thermopropionicum genome contains three [FeFe]-hydrogenases, one [NiFe]-hydrogenase and two formate dehydrogenases. One [FeFe]-hydrogenase (PTH_0668-0670) was shown to be down-regulated during syntrophic growth, while the other two [FeFe]-hydrogenases (PTH_1377-1379 and PTH_2010-2012) were up-regulated. The two formate dehydrogenases of P. thermopropionicum (I, PTH_1711-1714 and II, PTH_2645-2649) were both up-regulated during syntrophic growth [42]. According to TMHMM server v. 2.0 [43] formate dehydrogenase I of P. thermopropionicum has transmembrane helices. Therefore, it might play an essential role in the interspecies transfer of reducing equivalents in syntrophic growth.

The genome of D. kuznetsovii was screened for hydrogenase and formate dehydrogenase encoding gene clusters with BLAST analysis. Pfam search [44] was used to identify motifs in the amino acid sequences and the TMHMM Server v. 2.0 [43] was used to screen for transmembrane helices. The TatP 1.0 Server was used to screen for twin-arginine translocation (Tat) motifs in the N-terminus to predict protein localization in the cell [45]. The incorporation of selenocysteine (SeCys) was examined by RNA loop predictions with Mfold version 3.2 [46, 47]. The predicted RNA loop in the 50–100 bp region downstream of the UGA-codon was compared with the consensus loop described earlier [48].

Compared to P. thermopropionicum, D. kuznetsovii lacks membrane associated formate dehydrogenases and hydrogenases and also lacks [NiFe]-hydrogenase. This might explain why D. kuznetsovii cannot grow in syntrophic relation with methanogens. The genome of D. kuznetsovii indicates the presence of a confurcating selenocysteine-incorporated formate dehydrogenase (Desku_2987-2991), two trimeric confurcating [FeFe]-hydrogenases (Desku_2307-2309, Desku_2995-2997) and two [FeFe]-hydrogenases (Desku_0995, Desku_2934-2935) without NADH-binding sites (Figure 8). Several subunits of these enzymes are related to subunits of NADH dehydrogenase (complex I), including the NADH-binding proteins related to NuoF (Desku_2990, 2308 and 2996) and the electron transfer subunits related to NuoE (Desku_2991, 2935, and 2997) and to NuoG (Desku_2989). In three of the [FeFe]-hydrogenases this NuoG-like domain is fused with the catalytic subunit (Desku_2995, 2307 and 2934). Two of the multimeric hydrogenases are found next to [FeFe]-hydrogenases containing PAS-sensor domains (Desku_2932 and Desku_2994), suggesting they are involved in the regulation of the synthesis of those hydrogenases. All complexes are predicted to be cytoplasmic and not membrane bound.

Figure 8.
figure 8

Schematic representation of a putative confurcating formate dehydrogenase, two putative confurcating [FeFe]-hydrogenases and two ferredoxin re-oxidizing [FeFe]-hydrogenases in Desulfotomaculum kuznetsovii. Gene locus tag numbers and α-, β-, and γ-subunits are depicted. Moreover, predicted iron-sulfur clusters and metal-binding sites are indicated.

Apart from a possible involvement in the acetate oxidation pathway (Figure 4), it remains unclear for which purpose D. kuznetsovii uses its confurcating formate dehydrogenase and hydrogenases because our genome analysis indicates that pyruvate oxidation during propionate degradation generates formate instead of ferredoxin (Figure 6).

Vitamin synthesis

D. kuznetsovii is able to grow in medium without vitamins [1]. This indicates that D. kuznetsovii is able to synthesize all the vitamins that are required for its metabolism and that vitamin synthesis genes should be present in the genome. Vitamin B12 is essential for the methylmalonyl-CoA pathway and the acetyl-coA pathway. The biosynthesis of cobalamin (vitamin B12) is known to occur from uroporphyrinogen-III to adenosylcobalamin via two possible pathways, the aerobic and anaerobic pathway of the corrinoid ring [49,50]. The D. kuznetsovii genome contains all genes needed for the anaerobic pathway: cysGA (Desku_1520), cysGB (Desku_1460, Desku_1523), cbiA (Desku_1765, Desku_2368), cbiBCDEFGHJLPT (Desku_2369, 1459, 1468, 1467, 1464, 1463, 1462, 1461, 1465, 2370 and 1466, respectively), cobalt reductase (Desku_2757), btuR (Desku_0004, 1209), cobS (Desku_2367) and cobU (Desku_2371). Moreover, D. kuznetsovii has genes to convert glutamyl tRNA to uroporphyrinogen-III, hemABCDL (Desku_1522, 1518, 1521, 1520 and 1522, respectively). The genome also contains some unassigned cobalamin synthesis genes (P47K, Desku_0046, 0053; cbiM, Desku_2905), corrinoid transport proteins (Desku_0693, 702, 2237–2239, 2902–2904, 3025–3027) and, interestingly, two cobN genes (Desku_2189, 2227), genes involved in the aerobic pathway. It is unclear why D. kuznetsovii has these cobN genes, since all anaerobic pathway genes are present in the genome, and it is unclear if the products of these two genes are used for cobalamin synthesis by D. kuznetsovii.

Other vitamin synthesis genes present in the genome of D. kuznetsovii are genes involved in biotin synthesis (vitamin H) (Desku_1295-1297, 2246–2247, 2317), nicotinamide (vitamin B3) synthesis (Desku_0433, 0614, 0662, 0815, 1248, 1417, 1472, 1499, 1925, 1951, 3103, 3121, 3227, 3228, 3231, 3246, 3337), thiamin (vitamin B1) synthesis (0372, 0543, 0545, 2253, 2363, 2639), riboflavin (vitamin B2) synthesis (Desku_1244-1247), and pantothenate (vitamin B5) synthesis (Desku_3262). The genes involved in coenzyme A production from pantothenate are also present in the D. kuznetsovii genome (Desku_1254, 1307, 3145, 3200). Moreover, genes involved in the biosynthesis of pyridoxine (vitamin B6) via the deoxyxylulose 5-phosphate (DXP) independent route were found to be in the genome (Desku_0007, 0008). These genes code for two enzymes that facilitate the conversion of glutamine to the active form of vitamin B6, pyridoxal 5′-phosphate [51].

Menaquinone (vitamin K) and ubiquinone (coenzyme Q10) biosynthesis is important because of the electron transport function in the membranes. The genes that code for the biosynthesis enzymes from polyprenyldiphosphate to menaquinone and ubiquinone are present in the D. kuznetsovii genome (Desku_0124, 0126, 0629, 1551–1554, 1829, 2629 and 3525), except for the genes that code for a 2-polyprenyl-6-methoxyphenol 4-monooxygenase (UbiH) and 2-polyprenyl-3-methyl-6-methoxy-1,4-benzoquinone hydroxylase (UbiF). Additionally, three genes (Desku_1548-1550) could be identified as putative menaquinone biosynthesis genes and are part of a menaquinone biosynthesis gene cluster (Desku_1548-1554). The products of those three genes could be involved in the reactions of the missing UbiH and UbiF encoding genes.

Folate (vitamin B9) biosynthesis is also of great importance for D. kuznetsovii, because it is an essential part of the acetyl-CoA pathway. It is involved in the transfer of one-carbon compounds and can be biosynthesized from chorismate and guanosine triphosphate (GTP) [5255]. Both pathways use a dihydropteroate synthase to produce dihydropteroate. The genome of D. kuznetsovii contains the genes encoding the enzymes involved in the pathway from chorismate to dihydropteroate (Desku_0219, 2268–2269) and from GTP to dihydropteroate (Desku_0210, 0219–0221 and 1419). The gene encoding a phosphatase (Desku_0210) in the D. kuznetsovii genome is probably involved in the removal of phosphate groups from dihydropterine triphosphate as a substitute for an alkaline phosphatase encoding gene, which is not present in the genome. Additionally, the genome contains a bifunctional protein encoding gene (Desku_404) that is expected to be responsible for the production of dihydrofolate (DHF) and the addition of multiple glutamate moieties to DHF or tetrahydrofolate (THF). However, the D. kuznetsovii genome lacks the DHF reductase encoding gene, which is required to reduce DHF to THF. The DHF reductase encoding gene appears to be absent in many microorganisms [56]. Levin et al. (2004) propose that in Halobacterium salinarum a dihydrofolate synthase and a dihydropteroate synthase domain is able to replace the function of the DHF reductase. Additionally, the authors show that when using a BLAST search, homologs of polypeptides can be found in organisms that lack a DHF reductase [56]. However, BLAST results showed no homologous protein encoding gene in the genome of D. kuznetsovii (data not shown). How in D. kuznetsovii DHF is reduced to THF can currently not be deduced from the genome sequence.


  1. Nazina TN, Ivanova AE, Kanchaveli LP, Rozanova EP. A new sporeforming thermophilic methylotrophic sulfate-reducing bacterium, Desulfotomaculum kuznetsovii sp. nov. Mikrobiologiya 1988; 57:823–827.

    CAS  Google Scholar 

  2. Stackebrandt E, Sproer C, Rainey FA, Burghardt J, Pauker O, Hippe H. Phylogenetic analysis of the genus Desulfotomaculum: evidence for the mis-classification of Desulfotomaculum guttoideum and description of Desulfotomaculum orientis as Desulfosporosinus orientis gen. nov., comb. nov. Int J Syst Bacteriol 1997; 47:1134–1139. PubMed

    Article  CAS  PubMed  Google Scholar 

  3. Plugge CM, Zhang W, Scholten JC, Stams AJM. Metabolic flexibility of sulfate-reducing bacteria. Front Microbiol 2011; 2:81–87. PubMed

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  4. Field D, Garrity G, Gray T, Morrison N, Selengut J, Sterk P, Tatusova T, Thomson N, Allen MJ, Angiuoli SV, et al. The minimum information about a genome sequence (MIGS) specification. Nat Biotechnol 2008; 26:541–547. PubMed

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  5. Woese CR, Kandler O, Wheelis ML. Towards a natural system of organisms: proposal for the domains Archaea, Bacteria, and Eucarya. Proc Natl Acad Sci USA 1990; 87:4576–4579. PubMed

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  6. Gibbons NE, Murray RGE. Proposals Concerning the Higher Taxa of Bacteria. Int J Syst Bacteriol 1978; 28:1–6.

    Article  Google Scholar 

  7. Garrity GM, Holt JG. The Road Map to the Manual. In: Garrity GM, Boone DR, Castenholz RW (eds), Bergey’s Manual of Systematic Bacteriology, Second Edition, Volume 1, Springer, New York, 2001, p. 119–169.

    Chapter  Google Scholar 

  8. Murray RGE. The Higher Taxa, or, a Place for Everything…? In: Holt JG (ed), Bergey’s Manual of Systematic Bacteriology, First Edition, Volume 1, The Williams and Wilkins Co., Baltimore, 1984, p. 31–34.

    Google Scholar 

  9. List Editor. List of new names and new combinations previously effectively, but not validly, published. List no. 132. Int J Syst Evol Microbiol 2010; 60:469–472.

  10. Rainey FA. Class II. Clostridia class nov. In: De Vos P, Garrity G, Jones D, Krieg NR, Ludwig W, Rainey FA, Schleifer KH, Whitman WB (eds), Bergey’s Manual of Systematic Bacteriology, Second Edition, Volume 3, Springer-Verlag, New York, 2009, p. 736.

    Google Scholar 

  11. Skerman VBD, McGowan V, Sneath PHA. Approved Lists of Bacterial Names. Int J Syst Bacteriol 1980; 30:225–420.

    Article  Google Scholar 

  12. Prévot AR. In: Hauderoy P, Ehringer G, Guillot G, Magrou. J., Prévot AR, Rosset D, Urbain A (eds), Dictionnaire des Bactéries Pathogènes, Second Edition, Masson et Cie, Paris, 1953, p. 1–692.

    Google Scholar 

  13. Rogosa M. Peptococcaceae, a new family to include the Gram-positive, anaerobic cocci of the genera Peptococcus, Peptostreptococcus and Ruminococcus. Int J Syst Bacteriol 1971; 21:234–237.

    Article  Google Scholar 

  14. Campbell LL, Postgate JR. Classification of the spore-forming sulfate-reducing bacteria. Bacteriol Rev 1965; 29:359–363. PubMed

    PubMed Central  CAS  PubMed  Google Scholar 

  15. Campbell LL. Genus IV. Desulfotomaculum Campbell and Postgate 1965, 361. In: Buchanan RE, Gibbons NE (eds), Bergey’s Manual of Determinative Bacteriology, Eighth Edition, The Williams and Wilkins Co., Baltimore, 1974, p. 572–573.

    Google Scholar 

  16. Validation of the publication of new names and new combinations previously effectively published outside the IJSB. List no. 35. Int J Syst Bacteriol 1990; 40:470–471.

  17. Isaksen MF, Bak F, Jorgensen BB. Thermophilic sulfate-reducing bacteria in cold marine sediment. Microb Ecol 1994; 14:1–8.

    Article  CAS  Google Scholar 

  18. Goorissen HP, Boschker HT, Stams AJM, Hansen TA. Isolation of thermophilic Desulfotomaculum strains with methanol and sulfite from solfataric mud pools, and characterization of Desulfotomaculum solfataricum sp. nov. Int J Syst Evol Microbiol 2003; 53:1223–1229. PubMed

    Article  CAS  PubMed  Google Scholar 

  19. Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, et al. Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet 2000; 25:25–29. PubMed

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  20. Pagani I, Liolios K, Jansson J, Chen IM, Smirnova T, Nosrat B, Markowitz VM, Kyrpides NC. The Genomes OnLine Database (GOLD) v.4: status of genomic and metagenomic projects and their associated metadata. Nucleic Acids Res 2012; 40(Database issue):D571–D579. PubMed

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  21. JGI website.

  22. The Phred/Phrap/Consed software package.

  23. Zerbino DR, Birney E. Velvet: algorithms for de novo short read assembly using de Bruijn graphs. Genome Res 2008; 18:821–829. PubMed

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  24. Han C, Chain P. Finishing repeat regions automatically with Dupfinisher. In: H.R. A, H. V, editors2006 June 26–29, 2006. CSREA Press. p 141–6.

  25. Lapidus A, LaButti K, Foster B, Lowry S, Trong S, Goltsman E. POLISHER: An effective tool for using ultra short reads in microbial genome assembly and finishing. Marco Island, FL: AGBT; 2008.

    Google Scholar 

  26. Hyatt D, Chen GL, Locascio PF, Land ML, Larimer FW, Hauser LJ. Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinformatics 2010; 11:119. PubMed

    Article  PubMed Central  PubMed  Google Scholar 

  27. Pati A, Ivanova NN, Mikhailova N, Ovchinnikova G, Hooper SD, Lykidis A, Kyrpides NC. GenePRIMP: a gene prediction improvement pipeline for prokaryotic genomes. Nat Methods 2010; 7:455–457. PubMed

    Article  CAS  PubMed  Google Scholar 

  28. Markowitz VM, Mavromatis K, Ivanova NN, Chen IM, Chu K, Kyrpides NC. IMG ER: a system for microbial genome annotation expert review and curation. Bioinformatics 2009; 25:2271–2278. PubMed

    Article  CAS  PubMed  Google Scholar 

  29. Tourova TP, Kuznetsov BB, Novikova EV, Poltaraus AB, Nazina TN. Heterogeneity of the nucleotide sequences of the 16 S rRNA genes of the type strain of Desulfotomaculum kuznetsovii. Mikrobiologiya 2001; 70:788–795.

    Google Scholar 

  30. Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res 1997; 25:3389–3402. PubMed

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  31. Altschul SF, Wootton JC, Gertz EM, Agarwala R, Morgulis A, Schaffer AA, Yu YK. Protein database searches using compositionally adjusted substitution matrices. FEBS J 2005; 272:5101–5109. PubMed

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  32. Goevert D, Conrad R. Carbon isotope fractionation by sulfate-reducing bacteria using different pathways for the oxidation of acetate. Environ Sci Technol 2008; 42:7813–7817. PubMed

    Article  CAS  PubMed  Google Scholar 

  33. Spormann AM, Thauer RK. Anaerobic acetate oxidation to CO2 in Desulfotomaculum acetoxidans. Arch Microbiol 1988; 150:374–380.

    Article  CAS  Google Scholar 

  34. Spring S, Lapidus A, Schroder M, Gleim D, Sims D, Meincke L, Glavina Del Rio T, Tice H, Copeland A, Cheng JF, et al. Complete genome sequence of Desulfotomaculum acetoxidans type strain (5575). Stand Genomic Sci 2009; 1:242–253. PubMed

    Article  PubMed Central  PubMed  Google Scholar 

  35. Goorissen HP, Stams AJM, Hansen TA. Methanol dissimilation in Desulfotomaculum kuznetsovii PhD dissertation: Thermophilic methanol utilization by sulfate reducing bacteria 2002;Chapter 3:55–61.

  36. Pierce E, Xie G, Barabote RD, Saunders E, Han CS, Detter JC, Richardson P, Brettin TS, Das A, Ljungdahl LG, et al. The complete genome sequence of Moorella thermoacetica (f. Clostridium thermoaceticum). Environ Microbiol 2008; 10:2550–2573. PubMed

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  37. van der Meijden P, Heythuysen HJ, Pouwels A, Houwen F, van der Drift C, Vogels GD. Methyltransferases involved in methanol conversion by Methanosarcina barkeri. Arch Microbiol 1983; 134:238–242. PubMed

    Article  PubMed  Google Scholar 

  38. Imachi H, Sekiguchi Y, Kamagata Y, Hanada S, Ohashi A, Harada H. Pelotomaculum thermopropionicum gen. nov., sp. nov., an anaerobic, thermophilic, syntrophic propionate-oxidizing bacterium. Int J Syst Evol Microbiol 2002; 52:1729–1735. PubMed

    CAS  PubMed  Google Scholar 

  39. Kosaka T, Kato S, Shimoyama T, Ishii S, Abe T, Watanabe K. The genome of Pelotomaculum thermopropionicum reveals niche-associated evolution in anaerobic microbiota. Genome Res 2008; 18:442–448. PubMed

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  40. Schut GJ, Adams MW. The iron-hydrogenase of Thermotoga maritima utilizes ferredoxin and NADH synergistically: a new perspective on anaerobic hydrogen production. J Bacteriol 2009; 191:4451–4457. PubMed

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  41. Müller N, Worm P, Schink B, Stams AJM, Plugge CM. Syntrophic butyrate and propionate oxidation processes: from genomes to reaction mechanisms. Environmental Microbiology Reports 2010; 2:489–499.

    Article  PubMed  Google Scholar 

  42. Kato S, Kosaka T, Watanabe K. Substrate-dependent transcriptomic shifts in Pelotomaculum thermopropionicum grown in syntrophic co-culture with Methanothermobacter thermautotrophicus. Microb Biotechnol 2009; 2:575–584. PubMed

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  43. DTU. 2009 Center for biological sequence analysis: Technical University of Denmark.

  44. Sanger Institute CU. 2009.

  45. Bendtsen JD, Nielsen H, Widdick D, Palmer T, Brunak S. Prediction of twin-arginine signal peptides. BMC Bioinformatics 2005; 6:167. PubMed

    Article  PubMed Central  PubMed  Google Scholar 

  46. Mathews DH, Sabina J, Zuker M, Turner DH. Expanded sequence dependence of thermodynamic parameters improves prediction of RNA secondary structure. J Mol Biol 1999; 288:911–940. PubMed

    Article  CAS  PubMed  Google Scholar 

  47. Zuker M. Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Res 2003; 31:3406–3415. PubMed

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  48. Zhang Y, Gladyshev VN. An algorithm for identification of bacterial selenocysteine insertion sequence elements and selenoprotein genes. Bioinformatics 2005; 21:2580–2589. PubMed

    Article  CAS  PubMed  Google Scholar 

  49. Roth JR, Lawrence JG, Bobik TA. Cobalamin (co-enzyme B12): synthesis and biological significance. Annu Rev Microbiol 1996; 50:137–181. PubMed

    Article  CAS  PubMed  Google Scholar 

  50. Roessner CA, Santander PJ, Scott AI. Multiple biosynthetic pathways for vitamin B12: variations on a central theme. Vitam Horm 2001; 61:267–297. PubMed

    Article  CAS  PubMed  Google Scholar 

  51. Fitzpatrick TB, Amrhein N, Kappes B, Macheroux P, Tews I, Raschle T. Two independent routes of de novo vitamin B6 biosynthesis: not that different after all. Biochem J 2007; 407:1–13. PubMed

    Article  CAS  PubMed  Google Scholar 

  52. Bermingham A, Derrick JP. The folic acid biosynthesis pathway in bacteria: evaluation of potential for antibacterial drug discovery. Bioessays 2002; 24:637–648. PubMed

    Article  CAS  PubMed  Google Scholar 

  53. Dosselaere F, Vanderleyden J. A metabolic node in action: chorismate-utilizing enzymes in micro-organisms. Crit Rev Microbiol 2001; 27:75–131. PubMed

    Article  CAS  PubMed  Google Scholar 

  54. Hu SI, Drake HL, Wood HG. Synthesis of acetyl coenzyme A from carbon monoxide, methyltetrahydrofolate, and coenzyme A by enzymes from Clostridium thermoaceticum. J Bacteriol 1982; 149:440–448. PubMed

    PubMed Central  CAS  PubMed  Google Scholar 

  55. Stupperich E, Konle R. Corrinoid-Dependent Methyl Transfer Reactions Are Involved in Methanol and 3,4-Dimethoxybenzoate Metabolism by Sporomusa ovata. Appl Environ Microbiol 1993; 59:3110–3116. PubMed

    PubMed Central  CAS  PubMed  Google Scholar 

  56. Levin I, Giladi M, Altman-Price N, Ortenberg R, Mevarech M. An alternative pathway for reduced folate biosynthesis in bacteria and halophilic archaea. Mol Microbiol 2004; 54:1307–1318. PubMed

    Article  CAS  PubMed  Google Scholar 

Download references


The work conducted by the U.S. Department of Energy Joint Genome Institute is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231, and was also supported by grants CW-TOP 700.55.343 and ALW 819.02.014 of the Netherlands Science Foundation (NWO) and grant 323009 of the European Research Council.

Author information

Authors and Affiliations


Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Cite this article

Visser, M., Worm, P., Muyzer, G. et al. Genome analysis of Desulfotomaculum kuznetsovii strain 17T reveals a physiological similarity with Pelotomaculum thermopropionicum strain SIT. Stand in Genomic Sci 8, 69–87 (2013).

Download citation

  • Published:

  • Issue Date:

  • DOI: