Open Access

Veillonella, Firmicutes: Microbes disguised as Gram negatives

  • Tammi Vesth
  • , Aslı Ozen,
  • , Sandra C. Andersen
  • , Rolf Sommer Kaas
  • , Oksana Lukjancenko
  • , Jon Bohlin
  • , Intawat Nookaew
  • , Trudy M. Wassenaar
  • and David W. Ussery
Corresponding author

DOI: 10.4056/sigs.2981345

Received: 15 December 2013

Accepted: 15 December 2013

Published: 20 December 2013


The Firmicutes represent a major component of the intestinal microflora. The intestinal Firmicutes are a large, diverse group of organisms, many of which are poorly characterized due to their anaerobic growth requirements. Although most Firmicutes are Gram positive, members of the class Negativicutes, including the genus Veillonella, stain Gram negative. Veillonella are among the most abundant organisms of the oral and intestinal microflora of animals and humans, in spite of being strict anaerobes. In this work, the genomes of 24 Negativicutes, including eight Veillonella spp., are compared to 20 other Firmicutes genomes; a further 101 prokaryotic genomes were included, covering 26 phyla. Thus a total of 145 prokaryotic genomes were analyzed by various methods to investigate the apparent conflict of the Veillonella Gram stain and their taxonomic position within the Firmicutes. Comparison of the genome sequences confirms that the Negativicutes are distantly related to Clostridium spp., based on 16S rRNA, complete genomic DNA sequences, and a consensus tree based on conserved proteins. The genus Veillonella is relatively homogeneous: inter-genus pair-wise comparison identifies at least 1,350 shared proteins, although less than half of these are found in any given Clostridium genome. Only 27 proteins are found conserved in all analyzed prokaryote genomes. Veillonella has distinct metabolic properties, and significant similarities to genomes of Proteobacteria are not detected, with the exception of a shared LPS biosynthesis pathway. The clade within the class Negativicutes to which the genus Veillonella belongs exhibits unique properties, most of which are in common with Gram-positives and some with Gram negatives. They are only distantly related to Clostridia, but are even less closely related to Gram-negative species. Though the Negativicutes stain Gram-negative and possess two membranes, the genome and proteome analysis presented here confirm their place within the (mainly) Gram positive phylum of the Firmicutes. Further studies are required to unveil the evolutionary history of the Veillonella and other Negativicutes.


The genus Veillonella, belonging to Negativicutes, consists of anaerobic, non-fermentative, Gram-negative cocci, that are normally observed in pairs or short chains, and are non-sporulating and non-motile [1]. Veillonella spp. are abundant in the human microbiome and are found in the oral, respiratory, intestinal and genitourinary flora of humans and animals; they can make up as much as 10% of the bacterial community initially colonizing the enamel [2] and are found throughout the entire oral cavity [3], especially on the tongue dorsum and in saliva [4]. The importance of Veillonella spp. in human infections is uncertain, and they are generally considered to be of low virulence. Veillonella form biofilms, often with Streptococcus spp., and species of these genera have been found to be more abundant in the oral microflora of people with poor oral health [5]. Studies have shown that during formation of early dental plaque, the fraction of Veillonella spp. changes in mixed-microbial colonies with streptococci [6]. Thus, Veillonella spp. may play a role in caries formation as they utilize the lactic acid produced by the organisms conducive to caries [7]. Veillonella are also among the most common anaerobic species reported from pulmonary samples and are frequently recovered from cystic fibrosis cases [8]. The organisms are also abundant in the human gut flora, where their numbers were found to be higher in children with type I diabetes compared to healthy controls [9]. Currently, 12 species of Veillonella have been characterized [10,11] including V. parvula, V. atypica and V. dispar, which are found in the human oral cavity.

The Negativicutes are the only diderm (literally 'two skins') members of the phylum Firmicutes as they possess an inner and an outer membrane. Their placement within the Firmicutes has been widely accepted, and has been confirmed by 16S rRNA analysis [12]. However, their genomes have not been analyzed in detail to confirm their taxonomic position. This work presents a broad analysis of the Negativicutes with focus on the Veillonella spp. using comparative microbial genomics. A total of 24 genomes from the Negativicutes were compared to 121 genomes covering most of the taxonomic span of sequenced bacterial genomes. We investigated how the Negativicutes genomes compared to other bacterial genomes using three different and complementary approaches: 1) phylogenetic trees to visualize the relative distance of the Negativicutes genomes to other genomes; 2) amino acid composition, nucleotide tetramer frequency and metabolism analysis using 2-D clustering and heatmaps to compare genomes; and 3) proteomic comparison across the Negativicutes genomes.

Materials and Methods

Genome sequences used for analysis

The set of 145 genomes included in this study (24 Negativicutes genomes and 121 other prokaryotic genomes covering 26 phyla) are listed in Table 1.

Table 1

