Skip to main content

The complete genome sequence of the rumen bacterium Butyrivibrio hungatei MB2003


Butyrivibrio hungatei MB2003 was isolated from the plant-adherent fraction of rumen contents from a pasture-grazed New Zealand dairy cow, and was selected for genome sequencing in order to examine its ability to degrade plant polysaccharides. The genome of MB2003 is 3.39 Mb and consists of four replicons; a chromosome, a secondary chromosome or chromid, a megaplasmid and a small plasmid. The genome has an average G + C content of 39.7%, and encodes 2983 putative protein-coding genes. MB2003 is able to use a variety of monosaccharide substrates for growth, with acetate, butyrate and formate as the principal fermentation end-products, and the genes encoding these metabolic pathways have been identified. MB2003 is predicted to encode an extensive repertoire of CAZymes with 78 GHs, 7 CEs, 1 PL and 78 GTs. MB2003 is unable to grow on xylan or pectin, and its role in the rumen appears to be as a utilizer of monosaccharides, disaccharides and oligosaccharides made available by the degradative activities of other bacterial species.


Butyrivibrio are important rumen bacteria [1], and are among the small number of rumen genera capable of utilizing the complex plant structural polysaccharides xylan and pectin [2, 3]. They are classified as anaerobic, monotrichous, butyrate-producing, curved rods and have been isolated from the gastrointestinal tracts and feces of various ruminants, monogastric animals and humans [4, 5]. Butyrivibrio are metabolically versatile and are capable of growing on a range of carbohydrates, from simple mono- or oligosaccharides to complex plant polysaccharides such as pectins, mannans, starch and hemicelluloses [6]. Furthermore, xylans of diverse chemical and physical properties, from a range of forages are degraded by Butyrivibrio species [7]. Some Butyrivibrio species show strong proteolytic activity [8], and Butyrivibrio are thought to be the main butyrate producers in the rumen [9, 10]. The genus Butyrivibrio is classified within the family Lachnospiraceae , order Eubacteriales , and is phylogenetically diverse. The Butyrivibrio genus originally consisted of only one species, Butyrivibrio fibrisolvens [2]. In addition to phenotypic characterisations [11, 12], studies have utilized DNA-DNA hybridization [13, 14], 16S rRNA gene sequencing [15, 16] and 16S rRNA-based hybridization probes [17], to differentiate these organisms. To accommodate the observed diversity amongst the newly discovered bacterial strains, a new genus, Pseudobutyrivibrio , was described [18]. Four species are currently recognized: B. fibrisolvens , B. hungatei , B. proteoclasticus and B. crossotus [6], although B. crossotus is more distantly related to the other three. B. hungatei are common anaerobic rumen bacteria found in domestic and wild ruminants and the type strain is JK615T [19]. Butyrivibrio hungatei JK615T is non-proteolytic and non-fibrolytic, but is able to utilize oligo- and monosaccharides as substrates for growth. Gaining an insight into the role of these secondary degrader species in microbial plant polysaccharide breakdown is important for understanding rumen function. Here we present the complete genome sequence of Butyrivibrio hungatei MB2003, a strain isolated from a pasture-grazed dairy cow in New Zealand [20], and describe its comparison with genomes of closely related B. hungatei strains.

Organism information

Classification and features

MB2003 was isolated from the plant-adherent fraction of rumen contents from a New Zealand dairy cow grazing fresh forage [20, 21]. MB2003 cells are Gram positive, short rods, occurring singly or in pairs (Fig. 1). The morphological features of MB2003 cells were determined by electron microscopy of cells grown on RM02 medium [22], negatively stained with 1% phosphotungstic acid, mounted on Formvar-coated copper grids, and examined using a Philips model 201C electron microscope (Eindhoven, The Netherlands). MB2003 cells were observed to have a single polar flagellum (Fig. 2), although cells in growing cultures were non-motile. A phylogenetic analysis of the full-length 16S rRNA gene sequence placed MB2003 within the B. hungatei species, being 98% similar to the Butyrivibrio hungatei type strain JK615T [19] (Fig. 3). Additional characteristics of B. hungatei MB2003 are shown in Table 1.

Fig. 1

Morphology of B. hungatei MB2003. Micrograph of Gram stained B. hungatei MB2003 cells at 100 × magnification. Bar represents 10 μm

Fig. 2

Transmission electron micrograph of B. hungatei MB2003. Micrograph of negatively stained B. hungatei MB2003 cells at 10,000 × magnification

Fig. 3

Phylogenetic tree highlighting the relationship of B. hungatei MB2003 relative to the type strains of the other species within the genus Butyrivibrio. The evolutionary history was inferred using the Maximum Likelihood method based on the General Time Reversible model [55]. The tree with the highest log likelihood (−3712.3329) is shown. The percentage of replicate trees in which the associated taxa clustered together in the bootstrap test (10,000 replicates) is shown next to the branches [56]. Initial tree(s) for the heuristic search were obtained automatically by applying Neighbor-Join and BioNJ algorithms to a matrix of pairwise distances estimated using the Maximum Composite Likelihood (MCL) approach, and then selecting the topology with superior log likelihood value. A discrete Gamma distribution was used to model evolutionary rate differences among sites (5 categories (+G, parameter = 0.3950)). The tree is drawn to scale, with branch lengths measured in the number of substitutions per site. The analysis involved six nucleotide sequences. All positions with less than 95% site coverage were eliminated. There were a total of 1509 positions in the final dataset. Evolutionary analyses were conducted in MEGA6 [55]. GenBank accession numbers of the 16S rRNA gene sequences are shown in parentheses. Bar, 0.02 nucleotide substitutions per site. T, indicates type strain. All the type strains have their genome sequencing projects registered in the Genomes Online Database (GOLD) [57]

Table 1 Classification and general features of the rumen bacterium B. hungatei MB2003 in accordance with the MIGS recommendations [58]

