Research Article
1 April 2000

Phylogenetic Relationships of Butyrate-Producing Bacteria from the Human Gut

ABSTRACT

Butyrate is a preferred energy source for colonic epithelial cells and is thought to play an important role in maintaining colonic health in humans. In order to investigate the diversity and stability of butyrate-producing organisms of the colonic flora, anaerobic butyrate-producing bacteria were isolated from freshly voided human fecal samples from three healthy individuals: an infant, an adult omnivore, and an adult vegetarian. A second isolation was performed on the same three individuals 1 year later. Of a total of 313 bacterial isolates, 74 produced more than 2 mM butyrate in vitro. Butyrate-producing isolates were grouped by 16S ribosomal DNA (rDNA) PCR-restriction fragment length polymorphism analysis. The results indicate very little overlap between the predominant ribotypes of the three subjects; furthermore, the flora of each individual changed significantly between the two isolations. Complete sequences of 16S rDNAs were determined for 24 representative strains and subjected to phylogenetic analysis. Eighty percent of the butyrate-producing isolates fell within the XIVa cluster of gram-positive bacteria as defined by M. D. Collins et al. (Int. J. Syst. Bacteriol. 44:812–826, 1994) and A. Willems et al. (Int. J. Syst. Bacteriol. 46:195–199, 1996), with the most abundant group (10 of 24 or 42%) clustering with Eubacterium rectale, Eubacterium ramulus, and Roseburia cecicola. Fifty percent of the butyrate-producing isolates were net acetate consumers during growth, suggesting that they employ the butyryl coenzyme A-acetyl coenzyme A transferase pathway for butyrate production. In contrast, only 1% of the 239 non-butyrate-producing isolates consumed acetate.
The species diversity of the predominantly anaerobic bacterial communities from the human large bowel has been the subject of both conventional and molecular microbiological investigations (31, 32, 46, 49). These communities are believed to contribute to healthy gut function in a variety of ways, including protecting against pathogens and producing nutrients for the colonic mucosa (3, 11, 12, 15). We still know relatively little, however, about the contributions of individual anaerobic species to colonic fermentation and to the nutrition and health of the host.
Diet-derived substrates, particularly undigested fiber and starch reaching the large intestine, have major effects upon bacterial community structure and metabolism in the colon. Short-chain fatty acids (SCFA) formed by microbial fermentation have an important effect on colonic health (10, 44). Butyrate in particular has an important role in the metabolism and normal development of colonic epithelial cells and has been implicated in protection against cancer and ulcerative colitis (21). Butyrate is preferentially transported by gut epithelial cells (36), serves as a preferred energy source for colonocytes (37, 43), and has been shown to exert direct effects upon gene expression in mammalian cells through histone hyperacetylation and through interaction with butyrate response elements upstream of some genes (8, 45). Production of butyrate by mixed human fecal microflora in vitro is known to be strongly influenced by the growth substrate. Starch, for example, is strongly butyrogenic, whereas other polysaccharides such as pectin result in relatively less butyrate and more propionate and acetate (9). Thus, the relative production rates of these SCFA provide a potentially important link between diet and colonic health. Despite this, however, remarkably little is known about the physiologies, identities, and ecologies of the predominant species of butyrate-producing bacteria from the human large bowel. The most obvious explanation for the stimulation of butyrate synthesis by certain carbohydrates in vitro is direct selection for butyrate-producing components of the flora capable of utilizing the particular substrate. Greater knowledge of these bacteria should lead to a more mechanistic understanding of the effects of diet on butyrate production in the colon. The isolation and molecular characterization of butyrate-producing species from fecal samples of three human volunteers deliberately chosen to represent different ages and diets is described in this paper.

MATERIALS AND METHODS

Bacterial strains.

The following bacterial reference strains were used in restriction fragment length polymorphism (RFLP) analysis:Butyrivibrio fibrisolvens 1.230 (bovine rumen [41]), 2221 (bovine rumen, ATCC 19171 [6]), 16.4 (human feces [38]),Eubacterium multiforme (ATCC 25552), and E. limosum R5004 and Fusobacterium strains R5043, R7263, and R10012 (all from J. Brazier, University of Wales, Cardiff, United Kingdom).

Collection and processing of samples.