Genomes used in this study


   Name of organism and strain

   Strain designation

  Type strain?

   NCBI Taxon ID

    NCBI Project ID


   Acidobacterium capsulatum

   ATCC 51196





   “Korebacter versatiles”

   Ellin 345




   “Solibacter usitatus”





   Bifidobacterium bifidum






   Catenulispora acidiphila

   ID139908, DSM 44928





   Corynebacterium pseudotuberculosis






   Segniliparus rugosus

   ATCC BAA-974





   Streptomyces bingchenggensis


  Name not validly published




   Tropheryma whipplei






   Persephonella marina






   Sulfurihydrogenibium sp.


  No type strain available




   Thermocrinis albus

   HI 11/12, DSM 14484





   Bacteroides thetaiotaomicron






   Candidatus Sulcia muelleri





   Chitinophaga pinensis

   UQM 2034, DSM 2588





   Paludibacter propionicigenes

   WB4, DSM 17365





   Protochlamydia amoebophila






   Chlamydia trachomatis






   Chlamydophila pneumoniae






   Waddlia chondrophila

   WSU 86-1044





   Chlorobium chlorochromatii”


  Name not validly published




   Chlorobium tepidum






   Chloroflexus aggregans

   DSM 9485





   Dehalococcoides sp






   Herpetosiphon aurantiacus

   ATCC 23779





   Roseiflexus sp.


  No type strain available




   Anabaena variabilis 3

   ATCC 2941





   Cyanothece sp.

   PCC 7822





   Prochlorococcus marinus






   Synechocystis sp.






   Calditerrivibrio nitroreducens

   Yu37-1, DSM 19672





   Deferribacter desulfuricans

   SSM1, DSM 14783





   Denitrovibrio acetiphilus

   N2460, DSM 12809





   Oceanithermus profundus

   506, DSM 14977





   Thermus thermophilus






   Truepera radiovictrix

   RQ-24, DSM 17093





   Dictyoglomus turgidum

   DSM 6724





   Elusimicrobium minutum

   Pei 191





   Fibrobacter succinogenes






   Acetohalobium arabaticum

   Z-7288, DSM 5501





   Acidaminococcus fermentans

   VR4, DSM 20731





   Acidaminococcus sp.


  No type strain available




   Alkaliphilus oremlandii






   Bacillus subtilis subsp. subtilis






   Clostridium botulinum

   F Langeland





   Clostridium cellulolyticum






   Clostridium difficile

   630 (epidemic type X)





   Desulfotomaculum reducens”


  Name not validly published




   Dialister invisus

   DSM 15470





   Dialister micraerophilus

   Oral Taxon 843 DSM 19965





   Dialister micraerophilus






   Enterococcus faecalis






   Eubacterium cylindroides






   Eubacterium rectale

   A1-86, DSM 17629





   Exiguobacterium sibiricum






   Geobacillus kaustophilus






   Lactococcus lactis    cremoris






   Lysinibacillus sphaericus






   Megamonas hypermegale






   Megasphaera genomo sp.

   type 1 28L

  No type strain available




   Megasphaera micronuciformis






   Mitsuokella multacida

   A 405-1, DSM 20544





   Paenibacillus sp.






   Phascolarctobacterium sp.

   YIT 12067





   Selenomonas artemidis






   Selenomonas flueggei

   ATCC 43531





   Selenomonas noxia

   ATCC 43541





   Selenomonas sp.

   Oral Taxon 137 F0430

  No type strain available




   Selenomonas sp.

   Oral Taxon 149 67H29BP

  No type strain available




   Selenomonas sputigena

   DSM 20758





   Staphylococcus aureus aureus






   Streptococcus pneumoniae






   Thermoanaerobacter sp.


  Name not validly published




   Thermosinus carboxydivorans






   Turicibacter sp.

   PC909 702450 42765



   Veillonella atypica






   Veillonella atypica






   Veillonella dispar

   ATCC 17748





   Veillonella parvula

   ATCC 17745





   Veillonella parvula

   Te3, DSM 2008





   Veillonella sp.

   3 1 44

  Name not validly published




   Veillonella sp.

   6 1 27

  Name not validly published




   Veillonella sp.

   Oral Taxon 158 F0412

  Name not validly published




   Fusobacterium nucleatum nucleatum

   ATCC 25586





   Ilyobacter polytropus

   CuHBu1, DSM 2926





   Leptotrichia buccalis

   C-1013-b, DSM 1135





   Sebaldella termitidis

   NCTC 11300





   Streptobacillus moniliformis

   9901, DSM 12112





   Pirellula staleyi

   DSM 6068





   Planctomyces limnophilus

   Mu 290, DSM 3776





   Acinetobacter baumannii






   Alkalilimnicola ehrlichii






   Arcobacter nitrofigilis

   DSM 7299





   Burkholderia xenovorans

   (fungorum) LB400





   Campylobacter jejuni    doylei






   Candidatus Pelagibacter ubique

   SAR11 HTCC1062

  Name not validly published




   Candidatus Zinderia insecticola


  Name not validly published




   Cellvibrio japonicus






   Cupriavidus taiwanensis






   Escherichia coli

   K-12, MG1655





   Geobacter uraniireducens






   Hahella chejuensis

   KCTC 2396





   Haliangium ochraceum

   SMP-2, DSM 14365





   Helicobacter pylori






   Lawsonia intracellularis






   Magnetococcus sp.


  Name not validly published




   Methylobacterium nodulans






   Neisseria meningitidis






   Neorickettsia sennetsu






   Nitrosomonas eutropha

   C91 (C71)





   Photorhabdus luminescens laumondii






   Polynucleobacter necessarius






   Pseudomonas aeruginosa






   Pseudomonas fluorescens






   Pseudomonas stutzeri






   Salmonella enterica enterica

   PT4 P125109





   Shewanella oneidensis






   Sorangium cellulosum

   So ce56





   Stigmatella aurantiaca

   DW4 /3-1





   Sulfurospirillum deleyianum

   5175, DSM 6946





   Vibrio cholerae






   Borrelia turicatae






   Brachyspira murdochii

   56-150, DSM 12563





   Leptospira interrogans

   lai 56601





   Thermanaerovibrio acidaminovorans

   Su883, DSM 6589





   Acholeplasma laidlawii






   Candidatus Phytoplasma asteris

   yellows witches'-broom AY-WB 322098

  Name not validly published



   Candidatus Phytoplasma mali


  Name not validly published




   Mycoplasma genitalium






   Mycoplasma pneumoniae






   Ureaplasma parvum

   sv 3, ATCC 27815





   Fervidobacterium nodosum






   Kosmotoga olearia

   TBF 19.5.1





   Petrotoga mobilis






   Thermotoga naphthophila






   Akkermansia muciniphila

   ATCC BAA-835





   Opitutus terrae





   Sulfolobus solfataricus





   Thermosphaera aggregans

   M11TL, DSM 11486





   Halogeometricum borinquense

   PR3, DSM 11551





   Methanocella sp.


  Name not validly published




   Methanothermus fervidus

   V24S, DSM 2088





   Candidatus Korarchaeum cryptofilum


  Name not validly published




   “Nanoarchaeum equitans”


  Name not validly published



