- Short genome report
- Open Access
High quality draft genomic sequence of Flavobacterium enshiense DK69T and comparison among Flavobacterium genomes
Standards in Genomic Sciences volume 10, Article number: 92 (2015)
Flavobacterium enshiense DK69T is a Gram-negative, aerobic, rod-shaped, non-motile and non-flagellated bacterium that belongs to the family Flavobacteriaceae in the phylum Bacteroidetes. The high quality draft genome of strain DK69T was obtained and has a 3,375,260 bp genome size with a G + C content of 37.7 mol % and 2848 protein coding genes. In addition, we sequenced five more genomes of Flavobacterium type strains and performed a comparative genomic analysis among 12 Flavobacterium genomes. The results show some specific genes within the fish pathogenic Flavobacterium strains which provide information for further analysis the pathogenicity.
Flavobacterium enshiense DK69T (= CCTCC AB2011144T = KCTC 23775T ) is a type strain that belongs to the genus Flavobacterium of the family Flavobacteriaceae . In recent years, members of Flavobacterium were identified and widely distributed in soil, fresh water, marine water, sediment, microbial mat, and glaciers [2–5]. Some Flavobacterium strains are fish pathogens including Flavobacterium columnare ATCC 49512T causing columnaris disease , Flavobacterium psychrophilum JIP02/86T causing cold-water disease  and Flavobacterium branchiophilum FL-15T causing bacterial gill disease .
The common characters of Flavobacterium strains are Gram-negative, non-spore-forming, yellow-pigmented, rod-shaped, aerobic and with a low DNA G + C content (30–41 mol %) [2–12]. The Flavobacterium strains contained iso-C15:0 as the major fatty acid, phosphatidylethanolamine as the major polar lipid and menaquinone-6 as the major respiratory quinone [9–12].
In order to provide genome information of Flavobacterium species, we sequenced six Flavobacterium strains including F. enshiense DK69T , Flavobacterium beibuense F44-8T , Flavobacterium cauense R2A-7T , Flavobacterium rivuli WB 3.3-2T , Flavobacterium subsaxonicum WB 4.1-42T  and Flavobacterium suncheonense GH29-5T . In this study, we compared 12 genomes including the six strains that we sequenced and other six available Flavobacterium genomes in the NCBI, Flavobacterium indicum GPTSA100-9T , Flavobacterium frigoris PS1T , Flavobacterium sp. F52 , Flavobacterium columnare ATCC 49512T , Flavobacterium psychrophilum JIP02/86T and Flavobacterium branchiophilum FL-15T. Here, we present the description of the non-contiguous finished genomic sequencing of F. enshiense DK69T and the comparative genome analysis of the 12 Flavobacterium genomes.
Classification and features
F. enshiense DK69T is a Gram-negative, strictly aerobic, yellow-pigmented rod shaped bacterium isolated from soil collected at a pharmaceutical company in Enshi, Hubei province, China. The total soil C, N, P, S and Fe concentrations were 39.83, 3.34, 0.68, 0.36, 33.80 g kg−1, respectively, and the pH was 6.97 . A neighbor-joining phylogenetic tree based on the 16S rRNA gene sequences was built using MEGA 6  and showed that strain DK69T was clustered within a branch containing other species in the genus Flavobacterium (Fig. 1). In addition, the sequence of F. enshiense DK69T was compared with other sequenced strains of the family Flavobacteriaceae use BioLinux , and a total of 24 core protein sequences were obtained with 50 % identity and E-value exponent of e−10. A phylogenetic tree based on the 24 core protein sequences of the core genome (Fig. 2) is similar to the 16S rRNA gene based tree.
The colonies of F. enshiense DK69T are smooth with regular edges, circular, yellowish and about 1 mm in diameter after grown on R2A agar at 28 °C for 48 h. Growth occurs at 4–32 °C, pH 6.0–8.0 on R2A and TSA, but not on NA or LB media, and NaCl is not required . Cells are non-flagellated, non-spore-forming, non-motile, rod-shaped (Fig. 3). Oxidase- and catalase- positive. The DNA G + C content is 34.4 mol% . The general description of this strain is shown in Table 1.
The major cellular fatty acids of F. enshiense DK69T were iso-C15:0, iso-C17:1 ω9c, C15:0, iso-C17:0 3-OH and iso-C15:0 3-OH. The major polar lipids were phosphatidylethanolamine, one unidentified aminolipid and one unidentified lipid. F. enshiense DK69T contained menaquinone 6 as the major quinone .
Genome sequencing information
Genome project history
Genome of F. enshiense DK69T was sequenced by Majorbio Bio-pharm Technology Co., Ltd, Shanghai, China. The high-quality draft genome sequence was deposited in the National Center for Biotechnology Information. Contigs less than 200 bp were not included. The GenBank accession number is JRLZ00000000. The summary of the genome sequencing project information is shown in Table 2.
Growth conditions and genomic DNA preparation
F. enshiense DK69T was grown on R2A medium at 28 °C for 2 d with 160 rpm shaking. Cells in late-log-phase growth were harvested and lysed by EDTA, lysozyme, and detergent treatment, followed by proteinase K and RNase digestion. The DNA was extracted and purified using the QiAamp kit according to the manufacturer’s instruction (Qiagen, Germany). The quantity of DNA was measured by the NanoDrop Spectrophotometer to ensure that the DNA concentration is greater than 20 ng/μl, then 5 μg of DNA was sent to Majorbio (Shanghai, China) for sequencing.
Genome sequencing and assembly
The Illumina Hiseq2000 with the Paired-End library strategy was used to determine the whole-genome sequence of F. enshiense DK69T . TruSeq DNA Sample Preparation Kits are used to prepare DNA libraries with insert sizes of 300–500 bp for single, paired-end, and multiplexed sequencing. The protocol used 1 μg of DNA sheared by either sonication or nebulization . The genome raw data of F. enshiense DK69T generated 8,329,997 x 2 reads totaling 1,682,659,394 bp data with an average coverage of 498.4 x. Then SOAPdenovo v1.05  was used to perform the following steps to assemble the sequencing data: (1) removing the adapter sequences in the reads; (2) cutting the 5’ end bases without clear A, T, C and G; (3) trimming the quality read scores lower than 20; (4) removing the reads containing more than 10 % Ns; (5) removing the reads which the length were less than 25 bp. A total of 8,217,761 x 2 high quality reads totaling 1,645,393,073 bp data with an average coverage 487.4 × was generated. The assembled sequence contained 67 scaffolds with a genome size of 3.38 Mbp.
The annotation of the genomic sequences was completed using the NCBI Prokaryotic Genome Annotation Pipeline which was combined using Best-placed reference protein set and the gene caller GeneMarkS+. SignalP  and SOSUI  were used to predict signal peptides and transmembrane helices. The predicted CDSs were also used to search against the Pfam protein family database . The GenBank database  and the COG databases  BLASTP search were used to predict protein sequences.
The genome statistics are provided in Table 3 and Fig. 4. After genome annotation, the genome of F. enshiense DK69T was found to have a total length of 3,375,260 bp, a G + C content of 1,273,385 bp (37.7 mol %) and 74 contigs. From a total of 3,054 genes predicted, 2,848 genes are protein-coding genes, 50 are RNA genes, 57.9 % are assigned with putative functions and the remaining are annotated as hypothetical proteins or proteins of unknown functions. The distribution of genes into COGs functional categories is shown in Table 4.
Insights from the genome sequences
Profiles of metabolic network and pathway
The metabolic network and pathways of F. enshiense DK69T (Fig. 5) were predicted using the Kyoto Encyclopedia of Genes and Genomes . The metabolic network showed that F. enshiense DK69T possesses glycolysis, TCA cycle and pentose phosphate pathways and could utilize casein, tyrosine, sucrose and D-mannitol. The genome analysis results are in agreement with the phenotypes .
Comparison of the 12 Flavobacterium genomes
The genomic information of the 12 Flavobacterium genomes are summarized in Table 5. OrthoMCL  analysis was performed to identify the set of orthologs among the 12 Flavobacterium genomes. F. enshiense DK69T shared 1,190 genes with the other 11 Flavobacterium strains, and had 437 strain-specific genes which may contribute to the species-specific features (Fig. 6).
Three of the 12 Flavobacterium strains are fish pathogenic bacteria [6–8]. Using OrthoMCL  analysis, a total of ten proteins we found to be unique in the three fish-pathogenic species. Three of the putative proteins were reported to be related to the pathogenicity of pathogenic bacteria including polysaccharide deacetylase , ABC transporter ATPase and ABC transporter permease  (Table 6).
The genomic results of F. enshiense DK69T and related strains reveled useful information. (1) The genome based phylogenetic analysis results is in agreement with the 16S rRNA gene based one; (2) The genomic data are correlated with some phenotypes of strain DK69T; (3) Compared to the three fish pathogenic Flavobacterium strains, no pathogenic related genes was detected in the environmental strain DK69T which indicated its non-pathogenicity; and (4) Some specific genes were found within the three fish pathogenic Flavobacterium strains which provides information for further analysis the pathogenicity.
Dong K, Chen F, Du Y, Wang G. Flavobacterium enshiense sp. nov., isolated from soil, and emended descriptions of the genus Flavobacterium and Flavobacterium cauense, Flavobacterium saliperosum and Flavobacterium suncheonense. Int J Syst Evol Microbiol. 2013;63:886–92. doi:10.1099/ijs.0.039974-0.
Kim BY, Weon HY, Cousin S, Yoo SH, Kwon SW, Go SJ, et al. Flavobacterium daejeonense sp. nov. and Flavobacterium suncheonense sp. nov., isolated from greenhouse soils in Korea. Int J Syst Evol Microbiol. 2006;56:1645–9. doi:10.1099/ijs.0.64243-0.
Yoon JH, Kang SJ, Oh TK. Flavobacterium soli sp. nov., isolated from soil. Int J Syst Evol Microbiol. 2006;56:997–1000. doi:10.1099/ijs.0.64119-0.
Yoon JH, Kang SJ, Lee JS, Oh TK. Flavobacterium terrigena sp. nov., isolated from soil. Int J Syst Evol Microbiol. 2007;57:947–50. doi:10.1099/ijs.0.64776-0.
Bernardet JF, Bowman JP. Genus I. Flavobacterium Bergey et al. 1923. In: W. Whitman, editor. Bergey’s Manual of Systematic Bacteriology. 2nd edn, vol. 4. Baltimore: The Williams & Wilkins Co, Baltimore; 2011. p. 112–54.
Tekedar HC, Karsi A, Gillaspy AF, Dyer DW, Benton NR, Zaitshik J, et al. Genome sequence of the fish pathogen Flavobacterium columnare ATCC 49512. J Bacteriol. 2012;194(10):2763–4. doi:10.1128/JB.00281-12.
Duchaud E, Boussaha M, Loux V, Bernardet JF, Michel C, Kerouault B, et al. Complete genome sequence of the fish pathogen Flavobacterium psychrophilum. Nat Biotechnol. 2007;25(7):763–9. doi:10.1038/nbt1313.
Touchon M, Barbier P, Bernardet JF, Loux V, Vacherie B, Barbe V, et al. Complete genome sequence of the fish pathogen Flavobacterium branchiophilum. Appl Environ Microbiol. 2011;77(21):7656–62. doi:10.1128/AEM.05625-11.
Park M, Lu S, Ryu SH, Chung BS, Park W, Kim CJ, et al. Flavobacterium croceum sp. nov., isolated from activated sludge. Int J Syst Evol Microbiol. 2006;56:2443–7. doi:10.1099/ijs.0.64436-0.
Ryu SH, Park JH, Moon JC, Sung Y, Lee SS, Jeon CO. Flavobacterium resistens sp. nov., isolated from stream sediment. Int J Syst Evol Microbiol. 2008;58:2266–70. doi:10.1099/ijs.0.65720-0.
Sheu SY, Chiu TF, Young CC, Arun AB, Chen WM. Flavobacterium macrobrachii sp. nov., isolated from a freshwater shrimp culture pond. Int J Syst Evol Microbiol. 2011;61:1402–7. doi:10.1099/ijs.0.025403-0.
Xu M, Xin Y, Tian J, Dong K, Yu Y, Zhang J, et al. Flavobacterium sinopsychrotolerans sp. nov., isolated from a glacier. Int J Syst Evol Microbiol. 2011;61:20–4. doi:10.1099/ijs.0.014126-0.
Fu Y, Tang X, Lai Q, Zhang C, Zhong H, Li W, et al. Flavobacterium beibuense sp. nov., isolated from marine sediment. Int J Syst Evol Microbiol. 2011;61(Pt 1):205–9. doi:10.1099/ijs.0.018846-0.
Qu JH, Yuan HL, Li HF, Deng CP. Flavobacterium cauense sp. nov., isolated from sediment of a eutrophic lake. Int J Syst Evol Microbiol. 2009;59(Pt 11):2666–9. doi:10.1099/ijs.0.009688-0.
Ali Z, Cousin S, Frühling A, Brambilla E, Schumann P, Yang Y, et al. Flavobacterium rivuli sp. nov., Flavobacterium subsaxonicum sp. nov., Flavobacterium swingsii sp. nov. and Flavobacterium reichenbachii sp. nov., isolated from a hard water rivulet. Int J Syst Evol Microbiol. 2009;59(Pt 10):2610–7. doi:10.1099/ijs.0.008771-0.
Barbier P, Houel A, Loux V, Poulain J, Bernardet JF, Touchon M, et al. Complete genome sequence of Flavobacterium indicum GPSTA100-9 T, isolated from warm spring water. J Bacteriol. 2012;194(11):3024–5. doi:10.1128/JB.00420-12.
Van Trappen S, Vandecandelaere I, Mergaert J, Swings J. Flavobacterium degerlachei sp. nov., Flavobacterium frigoris sp. nov., and Flavobacterium micromati sp. nov., novel psychrophilic bacteria isolated from microbial mats in Antarctic lakes. Int J Syst Evol Microbiol. 2004;54(Pt 1):85–92. doi:10.1099/ijs.0.02857-0.
Kolton M, Green SJ, Harel YM, Sela N, Elad Y, Cytryn E. Draft genome sequence of Flavobacterium sp. strain F52, isolated from the rhizosphere of bell pepper (Capsicum annuum L. cv. Maccabi). J Bacteriol. 2012;194(19):5462–3. doi:10.1128/JB.01249-12.
Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30(12):2725–9. doi:10.1093/molbev/mst197.
Krampis K, Booth T, Chapman B, Tiwari B, Bicak M, Field D, et al. Cloud BioLinux: pre-configured and on-demand bioinformatics computing for the genomics community. BMC Bioinformatics. 2012;13:42. doi:10.1186/1471-2105-13-42.
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. doi:10.1038/nbt1360.
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:4576–9.
Gherna R, Woese CR. A partial phylogenetic analysis of the “flavobacter-bacteroides” phylum: basis for taxonomic restructuring. Syst Appl Microbiol. 1992;15:513–21. doi:10.1016/S0723-2020(11)80110-4.
Berbardet JF, Krieg NR, Staley JT, Brown DR, Hedlund BP, Paster BJ, et al. Class II. Flavobacteriia class. nov. In: Bergey’s Manual of Systematic Bacteriology, The Bacteroidetes, Spirochaetes, Tenericutes (Mollicutes), Acidobacteria, Fibrobacteres, Fusobacteria, Dictyglomi, Gemmatimonadetes, Lentisphaerae, Verrumicrobia, Chlamydiae, and Planctomycetes, vol. 4. 2nd ed. New York: Springer; 2011. p. 105.
Bernardet JF, Segers P, Vancanneyt M, Berthe F, Kersters K, Vandamme P. Cutting a Gordian knot: emended classification and description of the genus Flavobacterium, emended description of the family Flavobacteriaceae, and proposal of Flavobacterium hydatis nom. nov. (Basonym, Cytophaga aquatilis Strohl and Tait 1978). Int J Syst Bacteriol. 1996;46:128–48.
Bergey DH, Harrison FC, Breed RS, Hammer BW, Huntoon FM. Bergey’s Manual of Determinative Bacteriology. 1st ed. The Williams and Wilkins Co: Baltimore; 1923. p. 1–442.
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 Con-sortium. Nat Genet. 2000;25:25–9. doi:10.1038/75556.
Illumina official website [www.illumina.com.cn].
Xie Y, Wu G, Tang J, Luo R, Patterson J, Liu S, et al. SOAPdenovo-Trans: de novo transcriptome assembly with short RNA-Seq reads. Bioinformatics. 2014;30(12):1660–6. doi:10.1093/bioinformatics/btu077.
Petersen TN, Brunak S, von Heijne G, Nielsen H. SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat Methods. 2011;8:785–6. doi:10.1038/nmeth.1701.
Hirokawa T, Boon-Chieng S, Mitaku S. SOSUI: classification and secondary structure prediction system for membrane proteins. Bioinformatics. 1998;14:378–9. doi:10.1093/bioinformatics/14.4.378.
Finn RD, Bateman A, Clements J, Coggill P, Eberhardt RY, Eddy SR, et al. The Pfam protein families database. Nucleic Acids Res. 2014;Database Issue 42(D):222–30.
GenBank database. www.ncbi.nlm.nih.gov/genbank. Accessed 2 Oct 2014.
Tatusov RL, Galperin MY, Natale DA, Koonin EV. The COG database: a tool for genome-scale analysis of protein functions and evolution. Nucleic Acids Res. 2000;28:33–6. doi:10.1093/nar/28.1.33.
Kanehisa M, Goto S, Sato Y, Kawashima M, Furumichi M, Tanabe M. Data, information, knowledge and principle: back to metabolism in KEGG. Nucleic Acids Res. 2014;42:D199–205. doi:10.1093/nar/gkt1076.
Li L, Stoeckert Jr CJ, Roos DS. OrthoMCL: identification of ortholog groups for eukaryotic genomes. Genome Res. 2003;13(9):2178–89. doi:10.1101/gr.1224503.
Milani CJ, Aziz RK, Locke JB, Dahesh S, Nizet V, Buchanan JT. The novel polysaccharide deacetylase homologue Pdi contributes to virulence of the aquatic pathogen Streptococcus iniae. Microbiol. 2010;156(Pt 2):543–54. doi:10.1099/mic.0.028365-0.
Zhang M, Han X, Liu H, Tian M, Ding C, Song J, et al. Inactivation of the ABC transporter ATPase gene in Brucella abortus strain 2308 attenuated the virulence of the bacteria. Vet Microbiol. 2013;164(3–4):322–9. doi:10.1016/j.vetmic.2013.02.017.
Saitou N, Nei M. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987;4:406–25.
This work was supported by the National Natural Science Foundation of China (31470226).
The authors declare that they have no competing interests.
ZZ performed genome analysis the data and wrote the draft manuscript. CC and HD helped to analyze the data. GW organized the study and revised the manuscript. ML performed the comparative genomics analysis and revised the manuscript. All authors read and approved the final manuscript.
About this article
Cite this article
Zeng, Z., Chen, C., Du, H. et al. High quality draft genomic sequence of Flavobacterium enshiense DK69T and comparison among Flavobacterium genomes. Stand in Genomic Sci 10, 92 (2015). https://doi.org/10.1186/s40793-015-0084-z
- Flavobacterium enshiense
- Comparative genomics
- Genome sequence