Freshly voided fecal samples were collected from three healthy individuals, an infant (initially sampled at 11 months of age), a lacto-ovo-vegetarian adult (adult B, aged 46), and an omnivorous adult (adult A, aged 32), on two occasions approximately 1 year apart. None of the individuals had received antibiotics or other drugs in the months prior to sampling.
Fecal samples were placed in sterile universal bottles and processed within 30 min of collection. A subsample of feces (1 g) was aseptically transferred into a further sterile preweighed bottle. To the samples, 9 ml of M2GSC diluent (modified Hobson [30]), containing (per 100 ml) 1 g of casitone, 0.25 g of yeast extract, 0.4 g of NaHCO3, 0.2 g of glucose, 0.2 g of cellobiose, 0.2 g of soluble starch, 30 ml of clarified rumen fluid, 0.1 g of cysteine, 0.045 g of K2HPO4, 0.045 g of KH2PO4, 0.09 g of (NH4)2SO4, 0.09 g of NaCl, 0.009 g of MgSO4 · 7H2O, 0.009 g of CaCl2, and 0.1 mg of resazurin, was added aseptically while the bottle was being flushed with CO2 according to the anaerobic Hungate method (25). This diluent corresponded to the first 10-fold dilution, which was then mixed by vortexing for 3 min to form an evenly distributed suspension. The first dilution was subsequently diluted by 10-fold serial dilutions through to a 10−9 dilution.
Anaerobic roll tubes (5) were prepared in 16- by 125-mm Hungate tubes (25) sealed with butyl septum stoppers (Bellco Glass Inc., Vineland, N.J.) using medium M2GSC containing 0.75% agar (Difco) (30) and inoculated with 0.5-ml aliquots of appropriate serial dilutions. Undiluted aliquots (1 ml) were dispensed into 1.5-ml Eppendorf tubes and centrifuged (13,000 × g, 10 min) to pellet the bacteria, which were stored at −70°C prior to DNA extraction. Roll tubes were incubated at 37°C for 48 h prior to the picking of 30 to 80 colonies from each fecal sample. Cultures picked and grown in broths of M2GSC were used for examination of gram-stained smears (23, 24), for determination of fermentation products by capillary gas chromatography (35), and for extraction of DNA (see below).

DNA extraction.

DNA was extracted from bacterial pellets by following the method of Ausubel et al. (1). Bacterial pellets were resuspended in 567 μl of 1× Tris–EDTA buffer, pH 8.0. Three microliters of 20-mg/ml proteinase K and 30 μl of 10% (wt/vol) sodium dodecyl sulfate were added, and the solution was mixed and incubated for 1 h at 37°C. One hundred microliters of 5 M NaCl was added, followed by 80 μl of cetyltrimethylammonium bromide-NaCl solution (10% [wt/vol] cetyltrimethylammonium bromide, 0.7 M NaCl). The solutions were mixed and incubated for 10 min at 65°C. An equal volume of chloroform-isoamyl alcohol (24:1) was added, and the solution was mixed thoroughly and centrifuged for 10 min at 14,000 × g. The supernate was extracted twice with phenol-chloroform-isoamyl alcohol (25:24:1) and isopropanol precipitated. The DNA pellets were resuspended in sterile distilled water and stored at 4°C for immediate use or aliquoted at −20°C for long-term storage.

PCR amplification of 16S rDNA.

Fifty nanograms of total genomic DNA was used as a target for amplification of approximately 1,500 bp of 16S ribosomal DNA (rDNA) using the eubacterial primers fD1, 5′ AGAGTTTGATCCTGGCTCAG 3′ (Escherichia colipositions 8 to 27), and rP2, 5′ ACGGCTACCTTGTTACGACTT 3′ (E. coli positions 1494 to 1513) (47). PCR amplification was carried out as described previously (26).
Selected strains were sequenced. Sequencing was carried out using an automated ABI 377 sequencer. For the sequence reactions various universal primers (Table 1) and a Big Dye Ready Reaction DyeDeoxy Terminator Cycle Sequencing kit (Perkin-Elmer) were employed. All selected strains were sequenced in full (approximately 1,500 bases of 16S rDNA).
Table 1.
Table 1. Oligonucleotide primers used (positions are relative to those of E. coli) for 16S rDNA sequence analysis and ribotyping
PrimerSequencePositions relative to those of E. coli
fD1a5′ AGAGTTTGATCCTGGCTCAG 3′8–27
rP2a5′ ACGGCTACCTTGTTACGACTT 3′1513–1494
EBA-1b5′ GCCACATTGGGACTGAGACA 3′310–329
But-3b5′ GGAATTCCTAGTGTAGCGGT 3′673–692
Uniseq1c5′ ATGTTGGGTTAAGTCCCGCAAC 3′1082–1103
P3-Modd5′ ATTAGATACCCTDGTAGTCC 3′787–806
PC5d5′ TACCTTGTTACGACTT 3′1507–1487
a
Reference 47.
b
This work.
c
Reference 34.
d
Reference 49.

Nucleotide sequence accession numbers.

The 16S rDNA sequences of isolates used in the phylogenetic analysis have been deposited in the EMBL data library under accession numbers AJ70469 toAJ270492 . That of human B. fibrisolvens strain 16.4 wasAJ250365 (38), and that of rumen B. fibrisolvensstrain 1.230 was AJ270493 (46). The reference stains used in phylogenetic analysis were also from the EMBL database.

PCR-RFLP analysis.

PCR products were digested to completion with the enzymes HhaI, AluI, TaqI, andMspI and then analyzed by electrophoresis in a 2.0% (wt/vol) agarose gel (40).

Phylogenetic analysis.

The rDNA sequences corresponding toE. coli 16S rRNA bases 1 to 1500 (4) were compared directly with sequences in the EMBL and GenBank nonredundant nucleotide database using BLAST (19) and data of the Ribosomal Database Project (29). Database sequences with high similarity were then directly aligned over equalized lengths with the isolate sequences and used in phylogenetic analysis.
Sequences derived from previously cultured and described organisms of known phylogeny which corresponded to major subdivisions of the domainBacteria were included in the phylogenetic analysis of the fecal isolates. The sequence data approximating to E. colipositions 1 to 1500 were aligned using CLUSTAL V (22), and a phylogenetic tree was generated using software from the PHYLIP package (17). The DNADIST program analyzed distances using the Kimura-Nei correction (27) trees generated from distance matrices that employed the neighbor-joining method (39). Sequence data for distance matrices and analysis were subjected to bootstrap resampling (data resampled 100 times) using the SEQBOOT program, and consensus trees were generated by the CONSENSE program (16).