Strain MB2003 grew to highest optical density (OD) at pH values of 6.1 to 6.5 and at a temperature of 39 °C, conditions which are typical of its rumen environment. VFA production was determined from triplicate broth cultures grown overnight in RM02 medium with cellobiose as substrate and analysed for formate, acetate, propionate, n-butyrate, iso-valerate and lactate on a HP 6890 series GC (Hewlett Packard) with 2-ethylbutyric acid (Sigma-Aldrich, St. Louis, MO, USA) as the internal standard. To derivatize formic, lactic and succinic acids, samples were mixed with HCl ACS reagent (Sigma-Aldrich, St. Louis, MO, USA) and diethyl ether, with the addition of N-methyl-N-t-butyldimethylsilyltri-fluoroacetamide (MTBSTFA) (Sigma-Aldrich, St. Louis, MO, USA) [23]. Under these conditions MB2003 produced 16.4 mM formate, 3.6 mM acetate and 4.7 mM butyrate. MB2003 was able to grow in CO2-containing media with various soluble carbon sources and the semi-soluble inulin (all tested at 0.5% w/v final concentration). Growth on soluble substrates was assessed as an increase in culture density OD600nm compared to cultures without carbon source added, whereas total VFA production was used as an indicator of substrate utilization and growth for insoluble polymers (Table 2). All strains tested were net producers of formate, acetate and n-butyrate, which is characteristic of Butyrivibrio . Cellobiose and glucose supported the growth of MB2003, JK615T and B316T to high cell densities. Therefore, cellobiose was used to examine the growth of MB2003 over a 24 h period. The exponential phase of growth was between 4 and 8 h, with the maximum cell density reached at 8 to 10 h, and stationary phase between 10 to 24 h (Fig. 4).

Table 2 Carbon source utilization of the Butyrivibrio strains
Fig. 4

Culture density achieved in 24 h by MB2003 growing in media with cellobiose as the sole substrate. Points indicate means of three replicates, and the error bars represent +/−one standard error

Genome sequencing information

Genome project history

Butyrivibrio hungatei MB2003 was selected for genome sequencing as a NZ strain of B. hungatei . A summary of the genome project information is shown in Table 3 and in Additional file 1: Table S1.

Table 3 MB2003 genome project information

Growth conditions and genomic DNA preparation

MB2003 was grown in RM02 medium [22] with 10 mM glucose and 0.1% yeast extract but without rumen fluid. Culture purity was confirmed by Gram stain and sequencing of the 16S rRNA gene. Genomic DNA was extracted from freshly grown cells by a modification of the standard cell lysis method of Saito and Miura [24], using lysozyme, proteinase K and sodium dodecyl sulphate, followed by phenol-chloroform extraction, and purification using the Qiagen Genomic-Tip 500 Maxi kit (Qiagen, Hilden, Germany). Genomic DNA was precipitated by the addition of a 0.7 volume of isopropanol, and collected by centrifugation at 12,000×g for 10 min at room temperature. The supernatant was removed, and the DNA pellet was washed in 70% ethanol, re-dissolved in TE buffer (10 mM Tris-HCl, 1 mM EDTA, pH 7.5) and stored at −20 °C until required.

Genome sequencing and assembly

The complete genome sequence of MB2003 was determined by pyrosequencing 3 kb mate paired-end sequence libraries using the 454 GS FLX platform with Titanium chemistry (Macrogen, Korea). Pyrosequencing reads provided 234× coverage of the genome and were assembled using the Newbler assembler (version 2.7, Roche 454 Life Sciences, USA) which resulted in 31 contigs across 7 scaffolds. Gap closure was managed using the Staden package [25] and gaps were closed using additional Sanger sequencing by standard and inverse PCR techniques. In addition, MB2003 genomic DNA was sequenced using shotgun sequencing of 2 kb paired-end sequence libraries using the Illumina MiSeq platform (Macrogen, Korea) which provided 800-fold sequencing coverage. Illumina reads were analysed using the Galaxy web-based platform [26] and de novo assembly was performed using the Velvet assembler, version 3.0 [27]. The Velvet assembled MB2003 genome MiSeq sequences were combined with the Newbler assembly using the Staden package and Geneious, version 8.1 [28]. Genome assembly was confirmed by pulsed-field gel electrophoresis.

Genome annotation

Annotation of the MB2003 genome was performed as described previously [29]. The MB2003 genome sequence was prepared for NCBI submission using Sequin [30], and the adenine residue of the start codon of the chromosomal replication initiator protein DnaA1 (bhn_I0001, bhn_RS00450) gene was chosen as the first base for the MB2003 genome.

Genome properties

The genome of B. hungatei MB2003 consists of four replicons [21, 31]; a single chromosome (3,143,784 bp, %G + C 39.91), a chromid or secondary chromosome (BhuII, 91,776 bp, %G + C 37.71), a megaplasmid (pNP144, 144,470 bp, %G + C 36.86) and a plasmid (pNP6, 6284 bp, %G + C 35.71). The total size of the closed genome is 3,386,314 bp with an overall %G + C content of 39.71%. A total of 3064 genes were predicted, of which 2983 (97.4%) were protein-coding genes. A putative function was assigned to 2225 of the protein-coding genes, while 775 protein coding genes were annotated as hypothetical proteins. The MB2003 chromosome encodes 2758 genes, and BhuII, pNP144 and pNP6 encode 89, 147 and 6 genes, respectively. The properties and statistics of the MB2003 genome are summarized in Tables 4, 5 and 6. The nucleotide sequences of the MB2003 chromosome, chromid (BhuII), megaplasmid (pNP144) and plasmid (pNP6) have been deposited in Genbank under accession numbers CP017831, CP017830, CP017832 and CP017833. The genome atlas for B. hungatei MB2003 is shown in Fig. 5.

Table 4 Summary of MB2003 genome replicon features
Table 5 MB2003 genome statistics
Table 6 Number of genes associated with the general COG functional categories
Fig. 5