16S rRNA tree

For this analysis, 16S rRNA sequences were predicted from the whole genome sequences of the selected organisms, using the RNAmmer algorithm [13]. These sequences were aligned using the MAFFT program, with the iterative refinement algorithm using maximum iteration (1000) and default parameters for gap penalties [14]. A distance tree was constructed using MEGA5 [15] with the Neighbor-joining algorithm [16] and 1,000 bootstrap re-samplings. The taxa in the resulting tree were collapsed to phyla, except for the Negativicutes.

Composition Vector Tree (CV)

A Composition Vector Tree was constructed based on protein sequences of the 145 selected genomes using a webserver (available at Web Site) with the K parameter set at 6 [17]. The outcome from the program is a distance matrix based on amino acid sequence comparisons, which is then used to generate a phylogenetic tree with the neighbor-joining method. In the shown tree, the outgroup chosen was Methanothermus fervidus (an Archaea). After tree visualization with MEGA5, branches were collapsed wherever possible with the exception of the Negativicutes branch, which remained expanded.

Consensus tree of conserved genes

Using the list of universally conserved core genes, previously identified by Ciccarelli et al. [18], and an implementation of BLAST, a set of genes that was shared among all 145 genomes was identified. Proteins that had no match in at least one genome or showed poor E-value were eliminated. The 27 conserved core genes were extracted (Table 1) and a multiple alignment was produced using MUSCLE software [19]. A set of phylogenetic trees was constructed by PAUP [20] and a best-fit consensus tree was generated using Phylogeny Inference package (PHYLIP) as described elsewhere [21]. Bootstrap values were found after 27 re-samplings, which is equal to the number of gene families conserved in all the analyzed genomes.

DNA tetramer analysis and amino acid usage

A tetramer frequency heatmap was constructed from the observed ratios of tetra-nucleotide frequencies divided by estimated tetra-nucleotide frequencies for each genome [22]. The estimated tetra-nucleotides were computed from the genomes' base composition. The ratio of observed over expected frequency was used for hierarchical clustering using complete linkage and Euclidean distance, which was subsequently performed with respect to both strain and tetramer frequencies.

The amino acid heatmap is based on frequencies of deduced proteomic amino acids from each genome normalized with respect to the total number of amino acids in each genome. The amino acid frequencies for each genome were clustered using complete linkage and Euclidean distance with respect to both genomes and amino acids. The heatmap was made using the R package ggplot2 [23].

Comparison of metabolism potential

The protein sequences of Kyoto Encyclopedia of Genes and Genomes (KEGG) orthology categories [24] were downloaded and only the Bacterial sequences were considered. The Hidden Markov model (HMM) of each ortholog was generated using HMMER version 3 [25] based on the multiple alignment of each orthologous set of KEGG proteins, using MUSCLE software [19]. The 145 proteomes were queried against the HMMs to infer their ontology. A cutoff of 1×10−30 was used for statistical significance. A heatmap of each pathway and process derived from the database KEGG was illustrated based on normalized abundance of the enzymes present in each pathway. The heatmap and hierarchical clustering were performed in the software R [23].