RESULTS

Isolation of butyrate producers.

Bacteria were isolated from anaerobic roll tubes inoculated with freshly voided fecal samples from the three subjects. Each individual was sampled on two occasions, separated by approximately 1 year. A rumen fluid-based medium which has been shown to support growth of a wide range of anaerobic bacteria (14, 30) was used. A total of 313 colonies (20 to 30 from the first isolation and 80 from the second isolation) picked from the highest dilutions (10−8 or 10−9) were purified, grown in broth culture, and tested for the production of volatile fatty acids. Isolates producing detectable butyrate varied widely with respect to the quantity of butyrate produced in batch culture (Table 2). A net value of 2 mM butyrate was chosen as a cutoff, since this was considered the lowest value that could be distinguished unambiguously from butyrate concentrations present in uninoculated media. The proportion of isolates producing more than 2 mM butyrate ranged from 3 to 76% for the six samples analyzed. The largest proportion of high-concentration-butyrate-producing isolates was from the first infant sample, with around 50% of isolates producing >10 mM butyrate. The vegetarian adult fecal samples yielded the smallest proportion of butyrate producers and yielded no strains producing >5 mM butyrate. Net utilization of acetate (by >2 mM) present in the rumen fluid component of the medium was detected in 50% of butyrate-producing isolates. Comparison of the amounts of acetate utilized with butyrate produced gave a positive relationship, a correlationr2 value of 0.6, and a slope of 0.51 ± 0.071 (mean ± standard deviation) (Fig.1). Thus, by averaging across all acetate-consuming strains, approximately 1 mol of acetate initially present in the medium was consumed for every 2 mol of butyrate produced. In marked contrast, only 2 out of 239 non-butyrate-producing strains showed net acetate consumption (Table 2).
Table 2.
Table 2. Occurrence of butyrate-producing and acetate-consuming strains in six human fecal samples
IndividualcSampleButyrate production (mM) among human fecal isolates numbering:No. (%)a of strains showing net consumption of acetateb among strains producing:
0–2.02.1–5.05.1–10.010.1–15.0>15.0>2 mM butyrate<2 mM butyrate
Infant15054715/160/5
 26394407/172/63
Adult A115111005/120/15
 255885310/240/55
Adult B12520000/20/25
 27630000/30/76
 Total 2393318131037/74 (50)2/239 (1)
a
Number of strains relative to total numbers. N.B., most (90%) strains not showing acetate consumption produced more than 2 mM acetate.
b
Defined as >2 mM acetate, net.
c
Adult A is omnivorous; adult B is vegetarian.
Fig. 1.
Fig. 1. Summary of butyrate synthesis and acetate utilization in a selection of human colonic butyrate-producing isolates. Key to strains: 1, A1-86; 2, A1-189; 3, A1-816; 4, L2-6; 5, L2-21; 6, A2-204; 7, A2-207; 8, L2-15; 9, A1-87; 10, A1-815; 11, L2-10; 12, A2-225; 13, L1-81; 14, L1-83; 15, A2-194; 16, L2-7; 17, L1-92; 18, A2-215; 19, L1-93; 20, L2-39; 21, L2-61; 22, L1-911; 23, A2-223; 24, L1-872; 25, A2-173; 26, A2-165; 27, L1-910; 28, A2-227; 29, L1-94; 30, L1-9171; 31, L1-8151; 32, L1-97; 33, A2-181; 34, L1-91; 35, L1-82; 36, A2-183; and 37, L1-952.

Ribotypes of the butyrate producers based on PCR-RFLP analysis.

In order to group similar isolates, all 74 butyrate-producing isolates (30 from the first set of samples and 44 from the second set of fecal isolations, taken approximately 1 year later), were examined by 16S rDNA ribosomal profiling and compared with reference strains available within the laboratory (see Materials and Methods). The 16S rDNAs were amplified by PCR using universal primers (47), and the amplified material was cleaved initially with the restriction enzymes HhaI, AluI,TaqI, and MspI (33). The RFLP profiles obtained with the four restriction enzymes indicated that the enzymeAluI allowed the best discrimination between the isolates, distinguishing 18 different ribotypes. The collected results from ribotyping for all six fecal isolations are shown in Table3. It is apparent that samples from different individuals and from different sampling times show very little overlap with respect to the most abundant ribotypes present. Although none of the 18 ribotypes observed had a profile identical to that of any of the type strains, some of the isolates exhibited bands identical in size with those of type strains. One of the main identifiable bands is at 610 bp, which is characteristic ofEubacterium and Butyrivibrio species (Fig.2).
Table 3.
Table 3. Distribution of the 74 butyrate-producing strains isolated in this study in each of the 18 PCR-RFLP ribotypes obtained with the restriction enzyme AluI
IndividualSampleNo. of butyrate-producing strains from the indicated isolation in RFLP group:Total
123456789101112131415161718
Infant11411               16
 2  9   4211        17
Adult A1   12              12
 2   3  5   423421  24