Genome atlas for B. hungatei MB2003. The figure represents a circular view of the four replicons that make up the B. hungatei MB2003 genome. The key at the right describes the concentric circles within each replicon in the outermost to innermost direction. The diagram was created using GENEWIZ [66] and custom-developed software. The innermost circle 1 shows GC-skew; Circle 2 shows COG classification: predicted ORFs were analysed using the COG database and grouped into the five major categories: yellow, information storage and processing; red, cellular processes and signalling; green, metabolism; blue, poorly characterised; and uncoloured, ORFs with uncharacterized COGs or no COG assignment. Circle 3 shows transmembrane helices (TMH) and SignalP domains: the four categories represent, uncoloured, absent; red, TMH; blue, SignalP; and black, both TMH and SignalP present. Circle 4 shows ORF orientation: ORFs in sense orientation (ORF+) are shown in blue; ORFs oriented in antisense direction (ORF-) are shown in red. Circle 5 shows ribosomal machinery: tRNAs and rRNAs are shown as green or red lines, respectively. Clusters are represented as coloured boxes to maintain readability. Circle 6 shows G + C content deviation from the average: GC-content is shown in either green (low GC spike) or orange (high GC spike). A box filter was applied to visualize contiguous regions of low or high GC deviations. Circle 7 shows BLAST similarities: deduced amino acid sequences were compared against the nonredundant (nr) database using gapped BLASTP [67]. Regions in blue represent unique proteins in MB2003, whereas highly conserved features relative to sequences in the nr database are shown in red. The degree of colour saturation corresponds to the level of similarity. The predicted origin and terminus of DNA replication are indicated

Insights from the genome sequence

Comparison of the MB2003, B. hungatei JK615T, and B. proteoclasticus B316T genomes

A comparison of the B. hungatei MB2003 genome with the draft genome of B. hungatei JK615T [32] and the complete B. proteoclasticus B316T genome is shown in Table 7. The MB2003 genome is 8633 bp smaller than JK615T and contains 27 fewer protein-coding genes. Although several plasmid replication genes have been identified in the JK615T draft genome, the presence of extrachromosomal elements requires experimental validation.

Table 7 Genome statistics of MB2003, JK615T and B316T

A novel feature of both the MB2003 and B316T genomes is the presence of chromids or secondary chromosomes [33]. Chromids are replicons that have %G + C content similar to that of their main chromosome, but have plasmid-type maintenance and replication systems, are smaller than the chromosome, but are usually larger than any other plasmids present. Chromids contain genes essential for growth and maintenance of the organism along with several core genus-specific genes that can be found on the chromosome in other species of bacteria [33]. The Bhu II replicon has most of these characteristics and therefore has been designated as a chromid of MB2003. In B316T, almost 10% of the genes encoding enzymes that have a role in carbohydrate metabolism and transport are found on the chromid [29]. The Bhu II chromid of MB2003 also encodes genes with similar predicted functions (Table 9). Since the Bhu II chromid of MB2003 is smaller than the BPc2 chromid of B316T (186,325 bp), it is now the smallest chromid reported for bacteria. Comparison of MB2003, JK615T and B316T genomes based on COG category (Table 8) and synteny analysis (Fig. 6), show that these Butyrivibrio species and strains are genetically similar. Although the MB2003 and B316T genome sizes differ, the basic metabolism of these two rumen bacterial species is indicated to be similar.

Table 8 Comparison of MB2003, JK615T and B316T protein coding gene percentages to COG functional categories
Fig. 6

Genome synteny analysis. Alignment of the B. hungatei MB2003 genome against the draft genome of B. hungatei JK615T (a) and the complete genome of B. proteoclasticus B316T (b). Whenever the two sequences agree, a colored line or dot is plotted. If the two sequences were perfectly identical, a single line would go from the bottom left to the top right. Units displayed in base-pairs. Color codes: blue, forward sequence, red, reverse sequence

Butyrate production

For the production of butyrate and H2 from glucose, the MB2003 genome possesses a pyruvate:ferredoxin oxidoreductase gene, nifJ (bhn_I2528) required for pyruvate conversion to acetyl-CoA, as well as a butyryl-CoA dehydrogenase/electron transferring flavoprotein bcd-etfAB (bhn_I2225, bhn_I2221 and bhn_I2222) to generate ATP by classic substrate level phosphorylation. In addition, MB2003 possesses genes that encode all six subunits of the Rnf (rnfA, rnfB, rnfC, rnfD, rnfE, rnfG) and Ech (echA, echB, echC, echD, echE, echF) hydrogenases. These pathways involve the transmembrane ion pumps Ech [34] or Rnf [35,36,37,38], that generate a transmembrane proton and/or sodium electrochemical potential from redox cofactors for ATP synthesis by ETP [34, 36]. The MB2003 genome does not possess genes for PorABDG, a pyruvate ferredoxin oxidoreductase similar in function to NifJ or genes for EhaA-R, EhbA-P, HydA-C, MbhLKJ, or MvhADG/HdrABC similar in function to the Fd-dependent Ech hydrogenase. In addition, an alternative pathway exists where formate is predicted to be the end product and involves the decarboxylation of acetyl-CoA by a pyruvate formate lyase pflB (bhn_I0124) instead of NifJ. It has been proposed that Ech and Rnf work in concert with NifJ and Bcd-Etf complex to drive ATP synthesis by ETP during glucose fermentation to butyrate [34, 36, 39]. Interestingly, the vast majority of anaerobic prokaryotes appear to possess either an Ech or Rnf but not both [40, 41]. However, a recent analysis of rumen prokaryotic genomes identified Butyrivibrio and Pseudobutyrivibrio as a rare group of bacteria that possess genes for both Ech and Rnf [42]. These findings warrant further biochemical investigation to determine the activity of Ech and Rnf in Butyrivibrio .

