- Short genome report
- Open Access
Genome sequencing and description of Oerskovia enterophila VJag, an agar- and cellulose-degrading bacterium
Standards in Genomic Sciencesvolume 12, Article number: 30 (2017)
A nonmotile, Gram-positive bacterium that shows an elongated and branching cell shape was isolated from soil samples from the botanical garden of Ulm University, Ulm, Germany. Here, the isolation procedure, identification, genome sequencing and metabolic features of the strain are described. Phylogenetic analysis allowed to identify the isolated strain as Oerskovia enterophila. The genus Oerskovia belongs to the family Cellulomonadaceae within the order Actinomycetales. The length of cells of O. enterophila ranges from 1 μm to 15 μm, depending on the growth phase. In the exponential growth phase, cells show an elongated and branching shape, whereas cells break up to round or coccoid elements in the stationary growth phase. The 4,535,074 bp long genome consists of 85 contigs with 3918 protein-coding genes and 57 RNA genes. The isolated strain was shown to degrade numerous complex carbon sources such as cellulose, chitin, and starch, which can be found ubiquitously in nature. Moreover, analysis of the genomic sequence revealed the genetic potential to degrade these compounds.
Oerskovia enterophila was formerly characterized as Promicromonospora enterophila by Jàger et al. in 1983 . Later, P. enterophila was re-classified as O. enterophila by Stackebrandt et al. , since only spore-like elements and no real spores are formed. Furthermore, a phylogenetic tree based on the 16S rRNA gene sequences of strains of the genera Cellulomonas and Promicromonospora shows that O. enterophila did not cluster with the type species of Promicromonospora , Promicromonospora citrea , or Promicromonospora sukumoe [2, 3]. The genus Oerskovia was initially described in 1970 by Prauser et al.  and harbors currently four species with O. turbata as type species . Bacteria of the genus Oerskovia belong to the phylum Actinobacteria , which is one of the largest taxonomic units among the domain Bacteria . Bacteria belonging to Actinobacteria show a wide range of G + C-content, from 51% to more than 70% [5,6,7]. Actinobacteria are widely distributed in terrestrial as well as in aquatic habitats [8, 9]. In general, members of the class Actinobacteria show a high morphological variety, which is also true for species of the genera Oerskovia and Cellulomonas . Furthermore, members of the family Cellulomonadaceae are known for their ability to decompose plant-derived biopolymers such as starch, cellulose or chitin . Due to the close relationship of members of the genera Oerskovia and Cellulomonas [12, 13] it is likely that both share genetic features enabling them to degrade these biopolymers. To investigate the genetic potential for biopolymer degradation, the genome of the isolate was sequenced. Furthermore, a genome wide comparison of the isolated strain with other Oerskovia type strains was performed. Additionally, the isolated strain was aerobically grown on respective carbon sources to validate the functionality of the proposed degradation pathways.
In this contribution, the classification, the metabolic features, and the genome insights of the isolated strain are provided.
Classification and features
The isolated strains were identified as Oerskovia enterophila based on 16S rRNA gene sequence identities of more than 99% compared to the type strain of O. enterophila DSM 43852 . All subsequent analyses were performed using the strain designated as O. enterophila VJag. Information regarding the enrichment and isolation procedures as well as identification of Oerskovia strains are described in the Additional files 1 and 2: S1 and S2.
Investigations of the cell morphology of the isolated strain O. enterophila VJag (Table 1) using scanning electron microscopy revealed that cells show different morphologies in exponential and stationary growth stage. In the exponential growth phase, cells show extensive branches with an overall length up to 15 μm, whereas the cells are smaller and less branched in the stationary growth phase (Fig. 1). These different cell morphologies were also previously observed by Stackebrandt et al. .
The 16S rRNA gene sequence (OJAG_11220, LRIE01000058.1) of O. enterophila VJag was blasted  and used for subsequent phylogenetic analysis. Therefore, 16S rRNA reference sequences of 17 closely related type strains were aligned using MAFFT version 7.215 [16, 17] and was performed using EMBL-EBI web services. The length of the 17 references ranged from 1395 to 1612 bp and had average length of 1486 bp. The phylogenetic tree was reconstructed using the software MrBayes version 3.2.6 . The recommended settings in the manual for tree reconstruction use a generalized time reversible evolutionary model. The quick start instructions were followed to run Bayesian phylogenetic analysis. The run was stopped since the standard deviation of split frequencies was below 0.0042 after 1,000,000 generations.
The resulting phylogenetic tree is shown in Fig. 2. Described species of the genera Oerskovia and Cellulomonas belong to the same family of Cellulomonadaceae . On the other hand, Sanguibacter belongs to the family of Sanguibacteriaceae which is defined as a neighboring group to Cellulomonadaceae . Sanguibacter is the only described genus within the respective family with currently six species [20,21,22,23,24].
Genome sequencing information
Genome project history
The genome of O. enterophila VJag was sequenced to get insights in the genomic features and the metabolic potential of this strain. Furthermore, no genomes of members of this species were available at the time of writing. A draft sequence is available at NCBI for the species O. turbata NRRL B-8019 (JOFV00000000) . The complete genome of O. enterophila VJag has a size of 4,535,074 bp and consists of 85 contigs. In this contribution the version LRIE01000000 is described. The genome sequencing and gene annotation was performed by Goettingen Genomics Laboratory (Germany). The sequence can be found under the accession number LRIE00000000. Table 2 shows the project information according to MIGS specification .
Growth conditions and genomic DNA preparation
O. enterophila VJag was cultivated in 5 ml TSYE-medium (medium 92, DSMZ) at 28 °C overnight in an orbital shaker at 120 rpm for the isolation of genomic DNA. Genomic DNA was isolated using MasterPure Gram positive DNA Purification kit (Epicentre, Madison, WI, USA) according to the manufacturer’s instructions. DNA concentrations and purity were analyzed using the UV-Vis spectrophotometer NanoDrop 2000 (Thermo Fisher Scientific, Waltham, MA, USA). The genomic DNA yield was 2463 ng/μl. The DNA purity was determined using the UV absorbance ratio 260/208 nm and 260/230 nm and revealed ratios of 2.01 and 2.17, respectively.
Genome sequencing and assembly
A combined approach was used for the whole-genome sequencing of O. enterophila VJag using the 454 GS-FLX TitaniumXL system (titanium GS70 chemistry, Roche Life Science, Mannheim, Germany) and the Genome Analyzer II (Illumina, San Diego, CA). According to the manufacturer’s protocols, the shotgun libraries were prepared, which resulted in 97,681 reads for 454 shotgun sequencing (11.46 × coverage) and 4,756,630 112-bp paired end Illumina reads (68.28 × coverage). Illumina reads were trimmed using Trimmomatic 0.32  to remove sequences with quality scores lower than 20 (Illumina 1.9 encoding) and remaining adaptor sequences, respectively. The initial hybrid de novo assembly was performed using the MIRA 3.4  and Newbler 2.9 (Roche Life Science, Mannheim, Germany) software. The final assembly resulted in 85 contigs with an average coverage of 79.60, an N50 value of 96,617 bp and an N90 value of 28,097 bp, respectively.
The Prodigal software tool  was used for automatic gene prediction , rRNA and tRNA gene identification was performed using RNAmmer  and tRNAscan , respectively. The automatic gene-annotation was performed by using the IMG-ER system [32, 33]. The annotation was manually curated using the Swiss-Prot, TrEMBL, and InterPro databases .
The genome of O. enterophila VJag is 4,535,074 bp in length and has an average G + C content of 72.4% (Fig. 3). The genome sequence shows 3975 genes in total, 3918 are protein-coding genes, 57 are RNA genes, of which 6 code for rRNA. The remaining genes code for proteins with unknown function or hypothetical proteins. All statistics and properties are listed in Table 3, the number of protein-coding genes associated with general COG functional categories is shown in Table 4.
A circular representation of the O. enterophila VJag genome sequence and comparison to O. enterophila DFA-19T  and O. turbata NRRL B-8019 genome sequences is shown in Fig. 3. For O. enterophila VJag the genes encoded by the leading and the lagging strand (outer circles 1 and 2) are marked in COG colors in the artificial chromosome map. The third and fourth circle show the positions of rRNA genes and tRNA genes, respectively. The fifth and sixth circle show comparisons of genes present in the strains O. enterophila VJag and O. enterophila DFA-19T chromosome as well as O. enterophila VJag and O. turbata NRRL B-8019 chromosome, respectively. The red colored regions indicate high similarity, whereas yellow colored regions indicate low similarity (see color code, Fig. 3). The two inner most plots represent the GC content and the GC skew (circle 7-8). Furthermore, a pairwise ANI analysis of the VJag strain and type strain O. enterophila DFA-19  showed a similarity value of 99.36%, whereas a respective analysis of VJag strain and O. turbata NRRL B-8019 resulted in 89.31% similarity.
Insights from the genome sequence
Because of the close relationship to members of the genus Cellulomonas , O. enterophila VJag was expected to use cellulose as carbon source. According to the KEGG pathway, genes coding for enzymes probably responsible for the degradation of cellulose to cellobiose and β-D-glucose were found in O. enterophila VJag. Cellulose is one of the main components of plant material and is one of the most abundant biopolymers in the environment . Plate assays revealed that O. enterophila VJag is able to utilize cellulose [Additional file 3: Figure S1]. The used plates contained CMC as sole carbon source and Congo red to stain CMC. O. enterophila VJag hydrolyzed CMC to glucose whereby the Congo red was eluted, the red color got lost and resulted in formation of bright halos around cell spots.
A gene (OJAG_15690) encoding a cellulose 1,4- β -cellobiosidase is present in genome that converts cellulose to 1,4-β-D-glucan. 1,4-β-D-glucan would be further converted to β-D-glucose through the action of a β-glucosidase. The genome sequence of O. enterophila VJag comprises 13 genes encoding such β-glucosidases (OJAG_01470, OJAG_39370, OJAG_33570, OJAG_33160, OJAG_31620, OJAG_25090, OJAG_25070, OJAG_16840, OJAG_16640, OJAG_15960, OJAG_15000, OJAG_14990, OJAG_11840).
Furthermore, cellulose can be converted to cellobiose, using endoglucanases (encoded by OJAG_04410; OJAG_07660), and can also be converted to β-D-glucose through the action of a β-glucosidase.
Starch is also ubiquitous in nature as it accumulates in plants as storage compound . The genome sequence of O. enterophila VJag harbors genes coding for α-amylases (OJAG_12050; OJAG_09450) and a starch phosphorylase (OJAG_12070). Thus, starch is either converted to glycogen, dextrin, or amylose by O. enterophila VJag. Starch or glycogen could also be degraded to trehalose by respective enzymes (glycogen debranching enzyme encoded by OJAG_00790 or OJAG_12120). Subsequently, trehalose would be further converted to β-D-glucose-1-phosphate or D-glucose via an α-trehalose phosphorylase (encoded by OJAG_12210). Dextrin would be converted to α-D-glucose by an oligo-1,6-glucosidase (encoded by OJAG_08510). A plate assay using Jag-MM-agar plates containing starch (2% w/v) as carbon source showed that starch is utilized during growth of O. enterophila VJag [Additional file 3: Figure S2]. After incubation, starch was stained using Lugol’s solution and bright halos around cell spots showed starch consumption by O. enterophila VJag (see Additional file 3: Figure S2).
Another commonly occurring compound in natural environments besides cellulose and starch is chitin. Chitin is a major structural polymer of the cell walls of fungi and the exoskeletons of invertebrates . Numerous genes which encode enzymes for the degradation of chitin to chitobiose or N-acetylglucosamine were found in the genome sequence of O. enterophila VJag (OJAG_26450; OJAG_36940; OJAG_36950; OJAG_36990; OJAG_38030; OJAG_38450; OJAG_38460). Chitobiose could be converted to N-acetylglucosamine or N-acetylglucosamine-1-phosphate by hexosaminidases and β-N-acetylhexosaminidase (encoded by OJAG_07390; OJAG_13640; OJAG_13650; OJAG_33360; OJAG_35500; OJAG_09950; OJAG_30030; OJAG_15920). Furthermore, one gene (OJAG_13250) was found that encodes a glucosamine-1-phosphate N-acetyltransferase, which converts N-acetylglucosamine-1-phosphate to UDP-acetylglucosamine. This intermediate would be further transformed to N-acetylglucosamine enopyruvate by an UDP-N-acetylglucosamine-1-carboxyvinyltransferase (encoded by OJAG_15040; OJAG_22690). N-acetylglucosamine enopyruvate can subsequently be converted to N-acetylmuramic acid via an UDP-N-acetylmuramate dehydrogenase (encoded by OJAG_01210). N-acetylmuramic acid would be metabolized through the peptidoglycan biosynthesis pathway or the D-glutamin and D-glutamate metabolism (OJAG_14230; _ OJAG 14240).
Additionally, genes encoding enzymes for xylose degradation were found in the O. enterophila VJag genome sequence. D-xylose could be converted to D-xylulose by a xylose isomerase (encoded by OJAG_26770). Furthermore, D-xylulose would be phosphorylated to D-xylulose-5-phosphate via a xylulokinase (OJAG_26780). D-xylulose-5-phosphate would be converted to D-ribulose-5-phosphate by a ribulose-5-phosphate 3-epimerase (OJAG_00210), and then metabolized via the pentose phosphate pathway, or D-xylulose-5-phosphate would be converted to L-ribulose -5-posphate via a L-ribulose-5-phosphate 4-epimerase (OJAG_27380). This also fits into the overall picture since xylose is a main part of hemicellulose and makes up a part of plant materials .
The genome of O. enterophila VJag, which was isolated from forest soil, is described. Furthermore, the phylogenetic and phenotypic characteristics of the isolated strain are presented. It has been shown that the isolate belongs to the family of Cellulomonadaceae . Scanning electron micrographs confirmed the variable phenotype in exponential or stationary growth phase. Genome sequences analysis revealed that O. enterophila VJag has the genetic properties to degrade compounds typically abundant in forest soils. Plate assays showed that the isolated strain is able to use starch and cellulose as sole carbon and energy source. The genome sequence of O. enterophila VJag has been deposited at DDBJ/EMBL/GenBank and can be found under the accession number LRIE00000000. The version described in this paper is version LRIE01000000.
Jáger K, Márialigeti K, Hauck M, Barabás G. Promicromonospora enterophila sp. nov., a new species of monospore Actinomycetes. Int J Syst Bacteriol. 1983;33:525–31.
Stackebrandt E, Breymann S, Steiner U, Prauser H, Weiss N, Schumann P. Re-evaluation of the status of the genus Oerskovia, reclassification of Promicromonospora enterophila (Jáger et al. 1983) as Oerskovia enterophila comb. nov. and description of Oerskovia jenensis sp. nov. and Oerskovia paurometabola sp. nov. Int J Syst Evol Microbiol. 2002;52:1105–11.
Takahashi Y, Tanaka Y, Iwai Y, Omura S. Promicromonospora sukumoe sp. nov., a new species of the Actinomycetales. J Gen Appl Microbiol. 1987;33:507–19.
Prauser H, Lechevalier MP, Lechevalier H. Description of Oerskovia gen. n., to harbor Ørskov’s motile Nocardia. Appl Microbiol. 1970;19:534.
Goodfellow M. Phylum XXVI. Actinobacteria phyl. nov. In: Goodfellow M, Kämpfer P, Busse H-J, Trujillo ME, Suzuki K-I, Ludwig W, Whitman WB, editors. Bergey’s manual of systematic bacteriology, vol. 5, Part A, vol. 2nd. New York: Springer; 2012. p. 33.
Lightfield J, Fram NR, Ely B. Across bacterial phyla, distantly-related genomes with similar genomic GC content have similar patterns of amino acid usage. PLoS ONE. 2011;6(3):e17677.
Hill LR. An index to deoxyribonucleic acid base compositions of bacterial species. J Gen Microbiol. 1996;44:419–37.
Ghai R, McMahon KD, Rodriguez-Valera F. Breaking a paradigm: cosmopolitan and abundant freshwater Actinobacteria are low GC. Environ Microbiol Rep. 2012;4:29–35. _2.
Servin JA, Herbold CW, Skophammer RG, Lake JA. Evidence excluding the root of the tree of life from the Actinobacteria. Mol Biol Evol. 2008;25:1–4.
Embley TM, Stackebrandt E. The molecular phylogeny and systematics of the Actinomycetes. Annu Rev Microbiol. 1994;48:257–89.
Stackebrandt E, Schumann P, Prauser H. The family Cellulomonadaceae, The prokaryotes. New York: Springer; 2006. p. 983–1001.
Stackebrandt E, Rainey FA, Ward-Rainey NL. Proposal for a new hierarchic classification system, Actinobacteria classis nov. Int J Syst Bacteriol. 1997;47:479–91.
Stackebrandt E, Häringer M, Schleifer KH. Molecular genetic evidence for the transfer of Oerskovia species into the genus Cellulomonas. Arch Microbiol. 1980;127:179–85.
Jag V, Poehlein A, Bengelsdorf FR, Daniel R, Dürre P. Genome sequence of the facultative anaerobe Oerskovia enterophila DFA-19 (DSM 43852T). Genome Announc. 2016. doi:10.1128/genomeA.00973-16.
Johnson M, Zaretskaya I, Raytselis Y, Merezhunk Y, McGinnis S, Madden TL. NCBI BLAST: a better web interface. Nucleic Acids Res. 2008. doi:10.1093/nar/gkn201.
Katoh K, Toh H. Improved accuracy of multiple ncRNA alignment by incorporating structural information into a MAFFT-based framework. BMC Bioinformatics. 2008. doi:10.1186/1471-2105-9-212.
Katoh K, Standley M. MAFFT sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30:772–80.
Ronquist F, Teslenko M, van der Mark P, Ayres DL, Darling A, Höhna S, Larget B, Liu L, Suchard MA, Huelsenbeck JP. MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Syst Biol. 2012;61:539–42.
Stackebrandt E, Prauser H. Assignment of the genera Cellulomonas, Oerskovia, Promicromonospora and Jonesia to Cellulomonadaceae fam. nov. Syst Appl Microbiol. 1991;14:261–5.
Fernández-Garayzábal JF, Dominguez L, Pascual C, Jones D, Collins MD. Phenotypic and phylogenetic characterization of some unknown coryneform bacteria isolated from bovine blood and milk: description of Sanguibacter gen. nov. Lett Appl Microbiol. 1995;20:69–75.
Pascual C, Collins MD, Grimont PAD, Dominguez L, Fernandez-Garayzabal JF. Sanguibacter inulinus sp. nov. Int J Syst Bacteriol. 1996;46:811–3.
Huang Y, Dai X, He L, Wang YN, Wang BJ, Liu Z, et al. Sanguibacter marinus sp. nov., isolated from coastal sediment. Int J Syst Evol Microbiol. 2005;55:1755–8.
Hong SG, Lee YK, Yim JH, Chun J, Lee HK. Sanguibacter antarcticus sp. nov., isolated from Antarctic sea sand. Int J Syst Evol Microbiol. 2008;58:50–2.
Kim MK, Pulla RK, Kim SY, Yi TH, Soung NK, Yang DC. Sanguibacter soli sp. nov., isolated from soil of a ginseng field. Int J Syst Evol Microbiol. 2008;58:538–41.
Doroghazi JR, Albright JC, Goering AW, Ju K-S, Haines RR, Tchalukov KA, et al. A roadmap for natural product discovery based on large-scale genomics and metabolomics. Nat Chem Biol. 2014;10:963–8.
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:541–7.
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30:2114–20. http://dx.doi.org/10.1093/bioinformatics/btu170.
Chevreux B, Wetter T, Suhai S. Genome sequence assembly using trace signals and additional sequence information. Computer Science and Biology: Proceedings of the German Conference on Bioinformatics. 1999;99:45–56.
Hyatt D, Chen G-L, LoCascio P, Land M, Larimer F, Hauser L. Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinformatics. 2010;11(1):119. http://www.biomedcentral.com/1471-2105/11/119.
Lagesen K, Hallin P, Rødland EA, Stærfeldt HH, Rognes T, Ussery DW. RNAmmer: consistent and rapid annotation of ribosomal RNA genes. Nucleic Acids Res. 2007;35:3100–8.
Lowe TM, Eddy SR. tRNAscan-SE: A program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997;25:955–64.
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–8.
Markowitz VM, Chen IM, Palaniappan K, Chu K, Szeto E, Grechkin Y, et al. IMG: the integrated microbial genomes database and comparative analysis system. Nucleic Acids Res. 2012;40:D115–22.
Zdobnov EM, Apweiler R. InterProScan – an integration platform for the signature-recognition methods in InterPro. Bioinformatics. 2001;17:847–8.
Leschine SB. Cellulose degradation in anaerobic environments. Annu Rev Microbiol. 1995;49:399–426.
Britz SJ, Hungerford WE, Lee DR. Rhythms during extended dark periods determine rates of net subsequent light periods in leaves of Sorghum. Planta. 1987;171:339–45.
Reguera G, Leschine SB. Chitin degradation by cellulolytic anaerobes and facultative aerobes from soils and sediments. FEMS Microbiol Lett. 2001;204:367–74.
Loewus FA, Kelly S, Neufeld EF. Metabolism of myo-inositol in plants: conversion to pectin, hemicellulose, D-xylose, and sugar acids. Proc Natl Acad Sci USA. 1962;48:421–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–9.
Buchanan RE. Studies in the nomenclature and classification of bacteria. II. The primary subdivisions of the Schizomycetes. J Bacteriol. 1917;2:155–64.
Skerman VBD, McGowan V, Sneath PHA. Approved lists of bacterial names. Int J Syst Evol Microbiol. 1980;30:225–420.
Zhi XY, Li WJ, Stackebrandt E. An update of the structure and 16S rRNA gene sequence-based definition of higher ranks of the class Actinobacteria, with the proposal of two new suborders and four new families and emended descriptions of the existing higher taxa. Int J Syst Evol Microbiol. 2009;59:589–608.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, et al. Gene ontology: tool for the unification of biology. Nat Genet. 2000;25:25–9.
Financial support by Ulm University is gratefully acknowledged.
AP and RD planned the genome sequencing, AP did the genome sequencing and the genome annotations, FB and VJ planned the isolation procedure, VJ performed the strain isolation, growth experiments, and the genomic DNA preparation and prepared the electron micrographs and the phylogenetic tree. VJ and AP wrote the manuscript. PD conceived of the study, participated in its design and coordinated and helped to draft the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Enrichment, isolation and selection of bacterial strains; identification of isolated strains (S1). (DOCX 15 kb)
Detailed composition of Jag-MM agar and silica plates (S2). (DOCX 17 kb)
Jag-MM-silica plates with CMC and Congo red; Figure S2. Jag-MM-agar plates with starch. (ZIP 964 kb)