- Extended genome report
- Open Access
Potential probiotic-associated traits revealed from completed high quality genome sequence of Lactobacillus fermentum 3872
Standards in Genomic Sciencesvolume 12, Article number: 19 (2017)
The article provides an overview of the genomic features of Lactobacillus fermentum strain 3872. The genomic sequence reported here is one of three L. fermentum genome sequences completed to date. Comparative genomic analysis allowed the identification of genes that may be contributing to enhanced probiotic properties of this strain. In particular, the genes encoding putative mucus binding proteins, collagen-binding proteins, class III bacteriocin, as well as exopolysaccharide and prophage-related genes were identified. Genes related to bacterial aggregation and survival under harsh conditions in the gastrointestinal tract, along with the genes required for vitamin production were also found.
Probiotics are widely used for treatment of auto-immune conditions including allergic reactions, as well as metabolic disorders and are being applied as alternatives or additives to antibiotic treatment [1,2,3]. Probiotics may provide a beneficial effect by modulating the host immune system, via the release of antimicrobial substances, or through competitive exclusion of pathogenic bacteria . Various bacteria belonging to the Lactobacillus genus (including L. fermentum ) are commonly used as probiotics . The efficacy of these bacteria is not only species-specific, but also varies between the strains of the same species. Lactobacillus bacteria have a generally accepted as safe status. They are commonly found in various food products and are a part of the normal flora in animals and humans . However, some lactobacilli have been found to lower the intestinal barrier in vitro . L. fermentum 3872 has been patented in Russia along with a consortium of other Lactobacilli relating to their antimicrobial and probiotic uses . L. fermentum 3872 was sequenced in order to determine molecular modes of actions that may potentially be used against pathogenic bacteria that live in the same habitat as strain 3872, along with genes relating to its ability to survive harsh conditions of the GIT. Genomic data relating to the microflora of humans are also important for better understanding the role these bacteria play within its natural environment. With more high quality genomic data being made available a consortium of probiotics with similar modes of action may be utilised to effectively combat pathogenic bacteria. Currently, the genome sequence of L. fermentum strain 3872 reported here is one of only three complete genome sequences deposited in GenBank, with genome sequences of 16 more strains either being incomplete (draft) or containing ambiguities. For example, the genome of strain CECT 5716 (GenBank accession number CP002033) is shown in the GenBank as ‘complete’ and circular despite having a large number of ambiguities in the sequence. The aim of this study was to determine and characterise a complete genome sequence of this microorganism and to identify its specific genetic features.
Classification and features
Lactobacillus fermentum 3872 is a Gram-positive, rod-shaped (Fig. 1), facultative anaerobic bacteria  (Table 1). The strain is deposited under accession number VKM B-2793D at the All-Russian Collection of Microorganisms, Pushchino, Moscow Regions, Russia. Isolated from milk of a healthy woman. Identified as Lactobacillus fermentum in 2011 at the Institute of Engineering Immunology, Lyubuchany, Chekhov District, Moscow Regions, Russia. When grown in MRS agar L. fermentum 3872 forms medium sized, white colonies, that are round, smooth, and convex . L. fermentum 3872 was isolated from the milk of a healthy human female and has been found in infant and mother fecal matter along with vaginal secretions, indicating the strains ability to be present in different human ecological habitats . The bacterium has shown to be resistant to gastric and intestinal stresses, have high adhesion to human HeLa and buccal cells and has the ability to produce hydrogen peroxide and lactic acid, the release of which can be damaging to pathogenic bacteria . L. fermentum 3872 when present with a mixture of probiotics has been found to be a promising tool for the treatment of mastitis . L. fermentum 3872 belongs to the phylum firmicutes, among the circular genome sequences of L. fermentum the genome of strain 3872 appears to be most closely related to L. fermentum F6 (Fig. 2).
Genome sequencing information
Genome project history
Determination of a draft genome sequence of L. fermentum 3872 allowed the identification of a number of genes that may potentially be involved in probiotic activity, including a gene encoding a collagen-binding protein . The latter was subsequently found to be located on plasmid pLF3872, the sequence of which was reported in 2015 . In addition to the cbp gene, this plasmid, also contained a number of conjugation-related genes, as well as two toxin-antitoxin gene pairs required for stable maintenance of the plasmid within the bacterial cell . The current article conducts a detailed analysis of the recently completed chromosomal sequence of L. fermentum 3872, the assembly is of high quality due to the use of a hybrid sequencing approach along with a physical map of the genome described below. The article also conducts comparative analysis with other completed genome sequences belonging to the same species in order to determine targets for future probiotic experiments.
Growth conditions and genomic DNA preparation
L. fermentum 3872 was grown at 37 °C overnight on MRS agar plates under anaerobic conditions. DNA was isolated using Gentra Puregene Yeast/Bact Kit (Qiagen). For IonTorrent sequencing the NanoView photometer result indicated DNA concentration of 347ug/ul with DNA quality of A260/A280: 1.922 and A260/A230: 1.881. For Pacbio sequencing the NanoView photometer result for the extracted DNA was 314 ng/ul, A260/A280: 1.78 and A260/A230: 1.43, the Qubit DNA concentration result was 318 ng/ul. The DNA quality was also assessed by using agarose gel electrophoresis which indicated high concentration and good quality DNA (data not shown).
Genome sequencing and assembly
The complete circular genome sequence of L. fermentum 3872 was determined by employing a hybrid sequencing approach, including PacBio and IonTorrent PGM sequencing, as well as OpGen optical mapping. Long but high error and low coverage reads generated by PacBio were used as a scaffolding tool. PacBio sequencing was conducted using an RSII sequencing machine with P6/C4 sequencing chemistry and a single SMRT cell. HGAP and CELERA bioinformatics tools were used for the removal of low quality reads and generation of one large contig representing a circular 2.3 Mb chromosomal sequence of L. fermentum 3872. Short, but low error and high coverage reads produced by IonTorrent PGM using 314v2 chip and 400 bp kit were used for sequence verification and correction, which was essential for the low coverage areas. Three runs of IonTorrent sequencing were conducted producing 1,290,864 reads. Genome coverage by PacBio was 19.6 fold, as estimated by mapping of 4,902 reads between 500 and 21,671 bases long, with 4,871 of reads (99.37%) representing 99.07% nucleotides mapped. When combined with IonTorrent data, read mapping resulted in 413,661,861 bases (99.57%) mapped onto the assembly (2,330,492 nt) corresponding to 177.5 fold coverage (173.9 and 293.8 for chromosome and plasmid respectively). An optical map generated by OpGen technology was used for validation of the assembly, as well as for trimming and circularisation of the genome sequence. The genome information is summarised in Table 2 and Additional file 1: Table S1.
The genome sequence was annotated using PROKKA , BASys  and RAST  tools. In addition, the genome was annotated by NCBI GenBank annotation pipeline . Some annotation irregularities (such as e.g. truncated coding sequences) produced by these four annotation tools were identified and corrected using Geneious software .
The size of the L. fermentum 3872 genome (including the plasmid) is 2,330,492 bp. The G + C content of the circular chromosome (2,297,851 bp) is 55.6%. It contains 2328 genes, 2127 of which encode proteins and 128 are pseudogenes. There are 15 genes encoding rRNAs (23S, 16S and 5S) and 58 genes encoding tRNAs. The genome summary is presented in Tables 3, 4, 5.
Insights from the genome sequence
The circular view of the chromosome of L. fermentum 3872 was generated by using BRIGS software . The diagram indicates the leading (high G and low C region) and lagging (low G and high C region) strands of the L. fermentum 3872 chromosomal sequence (Fig. 3). Local GC skew deviations within the leading or lagging strand may indicate newly incorporated DNA, inversion or translocations . The diagram shows comparison of the genomic sequence of L. fermentum 3872 with those of L. fermentum CECT 5716, IFO 3956 and F6 strains.
L. fermentum 3872 contains genes required for the synthesis of such vitamins as B1, B2, B5, B7 and B9. These genes may play a crucial role in providing the natural hosts with essential vitamins. There are symporter encoding genes that allow the bacteria to survive acidic conditions of the stomach and thrive within the gastrointestinal tract. Among such genes are those encoding Na+/H+ (four copies), as well as gluconate/H+, sugar/H+, amino acid/H+, and glutamate/H+ symporters.
Survival of lactic acid bacteria within the gut is dependent on sugar metabolism and amino acid decarboxylation/deamination assisting in maintaining optimal pH levels . Among relevant genes of L. fermentum 3872 are those involved in arginine and proline metabolism (27 genes). There are also 14 genes involved in glutathione metabolism, which in Lactobacillus salivarius was found to be required for acid stress response . There is a gene encoding dTDP-glucose 4,6-dehydratase (Locus tag: N573_RS00605). In Lactobacillus plantarum this protein was found to be associated with gastric acid tolerance . There is a gene encoding Undecaprenyl-diphosphatase (EC 18.104.22.168) (Locus tag: N573_RS09665) with a possible role in bacitracin resistance by similarity to E. coli producing a similar protein . In other bacteria, such as Lactobacillus rhamnosus , the genes encoding DnaK (L. ferementum 3872 Locus tag: N573_RS04975) and GroEL (L. ferementum 3872 Locus tag: N573_RS01895) are known to play a role in heat and hyperosmotic shock tolerance. In addition, in Lactobacillus plantarum both genes are also implicated in mucin binding  potentially inhibiting adherence of pathogenic bacteria to the mucus layer. Furthermore, the GroEL of Lactobacillus johnsonii La1 was found to be a cell surface located protein capable of inducing aggregation of a gastric pathogen Helicobacter pylori in vitro . A gene (Locus tag: N573_RS03470) encoding a protein similar to Lactobacillus johnsonii La1 Translational Elongation Factor involved in bacterial adhesion to host cells , was also found. By similarity to function of similar genes found in Lactobacillus plantarum , L. fermentum 3872 genes encoding D-Lactate dehydrogenase (Locus tag: N573_RS11010) and 6-phosphogluconate dehydrogenase (Locus tag: N573_RS10960) are likely to promote bacterial adhesion to mucin and intestinal epithelial cells. There is a number of genes (e.g. loci N573_RS00495, N573_RS00500 and N573_RS00505, located to the same gene cluster) potentially involved in the biosynthesis of exopolysaccharides, which in other lactic acid bacteria were found to be important for bacterial survival and protection from toxic compounds .
Comparison of the complete chromosomal sequences of L. fermentum using LASTZ software  revealed a unique region of the L. fermentum 3872 genome (between positions 748,875 bp and 919,330 bp) (Fig. 4) This region contains genes encoding hypothetical proteins, enterolysin A (835,633 bp-836,847 bp) and ‘CAAX amino terminal protease self-immunity’ (838,683 bp-839,366 bp) protein, suggesting the bacterial ability to produce a bacteriocin. This was confirmed by running BAGEL3 bacteriocin prediction software , which identified a region (830,634 bp-840,633 bp) responsible for the biosynthesis of class III bacteriocin (Fig. 5d). No similarities were found for this region when using NCBI BlastN and the non-redundant database.
The region between 1,564,375 bp and 1,603,857 bp of the L. fermentum 3872 genome sequence contains inversions of respective parts of the genomes of L. fermentum strains F6, CECT 5716 and IFO 3956. This region also contains some prophage-related genes not found in the genomes of strains used for comparison. The region between 1,829,274 bp and 1,857,186 bp has a counterpart in L. gasseri ATCC 33323 (GenBank accession: CP000413) genome and may have been acquired via horizontal gene transfer (data not shown). The region between 2,212,692 bp and 2,237,160 bp has no matching sequences in the genomes of L. fermentum strains F6, CECT 5716 and IFO 3956, and contains conjugation and peptidoglycan hydrolase genes. NCBI BlastN analysis using the non-redundant database revealed high similarities to plasmid sequences, particularly with plasmid pPECL-5 from Pediococcus claussenii ATCC BAA-344 (e-value 0.0, query cover 55%). The other parts of this region contain the genes encoding transposases and an internalin J-like protein (InlJ, locus tag: N573_011130), containing an MucBP (mucin binding protein) domain [28, 29].
The genome of L. fermentum 3872 contains putative mucus binding protein-encoding gene also present in the genomes of strains F6, IFO 3956 and CECT 5716, but not in any other Lactobacillus genomes sequenced to date. Moreover, a gene, encoding a partial collagen-binding protein (Locus tag: N573_000435) is also found. This protein contains an LPXTG_anchor domain and a single B domain, but lacks the collagen-binding A domain . The gene encoding this protein was not found in any other L. fermentum strain.
The L. fermentum 3872 genome also contains an aggregation substance precursor protein encoding gene (Locus tag: N573_004020). The gene may potentially contribute to bacterial adhesion and aggregation . There are a number of exopolysaccharide production-related genes. In particular, epsH (Locus tag: N573_008790) predicted to be involved in biofilm formation, and may also contribute to protection against colitis . Remarkably, neither of these two genes (Locus tags: N573_004020, N573_008790) are present in the genomes of the strains used for comparison. An enolase encoding gene (Locus tag: N573_002185) present in L. fermentum 3872 may promote bacterial adhesion to collagen .
Comparative analysis of the genomes of L. fermentum strains 3872, F6 and IFO 3956 using Spine/AGent Pan-Core genome analysis tools with default parameters  allowed the identification of 428 unique ORFs of the L. fermentum 3872 genome with further 1650 ORFs representing core genes. One hundred and forty eight of the unique ORFs encode hypothetical proteins, with the other genes representing mobile elements, CRISPR-related and those involved in conjugal transfer. Among other genes were those encoding ABC transporters and those involved in bacteriocin biosynthesis, heavy metal resistance, and prophage-related genes.
PHAST software  allowed the identification of four prophage related regions (Fig. 5), each containing a phage attachment (ATT) site. A 34.5 kb region between 550,236 bp and 584,763 bp includes a number of genes encoding phage tail proteins, as well as transposases and integrases (Fig. 5a). Another 32 kb region (886,091 bp - 918,126 bp) also contains transposase, terminase and integrase encoding genes (Fig. 5d). A 39.4 kb region between 1,564,361 bp and 1,603,857 bp contains genes related to the biosynthesis of tail and head proteins, a protease, portal protein, terminase and integrase (Fig. 5b). A 30.2 kb region between 1,826,924 bp and 1,857,190 bp contains genes encoding a transposase, terminase, portal protein, capsid, head and recombinase. This region also contains an additional gene annotated as mucBP (Locus tag: N573_RS03620), which encodes amino acid protein containing 17 MucBP binding domain repeats. However, because of the absence of a cell wall anchor domain required for attachment, it is unlikely that this protein plays a role in adhesion (Fig. 5c). In addition, there are prophage-related genes (not identified by PHAST) adjacent to the bacteriocin encoding region. The prophage-related regions 550,236 bp - 584,763 bp, 1,564,361 bp - 1,603,857 bp and 1,826,924 bp - 1,857,190 have similarities in completely sequenced genomes of the species (Fig. 5a-b), whilst region 749,875 bp - 919,330 bp (containing prophage-related genes between 826,924 bp and 857,190 bp) is unique for strain 3872 (Fig. 5d).
Completion of the genome sequence of L. fermentum 3872 allowed the identification of various features that may contribute to probiotic properties of this bacterium, in addition to the already described CBP-encoding gene carried by the pLF3872 plasmid [9, 10]. Among these is a novel putative bacteriocin-encoding gene not found in any other genomes sequenced to date. Since a gene encoding a putative mucus-binding protein (Locus tag: N573_RS03620) suggests leucine as start codon, it remains to be verified whether the protein is actually expressed. There is a number of other genes (shared with other lactic acid bacteria) potentially required for bacterial attachment to host cells, survival in unfavourable conditions and resistance to toxic compounds. Despite the presence of some conserved features shared by all L. fermentum genomes, and a very high similarity between their sequences, the genome of strain 3872 has a large number of unique genes such as epsH, and a putative adhesion gene, inlJ. A gene that may promote bacterial aggregation has also been found. These genes could be a subject of further investigation. Conservation within the genome of as many as four large prophage-related gene clusters may also contribute to the lifestyle and probiotic properties of this microorganism. In particular, some bacteriocins produced by other bacteria resemble components of bacteriophages, and are encoded by prophage regions of the chromosomes . The bacteriophage-related gene products are being studied as alternatives to antibiotics due to their high potency and specificity, and thus may be of interest for further investigation [35, 37]. As L. fermentum 3872 was isolated from the milk of a healthy human female, the presence of multiple vitamin synthesising genes, along with the genes allowing the bacterium to thrive in the gut environment, would make L. fermentum an ideal candidate for probiotic studies. The ability of these bacteria to produce various adhesins may allow competitive exclusion of pathogenic microorganisms employing similar mechanisms of adhesion and interacting with the same host cell receptors. The presence of a novel bacteriocin-encoding gene may also contribute to beneficial properties of this strain.
Bacterial annotation system
Hierarchical genome assembly process
- MRS agar:
Agar developed by de Man, Rogosa and Sharpe
Rapid prokaryotic genome annotation tool
Bubnov RV, Spivak MY, Lazarenko LM, Bomba A, Boyko NV. Probiotics and immunity: provisional role for personalized diets and disease prevention. The EPMA Journal. 2015;6:14.
Yang G, Liu Z, Yang P. Treatment of allergic rhinitis with probiotics: An alternative approach. North Am J Med Sci. 2013;5:465–8.
Azcárate-Peril MA, Sikes M, Bruno-Bárcena JM. The intestinal microbiota, gastrointestinal environment and colorectal cancer: a putative role for probiotics in prevention of colorectal cancer? Am J Physiol. 2011;301:G401–24.
Collado MC, Isolauri E, Salminen S, Sanz Y. The impact of probiotic on gut health. Curr Drug Metab. 2009;10:68–78.
Jensen H, Dromtorp SM, Axelsson L, Grimmer S. Immunomodulation of Monocytes by Probiotic and Selected Lactic Acid Bacteria. Probiotics and Antimicrobial Proteins. 2015;7:14–23.
Tajabadi N, Mardan M, Saari N, Mustafa S, Bahreini R, Mohd Yazid AM. Identification of Lactobacillus plantarum, Lactobacillus pentosus and Lactobacillus fermentum from honey stomach of honeybee. Brazilian J Microbiol. 2013;44:717–22.
Anderson RC, Young W, Clerens S, Cookson AL, McCann MJ, Armstrong KM, et al. Human Oral Isolate Lactobacillus fermentum AGR1487 Reduces Intestinal Barrier Integrity by Increasing the Turnover of Microtubules in Caco-2 Cells. PLoS One. 2013;8:e78774.
Abramov VM, Khlebnikov VS, Pchelintsev SJ, Kosarev IV, Karlyshev AV, Vasilenko RN, Melnikov VG. Strain Lactobacillus fermentum having broad spectrum of antagonistic activity and probiotic lactobacterium consortium for manufacturing bacterial preparations” RU 2528862 C1, (20.09.14, Russia). Patent. Application Number: 2013118084/10, Application Date: 19.04.2013, Publication Number: 0002528862, Publication Date: 20.09.2014.
Karlyshev AV, Raju K, Abramov VM. Draft Genome Sequence of Lactobacillus fermentum Strain 3872. Genome Announcements. 2013;doi:10.1128/genomeA.01006-13 .
Lehri B, Seddon AM, Karlyshev AV. Lactobacillus fermentum 3872 genome sequencing reveals plasmid and chromosomal genes potentially involved in a probiotic activity. FEMS Microbiol Lett. 2015;doi:10.1093/femsle/fnv068
Seemann T. Prokka: rapid prokaryotic genome annotation. Bioinformatics. 2014;doi:10.1093/bioinformatics/btu153
Van Domselaar GH, Stothard P, Shrivastava S, Cruz JA, Guo A, Dong X, et al. BASys: a web server for automated bacterial genome annotation. Nucleic Acids Res. 2005;33:W455–9.
Overbeek R, Olson R, Pusch GD, Olsen GJ, Davis JJ, Disz T, et al. The SEED and the Rapid Annotation of microbial genomes using Subsystems Technology (RAST). Nucleic Acids Res. 2014;42:D206–14.
Angiuoli SV, Gussman A, Klimke W, Cochrane G, Field D, Garrity GM, et al. Toward an Online Repository of Standard Operating Procedures (SOPs) for (Meta) genomic Annotation. OMICS: J Integr Biol. 2008;12:137–41.
Kearse M, Moir R, Wilson A, Stones-Havas S, Cheung M, Sturrock S, et al. Geneious Basic: An integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics. 2012;28:1647–9.
Alikhan N, Petty NK, Ben Zakour NL, Beatson SA. BLAST Ring Image Generator (BRIG): simple prokaryote genome comparisons. BMC Genomics. 2011;12:1–10.
Grigoriev A. Analyzing genomes with cumulative skew diagrams. Nucleic Acids Res. 1998;26:2286–90.
Pessione E. Lactic acid bacteria contribution to gut microbiota complexity: lights and shadows. Front Cell Infect Microbiol. 2012;2:UNSP 86.
Lee K, Pi K, Kim EB, Rho B, Kang S, Lee HG, et al. Glutathione-mediated response to acid stress in the probiotic bacterium, Lactobacillus salivarius. Biotechnol Lett. 2010;32:969–72.
Hamon E, Horvatovich P, Marchioni E, Aoude-Werner D, Ennahar S. Investigation of potential markers of acid resistance in Lactobacillus plantarum by comparative proteomics. J Appl Microbiol. 2014;116:134–44.
El Ghachi M, Derbise A, Bouhss A, Mengin-Lecreulx D. Identification of multiple genes encoding membrane proteins with undecaprenyl pyrophosphate phosphatase (UppP) activity in Escherichia coli. J Biol Chem. 2005;280:18689–95.
Prasad J, McJarrow P, Gopal P. Heat and osmotic stress responses of probiotic Lactobacillus rhamnosus HN001 (DR20) in relation to viability after drying. Appl Environ Microbiol. 2003;69:917–25.
Izquierdo E, Horvatovich P, Marchioni E, Aoude-Werner D, Sanz Y, Ennahar S. 2-DE and MS analysis of key proteins in the adhesion of Lactobacillus plantarum, a first step toward early selection of probiotics based on bacterial biomarkers. Electrophoresis. 2009;30:949–56.
Bergonzelli G, Granato D, Pridmore R, Marvin-Guy L, Donnicola D, Corthesy-Theulaz I. GroEL of Lactobacillus johnsonii La1 (NCC 533) is cell surface associated: Potential role in interactions with the host and the gastric pathogen Helicobacter pylori. Infect Immun. 2006;74:425–34.
Granato D, Bergonzelli G, Pridmore R, Marvin L, Rouvet M, Corthesy-Theulaz I. Cell surface-associated elongation factor Tu mediates the attachment of Lactobacillus johnsonii NCC533 (La1) to human intestinal cells and mucins. Infect Immun. 2004;72:2160–9.
Harris RS. Improved pairwise alignment of genomic DNA. [PhD dissertation]. The Pennsylvania State University; 2007.
Van Heel AJ, de Jong A, Montalban-Lopez M, Kok J, Kuipers OP. BAGEL3: automated identification of genes encoding bacteriocins and (non-)bactericidal posttranslationally modified peptides. Nucleic Acids Res. 2013;41:W448–53.
Sabet C, Lecuit M, Cabanes D, Cossart P, Bierne H. LPXTG Protein InlJ, a Newly Identified Internalin Involved in Listeria monocytogenes Virulence. Infect Immun. 2005;73:6912–22.
Juge N. Microbial adhesins to gastrointestinal mucus. Trends Microbiol. 2012;20:30–9.
Deivanayagam CC, Rich RL, Carson M, Owens RT, Danthuluri S, Bice T, et al. Novel fold and assembly of the repetitive B region of the Staphylococcus aureus collagen-binding surface protein. Structure. 2000;8:67–78.
Suessmuth SD, Muscholl-Silberhorn A, Wirth R, Susa M, Marre R, Rozdzinski E. Aggregation Substance Promotes Adherence, Phagocytosis, and Intracellular Survival of Enterococcus faecalis within Human Macrophages and Suppresses Respiratory Burst. Infect Immun. 2000;68:4900–6.
Jones SE, Paynich ML, Kearns DB, Knight KL. Protection from intestinal inflammation by bacterial exopolysaccharides. J Immunol (Baltimore, Md: 1950). 2014;192:4813–20.
Salzillo M, Vastano V, Capri U, Muscariello L, Sacco M, Marasco R. Identification and characterization of enolase as a collagen-binding protein in Lactobacillus plantarum. J Basic Microbiol. 2015;55:890–7.
Ozer E, Allen J, Hauser A. Characterization of the core and accessory genomes of Pseudomonas aeruginosa using bioinformatic tools Spine and AGEnt. BMC Genomics. 2014;15:737.
Zhou Y, Liang Y, Lynch KH, Dennis JJ, Wishart DS. PHAST: A Fast Phage Search Tool. Nucleic Acids Research. 2011;doi:10.1093/nar/gkr485
Nakayama K, Takashima K, Ishihara H, Shinomiya T, Kageyama M, Kanaya S, et al. The R-type pyocin of Pseudomonas aeruginosa is related to P2 phage, and the F-type is related to lambda phage. Mol Microbiol. 2000;38:213–31.
Coates ARM, Hu Y. Novel approaches to developing new antibiotics for bacterial infections. Br J Pharmacol. 2007;152:1147–54.
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.
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.
Gibbons NE, Murray RGE. Proposals Concerning the Higher Taxa of Bacteria. Int J Syst Bacteriol. 1978;28:1–6.
Euzéby: Validation List No. 132. List of new names and new combinations previously effectively, but not validly, published. Int J Syst Evol Microbiol. 2010;60:469–72.
Ludwig W, Schleifer K-H, Whitman WB. Class I. Bacilli class nov. In: Bergey's Manual of Systematic Bacteriology. 2009;3:19–20.
Winslow CEA, Broadhurst J, Buchanan RE, Krumwiede C, Rogers LA, Smith GH. The Families and Genera of the Bacteria: Preliminary Report of the Committee of the Society of American Bacteriologists on Characterization and Classification of Bacterial Types. J Bacteriol. 1917;2:505–66.
Skerman VBD, McGowan V, Sneath PHA. Approved Lists of Bacterial Names. Int J Syst Bacteriol. 1980;30:225–420.
Beijerinck MW. Sur les ferments lactiques de l'industrie. Archives Néerlandaises des Sciences Exactes et Naturelles. 1901;6:212–43.
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:25–9.
Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32(5):1792–7.
Huelsenbeck JP, Ronquist F. MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics. 2001;17:754–5.
This work was not supported by any external funding.
BL participated in experimental design, genome sequencing and analysis, and in drafting the manuscript. AMS participated in experimental design, discussion of this study and in manuscript preparation. AVK conceived the study, participated in experimental design, genome sequencing and analysis, preparation and finalizing the manuscript. All authors read and approved the manuscript.
The authors declare that they have no competing interests, and hereby confirm that one of the authors (Andrey Karlyshev) is a named inventor on a patent mentioned in the Introduction section.
Associated MIGS record. (DOC 72 kb)