The MB2003 pathways for butyrate production presume the possession of a complete Embden-Meyerhof-Parnas glycolytic pathway. Enolase (eno, EC4.2.1.11), converts 2-phospho-D-glycerate to phosphoenolpyruvate in the second to last step of the EMP pathway. Previous work has shown that B316T lacks a detectable enolase [29], and the Methylglyoxal Shunt was proposed as a possible alternative to the EMP pathway. In this pathway the dihydroxyacetone phosphate is transformed to pyruvate via methylglyoxal and d-lactate dehydrogenase, encoded by ldhA [43]. The MB2003 genome possesses two methylglyoxal synthase genes, mgsA (bhn_I1328 and bhn_I1996), glyoxalase gene gloA (bhn_I1783) and an alternative l-lactate dehydrogenase, encoded by ldh (bhn_I0363). MB2003 has the same set of genes as B316T for the production of butyrate, formate, acetate and lactate, but also is the only B. hungatei reported to date that lacks a detectable enolase gene. Genome sequences from a wider range of B. hungatei and B. proteoclasticus strains are required to determine if these are common features in these organisms.

Polysaccharide degradation

The Carbohydrate-Active enZYmes database was used to identify glycoside hydrolases, glycosyl transferases, polysaccharide lyases, carbohydrate esterases and carbohydrate-binding protein module families within the MB2003 genome. MB2003 has a similar CAZyme profile to B316T [21, 31], and analysis of the functional domains of enzymes involved in the breakdown or synthesis of complex carbohydrates, has revealed the polysaccharide-degrading potential of this rumen bacterium.

Approximately 3% of the MB2003 genome (90 CDSs) is predicted to encode either secreted or intracellular proteins dedicated to polysaccharide degradation, similar to that found in B316T. The MB2003 genome is predicted to encode 19 secreted (16 GHs, two CEs and one CBP) and 65 intracellular (59 GHs, 5 CEs and one PL) proteins involved in polysaccharide breakdown (Table 9). The enzymatic profiles of MB2003 and JK615T are almost identical, as both possess the same genes encoding predicted secreted and intracellular CAZymes in their genomes (Table 9). Out of the 19 genes predicted to encode secreted polysaccharide degrading enzymes, only two, lysozyme lyc25B (bhn_III074) and feruloyl esterase est1A (bhn_III076), are encoded by the MB2003 chromid (Bhu II). MB2003 has no secreted enzyme larger than 1000 aa in size, with the average size secreted enzymes being 510 aa. The majority (59) of MB2003 genes involved in polysaccharide breakdown (excluding GTs), had corresponding homologues in B316T and JK615T. Three of the genes encoding intracellular proteins were found in the Bhu II chromid: a β-glucosidase bgl3A (bhn_III062), a β-galactosidase bga42A (bhn_III010) and a polysaccharide deacetylase est4A (bhn_III070). The analysis of the Pfam domains from the most abundant GH families (GH2, GH31, GH3, GH13 and GH43) showed they did not contain signal sequences and hence were predicted to be located intracellularly. Similarly, CAZymes with predicted roles in xylan and pectin degradation, the GH8, GH28, GH51, GH67, GH88, GH105, GH115, CE2 and CE10 families were also predicted to be intracellular. Of these, MB2003 contains CAZymes with homologues in B316T except for the α-L-arabinofuranosidase arf51C (bhn_I1509). These findings suggest that a variety of complex oligosaccharides resulting from extracellular hydrolysis are metabolized within the cell.

Table 9 Genes encoding predicted polysaccharide degrading enzymes in the MB2003 genome

Growth experiments showed MB2003 to be a metabolically versatile bacterium able to grow on a wide variety of monosaccharides, disaccharides and glycosides (Table 2). However, unlike B316T, MB2003 and JK615T were unable to utilize the insoluble substrates pectin and xylan for growth (Table 2). In addition, MB2003, JK615T and B316T are unable to degrade cellulose, however among these organisms, only B316T is able to utilize a range of other insoluble plant polysaccharides. The ability of B316T to breakdown pectin, starch and xylan is predicted to be based on nine large (>1000 aa) cell-associated proteins shown to be significantly up-regulated in B316T cells grown on xylan [44]. These are: α-amylase amy13A (bpr_I1087), arabinogalactan endo-1,4-β-galactosidase agn53A (bpr_I2041), carbohydrate esterase family 12 est12B (bpr_I1204), endo-1,3(4)-β-glucanase lic16A (bpr_I2326), pectate lyase pel1A (bpr_I2372), pectin methylesterase pme8B (bpr_I2473), xylosidase/arabinofuranosidase xsa43J (bpr_I2935), endo-1,4-β-xylanase xyn10B (bpr_I0026), and the cell wall binding domain-containing protein (bpr_I0264). These proteins contain multiple cell wall binding repeat domains (CW-binding domain, Pfam01473) at their C-termini that are predicted to anchor the protein to the peptidoglycan cell membrane. Among these secreted polysaccharidases, some contain single or combinations of catalytic activities: GH10 (endo-1, 4-β-xylanase, xyn10B), GH43 (xylosidase/arabinofuranosidase, xsa43J), PL1 (pectate lyase, pel1A), CE8 and PL9 (pectin methylesterase, pme8B) [45, 46]. Neither MB2003 nor JK615T contain any genes encoding CW-binding domains and are thus are markedly different from B316T.

A curious feature of MB2003 was the presence of a single large (983 aa) carbohydrate binding protein (CBP, bhn_I1848), also present in JK615T (EJ23DRAFT_00192). The domain structures of bhn_I1848 and EJ23DRAFT_00192 are unusual, containing six CBM6 (Pfam03422) domains towards the N-terminus and a single C-terminal CBM2a (Pfam00553) domain. In contrast, B316T encodes two CBPs (bpr_I0736 and bpr_I1599) where both contain two CBM2a domains, and bpr_I1599 also contains two CBM6 domains [29]. CBM6 non-catalytic modules characteristically bind xylose and are associated with xylanase activity with ligand specificity for xylan [47, 48]. CBM2 domains, are divided into two sub-families: 2a, that bind to crystalline cellulose even when associated with xylanases [49], and 2b, that bind to xylan [50]. Recent studies have shown that in discrete regions of plant cell walls, initial enzymatic attack of pectin increases the access of CBMs to cellulose [51], effectively loosening the polysaccharide interactions to expose the xylan and xyloglucan substrates [52, 53]. This initial stage in enzymatic saccharification of plant cell walls termed amorphogenesis [54], and is a possible role for such CBPs containing multiple non-catalytic domains. In the rumen, MB2003, B316T and JK615T may secrete these non-catalytic CBPs synergistically with polysaccharide-active enzymes as a mechanism to disrupt the interface between polysaccharides to enhance the rate and extent of plant cell wall degradation.


