- Short genome report
- Open access
- Published:
Working draft genome sequence of the mesophilic acetate oxidizing bacterium Syntrophaceticus schinkii strain Sp3
Standards in Genomic Sciences volume 10, Article number: 99 (2015)
Abstract
Syntrophaceticus schinkii strain Sp3 is a mesophilic syntrophic acetate oxidizing bacterium, belonging to the Clostridia class within the phylum Firmicutes, originally isolated from a mesophilic methanogenic digester. It has been shown to oxidize acetate in co-cultivation with hydrogenotrophic methanogens forming methane. The draft genome shows a total size of 3,196,921 bp, encoding 3,688 open reading frames, which includes 3,445 predicted protein-encoding genes and 55 RNA genes. Here, we are presenting assembly and annotation features as well as basic genomic properties of the type strain Sp3.
Introduction
During anaerobic degradation of organic material, acetate is formed as a main fermentation product, which is further converted to methane. Two mechanisms for methane formation from acetate have been described: The first one is carried out by aceticlastic methanogens converting acetate to methane and CO2 under low ammonia conditions [1]. The second mechanism, dominating under high ammonia conditions, occurs in two steps, and is performed by acetate-oxidizing bacteria oxidizing acetate to H2 (formate) and CO2 and a methanogenic partner using the hydrogen (formate) to reduce CO2 to methane [2–4]. Most fascinating on this syntrophic relationship is, that the overall reaction operates with a ΔG°´ of -36 kJ x mol−1 close to the thermodynamic equilibrium.
The number of isolated and characterized SAOB is restricted most likely due to their considerable differences in substrate utilization abilities and cultivation requirements. To date three mesophilic SAOB, namely Clostridium ultunense [5], Syntrophaceticus schinkii [6], “ Tepidanaerobacter acetatoxydans ” [7] and two thermophilic SAOB, namely Thermacetogenium phaeum [2] and Thermotoga lettingae [8] currently renamed to Pseudothermotoga lettingae have been isolated and characterized. Among those, two complete genome sequences of T. phaeum [9], “T. acetatoxydans” [10] and one draft genome sequence of C. ultunense [11] have been published, the later two by this working group. Here, we are presenting the draft genome sequence of the third mesophilic SAOB Syntrophaceticus schinkii strain Sp3. To date, strain Sp3 is the only isolated and characterized representative of the species S. schinkii and was recovered from an up flow anaerobic filter treating wastewater from a fishmeal factory [6]. This process was characterized by high ammonium concentration (6.4 g l−1 NH4 +). S. schinkii shows the least narrow substrate spectrum compared to all known SAOB, when growing heterotrophically [6]. The main end product formed is acetate, what allocates the species to the physiological group of acetogens.
Since the recovery of S. schinkii we found it at high abundance in all mesophilic large scale and lab scale biogas producing process we have investigated so far. Genome analysis and comparative genomics might help us to understand general features of syntrophy in particular energy conservation and electron transfer mechanisms during syntrophic acetate oxidation.
The present study summarizes genome sequencing, assembly and annotation as well as general genomic properties of the Syntrophaceticus schinkii strain Sp3 genome.
Organism information
Classification and features
Syntrophaceticus schinkii Sp3 (Fig. 1) is an obligate anaerobic, endospores forming bacterium, whose cells were found to be Gram variable with changing shapes dependent on the growth condition (Table 1, [6]). No flagella have been observed under any condition tested. It can grow up to 0.6 M NH4Cl in pure culture between 25 °C and 40 °C. A more detailed physiological description can be found in Westerholm et al. [6]. Minimum Information about the Genome Sequence (MIGS) of S. schinkii strain Sp3 is provided in Table 1 and Table S1 (Additional file 1).
Phylogentic analysis of the single 16S rRNA gene copy affiliates S. schinkii strain Sp3 to the Clostridia class within the phylum Firmicutes . The RDP Classifier ([12] 2015-08-05) confirmed further the affiliation to Thermoanaerobacteraceae as published by [6] in 2011 (Table 1). The comparison of the 16S rRNA gene sequence with the latest available databases from GenBank (2015-08-05) using NCBI BLAST [13] under default settings identified the thermophilic SAOB T. phaeum ( NR_074723.1) as the closest characterized relative sharing 92.12 % identity (Fig. 2). S. schinkii is only distantly related to the characterized mesophilic SAOB C. ultunense ( 82.54 % identity), and “T. acetatoxydans” (84.1 % identity) and the thermophilic P. lettingae (79.64 %). Although S. schinkii has been physiologically affiliated to the group of acetogens, Fig. 2 illustrates a distant relationship to this group, as represented by e.g. the model acetogen Moorella thermoacetica (89.15 % identity).
Genome sequencing information
Genome project history
Syntrophaceticus schinkii strain Sp3 was sequenced and annotated by the SLU-Global Bioinformatics Centre at the Swedish University of Agricultural Sciences, Uppsala, Sweden. The genome project is deposited in the Genomes OnLine Database [14] with GOLD id Gi0035837 and the working draft genome is deposited in the European Nucleotide Archive database with accession number ERP005192. The SAOB was selected for sequencing on the basis of environmental relevance to issues in global carbon cycling, alternative energy production, and biochemical importance. Table 2 contains the summary of project information.
Growth conditions and genomic DNA preparation
Since isolation by our research group, the strain has been kept in liquid cultures and a live culture and medium have been sent to DSMZ, (DSM21860). For DNA isolation batch cultures were grown in basal medium supplemented with 20 mM betaine as described by Westerholm et al. [6]. Cells were grown for 4 weeks at 37 °C without shaking and harvested at 5000 × g. DNA was isolated using the Blood & Tissue Kit from Qiagen (Hilden, Germany) according to the standard protocol recommended by the manufacturer.
Genome sequencing and assembly
The genome of Syntrophaceticus schinkii was sequenced at the SciLifeLab Uppsala, Sweden using Ion torrent PM systems with the mean length of 206 bp, longest read length 392 bp and a total of final library reads of 2,985,963 for single end reads. All general aspects of sequencing performed can be found at Scilifelab website [15]. The FastQC software package [16] was used for reads quality assessment. After preassembly quality checking, the reads were assembled with MIRA 4.0 and Newbler 2.8 assemblers. Possible miss-assemblies were corrected manually by using Tablet, a graphical viewer for visualization of assemblies and read mappings [17]. A comparison of two assemblies obtained from both of the assemblers was used to fill the gaps between contigs. The multiple genome alignment tool Mauve was used for this purpose [18]. The working draft genome sequence of S. schinkii Sp3 contains 3,196,921 bp based on the analysis done with the tools summarized above.
Genome annotation
Automated gene modeling was completed by MaGe [19] a bacterial genome annotation system. Genes were identified using Prodigal [20] and AMIGene [21] as part of MaGe genome annotation pipeline. The predicted CDSs were translated and used to search the NCBI non-redundant database, UniProt, TIGRFam, Pfam, PRIAM, KEGG, COG, and InterPro databases using BLASTP. Predicted coding sequences were subjected to manual analysis using MaGe web-based platform, which also provides functional information of proteins, and which was used to assess and correct genes predicted through the automated pipeline. The predicted functions were also further analyzed by the MaGe annotation system (Fig. 4).
Genome properties
The working draft genome comprises 301 contigs in 215 scaffolds with a total size of 3,196,921 bp and a calculated GC content of 46.59 %. The genome shows a protein coding density of 75.21 % with an average intergenic length of 230.2 bp. The genome encodes further 50 tRNA genes and 5 rRNA genes, more precisely three 5S genes, one 16S and one 23S rRNA gene (Table 3, Fig. 3).
The genome of S. schinkii genome contains 3,441 predicted protein-encoding genes, of which 2,099 (61 %) have been assigned tentative functions. The remaining 1,346 ORFs are hypothetical / unknown proteins. 2,586 (app. 75 %) of all predicted protein-encoding genes could be allocated to the 22 functional COGs. This is in the same range as described for other acetogenic bacteria such as Acetobacterium woodii WB1 and M. thermoacetica ATCC39073, acetate oxidizing sulfate reducers such as Desulfobacterium autotrophicum HRM2 and Desulfotomaculum kuznetsovii , and the SAOB P. lettingae TMO. Analysis of COGs revealed that ~28 % of all protein-encoding genes fall into four main categories: amino acid transport and metabolism (9.8 %), replication, recombination and repair (6.6 %), energy metabolism (5.9 %), and coenzyme transport and metabolism (4.9 %) (Table 4).
Insights from the genome sequence
Synteny-based analyses with all bacterial genomes present in the NCBI Reference Sequence database confirmed again that T. phaeum is the closest relative of S. schinkii having approximately 50 % of the total genome size in synteny (Fig. 4). A comparison of all inferred proteins of S. schinkii with all proteins collected in the NCBI RefSeq database revealed the highest number of orthologous (1788: 51.90 %) with T. phaeum . Both S. schinkii and T. phaeum , are known as syntrophic acetate oxidizing bacteria able to oxidize acetate in co-culture with a hydrogenotrophic methanogenic partner, but differ clearly in their substrate utilization patterns [2, 6] Moreover, in contrast to the thermophilic T. phaeum , S. schinkii possess mesophilic characteristics and cannot switch to a chemolithoautotrophic lifestyle.
The genome has been analyzed regarding general phenotypic features such as sporulation, oxygen tolerance, secreted and selenocystein-containing proteins and motility. The genome contains the master regulator Spo0A (SSCH_630004) needed for sporulation but lacks genes encoding the phosphorelays Spo0F and Spo0B as it has been observed in other clostridia. All the sporulation-specific sigma factors SigE (SSCH_460001), SigG (SSCH_1070017), and SigK (SSCH_700028) were predicted except for SigF. Two putative manganese containing catalases (SSCH_1760003, SSCH_2560004) and two putative rubrerythrin encoding genes (SSCH_590006, SSCH_180042) identified within the genome give reasons to believe, that this organism posses the ability to tolerate small amounts of oxygen. According to the observed immobility S. schinkii does not harbor any flagellum related genes including hook-associated proteins (FlgE, FlgK, FlgL), basal and hook proteins (FlgE), capping proteins (FliD), biosynthesis secretory proteins (FlhA, FlhB, FliF, FliH and FliI), flagella formation proteins, motor proteins (FliG and FliM) and the basal proteins (FlgC and FlgB).
Genes encoding key components of the selenocysteine-decoding (SelA, SelB, SelC, SelD) machinery are widely distributed in bacterial genomes. Also S. schinkii appears to have the ability to express selenocysteine proteins: The genome contains a single copy of the L-selenocysteinyl-tRNASec transferase (selA: SSCH_110005/6), monoselenophosphate synthase (selD: SSCH_970007), the selenocysteinyl-tRNA specific elongation factor (selB: SSCH_110004) and potential selenocysteine-specific tRNASec (selC: SSCH_tRNA31). We found two potential selenocysteine containing glycine/sarcosine/betaine reductase complexes encoded by the genome (SSCH_440002-8, SSCH_960012-15) consisting of selenoprotein subunit A, the substrate specific selenoprotein subunit B and acetyl phosphate forming subunit C. Since S. schinkii can only grow on betaine but not on glycine or sarcosine [6], this reductase complex might be specifically involved in betaine utilization. 57 CDSs were predicted to encode surface associated or secreted proteins identified by putative N-terminal signal peptides (signal peptide I and II).
Conclusions
Acetate oxidation under anoxic conditions is thermodynamically unfavorable and requires the metabolic cooperation of a partner organism in order to make endergonic reactions more exergonic through the efficient removal of the products. S. schinkii oxidizes acetate to hydrogen and/or formate, which is directly used by a hydrogenotrophic methanogen. Since the methanogenic partner has been isolated and sequenced S. schinkii appears to have great potential to serve as a model organism for studying methane producing syntrophic relationships. The working draft genome sequence presented here will open the door for understanding the preferred habitats, the metabolism behind different life styles, and the mechanisms initiating syntrophy. This knowledge will help us to trigger SAOB towards an efficient and stable hydrogen/biogas production in engineered anaerobic digestion processes suffering high ammonia release.
Abbreviations
- SAOB:
-
Syntrophic acetate-oxidizing bacteria
- DSMZ:
-
Deutsche Sammlung für Mikroorganismen und Zellkulturen
- MIRA:
-
Mimicking Intelligent Read Assembly
- MaGe:
-
Magnifying Genomes
- BLASTP:
-
Basic local alignment search tool for proteins
- NCBI:
-
National Center for Biotechnology Information
References
Ferry JG. Fermentation of Acetate. In: Ferry DJG, ed. Methanogenesis. Chapman & Hall Microbiology Series. Springer US; 1993:304–334. Available at: http://link.springer.com/chapter/10.1007/978-1-4615-2391-8_7. Accessed March 3, 2014.
Hattori S, Kamagata Y, Hanada S, Shoun H. Thermacetogenium phaeum gen. nov., sp. nov., a strictly anaerobic, thermophilic, syntrophic acetate-oxidizing bacterium. Int. J. Syst. Evol. Microbiol. 2000;50 Pt 4:1601–1609.
Petersen SP, Ahring BK. Acetate oxidation in a thermophilic anaerobic sewage-sludge digestor: the importance of non-aceticlastic methanogenesis from acetate. FEMS Microbiol. Lett. 1991;86:149–152.
Schnürer A, Svensson BH, Schink B. Enzyme activities in and energetics of acetate metabolism by the mesophilic syntrophically acetate-oxidizing anaerobe Clostridium ultunense. FEMS Microbiol. Lett. 1997;154:331–336.
Anna Schnürer BS. Clostridium ultunense sp. nov., a Mesophilic Bacterium Oxidizing Acetate in Syntrophic Association with a Hydrogenotrophic Methanogenic Bacterium. First Publ Int. J. Syst. Bacteriol. 1996:46(4);1145-1152.
Westerholm M, Roos S, Schnürer A. Syntrophaceticus schinkii gen. nov., sp. nov., an anaerobic, syntrophic acetate-oxidizing bacterium isolated from a mesophilic anaerobic filter. FEMS Microbiol. Lett. 2010;309:100–104.
Westerholm M, Roos S, Schnürer A. Tepidanaerobacter acetatoxydans sp. nov., an anaerobic, syntrophic acetate-oxidizing bacterium isolated from two ammonium-enriched mesophilic methanogenic processes. Syst. Appl. Microbiol. 2011;34:260–266.
Balk M, Weijma J, Stams AJM. Thermotoga lettingae sp. nov., a novel thermophilic, methanol-degrading bacterium isolated from a thermophilic anaerobic reactor. Int. J. Syst. Evol. Microbiol. 2002;52:1361–1368.
Oehler D, Poehlein A, Leimbach A, et al. Genome-guided analysis of physiological and morphological traits of the fermentative acetate oxidizer Thermacetogenium phaeum. BMC Genomics 2012;13:723.
Müller B, Manzoor S, Niazi A, Bongcam-Rudloff E, Schnürer A. Genome-guided analysis of physiological capacities of Tepidanaerobacter acetatoxydans provides insights into environmental adaptations and syntrophic acetate oxidation. PloS One 2015;10:e0121237.
Manzoor S, Müller B, Niazi A, Bongcam-Rudloff E, Schnürer A. Draft Genome Sequence of Clostridium ultunense Strain Esp, a Syntrophic Acetate-Oxidizing Bacterium. Genome Announc. 2013;1:e00107–13.
Anon. Classifier. Available at: https://rdp.cme.msu.edu/classifier/classifier.jsp. Accessed November 2, 2015.
Altschul SF, Madden TL, Schäffer AA, et al. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997;25:3389–3402.
Liolios K, Mavromatis K, Tavernarakis N, Kyrpides NC. The Genomes On Line Database (GOLD) in 2007: status of genomic and metagenomic projects and their associated metadata. Nucleic Acids Res. 2008;36:D475–D479.
Anon. SciLifeLab. Available at: https://www.scilifelab.se/. Accessed November 2, 2015.
Andrews S. FastQC A Quality Control tool for High Throughput Sequence Data. Available at: http://www.bioinformatics.babraham.ac.uk/projects/fastqc/
Milne I, Stephen G, Bayer M, et al. Using Tablet for visual exploration of second-generation sequencing data. Brief. Bioinform. 2013;14:193–202.
Darling AE, Mau B, Perna NT. ProgressiveMauve: Multiple Genome Alignment with Gene Gain, Loss and Rearrangement. PLoS ONE; 2010;5:e11147.
Vallenet D, Labarre L, Rouy Z, et al. MaGe: a microbial genome annotation system supported by synteny results. Nucleic Acids Res. 2006;34:53–65.
Hyatt D, Chen G-L, Locascio PF, Land ML, Larimer FW, Hauser LJ. Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinformatics 2010;11:119.
Bocs S, Cruveiller S, Vallenet D, Nuel G, Medigue C. AMIGene: Annotation of MIcrobial Genes. Nucleic Acids Res. 2003;31:3723–3726.
Field D, Garrity G, Gray T, et al. The minimum information about a genome sequence (MIGS) specification. Nat. Biotechnol. 2008;26:541–547.
Ludwig W, Schleifer K-H, Whitman WB. Revised road map to the phylumFirmicutes. In: Vos PD, Garrity GM, Jones D, et al., eds. Bergey’s Manual® of Systematic Bacteriology. Springer New York; 2009:1–13. Available at: http://link.springer.com/chapter/10.1007/978-0-387-68489-5_1. Accessed March 3, 2014.
Wolf M, Müller T, Dandekar T, Pollack JD. Phylogeny of Firmicutes with special reference to Mycoplasma (Mollicutes) as inferred from phosphoglycerate kinase amino acid sequence data. Int. J. Syst. Evol. Microbiol. 2004;54:871–875.
Gibbons NE, Murray RGE. Proposals Concerning the Higher Taxa of Bacteria. Int. J. Syst. Bacteriol. 1978;28:1–6.
Editor L. 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.
Rainey FA. Class II. Clostridia class nov. 2nd ed. (Vos PD, Garrity G, Jones D, et al., eds.). New York: Springer-Verlag; 2009.
Wiegel J. Order III. Thermoanaerobacterales ord. nov. 2nd ed. (Vos PD, Garrity G, Jones D, et al., eds.). New York: Springer-Verlag; 2009.
Wiegel J. Family I. Thermoanaerobacteraceae fam. nov. 2nd ed. (Vos PD, Garrity G, Jones D, et al., eds.). New York: Springer-Verlag; 2009.
Editor L. List of new names and new combinations previously effectively, but not validly, published. Int. J. Syst. Evol. Microbiol. 2011;61:1011–1013.
Ashburner M, Ball CA, Blake JA, et al. Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat. Genet. 2000;25:25–29.
Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32:1792–1797.
Felsenstein J. Evolutionary trees from DNA sequences: A maximum likelihood approach. J. Mol. Evol. 1981;17:368–376.
Tamura K, Dudley J, Nei M, Kumar S. MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) Software Version 4.0. Mol. Biol. Evol. 2007;24:1596–1599.
Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: Molecular Evolutionary Genetics Analysis Using Maximum Likelihood, Evolutionary Distance, and Maximum Parsimony Methods. Mol. Biol. Evol. 2011;28:2731–2739.
Felsenstein J. Confidence Limits on Phylogenies: An Approach Using the Bootstrap. Evolution 1985;39:783–791.
Acknowledgements
This work was supported by the Higher Education Commission, Pakistan, University of the Punjab, Lahore, Pakistan. Uppsala Genome Center performed sequencing supported by Science for Life Laboratory (Uppsala), the Swedish Bioinformatics Infrastructure for the Life Sciences supporting the SGBC bioinformatics platform at SLU, University of the Punjab, Lahore, Pakistan and Uppsala Multidisciplinary Center for Advanced Computational Science, Uppsala, Sweden. The contribution of SM and EB-R was supported by EU-COST action BM1006-SeqAhead. EB-R was also partially supported by EU FP7 ALLBIO project, grant number 289452, The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. We wish to thank Maria Westerholm, who isolated and characterized Syntrophaceticus schinkii strain Sp3 for providing the micrograph.
Author information
Authors and Affiliations
Corresponding author
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
SM, BM and AS contributed to the conception and design of this project. SM was involved in the acquisition and initial analysis of the data. SM and BM were involved in the interpretation of the data. SM prepared the first draft of the manuscript. EBR and AS provided financial support. All authors were involved in the critical revision of the manuscript and have given final approval of the version to be published and agree to be accountable for all aspects of the work.
Additional file
Additional file 1:
Table S1. Associated MIGS record for Syntrophaceticus schinkii strain Sp3. (DOCX 111 kb)
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Manzoor, S., Müller, B., Niazi, A. et al. Working draft genome sequence of the mesophilic acetate oxidizing bacterium Syntrophaceticus schinkii strain Sp3. Stand in Genomic Sci 10, 99 (2015). https://doi.org/10.1186/s40793-015-0092-z
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s40793-015-0092-z