Adult B1    11            2
 2             1  113
Fig. 2.
Fig. 2. PCR-RFLP profiles of a selection of human butyrate-producing isolates and standard strains digested with the restriction endonuclease AluI. (A) Lanes 1 to 14 contain human colonic isolates A2-165 (ribotype 7), A2-168 (ribotype 11), A2-166 (ribotype 11), A2-228 (ribotype 14), A2-175 (ribotype 14), A2-178 (ribotype 14), T2-132 (ribotype 14), A1-89 (ribotype 4), A2-171 (ribotype 4), A2-181 (ribotype 4), A2-183 (ribotype 4), T2-87 (ribotype 17), T2-145 (ribotype 18), and T1-815 (ribotype 5), respectively. (B) Lanes 1 to 12 contain human colonic isolates L2-6 (ribotype 7), L2-12 (ribotype 10), L2-10 (ribotype 9), L2-21 (ribotype 3b), L1-92 (ribotype 3a), L2-9 (ribotype 3b), L2-16 (ribotype 3b), L1-81 (ribotype 1), L2-7 (ribotype 8), L2-65 (ribotype 8), L1-83 (ribotype 2), and B. fibrisolvens 16.4, respectively. Lanes 13 and 14 contain rumen isolates B. fibrisolvens 1.230 and B. fibrisolvens 2221 (ATCC 19971), respectively. Size markers (1-kb ladder; Promega) are shown in lanes marked “M.” The 610-bp band in panel A, lanes 4 to 7, and in panel B, lanes 4, 7, and 8 and the 400-bp fragment in panel B, lane 11, represent partial products. Ribotype groups 6, 12, 13, 15, and 16 are not represented in this figure. TheirAluI banding profiles (in base pairs) are 610, 475, 210, 190, and 50 (ribotype 6); 610, 240, and 220 (ribotype 12); 475, 265, 220, and 190 (ribotype 13); 400, 265, 220, 190, and 127 (ribotype 15); and 797, 610, 265, and 190 (ribotype 16). The band pattern for ribotype 6 is similar to that in panel A, lanes 8 to 11, while that for ribotype 12 is similar to that in panel A, lanes 4 to 7. The band pattern for ribotype 13 is similar to that in panel A, lane 12, while that for ribotype 15 is similar to that in panel B, lane 13. The band pattern for ribotype 16 is similar to that in panel A, lanes 4 to 7.
The 25 strains of ribotypes 1, 2, and 3 were found uniquely in the flora from the infant, and ribotype 1 accounted for 14 of the 16 strains isolated from the first (preweaning) infant sample. Only two ribotypes (ribotypes 7 and 14) were recovered from more than one individual, and only four ribotypes (ribotypes 3, 4, 7, and 14) were recovered from more than one sample.

Phylogenetic comparison.

Complete 16S rDNA sequences were determined for 24 isolates representative of the most common ribotypes identified by RFLP analysis. These included strains producing large (>10 mM), medium (>5 and <10 mM), and small (<5 mM) amounts of butyrate in vitro (Table 4). These sequences were compared with available database sequences by a BLAST search (Table 4), and a phylogenetic tree was constructed (Fig.3). Twenty of the sequences were found to belong within cluster XIVa of gram-positive bacteria. A few isolates were shown to be closely related (>97% sequence identity) toEubacterium rectale (A1-86, L2-21, and T1-815),Eubacterium ventriosum (L2-12), or Fusobacterium prauznitzii (A2-165 and L2-6). The remaining 18 of the 24 sequenced isolates however shared 95% or less 16S rDNA sequence identity with their nearest relative.
Table 4.
Table 4. Relationships and fermentation products of the butyrate-producing strains shown in Fig. 3a
StrainRibotype groupSpecies identified% HomologybSCFA production (mM)H2 production (μ mol ml−1)
ButyrateAcetateLactate
L1-811UnidentifiedNA5.7−2.01.45.0
L1-8101UnidentifiedNA5.541.025.50.1
L1-9521UnidentifiedNA21.3−10.07.16.5
L1-931UnidentifiedNA9.8−10.65.7 NDc
L1-832UnidentifiedNA6.5−2.32.12.7
L1-923aUnidentifiedNA8.3−3.90.112.0
L2-213aEubacterium rectale983.6−3.6NDND
L2-503bUnidentifiedNA12.711.10.974.3
A1-864Eubacterium rectale984.7−2.126.53.2
A2-1834Roseburia cecicola9520.0−8.29.714.8
T1-8155Eubacterium rectale983.30.718.72.6
T1-8176Bifidobacterium infantis orlongum982.17.29.22.1
L2-67Fusobacterium prauznitzii983.3−4.8NDND
A2-1657Fusobacterium prauznitzii9814.4−7.42.9  0
L2-78Eubacteriumspp.977.9−2.2NDND
L2-109UnidentifiedNA4.9−4.3NDND
L2-1210Eubacterium ventriosum984.6−1.7NDND
A2-16611UnidentifiedNA9.72.4NDND
A2-20712Eubacteriumspp.964.6−2.001.7
A2-23113UnidentifiedNA5.65.2NDND
A2-17514UnidentifiedNA5.0−1.8NDND
T2-13214UnidentifiedNA2.247.66.70.1
A2-19415UnidentifiedNA7.3−3.57.1  0
T2-14518UnidentifiedNA2.71.0NDND
a
Data for SCFA (butyrate and acetate) and for production of lactate and H2 were from two separate experiments. None of the strains produced propionate, with the possible exception of L1-83 (3 mM).
b
Percent homology indicates the maximum base sequence identity obtained in database searches. NA, not applicable, indicates no value above 95% identity.
c
ND, not determined.
Fig. 3.
Fig. 3. Phylogenetic tree showing the relationships of 16S rDNA sequences from human butyrate-producing isolates falling within cluster XIVa of the Clostridium subphylum of low-G+C-content gram-positive bacteria (7, 48). The scale bar represents genetic distance (10 substitutions per 100 nucleotides). The tree was constructed using the neighbor-joining analysis of a distance matrix obtained from a multiple-sequence alignment. Bootstrap values (expressed as percentages of the value for 100 replications) are shown at branch points; values of 97% or more were considered significant. Sequences derived from the database are shown in italics (e.g.,B. fibrisolvens). Atopobium minutum is used as the outgroup sequence.
The largest group of 10 newly isolated strains clusters withEubacterium ramulus, Eubacterium rectale, andRoseburia cecicola. This group includes all sequenced strains belonging to ribotypes 1 and 4, together with one ribotype 3 strain and strains of ribotypes 5 and 15. Eight additional strains showed relationships to other Eubacterium species, includingEubacterium halii, Eubacterium uniforme,Eubacterium formicigenerans, and Eubacterium ventriosum, or to other members of the XIVa cluster, including two strains (L2-50 and A2-166) that are related to Coprococcus eutactus. The strains that lie outside the XIVa cluster were still related to gram-positive bacteria. These included A2-207 and the two strains (L2-6 and A2-165) belonging to ribotype 7 that are close toF. prauznitzii, as noted above. The remaining strain, T1-817, was related to Bifidobacterium longum, which is not expected to produce butyrate, but the amounts of butyrate produced by this strain were extremely small.