The B. hungatei MB2003 genome sequence adds valuable information regarding the polysaccharide-degrading potential present in the genus Butyrivibrio . Genomic comparisons revealed that B. hungatei MB2003 shows a high level of similarity with B. hungatei JK615T and B. proteoclasticus B316T type strains, including genes involved in production of butyrate, formate, acetate and lactate. While MB2003 and JK615T encode a large repertoire of enzymes predicted to metabolize insoluble polysaccharides such as xylan and pectin, they are unable to grow on these substrates and instead appear to be equipped to utilize mainly oligo- and monosaccharides as substrates for growth. Although MB2003 has similar phenotypic characteristics and occupies the same habitat as other Butyrivibrio species, its genome encodes fewer extracellular polysaccharide degrading enzymes, in particular, those that contain multiple cell wall binding repeat domains. The overall genome similarities, metabolic versatility and differences in the abundance of CAZymes observed in B. proteoclasticus and B. hungatei offers a new view of the genes required for polysaccharide degradation in the rumen. MB2003 appears to occupy a ruminal niche as a secondary degrader of oligosaccharides, in order to coexist with fibre-degrading organisms in this dynamic and competitive environment.



Base pair(s)


Carbohydrate-Active enZYmes


Carbohydrate-Binding Module(s)


Carbohydrate Esterase(s)


Escherichia coli hydrogenase-3-type hydrogenase


Energy-converting hydrogenase A


Energy-converting hydrogenase B




Electron transport phosphorylation


Glycoside Hydrolase(s)


Glycosyl Transferase(s)


Ferredoxin hydrogenase


Membrane-bound hydrogenase


Methyl viologen hydrogenase/heterodisulfide reductase


Pyruvate formate lyase


Polysaccharide Lyase(s)


Pyruvate ferredoxin oxidoreductase