Construction of BLAST matrix and proteome comparison

Reciprocal BLAST was performed between each genome pair. The program blastall version 2.2.25 was used for BLAST implementation using default settings (BLASTp, E-value set to 1×10−5 for non-homologs and 1×10−8 for homologs, without filtering). A hit was considered significant at a BLAST cutoff of 95% identity and 95% coverage (of the longest gene in comparison). The number of hits was then given as a percentage of the genes in the column representing the corresponding genome. The diagonal designates internal homologs, computed by blasting each genome with itself. To avoid including identical genes, the second highest scoring hits were used. Furthermore, we also performed homology reduction of the diagonal hits, using an implementation of the Hobohm algorithm [26].


Twenty-four Negativicutes genomes were compared to 121 other prokaryotic genomes covering 22 Bacterial and 4 Archaeal phyla. When available, at least two genomes were included for every phylum. The first analysis presented here is based on 16S rRNA alignments. A single 16S rRNA gene was extracted from each of the genomes and an alignment was produced spanning the maximum length of the gene. A phylogenetic tree was constructed based on this alignment, as shown in Figure 1. With the exception of the Negativicutes, branches of the tree were collapsed in those cases where the analyzed species within a phylum clustered together. With the exception of some Firmicutes, the analyzed genomes cluster according to their phylum, although the Deferribacteres phylum is mixed with the Proteobacteria phyla, and two members of Proteobacteria are not positioned with other members of their phylum (Lawsonia intracellularis and Magnetococcus). That most phyla could be collapsed is consistent with the weight of 16S rRNA similarities in currently accepted taxonomic descriptions of prokaryotes. The Firmicutes, however, show less consistency. Although most of the analyzed Firmicutes cluster together, two species are separated from the Firmicutes branch (Eubacterium cylindroides and Thermoanaerobacter sp., both members of Clostridia). The Negativicutes are positioned within the Firmicutes cluster, and this part of the tree is expanded in the figure for clarity. As can be seen, phylogeny of the 16S rRNA gene provides good resolution between the different genera of the analyzed Negativicutes. All Veillonella spp. are clustered within one branch of the Negativicutes. The Acidaminococcaceae (to which Phascolarctobacterium spp. also belong) are placed within the cluster of the Veillonellaceae, in accordance with their current classification [27]. The Acidaminococcaceae used to be recognized as a separate family within the Negativicutes, just like the Veillonellaceae, and during preparation of this contribution these two families were presented as such in the Taxonomy database at NCBI. Of note is the relatively close relationship between Negativicutes and two Clostridium species (C. botulinum and C. cellulolyticum), which does not cluster with other members of the Clostridium genus (Figure 1). That genus displays a high degree of variation and re-classification of some of the members of this genus is in progress (see for example [27]). That two members of the Clostridia are even placed outside the Firmicutes phylum is an indication of 16S rRNA gene sequence heterogeneity within this class.

Figure 1

Phylogenetic neighbor-joining tree based on 16S rRNA genes extracted from 145 genomes (24 Negativicutes and 121 prokaryotic genomes representing 26 phyla). Bootstrap values of 50 and higher are indicated. With the exception of the Negativicutes, branches where all organisms belong to the same phyla are collapsed and named by the phyla they represent. The green shading indicates the position of Firmicutes. The collapsed branch of the Bacilli, marked (1), contains Turicibacter sanguinis, a Firmicutes member of the Erysipelotrichales as well as Bacilli members. An uncollapsed tree is included in the supplementary material.

Next, all protein-coding genes of the analyzed genomes were compared and a composition vector tree (CVtree) was produced, based on amino acid sequences (Figure 2). The topology of the resulting tree is generally in accordance with the 16S rRNA tree shown in the previous figure. As indicated by the collapsed branches, the CVtree grouped most genomes according to their known taxonomic phyla, although not all Spirochaetes cluster together. In contrast to the 16S rRNA tree, in this protein tree all the Firmicutes cluster together, and are distinct from other phyla. The Negativicutes genomes, nested within the Firmicutes, again have the Acidaminococcaceae placed within the Veillonellaceae, while all Veillonella spp. are found in one cluster. All Clostridia, this time divided into two collapsed branches, are positioned as the closest relatives to Negativicutes. It is of interest that among the closest relatives to Firmicutes, based on this analysis, are the Fusobacteria and the Elusimicrobia; these are atypical diderm bacteria that produce lipopolysaccharides [28]. However, the spirochete, Brachyspira murdochii, does not possess two membranes, but is nevertheless grouped with atypical diderms. On the other hand while the Synergistetes are atypical diderm bacteria, they are placed elsewhere in the tree (Figure 2).

Figure 2

Phylogenetic tree based on composition vector analysis (CVtree) of all protein coding genes (amino acid sequences) derived from the analyzed genomes. Note that the branch lengths in this plot are artificial. The coloring is the same as in Figure 1 and branches have been collapsed. The Firmicutes branch Bacilli, marked (1), contains Turicibacter sanguinis. An uncollapsed tree is included in the supplementary material.