DISCUSSION

To our knowledge, this is the first study directed specifically at establishing the identities of the butyrate-producing bacteria in the human gut. In one respect the diversity of the butyrate-producing strains as revealed by phylogenetic analysis of 16S rDNA sequences was less than might have been anticipated. All isolates were related to gram-positive bacteria, 80% of which belonged to cluster XIVa, and none were related to gram-negative bacteria, since the few F. prauznitzii-like isolates reported here are in fact also related to gram-positive bacteria at the genotypic level. On the other hand, the RFLP analyses revealed extraordinary diversity and variability of the major strain types recovered between the different individuals and sampling times. Much more extensive work would be required to establish whether this diversity reflects differences in diet, age, or genotype between the individuals tested or simply large fluctuations in the predominant types within individuals between different times of sampling. While it is likely that the bacteria isolated include many of the predominant butyrate producers from the human colon, further work will be required to establish the occurrence of related strains in a larger number of individuals. This could be approached by using specific molecular probes to avoid possible cultural bias.
Phylogenetic analysis shows that the butyrate-producing strains isolated that belong to cluster XIVa are widely distributed within this cluster (Fig. 2). The majority of isolates appear to be broadly related to known species, but only 8 out of the 24 sequenced strains had greater than 95% 16S rDNA identity with type strains. This and the lack of recognized species type strains in some branches of the tree indicate that new uncharacterized phylogenetic groups are represented among our isolates. Of the three clusters, the Eubacterium rectale-R. cecicola group comprises the highest proportion of sequenced isolates (10 out of 24, 42%) while overall ribotypes 1 and 4, whose sequenced representatives fall within this group, account for 27 out of the 74 (36%) isolates examined. Within the main grouping theEubacterium rectale subgroup appears to contain isolates from all of the individuals sampled whereas the R. cecicolasubgroup was found mainly in infant samples. It would be of particular interest to establish in future work whether R. cecicola-related strains belonging to ribotype 1 are characteristic of preweaned infants.
A number of molecular approaches are now available for rapid bacterial strain typing. RFLP analyses based on hybridization with 16S rDNA probes have been widely used and depend on detection of variations in chromosomal sequences flanking multiple copies of rRNA target genes (2, 28). The approach used here is quite different in that it depends on PCR amplification of internal 16S rDNA fragments and cleavage with restriction enzymes (50). This approach proved very valuable here in revealing the diversity of butyrate-producing strains, and two ribotypes, 1 and 7, conformed well to branches within the phylogenetic tree obtained subsequently from sequence analysis. In general, however, it must be recognized that PCR-RFLP analysis based on a single enzyme will often group together genetically distant strains since extensive sequence changes can occur between strains without affecting sites for a particular restriction enzyme. The RFLP classes must therefore be seen as a useful comparative method for examining diversity rather than as a definitive phylogenetic typing approach, unless results from many enzymes are combined (33).
Fifty percent of butyrate-producing isolates examined here showed some net utilization of acetate in M2GSC medium. Such utilization is consistent with the operation of the butyryl coenzyme A-acetyl coenzyme A transferase route for butyrate synthesis (20). On average across the strains studied in vitro here, approximately 1 mol of acetate initially present in the medium was utilized for every 2 mol of butyrate produced. We were not able to estimate how much additional acetate may have been produced and utilized by each strain for butyrate synthesis. Some acetate may, of course, be used for biosynthetic reactions rather than butyrate synthesis. The high proportion of acetate-utilizing strains seen here suggests that a significant fraction of colonic acetate may be rerouted into butyrate in the human colon. Preliminary studies involving 13C-labeled acetate suggest that appreciable butyrate may be derived from exogenous acetate in mixed fecal fermentations (S. H. Duncan et al., unpublished data). Since for some of the strains isolated here growth was absolutely dependent on inclusion of acetate in the medium (S. H. Duncan et al., unpublished data), these strains would not have been recovered through isolations with media lacking acetate. The alternative butyrate kinase pathway is known to operate in manyClostridium spp. and in certain ruminal B. fibrisolvens strains, including the type strain 2221, and strains depending on this pathway would not be expected to utilize acetate (13, 20). It is also worth noting that only 1% of the non-butyrate-producing isolates examined showed net acetate utilization and that nearly all (95%) of the acetate utilizers recovered were butyrate producers.
Recent work by Diez-Gonzalez et al. (13) suggested that rumen butyrate-producing Butyrivibrio species could be divided into two distinct groups, based upon pathways that they utilized to produce butyrate and that these could be distinguished by production or consumption of acetate. Rumen Butyrivibriostrains have also been classified as lactate producing or non-lactate producing (42). None of the human isolates obtained in this study fell within either group of Butyrivibrio strains. However, the present work does confirm that one isolate, B. fibrisolvens 16.4, obtained previously from fermentor studies with the human fecal flora (38), is related to the group of ruminal B. fibrisolvens strains that includes B. fibrisolvens 1.230 and 2223.
We found no simple correlation between the phylogenetic position of the butyrate-producing isolates examined here and their metabolic behavior. Thus, for example, ribotype 1 strains included one acetate producer (L1-810) in addition to the predominant acetate-utilizing strains. In addition, we were able to detect significant production of lactate both in acetate-producing strains (e.g., L1-810) and in acetate-utilizing strains (e.g., L1-952) from human feces. On the other hand it is worth noting that ribotypes 1 and 7 contained particularly high proportions of strains that produced >10 mM butyrate in vitro (11 out of 14 [78%] and 6 out of 9 [67%], respectively, compared to 24 out of 74 [32%] butyrate producers over all ribotypes).
Recent investigations by fluorescent in situ hybridization with group-specific probes (18) showed that the highest proportion of human fecal organisms detected fell within theClostridium coccoides-Eubacterium rectale group (7.2 × 1010 cells/g [dry weight] of feces), which forms part of clostridial cluster XIVa. We have shown here that most of the butyrate-producing isolates in this study fall into cluster XIVa and that they are related to the Clostridium coccoides-Eubacterium rectale group. Further understanding of the phylogeny and physiology of this group of organisms will be crucial in understanding the role of the anaerobic microflora in colonic metabolism and gut health.