Rhodobacter nitrogen fixation


  1. 1.

    Henderson G, Cox F, Ganesh S, Jonker A, Young W, Collaborators GRC, Janssen PH. Rumen microbial community composition varies with diet and host, but a core microbiome is found across a wide geographical range. Sci Rep. 2015;5(14567):1–13.

    Google Scholar 

  2. 2.

    Bryant MP, Small N. The anaerobic monotrichous butyric acid-producing curved rod-shaped bacteria of the rumen. J Bacteriol. 1956;72(1):16–21.

    CAS  PubMed  PubMed Central  Google Scholar 

  3. 3.

    Hungate RE. The rumen and its microbes. New York: Academic Press; 1966.

    Google Scholar 

  4. 4.

    Willems A, Collins MD. Butyrivibrio. Bryant and Small 1956b, 18, emend. Moore, Johnson and Holdeman 1976, 241AL. In: De Vos P, Garrity G, Jones D, Krieg N, Ludwig W, Rainey F, Schleifer K-H, Whitman W, editors. Bergey's manual of systematic bacteriology. Volume 3. 2nd ed. New York: Springer; 2009. p. 927–37.

    Google Scholar 

  5. 5.

    Hespell RB. The genera Butyrivibrio, Lachnospira, and Roseburia. In: Balows A, Trüper HG, Dworkin M, Harder W, Schleifer K-H, editors. The prokaryotes: a handbook on the biology of bacteria: ecophysiology, isolation, identification, applications. New York: Springer; 1992. p. 2023–33.

    Google Scholar 

  6. 6.

    Moon CD, Pacheco DM, Kelly WJ, Leahy SC, Li D, Kopecny J, Attwood GT. Reclassification of Clostridium proteoclasticum as Butyrivibrio proteoclasticus comb. nov., a butyrate-producing ruminal bacterium. Int J Syst Evol Microbiol. 2008;58(9):2041–5.

    CAS  Article  PubMed  Google Scholar 

  7. 7.

    Hespell RB, Cotta MA. Degradation and utilization by Butyrivibrio fibrisolvens H17c of xylans with different chemical and physical properties. Appl Environ Microbiol. 1995;61(8):3042–50.

    CAS  PubMed  PubMed Central  Google Scholar 

  8. 8.

    Attwood GT, Reilly K. Characterization of proteolytic activities of rumen bacterial isolates from forage-fed cattle. J Appl Bacteriol. 1996;81(5):545–52.

    CAS  PubMed  Google Scholar 

  9. 9.

    Diez-Gonzalez F, Bond DR, Jennings E, Russell JB. Alternative schemes of butyrate production in Butyrivibrio fibrisolvens and their relationship to acetate utilization, lactate production, and phylogeny. Arch Microbiol. 1999;171(5):324–30.

    CAS  Article  PubMed  Google Scholar 

  10. 10.

    Paillard D, McKain N, Rincon MT, Shingfield KJ, Givens DI, Wallace RJ. Quantification of ruminal Clostridium proteoclasticum by real-time PCR using a molecular beacon approach. J Appl Microbiol. 2007;103(4):1251–61.

    CAS  Article  PubMed  Google Scholar 

  11. 11.

    Shane BS, Gouws L, Kistner A. Cellulolytic bacteria occurring in the rumen of sheep conditioned to low-protein teff hay. J Gen Microbiol. 1969;55(3):445–57.

    CAS  Article  PubMed  Google Scholar 

  12. 12.

    Hazlewood GP, Theodorou MK, Hutchings A, Jordan DJ, Galfre G. Preparation and characterization of monoclonal antibodies to a Butyrivibrio sp. and their potential use in the identification of rumen Butyrivibrios, using an enzyme-linked immunosorbent assay. J Gen Microbiol. 1986;132(1):43–52.

    Google Scholar 

  13. 13.

    Mannarelli BM. Deoxyribonucleic acid relatedness among strains of the species Butyrivibrio fibrisolvens. Int J Syst Bacteriol. 1988;38(4):340–7.

    Article  Google Scholar 

  14. 14.

    Mannarelli BM, Stack RJ, Lee D, Ericsson L. Taxonomic relatedness of Butyrivibrio, Lachnospira, Roseburia, and Eubacterium species as determined by DNA hybridization and extracellular-polysaccharide analysis. Int J Syst Bacteriol. 1990;40(4):370–8.

    CAS  Article  Google Scholar 

  15. 15.

    Forster R, Teather R, Gong J, Deng S-J. 16S rDNA analysis of Butyrivibrio fibrisolvens: phylogenetic position and relation to butyrate-producing anaerobic bacteria from the rumen of white-tailed deer. Lett Appl Microbiol. 1996;23(4):218–22.

    CAS  Article  PubMed  Google Scholar 

  16. 16.

    Willems A, Amat-Marco M, Collins MD. Phylogenetic analysis of Butyrivibrio strains reveals three distinct groups of species within the Clostridium subphylum of the gram-positive bacteria. Int J Syst Bacteriol. 1996;46(1):195–9.

    CAS  Article  PubMed  Google Scholar 

  17. 17.

    Forster RJ, Gong J, Teather RM. Group-specific 16S rRNA hybridization probes for determinative and community structure studies of Butyrivibrio fibrisolvens in the rumen. Appl Environ Microbiol. 1997;63(4):1256–60.

    CAS  PubMed  PubMed Central  Google Scholar 

  18. 18.

    Van Gylswyk NO, Hippe H, Rainey FA. Pseudobutyrivibrio ruminis gen. Nov., sp. nov., a butyrate-producing bacterium from the rumen that closely resembles Butyrivibrio fibrisolvens in phenotype. Int J Syst Bacteriol. 1996;46(2):559–63.

    Article  Google Scholar 

  19. 19.

    Kopecny J, Zorec M, Mrazek J, Kobayashi Y, Marinsek-Logar R. Butyrivibrio hungatei sp. nov. and Pseudobutyrivibrio xylanivorans sp. nov., butyrate-producing bacteria from the rumen. Int J Syst Evol Microbiol. 2003;53(1):201–9.

    CAS  Article  PubMed  Google Scholar 

  20. 20.

    Noel S. Cultivation and community composition analysis of plant-adherent rumen bacteria. PhD thesis. Massey University, Institute of Molecular Biosciences; 2013.

  21. 21.

    Palevich N. Genome sequencing of rumen bacteria involved in lignocellulose digestion. MSc thesis: Massey University, Institute of Molecular Biosciences; 2011.

    Google Scholar 

  22. 22.

    Kenters N, Henderson G, Jeyanathan J, Kittelmann S, Janssen PH. Isolation of previously uncultured rumen bacteria by dilution to extinction using a new liquid culture medium. J Microbiol Methods. 2011;84(1):52–60.

    Article  PubMed  Google Scholar 

  23. 23.

    Richardson A, Calder A, Stewart C, Smith A. Simultaneous determination of volatile and non-volatile acidic fermentation products of anaerobes by capillary gas chromatography. Lett Appl Microbiol. 1989;9(1):5–8.

    CAS  Article  Google Scholar 

  24. 24.

    Saito H, Miura K-I. Preparation of transforming deoxyribonucleic acid by phenol treatment. Biochim Biophys Acta. 1963;72:619–29.

    CAS  Article  PubMed  Google Scholar 

  25. 25.

    Staden R, Beal KF, Bonfield JK. The staden package. In: Misener S, Krawetz SA, editors. Bioinformatics methods and protocols, vol. Volume 132. New York: Humana Press; 1999. p. 115–30.

    Chapter  Google Scholar 

  26. 26.

    Goecks J, Nekrutenko A, Taylor J. Galaxy: a comprehensive approach for supporting accessible, reproducible, and transparent computational research in the life sciences. Genome Biol. 2010;11(8):R86.

    Article  PubMed  PubMed Central  Google Scholar 

  27. 27.

    Zerbino DR, Birney E. Velvet: algorithms for de novo short read assembly using De Bruijn graphs. Genome Res. 2008;18(5):821–9.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  28. 28.

    Kearse M, Moir R, Wilson A, Stones-Havas S, Cheung M, Sturrock S, Buxton S, Cooper A, Markowitz S, Duran C. Geneious basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics. 2012;28(12):1647–9.

    Article  PubMed  PubMed Central  Google Scholar 

  29. 29.

    Kelly WJ, Leahy SC, Altermann E, Yeoman CJ, Dunne JC, Kong Z, Pacheco DM, Li D, Noel SJ, Moon CD, Cookson AL, Attwood GT. The glycobiome of the rumen bacterium Butyrivibrio proteoclasticus B316T highlights adaptation to a polysaccharide-rich environment. PLoS One. 2010;5(8):e11942.

    Article  PubMed  PubMed Central  Google Scholar 

  30. 30.

    Benson DA, Clark K, Karsch-Mizrachi I, Lipman DJ, Ostell J, Sayers EW. GenBank. Nucleic Acids Res. 2015;43(D1):D30–5.

    CAS  Article  PubMed  Google Scholar 

  31. 31.

    Palevich N. Comparative genomics of Butyrivibrio and Pseudobutyrivibrio from the rumen. PhD thesis. Massey University, Institute of Fundamental Sciences; 2016.

  32. 32.

    Kyrpides NC, Hugenholtz P, Eisen JA, Woyke T, Göker M, Parker CT, Amann R, Beck BJ, Chain PS, Chun J. Genomic encyclopedia of bacteria and archaea: sequencing a myriad of type strains. PLoS Biol. 2014;12(8):e1001920.

    Article  PubMed  PubMed Central  Google Scholar 

  33. 33.

    Harrison PW, Lower RP, Kim NK, Young JP. Introducing the bacterial 'chromid': not a chromosome, not a plasmid. Trends Microbiol. 2010;18(4):141–8.

    CAS  Article  PubMed  Google Scholar 

  34. 34.

    Welte C, Krätzer C, Deppenmeier U. Involvement of Ech hydrogenase in energy conservation of Methanosarcina mazei. FEBS J. 2010;277(16):3396–403.

    CAS  Article  PubMed  Google Scholar 

  35. 35.

    Biegel E, Müller V. Bacterial Na+-translocating ferredoxin: NAD+ oxidoreductase. Proc Natl Acad Sci. 2010;107(42):18138–42.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  36. 36.

    Buckel W, Thauer RK. Energy conservation via electron bifurcating ferredoxin reduction and proton/Na+ translocating ferredoxin oxidation. Biochim Biophys Acta Bioenerg. 2013;1827(2):94–113.

    CAS  Article  Google Scholar 

  37. 37.

    Schlegel K, Welte C, Deppenmeier U, Müller V. Electron transport during aceticlastic methanogenesis by Methanosarcina acetivorans involves a sodium-translocating Rnf complex. FEBS J. 2012;279(24):4444–52.

    CAS  Article  PubMed  Google Scholar 

  38. 38.

    Tremblay P-L, Zhang T, Dar SA, Leang C, Lovley DR. The Rnf complex of Clostridium ljungdahlii is a proton-translocating ferredoxin: NAD+ oxidoreductase essential for autotrophic growth. MBio. 2013;4(1):e00406–12.

    CAS  Google Scholar 

  39. 39.

    Gutekunst K, Chen X, Schreiber K, Kaspar U, Makam S, Appel J. The bidirectional NiFe-hydrogenase in Synechocystis sp. PCC 6803 is reduced by flavodoxin and ferredoxin and is essential under mixotrophic, nitrate-limiting conditions. J Biol Chem. 2014;289(4):1930–7.

    CAS  Article  PubMed  Google Scholar 

  40. 40.

    Pereira IAC, Ramos A, Grein F, Marques M, Da Silva S, Venceslau S. A comparative genomic analysis of energy metabolism in sulfate reducing bacteria and archaea. Front Microbiol. 2011;2(69).

  41. 41.

    Weghoff MC, Bertsch J, Müller V. A novel mode of lactate metabolism in strictly anaerobic bacteria. Environ Microbiol. 2015;17(3):670–7.

    CAS  Article  PubMed  Google Scholar 

  42. 42.

    Hackmann TJ, Firkins JL. Electron transport phosphorylation in rumen butyrivibrios: unprecedented ATP yield for glucose fermentation to butyrate. Front Microbiol. 2015;6(622).

  43. 43.

    Cooper R. Metabolism of methylglyoxal in microorganisms. Annu Rev Microbiol. 1984;38(1):49–68.

    CAS  Article  PubMed  Google Scholar 

  44. 44.

    Leahy SC, Kelly WJ, Altermann E, Ronimus RS, Yeoman CJ, Pacheco DM, Li D, Kong Z, McTavish S, Sang C. The genome sequence of the rumen methanogen Methanobrevibacter ruminantium reveals new possibilities for controlling ruminant methane emissions. PLoS One. 2010;5(1):e8926.

    Article  PubMed  PubMed Central  Google Scholar 

  45. 45.

    Dunne JC, Li D, Kelly WJ, Leahy SC, Bond JJ, Attwood GT, Jordan TW. Extracellular polysaccharide-degrading proteome of Butyrivibrio proteoclasticus. J Proteome Res. 2011;11(1):131–42.

    Article  PubMed  Google Scholar 

  46. 46.

    Bond JJ, Dunne JC, Kwan FY, Li D, Zhang K, Leahy SC, Kelly WJ, Attwood GT, Jordan TW. Carbohydrate transporting membrane proteins of the rumen bacterium, Butyrivibrio proteoclasticus. J Proteome. 2012;75(11):3138–44.

    CAS  Article  Google Scholar 

  47. 47.

    McCartney L, Gilbert HJ, Bolam DN, Boraston AB, Knox JP. Glycoside hydrolase carbohydrate-binding modules as molecular probes for the analysis of plant cell wall polymers. Anal Biochem. 2004;326(1):49–54.

    CAS  Article  PubMed  Google Scholar 

  48. 48.

    Gilbert HJ, Knox JP, Boraston AB. Advances in understanding the molecular basis of plant cell wall polysaccharide recognition by carbohydrate-binding modules. Curr Opin Struct Biol. 2013;23(5):669–77.

    CAS  Article  PubMed  Google Scholar 

  49. 49.

    Gilkes NR, Jervis E, Henrissat B, Tekant B, Miller R, Warren R, Kilburn D. The adsorption of a bacterial cellulase and its two isolated domains to crystalline cellulose. J Biol Chem. 1992;267(10):6743–9.

    CAS  PubMed  Google Scholar 

  50. 50.

    Simpson PJ, Bolam DN, Cooper A, Ciruela A, Hazlewood GP, Gilbert HJ, Williamson MP. A family IIb xylan-binding domain has a similar secondary structure to a homologous family IIa cellulose-binding domain but different ligand specificity. Structure. 1999;7(7):853–64.

    CAS  Article  PubMed  Google Scholar 

  51. 51.

    Blake AW, McCartney L, Flint JE, Bolam DN, Boraston AB, Gilbert HJ, Knox JP. Understanding the biological rationale for the diversity of cellulose-directed carbohydrate-binding modules in prokaryotic enzymes. J Biol Chem. 2006;281(39):29321–9.

    CAS  Article  PubMed  Google Scholar 

  52. 52.

    Marcus SE, Verhertbruggen Y, Hervé C, Ordaz-Ortiz JJ, Farkas V, Pedersen HL, Willats WG, Knox JP. Pectic homogalacturonan masks abundant sets of xyloglucan epitopes in plant cell walls. BMC Plant Biol. 2008;8(60):1–12.

    Google Scholar 

  53. 53.

    Hervé C, Rogowski A, Gilbert HJ, Paul Knox J. Enzymatic treatments reveal differential capacities for xylan recognition and degradation in primary and secondary plant cell walls. Plant J. 2009;58(3):413–22.

    Article  PubMed  Google Scholar 

  54. 54.

    Coughlan MP. The properties of fungal and bacterial cellulases with comment on their production and application. Biotechnol Genet Eng Rev. 1985;3(1):39–110.

    CAS  Article  Google Scholar 

  55. 55.

    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.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  56. 56.

    Felsenstein J. Confidence limits on phylogenies: an approach using the bootstrap. Evolution. 1985:783–91.

  57. 57.

    Reddy T, Thomas AD, Stamatis D, Bertsch J, Isbandi M, Jansson J, Mallajosyula J, Pagani I, Lobos EA, Kyrpides NC. The genomes OnLine database (GOLD) v. 5: a metadata management system based on a four level (meta) genome project classification. Nucleic Acids Res. 2015;43(D1):D1099–106.

    CAS  Article  PubMed  Google Scholar 

  58. 58.

    Field D, Garrity G, Gray T, Morrison N, Selengut J, Sterk P, Tatusova T, Thomson N, Allen MJ, Angiuoli SV. The minimum information about a genome sequence (MIGS) specification. Nat Biotechnol. 2008;26(5):541–7.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  59. 59.

    Woese CR, Kandler O, Wheelis ML. Towards a natural system of organisms: proposal for the domains Archaea, bacteria, and Eucarya. Proc Natl Acad Sci. 1990;87(12):4576–9.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  60. 60.

    Rainey F. Phylum SK-H, XIII. Firmicutes gibbons and Murray 1978, 5 (Fimacutes [sic] gibbons and Murray 1978, 5). In: De Vos P, Garrity G, Jones D, Krieg N, Ludwig W, Rainey F, Schleifer K-H, Whitman W, editors. Bergey’s manual of systematic bacteriology, vol. Volume 3. 2nd ed. New York: Springer; 2009. p. 19–1317.

    Google Scholar 

  61. 61.

    Ludwig W, Schleifer K-H, Whitman W. Revised road map to the phylum Firmicutes. In: De Vos P, Garrity G, Jones D, Krieg N, Ludwig W, Rainey F, Schleifer K-H, Whitman W, editors. Bergey’s manual of systematic bacteriology. Volume 3. 2nd ed. New York: Springer; 2009. p. 1–13.

    Google Scholar 

  62. 62.

    Rainey F, Class II. Clostridia class. nov. In: De Vos P, Garrity G, Jones D, Krieg N, Ludwig W, Rainey F, Schleifer K-H, Whitman W, editors. Bergey’s manual of systematic bacteriology, vol. Volume 3. 2nd ed. New York: Springer; 2009. p. 736.

    Google Scholar 

  63. 63.

    Rainey F, Order I. Clostridiales Prévot 1953. In: De Vos P, Garrity G, Jones D, Krieg N, Ludwig W, Rainey F, Schleifer K-H, Whitman W, editors. Bergey’s manual of systematic bacteriology, vol. Volume 3. 2nd ed. New York: Springer; 2009. p. 736.

    Google Scholar 

  64. 64.

    Rainey F, Family V. Lachnospiraceae fam. nov. In: De Vos P, Garrity G, Jones D, Krieg N, Ludwig W, Rainey F, Schleifer K-H, Whitman W, editors. Bergey’s manual of systematic bacteriology, vol. Volume 3. 2nd ed. New York: Springer; 2009. p. 921–68.

    Google Scholar 

  65. 65.

    Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT. Gene ontology: tool for the unification of biology. Nat Genet. 2000;25(1):25–9.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  66. 66.

    Jensen LJ, Friis C, Ussery DW. Three views of microbial genomes. Res Microbiol. 1999;150(9):773–7.

    CAS  Article  PubMed  Google Scholar 

  67. 67.

    Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ, Gapped BLAST. PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997;25(17):3389–402.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  68. 68.

    Cantarel BL, Coutinho PM, Rancurel C, Bernard T, Lombard V, Henrissat B. The carbohydrate-active EnZymes database (CAZy): an expert resource for glycogenomics. Nucleic Acids Res. 2009;37(D1):D233–8.

    CAS  Article  PubMed  Google Scholar 

Download references


The MB2003 genome sequencing project was funded by the New Zealand Ministry of Business, Innovation and Employment New Economy Research Fund programme: Accessing the uncultured rumen microbiome, contract number C10X0803. Electron microscopy was conducted with the assistance of the Manawatū Microscopy and Imaging Centre at Massey University, Palmerston North, New Zealand.

Author information




NP, WJK, GTA conceived and designed the experiments. NP performed the sequencing and assembly experiments. NP, WJK, SCL, EA performed the genome annotation and comparative studies. NP performed the bacterial growth studies and polysaccharide utilization profiling. NP, WJK, GTA wrote the manuscript. All authors commented on the manuscript before submission. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Nikola Palevich.

Ethics declarations

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional file

Additional file 1: Table S1.

Associated MIGS record for B. hungatei MB2003, which links to the SIGS supplementary content website. (DOCX 17 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Palevich, N., Kelly, W.J., Leahy, S.C. et al. The complete genome sequence of the rumen bacterium Butyrivibrio hungatei MB2003. Stand in Genomic Sci 12, 72 (2017).

Download citation


  • Rumen
  • Bacteria
  • Hemicellulose
  • Pectin
  • Degradation
  • Butyrivibrio
  • Genome