A third analysis was based on a subset of proteins found conserved amongst all analyzed genomes. These conserved proteins were selected based on a protein BLAST (a cutoff of 50% identity and 50% coverage of the query length was used) and single linkage clustering. The analysis identified 29 genes that are shared among all 145 genomes [Table 2]. A consensus tree was constructed based on these 29 conserved proteins (Figure 3). The results confirm the global observations of the other two phylogenetic analyses: the Negativicutes cluster together and are most closely related to Clostridia (in this case the most closely related species are Desulfotomaculum reducens and Acetohalobium arabaticum). As before, the Acidaminococcaceae cluster together but within the Veillonellaceae. The position of Turicibacter sanguinis within the Bacilli group of Firmicutes is consistent with the other two trees but contrasts with its taxonomic description at NCBI as a member of the Erysipelotrichia.

Table 2

Universally conserved COGs


    Average length (aa)




    Predicted GTPase, probable translation factor



    Phenylalanine-tRNA synthethase alpha subunit



    Ribosomal protein S12



    Ribosomal protein S7



    Ribosomal protein S2



    Ribosomal protein L11



    Ribosomal protein L1



    Ribosomal protein L3



    Ribosomal protein L22



    Ribosomal protein S3



    Ribosomal protein L14



    Ribosomal protein L5



    Ribosomal protein S8



    Ribosomal protein L6P/L9E



    Ribosomal protein S5



    Ribosomal protein S11



    Ribosomal protein L13



    Ribosomal protein S9



    Seryl-tRNA synthetase



    Ribosomal protein S15P/S13E



    Ribosomal protein S17



    Ribosomal protein L16/L10E



    Ribosomal protein L15



    Preprotein translocase subunit SecY



    DNA-directed RNA polymerase, alpha subunit



    Ribosomal protein L18



    Leucyl-tRNA synthetase



    Ribosomal protein S4 and related proteins



    Metal-dependent proteases with chaperone activity

Figure 3

Consensus tree based on the phylogenetic trees of 27 genes conserved in all 145 genomes. The collapsed branch of the Bacilli, marked (1), contains Turicibacter sanguinis. An uncollapsed tree is available as a supplemental figure.

In conclusion, based on three independent phylogenetic analyses, the closest relatives to the Negativicutes seem to be the Clostridiaceae. The observed clustering of species within the Negativicutes is consistent with their assigned taxonomy. Furthermore, these analyses show that Veillonella spp. form a distinct branch, most different from the other Negativicutes, while the recent change of status of the Acidaminococcaceae (they are no longer a separate family) is confirmed by these analyses.

Apart from comparing proteins and genes, genomes can also be compared based on nucleotide composition irrespective of their coding capacity. For instance, the frequency of nucleotide combinations can reveal similarities between genomes that are independent of protein-coding information. We compared the frequency of tetranucleotides for all 145 genomes. The observed frequency of all 64 tetranucleotide combinations was extracted for each genome and these frequencies were divided by the theoretically calculated, expected frequencies (corrected for differences in base composition). This ratio, which could be interpreted as a genomic signature, was expected to reflect taxonomic divisions [29]. However, although the analysis identified a high similarity in tetranucleotide frequency for all of the analyzed Veillonella genomes, most of the clustering observed was not in accordance with known taxonomic relationships. Not only were Negativicutes other than Veillonella separated from each other and strewn across the phyla, but also several other Firmicutes were distributed over various branches (data shown as supplementary material). In fact, for most of the analyzed genomes, members of identical phyla did not cluster together and even the Archaea were mixed with Bacteria, although some closely related species were indeed clustered. This may explain why all Veillonella genomes grouped together. Several organisms with similar tetranucleotide frequencies did not share a common ecological niche, in contrast to previously reported observations (reviewed in [30]). Neither was the obtained clustering dictated by GC-content. The conclusion from this analysis was that tetranucleotide analysis is only taxonomically informative for closely related genomes.

We also compared whole-genome amino acid frequencies in each of the deduced proteomes. Although the results are slightly more in agreement with known taxonomy as compared with the genomic signatures discussed above, this analysis does not cluster organisms according to their phyla, and again some Archaea are mixed with Bacteria. The relevant part of the heatmap based on amino acid frequency is shown in Figure 4. All Veillonella genomes cluster together within the Negativicutes, with the exception of two of the three Dialister genomes, which are found most closely related to Clostridium species (See supplemental information for a version of this figure showing all the genomes). The major Negativicutes cluster also contains a Geobacillus (which is a Gram-positive Firmicutes) and a methanogenic Archaean. Interestingly, the closest relatives to this cluster are not Clostridia, as the previous phylogenetic trees suggest, but a number of Proteobacteria. It is striking that the amino acid frequency analysis detects similarities to Proteobacteria, with which the Negativicutes have their two membranes in common.

Figure 4

