Previous Article | Next Article ![]()
Journal of Bacteriology, April 2003, p. 2571-2581, Vol. 185, No. 8
0021-9193/03/$08.00+0 DOI: 10.1128/JB.185.8.2571-2581.2003
Copyright © 2003, American Society for Microbiology. All Rights Reserved.
Teagasc, Dairy Products Research Centre, Moorepark, Fermoy, County Cork,1 Department of Microbiology, University College Cork, Cork, Ireland2
Received 20 November 2002/ Accepted 3 February 2003
|
|
|---|
|
|
|---|
At present, the genus Bifidobacterium includes 34 species (5, 16). Approximately 20 species have been identified from fecal sources, of which only 10 appear to have been recovered from GIT samples (5, 24, 29). In addition, three species, Bifidobacterium ruminatium and Bifidobacterium merycicum from bovine rumen (4) and Bifidobacterium gallinarum from chicken cecum (51), do not appear to have been recovered from fecal samples. Therefore, it seems likely that the examination of GIT samples will identify more Bifidobacterium species. Indeed, a potential new species, based on 16S ribosomal DNA (rDNA) sequence comparisons, was isolated from the porcine intestine (24), and recently, bifidobacteria were recovered from the chicken crop, a habitat not considered suitable for the growth of strict anaerobes (31).
Bifidobacteria have been cultured from the cecal contents of chickens (31, 35, 51), rats (43), rabbits (20, 35, 46), pigs (37), and humans (21). Isolates were confirmed as belonging to the genus Bifidobacterium in all but the rat and pig studies, through the detection of fructose-6-phosphate phosphoketolase (F-6-PPK) (41). However, only isolates from chicken and rabbit cecum were characterized to the species level. In the former a new species, B. gallinarum, was proposed (51), and in the latter two subspecies, Bifidobacterium pseudolongum subsp. pseudolongum and Bifidobacterium pseudolongum subsp. globosum, were identified (20, 46). Although none of the studies included an assessment of genomic diversity within the bifidobacterial population, the use of the restriction enzyme XbaI followed by pulsed-field gel electrophoresis (PFGE) was reported to be an effective method for discriminating human fecal Bifidobacterium isolates (17, 23).
In addition, a variety of molecularly based techniques have been used to establish the genetic relatedness and species identity of Bifidobacterium isolates. Randomly amplified polymorphic DNA (RAPD) PCR was used to group bifidobacteria recovered from rat intestines (13), human feces (50), and human GIT samples (22). Whole-cell protein profiling (WCPP) (3, 54), 16S rDNA sequencing (11, 16, 26), and DNA-DNA hybridization (39, 47) have been used to identify Bifidobacterium species.
Bifidobacterium strains are reported to vary in their ability to adhere to cultured epithelial cells, and this has been attributed to specific protein, polysaccharide, lipoteichoic acid, hydrophobic surface, and autoaggregation factors (2, 6, 9, 10, 19, 30). In addition, studies on the general microflora of the epithelium and contents of the porcine GIT have suggested that distinct bacterial populations exist within each region (1, 44). However, to our knowledge no study has determined the genomic diversity of bifidobacteria isolated from the contents and epithelium of the same GIT sample.
The present study had two aims: to determine the level of genomic diversity of Bifidobacterium isolates recovered from different fractions of a porcine cecum, using PFGE; and to establish genetic relatedness and species identity of isolates based on PFGE, cell morphology, RAPD PCR, WCPP, 16S rDNA sequencing, and DNA-DNA hybridization.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Bifidobacterium strains used in this study
|
8 cm2, was aseptically removed, washed with MRD to remove loosely adhering contents, weighed, and diluted accordingly. To obtain weakly adherent bacteria, the tissue sample was gently shaken for 5 min before serial dilution. To obtain strongly adherent bacteria, the same piece of tissue was reweighed, resuspended with fresh MRD, and vigorously shaken using a Stomacher 400 (Seward, London, United Kingdom), set at high speed, for 5 min before serial dilution (sample III). Diluents from each cecal sample were pour plated with mMRS with and without 50 mg of mupirocin/liter and incubated for 2 days at 37°C under anaerobic conditions. From each sample, 50 to 60 colonies were randomly picked, inoculated into mMRS broth, and incubated as before for 48 h, from which 1-ml samples were taken for the preparation of frozen stocks, detection of F-6-PPK activity, and PFGE preparations. Selected isolates were deposited in the Dairy Products Centre (DPC) Collection at Teagasc, Moorepark, Fermoy, County Cork, Ireland. F-6-PPK activity. A protocol for detecting F-6-PPK activity in the intracellular extracts from 1 ml of mMRS broth cultures using Triton X-100 detergent-mediated cell lysis was employed based on protocols previously described (7, 28). Bifidobacterium thermophilum NCIMB 702553 and Lactobacillus rhamnosus ATCC 7469 were used as positive and negative controls, respectively.
PFGE. Preparation of high-molecular-weight DNA from mMRS broth cultures was as previously described (45) for pediococci except that 40 µg of mutanolysin/ml was added to the cell lysis buffer and the restriction enzymes XbaI and SpeI were used (New England Biolabs, Beverly, Mass.). DNA fragments were resolved using a CHEF-DR III pulsed-field system (Bio-Rad Laboratories, Richmond, Calif.) at 6 V/cm for 18 h with a 1- to 15-s linear ramp pulse time. A 0.5x Tris base-borate-EDTA buffer was maintained at 14°C during electrophoresis. Molecular size markers were included in each gel (#N0350S; New England Biolabs). Gels were stained in distilled water containing 0.5 µg of ethidium bromide/ml for 30 min and destained for 60 min in distilled water.
RAPD PCR. RAPD PCR was performed on DNA from selected pig cecal isolates. Details on the isolation of DNA from mMRS broth cultures, random primers P1 and P2, RAPD PCR conditions, and RAPD pattern comparisons using Gelcompar software were as previously described (45).
Whole-cell protein profile. Whole-cell protein extraction, sodium dodecyl sulfate-polyacrylamide gel electrophoresis, and computer-based comparisons of the resulting protein profiles were performed by BCCM/LMG (Laboratory for Microbiology, University of Ghent, Ghent, Belgium) as previously described (32). The normalized and digitized protein patterns were numerically analyzed and clustered with the reference profiles of 33 Bifidobacterium type strains.
16S rDNA sequencing. Two 16S rDNA primers, CO1, for the 5' end (5'AGTTTGATCCTGGCTCAG3'), and CO2, for the 3' end (5'TACCTTGTTACGACT3'), were used to generate an approximately 1.5-kb 16S rDNA product. The PCRs were performed in a 50-µl (each) solution containing 0.25 µM of each deoxynucleoside triphosphate (supplied as a 50x dNTP Mastermix; Bioline, USA Inc., Randolph, Mass.), 10 mM (NH4)2SO4, 67 mM Tris-HCl (pH 8.8), 0.1% (vol/vol) Tween-20 (supplied as a 10x NH4 reaction buffer; Bioline), 5 mM MgCl2 (supplied as a 50 mM stock; Bioline), 1 µM (each) primer, 2 µl of DNA, and 1 U of BIOTAQ DNA polymerase (Bioline). The mixture was denatured at 94°C for 5 min in a DNA thermal cycler (Hybaid Ltd., Middlesex, United Kingdom) followed by 30 amplification cycles, each consisting of three 1-min stages at 94°C (denaturation), 60°C (annealing), and 72°C (extension).
PCR amplicons were separated through a 1.5% (wt/vol) agarose gel using a 1x Tris base-acetate-EDTA running buffer containing 0.5 µg of ethidium bromide/ml and gel electrophoresis conditions of 100 V for 2 h. DNA was visualized by UV trans-illumination and sized using a 100-bp ladder (Amersham-Pharmacia-Biotech, Uppsala, Sweden) as a molecular weight standard. The
1.5-kb fragment was extracted from the gel using a QIAquick gel extraction kit (Qiagen, Valencia, Calif.). The recovered rDNA was sequenced using two forward primers: CO1 (5' AGT TTG ATC CTG GCT CAG 3') and MG3f (5' CTA CGG GAG GCA G 3') and three reverse primers, CO2 (5' TAC CTT GTT ACG ACT T 3'), 782R (5' ACC AGG GTA TCT AAT CCT GT 3'), and 765R (5' CTG TTT GCT CCC CAC GCT TTC 3') by MWG sequencing service (Ebersberg, Germany). Overlapping contigs and a consensus sequence were established using Seqmanager (DNAstar).
Phylogenetic analysis. The analysis of 16S rDNA sequences was as previously described (8). In brief, the 16S rDNA sequences obtained in the present study and those for each Bifidobacterium type strain (obtained from the National Center for Biotechnology Information database [www.ncbi.nlm.nih.gov]) were aligned using CLUSTAL X. Using the PHYLIP package, version 3.6 (obtained from www.evolutionary.genetics.washington.edu/phylip.html), the evolutionary distance matrices were generated with the Cantor and Juke coefficient by the DNADIST program. A phylogenetic tree was constructed according to the neighbor-joining method, using the program NEIGHBOR. The stability of the grouping was estimated by bootstrap analysis (1,000 replicons) using the programs SEQBOOT, DNADIST, NEIGHBOR, and CONSENSE. The phylogenetic tree was viewed using TREEVIEW software. The root of the unrooted tree based on the neighbor-joining method was estimated by using the outgroup Eubacterium combesii. Gardnerella vaginalis, Actinobaculum suis, Arcanobacterium haemolyticum, Mobiluncus curtisii, and Actinomyces bovis were also included as outgroups (16).
DNA-DNA hybridization reactions. For selected isolates, DNA-DNA hybridization reactions to Bifidobacterium type strains were performed by BCCM/LMG (Laboratory for Microbiology, University of Ghent, Ghent, Belgium) as previously described (12). Hybridization reactions with Bifidobacterium inopinatum and Bifidobacterium denticolens were performed at 44°C; all others were performed at 47°C.
Nucleotide sequence accession numbers. The GenBank accession numbers for PFGE types A through G are AY174103, AY174104, AY174105, AY174106, AY174107, AY174108, and AY174109, respectively.
|
|
|---|
0.1% of the total microflora enumerated. All of the 160 selected isolates tested positive for F-6-PPK activity, indicating their bifidobacterium status (41).
PFGE strain discrimination.
Previous studies have reported that PFGE following genomic digestion with the restriction enzyme XbaI can effectively discriminate strains from a number of Bifidobacterium species (29, 38). However, no strains of animal origin were included. In the present study, the discrimination of 38 animal-derived strains, including 11 type strains (Table 1), was examined based on genomic digests with XbaI and PFGE. The 38 strains produced 35 distinct macro-restriction patterns. The PFGE patterns for three Bifidobacterium suis strains (ATCC 27531, ATCC 27532, and ATCC 27533) were identical, as were those for two Bifidobacterium magnum strains (JCM 7120 and DSMZ 20220). In each case, following genomic digestion with the restriction enzyme SpeI, identical PFGE patterns within each group of strains were again observed, suggesting that they represented only two strains. In addition, Bifidobacterium thermophilum ATCC 25867 and Bifidobacterium saeculare DSMZ 6532 differed from their respective type strains by two extra fragments (at
40 and 80 kb) and one extra fragment (at
170 kb), respectively.
Given the effective level of strain discrimination achieved by XbaI and PFGE, 160 pig cecal isolates were examined. A total of 15 distinct macro-restriction patterns were observed (Fig. 1 and Table 2). Examples of different isolates having the same PFGE patterns are shown in Fig. 1 (see panel A, lanes 4 and 10, panel B, lanes 2 and 8, panel A, lane 5, and panel B, lane 7). In some instances, apparent differences were not recorded because they were considered to be an artifact of the weak staining associated with the low DNA content of fragments below 20 kb (see Fig. 1A, lane 8 and Fig. 1B, lane 11).
![]() View larger version (75K): [in a new window] |
FIG. 1. PFGE macro-restriction patterns following genomic DNA digests with the restriction enzyme XbaI. (A) Lanes 2 through 10 represent PFGE types G, F, E, D, C, B, and, PFGE subtype Ea, A, and PFGE type E, respectively. (B) Lanes 2 through 12 represent PFGE types and PFGE subtypes G, Caa, Ca, Aa, Ga, F, G, Cb, Ab, A, and Aa, respectively. Lanes 1, molecular size markers, lambda HindIII fragments with lambda concatemers.
|
|
View this table: [in a new window] |
TABLE 2. Frequency of PFGE types and subtypes representing 160 Bifidobacterium pig cecal isolates
|
280-kb fragment would then produce two smaller fragments at
170 and
110 kb, equal in size to the missing type A fragment. A single fragment of
85 kb in type A and an additional fragment of
60 kb, which comigrated with an existing fragment resulting in a fragment with a greater staining intensity, in subtype Ab, distinguished the two PFGE patterns (Fig. 1B, lanes 10 and 11, respectively). The difference could be explained by a
25-kb deletion in the
85-kb fragment. If the difference was explained by a
25-kb insertion in the subtype Ab fragment at
60 kb, a further mutation to remove an XbaI site would be required to produce subtype Aa.
The occurrence of the subtypes among the seven types was not equally distributed, with types A and C accounting for 39 and 44% of the subtypes, respectively. Both of the type A subtypes appeared to result from chromosomal mutations, unlike the majority of type C subtypes, which appeared to have putative plasmid differences (Table 2). Assuming each of the putative plasmids had a single XbaI site, the following sizes were estimated:
30 kb for E and Ea (Fig. 1A, lanes 10 and 9, respectively),
45 kb for G and Ga (Fig. 1B, lanes 8 and 6, respectively),
75 kb for C and Ca or Caa (Fig. 1A, lane 6, and 1B, lanes 4 and 3, respectively) and
300 kb for C and Cb (Fig. 1A, lane 6, and 1B, lane 9, respectively). The extreme size of the latter putative plasmid prompted analysis with a second restriction enzyme, SpeI, which also yielded a single fragment difference at
300 kb, suggesting that the fragment did in fact represent a single megaplasmid.
Overall, the relative distribution of the seven types was maintained within the three cecal samples analyzed (Table 2). This suggested that if bifidobacteria with specific abilities to adhere and colonize the epithelial layer do exist, these are also present at a similar relative frequency within the cecal contents, although in the latter sample, cell numbers were
200-fold higher. Type E was the most dominant isolate, accounting for 48, 51, and 35% of the bifidobacteria isolated from the two cecal epithelial samples and contents, respectively.
Relatedness of PFGE types. Using five different methods of analysis, the relatedness between the seven PFGE types that appeared to represent unrelated macro-restriction patterns were determined. Results for the following sections are summarized in Table 3.
|
View this table: [in a new window] |
TABLE 3. Relatedness of PFGE types A through G
|
![]() View larger version (55K): [in a new window] |
FIG. 2. Dendrograms based on computer comparisons of RAPD PCR patterns produced by random primers P1 and P2 with Bifidobacterium isolates (PFGE types A through G) recovered from a pig cecum and animal-derived culture collection strains (Table 1). The dendrograms were constructed using the Pearson product moment correlation coefficient (r) (similarity, 100 x r) and unweighted pair group algorithm method with arithmetic averages. Superscript T denotes type strain.
|
0.7 to 1 µm wide and
2 to 3 µm long, represented by type A and subtypes Aa and Ab only (Fig. 3A). Group 2 cells were also irregular rods, arranged singly or in pairs, with occasional bifurcated extremities,
1 to 1.2 µm wide and
2.5 to 5 µm long, with a tendency to form large aggregates (not shown), represented by types C, D, and G and all their related subtypes (Fig. 3C, D, and G). Group 3 cells were irregular short rods, arranged singly or in pairs,
0.6 to 0.9 µm wide and
0.8 to 1.5 µm long, represented by types B and E and subtype Ea (Fig. 3B and E). Group 4 cells were also short rods but were generally more globular than group 3 cells, arranged singly or in pairs,
0.7 to 1 µm wide and
1 to 1.5 µm long, represented by type F (Fig. 3F). Therefore, with the exception of type C, the morphology groups concurred with the genetic relatedness observed following RAPD PCR (Table 3).
![]() View larger version (68K): [in a new window] |
FIG. 3. Phase contrast cell morphology of Bifidobacterium isolates recovered from a pig cecum. Panels A through G represented PFGE types A through G, respectively. Four morphology groups, 1 through 4, were considered to represent type A, B and E, C, D, and G, and F, respectively (Table 3). Each subtype segregated with its related type. Scale bars in the main and inset panels are 10 and 2 µm, respectively.
|
(iv) 16S rDNA sequencing. The comparison of Bifidobacterium 16S rDNA sequences has been shown to group all species within the genus and identify isolates to specific species or groups of species (11, 16, 26). In general, when 16S rDNA similarity values exceed 97%, the strains are considered to belong to the same species (47). Therefore, a tree based on comparisons between almost complete 16S rDNA sequences (>1,400 bases) for types A through G and for each Bifidobacterium type strain was constructed (Fig. 4). The seven PFGE types segregated in good agreement with the groupings based on cell morphology, WCPP, and RAPD PCR (Table 3). PFGE Type G shared 100% similarity with the B. pseudolongum subsp. pseudolongum type strain, and PFGE types C and D shared 99.93 and 100% similarity with the B. pseudolongum subsp. globosum type strain, respectively. However, in general agreement with WCPP, the assignments appeared to contradict the RAPD PCR-based data, which grouped types D and G but showed type C to be distinct. However, the 16S rDNA sequences for the two subspecies are themselves 99.14% similar, indicating the difficulty in distinguishing the two subspecies based on 16S rDNA sequence comparisons (26). Previously strains considered as intermediates between the two subspecies were proposed (52), and type D could represent an intermediate strain. Type A was identified as B. minimum since it was 99.57% similar to the species type strain. The inclusion of type F changed the position of B. minimum within the phylogenetic tree, replacing the three insect-derived Bifidobacterium species as the closest neighbor (26). However, type F was only 95.08 and 94.85% similar to type A and B. minimum, respectively, and was considered to likely represent a new Bifidobacterium species. Types B and E were 99.72% similar to each other and were considered to be highly related and probably from the same species. However, they were only 91.21 and 91.36% similar, respectively, to their nearest neighbor, Bifidobacterium denticolens, and were considered to represent a new Bifidobacterium species.
![]() View larger version (40K): [in a new window] |
FIG. 4. Phylogenetic tree based on 16S rDNA sequences for each Bifidobacterium species type strain and seven PFGE types, termed A through G, representing bifidobacteria isolated from a pig cecum. Database accession numbers are given in parenthesis. The tree was rooted with E. combesii and constructed using the neighbor-joining method. Bootstrap values, expressed as percentages of 1,000 replicons, are given at each branch point where above 50%.
|
|
View this table: [in a new window] |
TABLE 4. DNA-DNA hybridization homologies between selected PFGE types and Bifidobacterium type strainsa
|
|
|
|---|
PFGE has previously been found to discriminate between Bifidobacterium strains (29, 38), a finding extended in the present study for a wide range of animal derived strains. In addition, PFGE was previously shown to be an effective technique for determining the genomic diversity of bifidobacteria isolated from human feces (17, 23). In these studies, the PFGE macro-restriction patterns for 10 Bifidobacterium isolates recovered at different times from a number of individuals were compared. Each PFGE type was considered to represent a distinct strain, and a complex population was arbitrarily assigned to subjects with four or more strains, while samples with fewer patterns were considered to be simple. Both kinds of population were observed with equal frequency. Following this convention our findings suggest that a simple bifidobacteria population, dominated by 15 Bifidobacterium strains, was maintained within a pig cecum.
A variety of studies have reported that Bifidobacterium strains vary in their ability to adhere to cultured epithelial cells (2, 6, 9, 10, 19, 30) and that the general microflora of the epithelium and contents of the porcine GIT are distinct (1, 44). However, in our study of bifidobacteria, although the total numbers of bifidobacteria associated with the cecal contents and epithelium differed considerably, distinct populations were not found. The data suggest that in vivo a relationship between "attached," "loosely unattached," and "unattached" cells exists. However, it is not clear if the bifidobacteria proliferate in the lumen and attach in a density-dependent manner or if proliferation occurs when they are attached to the epithelium with the lumen fraction represented an accumulation of detached bacteria.
Representatives for each PFGE type were examined by cell morphology, RAPD PCR, WCPP, 16S rDNA sequencing, and DNA-DNA hybridization to determine if a relatedness existed between the apparently unrelated macro-restriction patterns. Data from each technique were mostly complementary, and overall the seven types appeared to represent four distinct groups.
One group was identified as B. minimum, a species previously isolated only from sewage (5, 40). Another group appeared to include strains from B. pseudolongum subsp. pseudolongum and B. pseudolongum subsp. globosum (52). This study appears to be the first to report the recovery of either subspecies from a pig cecum, although they have previously been isolated from pig feces (25, 53), pork meat (14), and rabbit cecum (20, 46).
The two remaining groups were not assigned to any known Bifidobacterium species, although 16S rDNA sequencing did place them within the Bifidobacterium genus. In addition, both groups were positive for F-6-PPK activity. Based on the presented data and phenotypic analysis (unpublished data), PFGE types B, E, and subtype Ea strains were provisionally identified as a new species, termed Bifidobacterium aerophilum sp. nov. Similarly, the PFGE type F strain was provisionally identified as a new species termed Bifidobacterium psychroaerophilum sp. nov. The identification of additional species supports a recent study in which 16S rDNA sequences for 14 of 20 undefined isolates from the pig cecum were less than 95% similar to any database sequence (33).
A number of PFGE macro-restriction patterns differed by a single DNA fragment. It is possible that the extra band seen in these patterns represents an insertion, around or into a XbaI site, of a fragment with two restriction sites separated by a distance equal to the size of the insertion. However, it was considered more likely that the fragment represents the presence of a putative plasmid. In support of this, the estimated sizes for the majority of the suspected plasmids were within the range previously noted for the Bifidobacterium genus (41). However, two putative megaplasmids were also observed which have not previously been reported for the genus. Interestingly, most of the plasmids appeared to be associated with PFGE type C, which was identified as B. pseudolongum subsp. globosum, one of only four species in a survey of 25 Bifidobacterium species reported to contained plasmids (42). Work is ongoing to verify the PFGE interpretations and characterize the possible plasmids.
One interesting aspect of this study is the discovery of strains which from macro-restriction pattern analysis appear to show genomic rearrangements through reciprocal events involving either deletions or insertions. It is unclear if this genomic plasticity reflects a high rate of genomic change or a selection for very similar strains within the cecum. Interestingly, a recent study found that the commensal intestinal pathogen Bacteroides fragilis modified the structure of exopolysaccharides through chromosomal inversions within genetic promoter regions (18). The changes were considered to be part of a defense mechanism against the host immune response and allow for the long-term colonization of the GIT. It is tempting to suggest that the genomic plasticity observed in the present study may also have some bearing on long-term colonization of bifidobacteria.
In conclusion, the study presents a detailed examination of genomic diversity and relatedness of the bifidobacteria isolated from a pig cecum. The results suggest that the same relatively simple bifidobacteria population exists in the lumen and epithelium. Moreover, two potentially new species and evidence for genomic plasticity within strains were reported.
We are grateful to Tomotari Mitsuoka for providing B. thermophilum biotype information, Noelle Brennan for her assistance with 16S rDNA sequence comparisons and the determination of phylogenetic differences, and Shahenda Mohamed Mahmoud Al Abi for technical assistance.
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»