ACKNOWLEDGMENTS

We thank Tony Richardson and Peter Dewey for their contributions to the SCFA analyses and Moira Johnston for DNA sequencing.
This work was supported by SERAD (Scottish Executive Rural Affairs Department), by a BORC (Boyd Orr Research Consortium) Ph.D. studentship to A.B., and by a SERAD flexible fund grant.

REFERENCES

1.
Ausubel F. M., Brent R., Kingston R. E., Moore D. M., Seidman J. G., Smith J. A., and Struhl K. Current protocols in molecular biology I 1994 2.4 J. Wiley and Sons, Inc. New York, N.Y
2.
Avgustin G., Wright F., and Flint H. J. Genetic diversity and phylogenetic relationships among strains of Prevotella (Bacteroides) ruminicola from the rumen. Int. J. Syst. Bacteriol. 47 1994 284 -288
3.
Bengmark S. Econutrition and health maintenance—a new concept to prevent GI inflammation, ulceration and sepsis. Clin. Nutr. 15 1996 1 -10
4.
Brosius J., Palmer M. L., Kennedy P. J., and Noller H. F. Complete nucleotide sequence of a 16S ribosomal gene from Escherichia coli. Proc. Natl. Acad. Sci. USA 75 1978 4801 -4805
5.
Bryant M. P. Commentary on the Hungate technique for cultivation of anaerobic bacteria. Am. J. Clin. Nutr. 25 1972 1324 -1328
6.
Bryant M. P. and Small N. The anaerobic monotrichous butyric acid-producing curved rod-shaped bacteria of the rumen. J. Bacteriol. 72 1956 16 -21
7.
Collins M. D., Lawson P. A., Willems A., Cordoba J. J., Fernandez-Garayzabal J., Garcia P., Cai J., Hippe H., and Farrow J. A. E. The phylogeny of the genus Clostridium: proposal of five new genera and eleven new species combinations. Int. J. Syst. Bacteriol. 44 1994 812 -826
8.
Csordas, A. 1995. Toxicology of butyrate and short-chain fatty acids, p. 105–125. In M. J. (ed.), Role of gut bacteria in human toxicology and pharmacology. Taylor and Francis, London, United Kingdom.
9.
Cummings J. H., Pomare E. W., Branch W. J., Naylor C. P. E., and MacFarlane G. T. Short-chain fatty acids in human large intestine, portal, hepatic and venous blood. Gut 28 1987 1221 -1227
10.
Cummings J. H. and MacFarlane G. T. The colonic flora, fermentation and large bowel digestive function The large intestine: physiology, pathophysiology, and disease. Philips S. F., Pemberton J. H., and Shorter R. G. 1991 51 -91 Raven Press New York, N.Y
11.
Cummings J. H. and MacFarlane G. T. Colonic microflora: nutrition and health. Nutrition 13 1997 476 -478
12.
Cummings J. H. and MacFarlane G. T. Role of intestinal bacteria in nutrient metabolism. Clin. Nutr. 16 1997 3 -11
13.
Diez-Gonzalez F., Bond D. R., Jennings E., and Russell J. B. Alternative schemes of butyrate production in Butyrivibrio fibrisolvens and their relationship to acetate utilization, lactate production and phylogeny. Arch. Microbiol. 171 1999 324 -330
14.
Duncan S. H., Flint H. J., and Stewart C. S. Inhibitory activity of gut bacteria against Escherichia coli O157 mediated by dietary plant metabolites. FEMS Microbiol. Lett. 164 1998 283 -288
15.
Engelhardt W., Bartels J., Kirschberger S., Meyer zu Duttingdof H. D., and Busche R. Role of short-chain fatty acids in the hind gut. Vet. Q. 20 1998 S52 -S59
16.
Felsenstein J. Confidence limits on phylogenies: an approach using the bootstrap. Evolution 39 1985 783 -791
17.
Felsenstein J. PHYLIP—Phylogeny Inference Package (version 3.2). Cladistics 5 1989 164 -166
18.
Franks A. H., Harmsen H. J. M., Raangs G. C., Jansen G. J., and Welling G. W. Variations of bacterial populations in human feces quantified by fluorescent in situ hybridization with group-specific 16S rRNA-targeted oligonucleotide probes. Appl. Environ. Microbiol. 64 1998 3336 -3345
19.
Gish W. and States D. J. Identification of protein coding regions by database similarity search. Nat. Genet. 3 1993 266 -272
20.
Gottschalk G. Bacterial fermentations Bacterial metabolism. Starr M. P. 1979 167 -224 Springer-Verlag New York, N.Y
21.
Hague A., Singh B., and Paraskeva C. Butyrate acts as a survival factor for colonic epithelial cells: further fuel for the in vivo versus in vitro debate. Gastroenterology 112 1997 1036 -1040
22.
Higgins D. G., Bleasby A. J., and Fuchs R. CLUSTALV: improved software for multiple sequence alignment. Comput. Appl. Biosci. 8 1992 189 -191
23.
Hobson P. N. Rumen bacteria. Methods Microbiol. 3B 1969 133 -149
24.
Holdeman L. V., Cato E. P., and Moore W. E. C. Anaerobe laboratory manual 4th ed. 1977 Virginia Polytechnic Institute and State University Blacksburg
25.
Hungate R. E. The rumen and its microbes. 1966 Academic Press New York, N.Y
26.
Innis M. A., Gelfand D. H., Sninsky J. J., and White T. J. PCR products: a guide to methods and applications. 1990 Academic Press, Inc. London, United Kingdom
27.
Kimura M. A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J. Mol. Evol. 16 1980 111 -120
28.
Lawson P. A., Gharbia S. E., Shah H. N., and Clark D. R. Intrageneric relationships of members of the genus Fusobacterium as determined by reverse transcriptase sequencing of small-subunit rRNA. Int. J. Syst. Bacteriol. 41 1991 347 -354
29.
Maidak B. L., Larsen N., McCaughey M. J., Overbeek R., Olsen G. J., Forgel K., Blandy J., and Woese C. R. The Ribosomal Database Project. Nucleic Acids Res. 22 1994 3485 -3487
30.
Miyazaki K., Martin J. C., Marinsek-Logar R., and Flint H. J. Degradation and utilization of xylans by the rumen anaerobe Prevotella bryantii (formerly P. ruminicola subsp. brevis) B14. Anaerobe 3 1997 373 -381
31.
Moore W. E. C. and Holdeman L. V. Human fecal flora: the normal flora of 20 Japanese-Hawaiians. Appl. Environ. Microbiol. 27 1974 961 -979
32.
Moore W. E. C. and Moore L. H. Intestinal floras of populations that have a high risk of colon cancer. Appl. Environ. Microbiol. 61 1995 3202 -3207
33.
Moyer C. L., Dobbs F. C., and Carl D. M. Estimation of diversity and community structure through restriction fragment length polymorphism distribution analysis of bacterial 16S rDNA genes from a microbial mat at an active, hydrothermal vent system, Loihi Seamount, Hawaii. Appl. Environ. Microbiol. 60 1996 871 -879
34.
Pryde S. E., Richardson A. J., Stewart C. S., and Flint H. J. Molecular analysis of the microbial diversity present in the colonic wall, colonic lumen, and cecal lumen of a pig. Appl. Environ. Microbiol. 65 1999 5372 -5377
35.
Richardson A. J., Calder A. G., Stewart C. S., and Smith A. Simultaneous determination of volatile and non-volatile acidic fermentation products of anaerobes by capillary gas chromatography. Lett. Appl. Microbiol. 9 1989 5 -8
36.
Ritzhaupt A., Ellis A., Hosie K. B., and Shirazi-Beechey S. P. The characterization of butyrate transport across pig and human colonic luminal membrane. J. Physiol. 507 1998 819 -830
37.
Roediger W. E. W. The colonic epithelium in ulcerative colitis: an energy deficient disease? Lancet ii 1980 712 -715
38.
Rumney C., Duncan S. H., Henderson C., and Stewart C. S. Isolation and characteristics of a wheatbran degrading butyrivibrio from human faeces. Lett. Appl. Microbiol. 20 1995 232 -236
39.
Saitou N. and Nei M. The neighbor joining method: a new method for constructing phylogenetic trees. Mol. Biol. Evol. 4 1987 406 -425
40.
Sambrook J., Fritsch E. F., and Maniatis T. Molecular cloning: a laboratory manual 2nd ed. I 1989 5.30 -5.32 Cold Spring Harbor Laboratory Press Cold Spring Harbor, N.Y
41.
Scott K. P., Barbosa T. M., Forbes K. J., and Flint H. J. High-frequency transfer of naturally occurring chromosomal tetracycline resistance element in the ruminal anaerobe Butyrivibrio fibrisolvens. Appl. Environ. Microbiol. 6 1997 3405 -3411
42.
Shane B. S., Gouws L., and Kistner A. Cellulolytic bacteria occurring in the rumen of sheep conditioned to low-protein teff hay. J. Gen. Microbiol. 55 1969 445 -457
43.
Sheppach W., Bartram H. P., and Richter F. Role of short-chain fatty acids in the prevention of colorectal cancer. Eur. J. Cancer 31A 1998 107 -1080
44.
Szylit O. and Andrieux C. Physiological and pathophysiological effects of carbohydrate fermentation. World Rev. Nutr. Diet. 74 1993 88 -122
45.
Tran C. P., Familari M., Parker L. M., Whitehead R. H., and Giraud A. S. Short-chain fatty acids inhibit intestinal trefoil factor gene expression in colon cancer cells. Am. J. Physiol.-Gastrointest. Liver Physiol. 38 1998 G85 -G93
46.
Wang R. F., Cao W. W., and Cernaglia C. E. PCR detection and quantitation of predominant and anaerobic bacteria in human and animal fecal samples. Appl. Environ. Microbiol. 62 1996 1242 -1247
47.
Weisburg W. G., Barns S. M., Pelletier D. A., and Lane D. J. 16S ribosomal DNA for phylogenetic study. J. Bacteriol. 173 1991 697 -703
48.
Willems A., Amat-Marco M., and Collins M. D. Phylogenetic analysis of Butyrivibrio strains reveals three distinct groups of species within the Clostridium subphylum of the gram-positive bacteria. Int. J. Syst. Bacteriol. 46 1996 195 -199
49.
Wilson K. H. and Blitchington R. B. Human colonic biota studied by ribosomal DNA sequence analysis. Appl. Environ. Microbiol. 62 1996 2273 -2278
50.
Wood J., Scott K. P., Avgustin G., Newbold C. J., and Flint H. J. Estimation of the relative abundance of different Prevotella ribotypes in gut samples by restriction enzyme profiling of PCR-amplified 16S rRNA gene sequences. Appl. Environ. Microbiol. 64 1998 3683 -3689