A zoomed heatmap of the amino acid frequency found in the deduced proteomes of all 145 genomes. A fragment of the heatmap is shown, presenting the cluster in which all but two Negativicutes are found. The remaining two, both Dialister microaerophilus genomes, are positioned elsewhere in the tree, closest to Clostridium cellulolyticum (not shown in this zoom). The color scale indicates highly underrepresented (orange) to highly overrepresented amino acid frequency (magentum). The full figure is available as supplementary information.

The metabolic properties encoded by the genomes were analyzed next, based on KEGG comparisons [24]. The results are again visualized in a heatmap (Figure 5). We hypothesized that this analysis could identify similarities based on niche adaptation. For simplicity, only a selected number of phyla are shown: apart from the Firmicutes, genomes are included that represent Bacteroidetes and Proteobacteria (both of which contain members frequently found in the oral or gut microbiome), while Cyanobacteria are included as representatives of a phylum that occupy an environmental niche. Since the genomes are compared based on predicted proteomes, their annotation was standardized in order to reduce artificial variation caused by gene annotation differences. As can be seen in Figure 5, the Veillonella genomes all cluster together at the right-hand side of the plot, within a larger cluster containing most of the other Negativicutes and some Firmicutes. The three Dialister species are placed outside the Negativicutes cluster. The other Firmicutes that are found combined with the Negativicutes, based on their metabolic potential, are Clostridium cellulolyticum, Eubacterium rectale, Lactococcus lactis, Streptococcus pneumoniae and Turicibacter sanguinis. These are all common members of the oral or intestine microbiome. As expected, the metabolic pathway for lipopolysaccharide biosynthesis is shared between the Negativicutes and other Gram-negative species, as indicated by the arrows in Figure 5. Interestingly, the Cyanobacteria form a small cluster within, not outside the tree, together with a Haliangium and a Sorangium species as their closest neighbors (both are social Myxococcales belonging to the Deltaproteobacteria). The exclusive ability of carbon fixation by Cyanobacteria is apparent from the dark red square in the block 'energy'. The lanes of Veillonella in Figure 5 are dominated by light colors, indicative of medium metabolic potential; that is, in contrast to some genomes where most of the pathways are present (dark red for Proteobacteria for example) or missing (dark green for other Negativicutes), the Veillonella genomes have partial pathways (based on knowledge primarily from aerobic genomes). There is no reason to believe that the Veillonella genomes should have less metabolic potential than other Negativicutes. Indeed, it is likely that the differences in metabolic potential of Veillonella are truly reflective of alternative capabilities for these bacteria.

Figure 5

Heatmap of metabolism potential, based on Kyoto Encyclopedia of Genes and Genomes ontology (KEGG). The green color in the heatmap indicates weak metabolic potential, while red signals strong potential. The arrows to the right indicate the scores for lipopolysaccharide biosynthesis. A version summarizing the metabolism pathways and showing the species legend is available as supplementary material.

