- Short genome report
- Open Access
Complete genome sequences of Geobacillus sp. WCH70, a thermophilic strain isolated from wood compost
Standards in Genomic Sciencesvolume 11, Article number: 33 (2016)
Geobacillus sp. WCH70 was one of several thermophilic organisms isolated from hot composts in the Middleton, WI area. Comparison of 16 S rRNA sequences showed the strain may be a new species, and is most closely related to G. galactosidasius and G. toebii. The genome was sequenced, assembled, and annotated by the DOE Joint Genome Institute and deposited at the NCBI in December 2009 (CP001638). The genome of Geobacillus species WCH70 consists of one circular chromosome of 3,893,306 bp with an average G + C content of 43 %, and two circular plasmids of 33,899 and 10,287 bp with an average G + C content of 40 %. Among sequenced organisms, Geobacillus sp. WCH70 shares highest Average Nucleotide Identity (86 %) with G. thermoglucosidasius strains, as well as similar genome organization. Geobacillus sp. WCH70 appears to be a highly adaptable organism, with an exceptionally high 125 annotated transposons in the genome. The organism also possesses four predicted restriction-modification systems not found in other Geobacillus species.
Originally classified as members of the genus Bacillus , Geobacillus species were reclassified into a separate genus based on properties such as 16S rRNA gene sequence analysis, lipid and fatty acid analysis, phenotypic characterization, and DNA–DNA hybridization experiments . Geobacillus species have been isolated from high-temperature oilfields , a corroded pipeline in an extremely deep well , American [4, 5] African  and Russian  hot springs, marine vents , and the Mariana Trench . In addition to these extreme environments, Geobacillus species are commonly found in composting materials . Geobacillus. sp. WSUCF1 , G. galactosidasius  and G. toebii  were isolated from high-temperature composts. The ability of Geobacillus species to thrive in these varied and often hostile environments suggests that these species possess enzymes suitable for applications in challenging industrial environments [14, 15]. As part of a program to identify organisms, we isolated Geobacillus species from a variety of composts in Middleton, WI. We report here the isolation and genome sequence of Geobacillus sp. WCH70, isolated from high-temperature wood compost.
Classification and features
Geobacillus sp. WCH70 is a novel thermophilic species isolated from a hot wood compost pile (~70 °C) in Middleton, WI (43.097090° latitude and -89.504730° longitude). The organism was isolated from a piece of decaying wood by enrichment and plating on YTP-2 medium (YTP-2 media contains (per liter) 2.0 g yeast extract, 2.0 g tryptone, 2.0 g sodium pyruvate, 1.0 g KCl, 2.0 g KNO3, 2.0 g Na2HPO4.7H2O, 0.1 g MgSO4, 0.03 g CaCl2, and 2.0 ml clarified tomato juice) at 70 °C. The culture is available from the Bacillus Genetic Stock Center. Cultures are routinely grown on tryptic soy broth without glucose (Difco) media and maintained on TSB agar plates. C5-6 Technologies, Lucigen, and the Joint Genome Institute have placed no restrictions on the use of the culture or sequence data. Geobacillus sp. WCH70 is a Gram-positive, rod-shaped facultative anaerobe (Table 1), with optimum growth temperature of 70 °C and maximum growth temperature of 80 °C. Geobacillus sp. WCH70 appears to grow as a mixture of single cells and large clumps in liquid culture (Fig. 1).
A phylogenetic tree was constructed to identify the relationship of Geobacillus sp. WCH70 to other members of the Geobacillus family (Fig. 2). The phylogeny of Geobacillus sp. WCH70 was determined using one of the ten16S rRNA gene sequence (genome coordinates 10256 through 11801), as well as those of the type strains of all validly described Geobacillus spp. The 16S rRNA gene sequences were aligned using MUSCLE , pairwise distances were estimated using the Maximum Composite Likelihood (MCL) approach, and initial trees for heuristic search were obtained automatically by applying the Neighbour-Joining method in MEGA 5 . The alignment and heuristic trees were then used to infer the phylogeny using the Maximum Likelihood method based on Tamura-Nei . Comparison of 16 S rRNA sequences shows Geobacillus sp. WCH70 clades with other 42 to 45 % G + C content species including G. thermoglucosidasius , G. caldoxylolyticus, G. galactosidasius and G. toebii and is most closely related to G. galactosidasius and G. toebii . Bootstrap analysis indicates that G. galactosidasius and G. toebii are more closely related to each other than to Geobacillus sp. WCH70, suggesting Geobacillus sp. WCH70 may be a new Geobacillus sp. Essentially identical trees were obtained when the other nine Geobacillus sp. WCH70 16S rRNA gene sequences were used to generate phylogenetic trees.
Genome sequencing information
Genome project history
Geobacillus sp. WCH70 was selected for sequencing on the basis of its biotechnological potential as part of the U.S. Department of Energy’s Genomic Science program (formerly Genomics:GTL). The genome sequence is deposited in the Genomes On Line Database [19, 20] (GOLD ID = Ga0028898), and in GenBank (NCBI Reference Sequence = CP001638.1). Sequencing, finishing and annotation were performed by the DOE JGI. A summary of the project information and its association with MIGS identifiers is shown in Table 2.
Growth conditions and genomic DNA preparation
For preparation of genomic DNA, one liter cultures of Geobacillus sp. WCH70 were grown from a single colony in YTP-2 medium at 70 °C in flasks agitated at 200 rpm and collected by centrifugation. Culture stocks were maintained on YTP-2 agar plates grown at 70 °C. The cell concentrate was lysed using a combination of SDS and proteinase K, and genomic DNA was isolated using a phenol/chloroform extraction. The genomic DNA was precipitated, and treated with RNase to remove residual contaminating RNA. The purity and concentration of the recovered DNA was determined by gel electrophoresis in 0.7 % agarose containing ethidium bromide. Low and high molecular weight lambda DNA ladders were used as standards. The purity,and quantity of the recovered DNA was also independently confirmed by the JGI as suitable for sequencing prior to initiation of the project.
Genome sequencing and assembly
The genome of Geobacillus sp. WCH70 was sequenced at the JGI using a combination of Sanger and 454 technologies . Two Sanger libraries with average insert size of 8 Kb and 40 Kb (fosmid) were generated for this genome. In addition to Sanger sequencing, 454 pyrosequencing was done to a depth of 20x coverage. Draft assemblies were based on 52,102 total reads. All three libraries provided 12.7x coverage of the genome. The Phred/Phrap/Consed software package was used for sequence assembly and quality assessment [22–24] in the following finishing process. After the shotgun stage, reads were assembled with parallel phrap (High Performance Software, LLC). Possible mis-assemblies were corrected with gapResolutioin (Cliff Han, unpublished), Dupfinisher, 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. A total of 2,285 additional reactions were necessary to close gaps and to raise the quality of the finished sequence. The completed genome sequences of Geobacillus contains 56,142 reads, achieving an average of 13-fold sequence coverage per base with an error rate less than 1 in 100,000.
Genes were identified using Prodigal  as part of the Oak Ridge National Laboratory genome annotation pipeline, followed by a round of manual curation using the JGI GenePRIMP pipeline . The predicted CDSs were translated and used to search the National Center for Biotechnology Information nonredundant database, UniProt, TIGRFam, Pfam, PRIAM, KEGG, COG, and InterPro databases. These data sources were combined to assert a product description for each predicted protein. Non-coding genes and miscellaneous features were predicted using tRNAscan-SE , RNAMMer , Rfam , TMHMM , and signalP .
The genome of Geobacillus sp. WCH70 consists of one circular chromosome (Table 3 and Fig. 3) of 3,464,618 bp and an average G + C content of 43 % and two circular plasmids of 33,899 and 10,287 bp and an average G + C content of 40 % (Table 4). There are 92 tRNA genes and 28 rRNA genes. There are 3,477 predicted protein-coding regions and 309 pseudogenes in the genome. A total of 2,373 genes (66.0 %) have been assigned a predicted function while the rest have been designated as hypothetical proteins (Table 4). The numbers of genes assigned to each COG functional category are listed in Table 5. About 39 % of the annotated genes were not assigned to a COG or have an unknown function.
Insights from the genome sequence
The genome sequence of Geobacillus sp. WCH70 was compared to Geobacillus species with sequenced genomes. The lack of genome sequence information for G. galactosidasius and G. toebii prevents direct comparisons with these two organisms that are most closely related to Geobacillus sp. WCH70 based on rRNA gene sequences. Geobacillus sp. WCH70 Average Nucleotide Identity values  were 86.5 to 86.7 % to five G. thermoglucosidasius strains, 85.2 % to G. stearothermophilus NUB3621, and 84.9 and 85.0 % for two G. caldoxylosilyticus strains. ANI values ranged from 75.3 to 76.3 % for 20 other Geobacillus strains including G. stearothermophilus ATCC 7953, G. thermodenitrificans DSM 465, G. subterraneus PSS2, and G. kaustophilus HTA426. These values mirror the relationships of Geobacillus sp. WCH70 to other species seen in the phylogenetic tree based on rRNA. In addition to being closely related to G. thermoglucosidasius strains based on these two criteria, synteny plots reveal highly similar genome organizations in Geobacillus sp. WCH70 and G. thermoglucosidasius C56-YS93 (Fig. 4).
Geobacillus sp. WCH70 possesses a number of unusual features when compared to other Geobacillus species. A major feature of Geobacillus sp. WCH70 is the presence of 125 insertion-sequence (IS) elements predicted to code for transposons, significantly more elements than are found in sequenced strains of either G. thermoglucosidasius or G. caldoxylosilyticus (Table 6). In addition to these IS elements, Geobacillus sp. WCH70 possesses four predicted restriction-modification gene clusters not found in other Geobacillus species. Genes GWCH70_1298 through GWCH70_1302 code for a predicted Type I restriction system most closely related to a system in B. cereus VD021, while GWCH70_2032 through GWCH70_2034 and GWCH70_3440 through GWCH70_3444 code for predicted Type I restriction systems most closely related to systems in B. coagulans XZL4. Genes GWCH70_2067 through GWCH70_2069 code for a predicted Type III restriction system most closely related to a system in Thermincola ferriacetica DSM 14005™. Genes GWCH70_1385 and GWCH70_1386 code for restriction system proteins most closely related to proteins in Streptosporangium roseum DSM 43021™. These restriction systems may facilitate transfer of DNA to and from other organisms in the compost microbiome.
Surprisingly, the genome of Geobacillus sp. WCH70 is lacking many of the predicted polysaccharide degradation clusters seen in other Geobacillus species , including the metabolic cluster for degrading hemicellulose . The organism may utilize starch and other alpha-glucans based on the presence of a eleven-gene cluster GWCH70_0695 through GWCH70_0704 that is predicted to code for two, three-gene ABC carbohydrate transport systems, three alpha-amylase catalytic regions, an alpha-glucosidase, and a LacI family transcriptional regulator.
Geobacillus sp. WCH70 is a thermophilic gram-positive, spore-forming organism isolated from hot wood compost in the Middleton, WI area. Comparison of 16 S rRNA sequences showed the strain may be a new species, and is most closely related to G. galactosidasius and G. toebii . The genome of Geobacillus has an average G + C content of 43 %, similar to that reported for G. toebii (43.9 %) . G. galactosidasius is reported to possess a 53.5 % average G + C content  significantly higher than the value for Geobacillus sp. WCH70. Six G. thermoglucosidasius strains have 43.8 to 44.0 % average G + C content based on genomic sequence , similar to the value obtained for Geobacillus sp. WCH70. These G + C content values are lower than the 53 to 54 % obtained using chemical analyses [1, 12]. Genomic sequencing of G. galactosidasius and G. toebii is necessary to clarify the relationships among Geobacillus sp. WCH70, G. galactosidasius and G. toebii , and G. thermoglucosidasius .
The presence of 125 insertion-sequence (IS) elements predicted to code for transposons along with multiple restriction-modification systems suggests Geobacillus sp. WCH70 possesses a highly mutable chromosome, able to add or delete non-essential genes and gene clusters depending on the environmental conditions. Genomic sequencing of other Geobacillus species may help clarify if this mutability is a common element in other organisms in composts, or unique to Geobacillus sp. WCH70.
integrate microbial genomes database
joint genome institute
Nazina TN, Tourova TP, Poltaraus AB, Novikova EV, Grigoryan AA, Ivanova AE, et al. Taxonomic study of aerobic thermophilic bacilli: descriptions of Geobacillus subterraneus gen. nov., sp. nov. and Geobacillus uzenensis sp. nov. from petroleum reservoirs and transfer of Bacillus stearothermophilus, Bacillus thermocatenulatus, Bacillus thermoleovorans, Bacillus kaustophilus, Bacillus thermodenitrificans to Geobacillus as the new combinations G. stearothermophilus, G. thermocatenulatus, G. thermoleovorans, G. kaustophilus, G. thermoglucosidasius and G. thermodenitrificans. Int J Syst Evol Microbiol. 2001;51(Pt 2):433–46.
Kuisiene N, Raugalas J, Chitavichius D. Geobacillus lituanicus sp. nov. Int J Syst Evol Microbiol. 2004;54(Pt 6):1991–5.
Popova NA, Nikolaev Iu A, Turova TP, Lysenko AM, Osipov GA, Verkhovtseva NV, et al. Geobacillus uralicus, a new species of thermophilic bacteria. Mikrobiologiia. 2002;71(3):391–8.
Brumm P, Land M, Hauser LJ, Jeffries C, Chang YJ, Mead D. Complete genome sequence of Geobacillus strain Y4.1MC1, a novel CO-utilizing Geobacillus thermoglucosidasius strain isolated from Bath Hot Spring in Yellowstone National Park. Bioenerg Res. 2015;8(3):1039–1045.
Brumm P, De Maayer P, Cowan DA, MEAD DA. Genomic analysis of six new Geobacillus strains reveals highly conserved carbohydrate degradation architectures and strategies. Front Microbiol. 2015;6:430.
Hawumba JF, Theron J, Brozel VS. Thermophilic protease-producing Geobacillus from Buranga hot springs in Western Uganda. Curr Microbiol. 2002;45(2):144–50.
Nazina TN, Lebedeva EV, Poltaraus AB, Tourova TP, Grigoryan AA, Sokolova D, et al. Geobacillus gargensis sp. nov., a novel thermophile from a hot spring, and the reclassification of Bacillus vulcani as Geobacillus vulcani comb. nov. Int J Syst Evol Microbiol. 2004;54(Pt 6):2019–24.
Maugeri TL, Gugliandolo C, Caccamo D, Stackebrandt E. Three novel halotolerant and thermophilic Geobacillus strains from shallow marine vents. Syst Appl Microbiol. 2002;25(3):450–5.
Takami H, Nishi S, Lu J, Shimamura S, Takaki Y. Genomic characterization of thermophilic Geobacillus species isolated from the deepest sea mud of the Mariana Trench. Extremophiles. 2004;8(5):351–6.
Li R, Li L, Huang R, Sun Y, Mei X, Shen B, et al. Variations of culturable thermophilic microbe numbers and bacterial communities during the thermophilic phase of composting. World J Microbiol Biotechnol. 2014;30(6):1737–46.
Bhalla A, Kainth AS, Sani RK. Draft Genome sequence of lignocellulose-degrading thermophilic bacterium Geobacillus sp. Strain WSUCF1. Genome Announc. 2013;1(4):e00595-13.
Poli A, Laezza G, Gul-Guven R, Orlando P, Nicolaus B. Geobacillus galactosidasius sp. nov., a new thermophilic galactosidase-producing bacterium isolated from compost. Syst Appl Microbiol. 2011;34(6):419–23.
Sung MH, Kim H, Bae JW, Rhee SK, Jeon CO, Kim K, et al. Geobacillus toebii sp. nov., a novel thermophilic bacterium isolated from hay compost. Int J Syst Evol Microbiol. 2002;52(Pt 6):2251–5.
de Miguel BT, Barros-Velazquez J, Villa TG. Industrial applications of hyperthermophilic enzymes: a review. Protein Pept Lett. 2006;13(7):645–51.
Bergquist PL, Morgan HW, Saul D. Selected enzymes from extreme thermophiles with applications in biotechnology. Curr Biotechnol. 2014;3:45–59.
Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32(5):1792–7.
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(10):2731–9.
Tamura K, Nei M. Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Mol Biol Evol. 1993;10(3):512–26.
Liolios K, Tavernarakis N, Hugenholtz P, Kyrpides NC. The Genomes On Line Database (GOLD) v. 2: a monitor of genome projects worldwide. Nucleic Acids Res. 2006;34(Database issue):D332–4.
Liolios K, Chen IM, Mavromatis K, Tavernarakis N, Hugenholtz P, Markowitz VM, et al. The Genomes On Line Database (GOLD) in 2009: status of genomic and metagenomic projects and their associated metadata. Nucleic Acids Res. 2010;38(Database issue):D346–54.
Margulies M, Egholm M, Altman WE, Attiya S, Bader JS, Bemben LA, et al. Genome sequencing in microfabricated high-density picolitre reactors. Nature. 2005;437(7057):376–80.
Ewing B, Green P. Base-calling of automated sequencer traces using phred. II. Error probabilities. Genome Res. 1998;8(3):186–94.
Ewing B, Hillier L, Wendl MC, Green P. Base-calling of automated sequencer traces using phred. I. Accuracy assessment. Genome Res. 1998;8(3):175–85.
Gordon D, Abajian C, Green P. Consed: a graphical tool for sequence finishing. Genome Res. 1998;8(3):195–202.
Hyatt D, Chen GL, Locascio PF, Land ML, Larimer FW, Hauser LJ. Prodigal Prokaryotic Dynamic Programming Genefinding Algorithm. BMC Bioinformatics. 2010;11:119.
Lowe TM, Eddy SR. tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997;25(5):955–64.
Lagesen K, Hallin P, Rodland EA, Staerfeldt HH, Rognes T, Ussery DW. RNAmmer: consistent and rapid annotation of ribosomal RNA genes. Nucleic Acids Res. 2007;35(9):3100–8.
Griffiths-Jones S, Bateman A, Marshall M, Khanna A, Eddy SR. Rfam: an RNA family database. Nucleic Acids Res. 2003;31(1):439–41.
Krogh A, Larsson B, von Heijne G, Sonnhammer EL. Predicting transmembrane protein topology with a hidden Markov model: application to complete genomes. J Mol Biol. 2001;305(3):567–80.
Kim M, Oh HS, Park SC, Chun J. Towards a taxonomic coherence between average nucleotide identity and 16S rRNA gene sequence similarity for species demarcation of prokaryotes. Int J Syst Evol Microbiol. 2014;64(Pt 2):346–51.
De Maayer P, Brumm PJ, Mead DA, Cowan DA. Comparative analysis of the Geobacillus hemicellulose utilization locus reveals a highly variable target for improved hemicellulolysis. BMC Genomics. 2014;15(1):836.
Markowitz VM, Chen IM, Palaniappan K, Chu K, Szeto E, Pillay M, et al. IMG 4 version of the integrated microbial genomes comparative analysis system. Nucleic Acids Res. 2014;42(Database issue):D560–7.
Field D, Garrity G, Gray T, Morrison N, Selengut J, Sterk P, et al. The minimum information about a genome sequence (MIGS) specification. Nat Biotechnol. 2008;26(5):541–7.
Woese CR, Kandler O, Wheelis ML. Towards a natural system of organisms: proposal for the domains Archaea, Bacteria, and Eucarya. Proc Natl Acad Sci U S A. 1990;87(12):4576–9.
Gibbons NE, Murray RGE. Proposals Concerning the Higher Taxa of Bacteria. Int J Syst Bacteriol. 1978;28:1–6.
Ludwig W, Schleifer KH, Whitman WB. Class I. Bacilli class nov. Bergey’s Manual of Systematic Bacteriology. 2009;3:19–20.
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–72.
Prevot AR. Dictionnaire des Bactéries Pathogènes. 1953. p. 1–692.
Skerman VBD, McGowan V, Sneath PHA. Approved Lists of Bacterial Names. Int J Syst Bacteriol. 1980;30:225–420.
Fischer A. Untersuchungen über bakterien. Jahrbücher für Wissenschaftliche Botanik. 1895;27:1–163.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, et al. Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000;25(1):25–9.
Felsenstein J. Confidence limits on phylogenies: An approach using the bootstrap. Evolution. 1985;10:512–26.
This work was funded by the DOE Great Lakes Bioenergy Research Center (DOE Office of Science BER DE-FC02-07ER64494). Sequencing work was performed under the auspices of the US Department of Energy’s Office of Science, Biological and Environmental Research Program, and by the University of California, Lawrence Berkeley National Laboratory under contract No. DE-AC02-05CH11231, Lawrence Livermore National Laboratory under Contract No. DE-AC52-07NA27344, and Los Alamos National Laboratory under contract No. DE-AC02-06NA25396.
Great Lakes Bioenergy Research Center, C5•6 Technologies, and Lucigen Corporation provided support in the form of salaries for authors PB, & DM, but did not have any additional role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript. The specific roles of these authors are articulated in the ‘author contributions’ section.
PJB isolated and characterized the organism and wrote the manuscript, DAM managed the DNA preparation and submission to JGI for sequencing, and MLL did the genome annotation and document editing. All authors read and approved the final manuscript.