Information & Contributors

Information

Published In

cover image Applied and Environmental Microbiology
Applied and Environmental Microbiology
Volume 66Number 41 April 2000
Pages: 1654 - 1661
PubMed: 10742256

History

Received: 18 November 1999
Accepted: 28 January 2000
Published online: 1 April 2000

Permissions

Request permissions for this article.

Contributors

Authors

Adela Barcenilla
Rowett Research Institute, Bucksburn, Aberdeen AB21 9SB,1 and
Susan E. Pryde
Rowett Research Institute, Bucksburn, Aberdeen AB21 9SB,1 and
Jennifer C. Martin
Rowett Research Institute, Bucksburn, Aberdeen AB21 9SB,1 and
Sylvia H. Duncan
Rowett Research Institute, Bucksburn, Aberdeen AB21 9SB,1 and
Colin S. Stewart
Rowett Research Institute, Bucksburn, Aberdeen AB21 9SB,1 and
Colin Henderson
Robert Gordon University, Kepplestone Campus, Aberdeen AB9 2PG,2 United Kingdom
Harry J. Flint
Rowett Research Institute, Bucksburn, Aberdeen AB21 9SB,1 and

Metrics & Citations

Metrics

Note:

  • For recently published articles, the TOTAL download count will appear as zero until a new month starts.
  • There is a 3- to 4-day delay in article usage, so article usage will not appear immediately after publication.
  • Citation counts come from the Crossref Cited by service.

Citations

If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. For an editable text file, please select Medlars format which will download as a .txt file. Simply select your manager software from the list below and click Download.

View Options

Figures

Tables

Media

Share

Share

Share the article link

Share with email

Email a colleague

Share on social media

American Society for Microbiology ("ASM") is committed to maintaining your confidence and trust with respect to the information we collect from you on websites owned and operated by ASM ("ASM Web Sites") and other sources. This Privacy Policy sets forth the information we collect about you, how we use this information and the choices you have about how we use such information.
FIND OUT MORE about the privacy policy