The partial 60 kDa heat-shock protein (HSP60) genes of 36 Bifidobacterium strains representing 30 different Bifidobacterium species and subspecies and of the type strain of Gardnerella vaginalis were cloned and sequenced using a pair of universal degenerate HSP60 PCR primers. The HSP60 DNA sequence similarities were determined for the taxa at various ranks as follows: 99.4-100% within the same species, 96% at the subspecies level, and 73-96% (mean 85%) at the interspecies level (and 98% in the case of two groups of closely related species, Bifidobacterium animalis and Bifidobacterium lactis, Bifidobacterium infantis, Bifidobacterium longum and Bifidobacterium suis, whose 165 rRNA sequence similarities are all above 99%). The HSP60 DNA sequence similarities between different Bifidobacterium species and G. vaginalis, a closely related bacterium according to 16S rRNA analysis, ranged from 71 to 79% (mean 75%). Although the topology of the phylogenetic tree constructed using the HSP60 sequences determined was basically similar to that for 16S rRNA, it seemed to be more clear-cut for species delineation, and the clustering was better correlated with the DNA base composition (mol% G+C) than that of the 16S rRNA tree. In the HSP60 phylogenetic tree, all of the high-G+C (55-67 mol%) bifidobacteria were grouped into one cluster, whereas the low-G+C species Bifidobacterium inopinatum (45 mol %) formed a separate cluster with G. vaginalis (42 mol%) and Bifidobacterium denticolens (55 mol%); a Bifidobacterium species of intermediate G+C content formed another cluster between the two. This study demonstrates that the highly conserved and ubiquitous HSP60 gene is an accurate and convenient tool for phylogenetic analysis of the genus Bifidobacterium.

The lack of a simple and rapid identification system for Bifidobacterium species makes them difficult to use in industrial applications. To obtain valuable discriminating factor, we studied different strains, and human isolates by two molecular taxonomy methods. First method was based on chrono-differentiation. A metabolic gene (pyruvate kinase) was chosen to be used as a systematic discriminating factor. A comparison of about 40 pyruvate kinase protein sequences allowed us to synthesize two oligonucleotides that were able to amplify a fragment of this corresponding gene in our strains. Based on these partial pyruvate kinase gene sequences, several clusters could be identified. The second method used in this study was based on 16S rRNA sequences analysis. We compared sequences present in GenBank database, and this allowed to separate bifidobacteria species into different clusters. They were different from those obtained with partial pyruvate kinase gene sequences analysis. So, by combining both methods, we were able to identify our isolates, when only 10% of them could be strictly identified using the 16S rRNA method. Moreover, pyruvate kinase analysis allowed to differentiate very ambivalent groups such as B. animalis/B. lactis or B. infantis/B. longum, but created different clusters for B. infantis species group, questioning on the homogeneity of this species.

The bacterial heat shock response is characterized by the elevated expression of a number of chaperone complexes and transcriptional regulators, including the DnaJ and the HrcA proteins. Genome analysis of Bifidobacterium breve UCC 2003 revealed a second copy of a dnaJ gene, named dnaJ2, which is flanked by the hrcA gene in a genetic constellation that appears to be unique to the actinobacteria. Phylogenetic analysis using 53 bacterial dnaJ sequences, including both dnaJ1 and dnaJ2 sequences, suggests that these genes have followed a different evolutionary development. Furthermore, the B. breve UCC 2003 dnaJ2 gene seems to be regulated in a manner that is different from that of the previously characterized dnaJ1 gene. The dnaJ2 gene, which was shown to be part of a 2.3-kb bicistronic operon with hrcA, was induced by osmotic shock but not significantly by heat stress. This induction pattern is unlike those of other characterized dnaJ genes and may be indicative of a unique stress adaptation strategy by this commensal microorganism.

The incorporation and delivery of bifidobacterial strains as probiotic components in many food preparations expose these microorganisms to a multitude of environmental insults, including heat and osmotic stresses. We characterized the dnaK gene region of Bifidobacterium breve UCC 2003. Sequence analysis of the dnaK locus revealed four genes with the organization dnaK-grpE-dnaJ-ORF1, whose deduced protein products display significant similarity to corresponding chaperones found in other bacteria. Northern hybridization and real-time LightCycler PCR analysis revealed that the transcription of the dnaK operon was strongly induced by osmotic shock but was not induced significantly by heat stress. A 4.4-kb polycistronic mRNA, which represented the transcript of the complete dnaK gene region, was detected. Many other small transcripts, which were assumed to have resulted from intensive processing or degradation of this polycistronic mRNA, were identified. The transcription start site of the dnaK operon was determined by primer extension. Phylogenetic analysis of the available bifidobacterial grpE and dnaK genes suggested that the evolutionary development of these genes has been similar. The phylogeny derived from the various bifidobacterial grpE and dnaK sequences is consistent with that derived from 16S rRNA. The use of these genes in bifidobacterial species as an alternative or complement to the 16S rRNA gene marker provides sequence signatures that allow a high level of discrimination between closely related species of this genus.