It was further investigated how conserved the predicted proteomes are within the Negativicutes. As a quantitative measure for homology, shared protein-coding genes were identified by pairwise BLASTP comparison and expressed as a percentage of the combined proteomes. The results are shown in a matrix (Figure 6). In addition to the proteomes of the 24 Negativicutes, the comparison includes Clostridium botulinum, Cl. cellulolyticum and Desulfotomaculum reducens, as these Firmicutes were shown to share characteristics with Negativicutes in previous analyses (cf. Figures 1 and 3). The proteome of E. coli K12 is included as an example of a Gram-negative intestinal bacterium. The BLAST matrix was constructed using reciprocal best BLAST hits to determine the presence of shared protein family between two genomes. Inspection of Figure 6 shows that the genus Veillonella is relatively homogeneous; any two members of this genus share between 67% and 90% homology (1,357 to 1,682 protein families), irrespective of the species. The genus Selenomonas is more heterogeneous, with pairwise homology varying from 42% to 82% between any two species (980 to 1659 protein families). The three proteomes of Dialister spp., covering two species, share between 40% and 84% homology. The highest homologous fraction identified between two members of different genera within the Negativicutes is 43% (Mitsuokella multacida compared to Selenomonas sputigena, whereas the lowest homology is 15% (Dialister spp. compared to Thermosinus carboxydivorans). Negativicutes share between 9% and 33% homology with the analyzed Firmicutes, whereas slightly lower homology is detected with E. coli (between 7% and 24%).

Figure 6

Proteome comparison represented by a BLAST matrix, based on 24 Negativicutes genomes with reciprocal best hits. The genomes of Clostridium botulinum, Cl. cellulolyticum, Desulfotomaculum reducens and E. coli are added for comparison. Inter-genus comparisons are indicated by black squares. A version reporting the numerical values of homology percentages is available as supplementary information.

Finally, we assessed the gene pool conserved within all analyzed Negativicutes. Using the same cutoff for protein BLAST comparison as before, a core-genome is identified that contains about 300 conserved protein families (data not shown). This is a relatively low number of conserved proteins, reflective of the extensive genetic heterogeneity within this bacterial class.


The availability of complete sequences for a large and diverse set of Bacterial genomes has helped in exploring the conundrum of the genus Veillonella, a genus within the Negativicutes class, all of which are Gram negative Firmicutes. The 16S rRNA tree shown as Figure 1 illustrates how “close” the Negativicutes are to other Firmicutes. The closest Gram positive Clostridium species are actually quite distant to Veillonella and other Negativicutes genomes, as can be seen in the low fraction of shared protein families in Figure 6. The Gram-negative Firmicutes are even more distant to other Gram negatives, such as Proteobacteria (e.g., E. coli). It should be noted that the family Clostridiaceae is a largely diverse group with many members being re-classified [27]. It is therefore possible that the taxonomic description of some Clostridium genomes may change in future. However, our analyses did not identify one single Gram-positive Firmicutes (Clostrida or others) that consistently was identified as most closely related to Veillonella. As seen from three types of phylogenetic analysis, the Negativicutes class genomes form a distinct cluster within the Firmicutes, and the Veillonella genus forms a relatively homogeneous group of species within the Negativicutes, with relatively conserved metabolic properties (Figure 5). In comparison, the Selenomonas genus is more heterogeneous, at least based on their total gene comparison, as illustrated in Figure 6.

In contrast to expectations, relatively little homology between Negativicutes and other Gram-negative genomes was detected in our analyses. Neither gene-dependent phylogenetic analysis, nor gene-independent DNA tetramer analysis identified a significant commonness between Negativicutes and, say, Proteobacteria. Only whole-genome frequency analysis of amino acid usage identified some similarity to a few Proteobacteria, and this might be more reflective of environment the organism is adapted to, and not phylogeny. Using KEGG pathways for metabolic comparison of the proteomes we found few pathways in common, with the exception of a shared lipopolysaccharide biosynthesis pathway. From all analyses combined, it is clear that the taxonomic placement of Negativicutes within the Firmicutes reflects their genetic and genomic characteristics, although the proteins encoded by the Negativicutes genomes are quite distinct from their Gram-positive cousins. It could be speculated that the double membrane of the Negativicutes evolved in a lineage that used to be a single-membrane (Gram-positive) Firmicute. Whether this event co-evolved independently of the formation of other Gram-negative phyla, or was the result of lateral gene transfer, cannot be stated for certain at present; estimations of horizontally transferred regions in Veillonella parvula DSM 2008, the only fully assembled Veillonella genome available, using the least conservative method on the Islandviewer web-site [31], revealed that only 2% of the genome is of foreign origin. In comparison, 9% of the E. coli K-12 subsp. MG1655 genome was predicted as horizontally transferred. Further analyses are therefore needed to assess this in more detail.

Author’s contributions

Tammi Vesth was a main contributor to the writing of the manuscript and to the organization of the work. Trudy Wassenaar helped considerably in editing and improving the manuscript. Individual contributions: Asli Ozen (16s rRNA and CV tree), Oksana Lukjancenko (consensus tree), Sandra Andersen (initial investigations and background research, early version of the manuscript), Rolf Sommer Kaas (BLAST matrix), Jon Bohlin (tetramer and amino acid usage heatmaps), Intawat Nookaew (metabolism heatmaps). David Ussery provided the original idea for this manuscript, suggested the figures, helped in early drafts of the manuscript, and supervised the project.



This research was supported by grants from the Danish Research Council, and in part by a grant 09-067103/DSF from the Danish Council for Strategic Research.

This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


  1. Delwiche EA, Pestka JJ and Tortorello ML. The veillonellae: gram-negative cocci with a unique physiology. Annu Rev Microbiol. 1985; 39:175-193 View ArticlePubMed
  2. Diaz PI, Chalmers NI, Rickard AH, Kong C, Milburn CL, Palmer RJ and Kolenbrander PE. Molecular Characterization of Subject-Specific Oral Microflora during Initial Colonization of Enamel. Appl Environ Microbiol. 2006; 72:2837-2848 View ArticlePubMed
  3. Aas JA, Paster BJ, Stokes LN, Olsen I and Dewhirst FE. Defining the Normal Bacterial Flora of the Oral Cavity. J Clin Microbiol. 2005; 43:5721-5732 View ArticlePubMed
  4. Mager DL, Ximenez-Fyvie LA, Haffajee AD and Socransky SS. Distribution of selected bacterial species on intraoral surfaces. J Clin Periodontol. 2003; 30:644-654 View ArticlePubMed
  5. Olson JC, Cuff CF, Lukomski S, Lukomska E, Canizales Y, Wu B, Crout RJ, Thomas JG, McNeil DW and Weyant RJ. Use of 16S ribosomal RNA gene analyses to characterize the bacterial signature associated with poor oral health in West Virginia. BMC Oral Health. 2011; 11:7 View ArticlePubMed
  6. Chalmers NI, Palmer RJ, Cisar JO and Kolenbrander PE. Characterization of a Streptococcus sp.-Veillonella sp. Community Micromanipulated from Dental Plaque. J Bacteriol. 2008; 190:8145-8154 View ArticlePubMed
  7. Leuckfeld I, Paster BJ, Kristoffersen AK and Olsen I. Diversity of Veillonella spp. from subgingival plaque by polyphasic approach. APMIS. 2010; 118:230-242 View ArticlePubMed
  8. Tunney MM, Field TR, Moriarty TF, Patrick S, Doering G, Muhlebach MS, Wolfgang MC, Boucher R, Gilpin DF, McDowell A and Elborn JS. Detection of anaerobic bacteria in high numbers in sputum from patients with cystic fibrosis. Am J Respir Crit Care Med. 2008; 177:995-1001 View ArticlePubMed
  9. Murri M, Leiva I, Gomez-Zumaquero JM, Tinahones FJ, Cardona F, Soriguer F and Quepo-Ortuño MI. Gut microbiota in children with type I diabetes differs from that in healthy children: a case-control study. BMC Med. 2013; 11:46 View ArticlePubMed
  10. Kolenbrander PE, Moore LVH. The genus Veillonella in:H.G. Balows, M. Trüper, W. Dworkin, W. Harder, K.H. Schleifer (Eds.), The prokaryotes (2nd ed.), Springer, New York (1992), pp. 2034–2047.
  11. Mashima I, Nakazawa F. Identification of Veillonella tobetsuensis in tongue biofilm by using a species-specific primer pair. Anaerobe 2013; 22:77_81.
  12. De Vos P, Garrity GM, Jones D, Krieg NR, Ludwig W, Rainey FA, Schleifer KH, Whitman WB. Volume 3: The Firmicutes In Bergey’s Manual of Systematic Bacteriology, Springer 2009.
  13. Lagesen K, Hallin P, Rødland EA, Staerfeldt HH, Rognes T and Ussery DW. RNAmmer: consistent and rapid annotation of ribosomal RNA genes. Nucleic Acids Res. 2007; 35:3100-3108 View ArticlePubMed
  14. Katoh K and Toh H. Parallelization of the MAFFT multiple sequence alignment program. Bioinformatics. 2010; 26:1899-1900 View ArticlePubMed
  15. Tamura K, Peterson D, Peterson N, Stecher G, Nei M and Kumar S. MEGA5: Molecular Evolutionary Genetics Analysis using Maximum Likelihood, Evolutionary Distance, and Maximum Parsimony Methods. Mol Biol Evol. 2011; 28:2731-2739 View ArticlePubMed
  16. Saitou N and Nei M. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987; 4:406-425PubMed
  17. Xu Z and Hao B. CVTree update: a newly designed phylogenetic study platform using composition vectors and whole genomes. Nucleic Acids Res. 2009; 37(suppl 2):W174-W178 View ArticlePubMed
  18. Ciccarelli FD, Doerks T, von Mering C, Creevey CJ, Snel B and Bork P. Toward automatic reconstruction of a highly resolved tree of life. Science. 2006; 311:1283-1287 View ArticlePubMed
  19. Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004; 32:1792-1797 View ArticlePubMed
  20. Fink WL. Microcomputers and phylogenetic analysis. Science. 1986; 234:1135-1139 View ArticlePubMed
  21. Retief JD. Phylogenetic analysis using PHYLIP. Methods Mol Biol. 2000; 132:243-258PubMed
  22. Karlin S and Burge C. Dinucleotide relative abundance extremes: a genomic signature. Trends Genet. 1995; 11:283-290 View ArticlePubMed
  23. Wickham H: ggplot2: Elegant Graphics for Data Analysis (Use R!). Springer New York, 2009. ISBN-10: 0387981403 | ISBN-13: 978-0387981406
  24. Kanehisa M, Goto S, Kawashima S, Okuno Y and Hattori M. The KEGG resource for deciphering the genome. Nucleic Acids Res. 2004; 32:277-280 View ArticlePubMed
  25. Eddy SR. BIOINFORMATICS REVIEW Profile hidden Markov models. [PubMed]. Bioinformatics. 1998; 14:755-763 View ArticlePubMed
  26. Hobohm U and Sander C. Enlarged representative set of protein structures. Protein Sci. 1994; 3:522-524 View ArticlePubMed
  27. Ludwig W, Schleifer KH, Whitman W. Revised road map to the phylum Firmicutes. In Bergey’s Manual of Systematic Bacteriology, Springer 2009:1-13.
  28. Gupta RS. Origin of diderm (Gram-negative) bacteria: antibiotic selection pressure rather than endosymbiosis likely led to the evolution of bacterial cells with two membranes. Antonie van Leeuwenhoek. 2011; 100:171-182 View ArticlePubMed
  29. Pride DT, Meinersmann RJ, Wassenaar TM and Blaser MJ. Evolutionary implications of microbial genome tetranucleotide frequency biases. Genome Res. 2003; 13:145-158 View ArticlePubMed
  30. Dutta C and Paul S. Microbial lifestyle and genome signatures. Curr Genomics. 2012; 13:153-162 View ArticlePubMed
  31. Langille MG and Brinkman FS. IslandViewer: an integrated interface for computational identification and visualization of genomic islands. Bioinformatics. 2009; 25:664-665 View ArticlePubMed