The bacterial heat shock response is characterized by the elevated expression of a number of chaperone complexes, including the GroEL and GroES proteins. The groES and groEL genes are highly conserved among eubacteria and are typically arranged as an operon. Genome analysis of Bifidobacterium breve UCC 2003 revealed that the groES and groEL genes are located in different chromosomal regions. The heat inducibility of the groEL and groES genes of B. breve UCC 2003 was verified by slot blot analysis. Northern blot analyses showed that the cspA gene is cotranscribed with the groEL gene, while the groES gene is transcribed as a monocistronic unit. The transcription initiation sites of these two mRNAs were determined by primer extension. Sequence and transcriptional analyses of the region flanking the groEL and groES genes of various bifidobacteria revealed similar groEL-cspA and groES gene units, suggesting a novel genetic organization of these chaperones. Phylogenetic analysis of the available bifidobacterial groES and groEL genes suggested that these genes evolved differently. Discrepancies in the phylogenetic positioning of groES-based trees make this gene an unreliable molecular marker. On the other hand, the bifidobacterial groEL gene sequences can be used as an alternative to current methods for tracing Bifidobacterium species, particularly because they allow a high level of discrimination between closely related species of this genus.

The atp operon is highly conserved among eubacteria, and it has been considered a molecular marker as an alternative to the 16S rRNA gene. PCR primers were designed from the consensus sequences of the atpD gene to amplify partial atpD sequences from 12 Bifidobacterium species and nine Lactobacillus species. All PCR products were sequenced and aligned with other atpD sequences retrieved from public databases. Genes encoding the subunits of the F(1)F(0)-ATPase of Bifidobacterium lactis DSM 10140 (atpBEFHAGDC) were cloned and sequenced. The deduced amino acid sequences of these subunits showed significant homology with the sequences of other organisms. We identified specific sequence signatures for the genus Bifidobacterium and for the closely related taxa Bifidobacterium lactis and Bifidobacterium animalis and Lactobacillus gasseri and Lactobacillus johnsonii, which could provide an alternative to current methods for identification of lactic acid bacterial species. Northern blot analysis showed that there was a transcript at approximately 7.3 kb, which corresponded to the size of the atp operon, and a transcript at 4.5 kb, which corresponded to the atpC, atpD, atpG, and atpA genes. The transcription initiation sites of these two mRNAs were mapped by primer extension, and the results revealed no consensus promoter sequences. Phylogenetic analysis of the atpD genes demonstrated that the Lactobacillus atpD gene clustered with the genera Listeria, Lactococcus, Streptococcus, and Enterococcus and that the higher G+C content and highly biased codon usage with respect to the genome average support the hypothesis that there was probably horizontal gene transfer. The acid inducibility of the atp operon of B. lactis DSM 10140 was verified by slot blot hybridization by using RNA isolated from acid-treated cultures of B. lactis DSM 10140. The rapid increase in the level of atp operon transcripts upon exposure to low pH suggested that the ATPase complex of B. lactis DSM 10140 was regulated at the level of transcription and not at the enzyme assembly step.

We analyzed the tuf gene, encoding elongation factor Tu, from 33 strains representing 17 Lactobacillus species and 8 Bifidobacterium species. The tuf sequences were aligned and used to infer phylogenesis among species of lactobacilli and bifidobacteria. We demonstrated that the synonymous substitution affecting this gene renders elongation factor Tu a reliable molecular clock for investigating evolutionary distances of lactobacilli and bifidobacteria. In fact, the phylogeny generated by these tuf sequences is consistent with that derived from 16S rRNA analysis. The investigation of a multiple alignment of tuf sequences revealed regions conserved among strains belonging to the same species but distinct from those of other species. PCR primers complementary to these regions allowed species-specific identification of closely related species, such as Lactobacillus casei group members. These tuf gene-based assays developed in this study provide an alternative to present methods for the identification for lactic acid bacterial species. Since a variable number of tuf genes have been described for bacteria, the presence of multiple genes was examined. Southern analysis revealed one tuf gene in the genomes of lactobacilli and bifidobacteria, but the tuf gene was arranged differently in the genomes of these two taxa. Our results revealed that the tuf gene in bifidobacteria is flanked by the same gene constellation as the str operon, as originally reported for Escherichia coli. In contrast, bioinformatic and transcriptional analyses of the DNA region flanking the tuf gene in four Lactobacillus species indicated the same four-gene unit and suggested a novel tuf operon specific for the genus Lactobacillus.

This study analyzes the application of degenerative primers for the screening of bile salt hydrolase-encoding genes (bsh) in various intestinal bifidobacteria. In the first stage, the design and evaluation of the universal PCR primers for amplifying the partial coding sequence of bile salt hydrolase in bifidobacteria were performed. The amplified bsh gene fragments were sequenced and the obtained sequences were compared to the bsh genes present in GenBank. The determined results showed the utility of the designed PCR primers for the amplification of partial gene encoding bile salt hydrolase in different intestinal bifidobacteria. Moreover, sequence analysis revealed that bile salt hydrolase-encoding genes may be used as valuable molecular markers for phylogenetic studies and identification of even closely related members of the genus Bifidobacterium.