- Methodology article
- Open Access
Identification of species belonging to the Bifidobacterium genus by PCR-RFLP analysis of a hsp60 gene fragment
BMC Microbiology volume 13, Article number: 149 (2013)
Bifidobacterium represents one of the largest genus within the Actinobacteria, and includes at present 32 species. These species share a high sequence homology of 16S rDNA and several molecular techniques already applied to discriminate among them give ambiguous results.
The slightly higher variability of the hsp60 gene sequences with respect to the 16S rRNA sequences offers better opportunities to design or develop molecular assays, allowing identification and differentiation of closely related species. hsp60 can be considered an excellent additional marker for inferring the taxonomy of the members of Bifidobacterium genus.
This work illustrates a simple and cheap molecular tool for the identification of Bifidobacterium species. The hsp60 universal primers were used in a simple PCR procedure for the direct amplification of 590 bp of the hsp60 sequence. The in silico restriction analysis of bifidobacterial hsp60 partial sequences allowed the identification of a single endonuclease (HaeIII) able to provide different PCR-restriction fragment length polymorphism (RFLP) patterns in the Bifidobacterium spp. type strains evaluated. The electrophoretic analyses allowed to confirm the different RFLP patterns.
The developed PCR-RFLP technique resulted in efficient discrimination of the tested species and subspecies and allowed the construction of a dichotomous key in order to differentiate the most widely distributed Bifidobacterium species as well as the subspecies belonging to B. pseudolongum and B. animalis.
Members of the genus Bifidobacterium are Gram-positive, obligate anaerobic, non-motile, non-spore forming bacteria, and are the most important constituents of human and animal intestinal microbiota[2, 3]. Recently, news species of bifidobacteria have been described[4–6] and now more than 30 species have been included in this genus.
Bifidobacterium spp. can be detected in various ecological environments, such as intestines of different vertebrates and invertebrates, dairy products, dental caries and sewage. Considering the increasing application of Bifidobacterium spp. as protective and probiotic cultures[7–9], and the fast enlargement of the genus, easy identification tools to discriminate new isolates are essential. Moreover, their correct taxonomic identification is of outmost importance for their use as probiotics. Conventional identification and classification of Bifidobacterium species have been based on phenotypic and biochemical features, such as cell morphology, carbohydrate fermentation profiles, and polyacrylamide gel electrophoresis analysis of soluble cellular proteins. In the last years several molecular techniques have been proposed in order to identify bifidobacteria. Most available bifidobacterial identification tools are based on 16S rRNA gene sequence analysis, such as ARDRA[11, 12], DGGE and PCR with the use of species-specific primers[14–16]. However, 16S rDNA of Bifidobacterium spp. has a high similarity, ranging from 87.7 to 99.5% and bifidobacterial closely related species (e.g. B. catenulatum and B. pseudocatenulatum) or subspecies (e.g. B. longum and B. animalis subspecies) even possess identical 16S rRNA gene sequences[17, 18]. For this reason different molecular approaches have been tested based on repetitive genome sequences amplification, such as ERIC-PCR[19, 20], BOX-PCR[21, 22] or RAPD fingerprinting analysis. These fingerprinting methods have the disadvantage of a low reproducibility, and they need strict standardization of PCR conditions. The use of different polymerases, DNA/primer ratios or different annealing temperatures may lead to a discrepancy in the results obtained in different laboratories.
In recent years alternative molecular markers have been proposed for bifidobacteria identification (e.g. hsp60, recA, tuf, atpD, dnaK) and Ventura et al. developed a multilocus approach, based on sequencing results, for the analysis of bifidobacteria evolution. The hsp60 gene, coding for a highly conserved 60 kDa heat-shock-protein (a chaperonin), has been evaluated for phylogenetic analysis in bifidobacteria by Jian et al.. The sequence comparison of this gene has been already used for species identification and phylogenetic analysis of other genera (e.g. Staphylococcus, Lactobacillus) and enteric pathogens[26–28]. A chaperonin database (cpnDB) is available on line, collecting bacterial and eukaryotic sequences (http://www.cpndb.ca/cpnDB/home.php).
The purpose of this study is the development of a rapid, reproducible and easy-to-handle molecular tool for the identification of Bifidobacterium species isolated from various environments. The protocol is based on the restriction endonuclease analysis of the PCR-amplified hsp60 partial gene sequence (hsp60 PCR-RFLP) with the use of a single restriction enzyme and has been tested on the 30 most widely distributed Bifidobacterium species and subspecies. A diagnostic dichotomous key to speed up profile interpretation has also been proposed.
Bacterial strains and culture conditions
The type strains used to develop the technique are listed in Table 1, whereas the strains used to validate the method are reported in Table 2. The strains, belonging to BUSCoB (Bologna University Scardovi Collection of Bifidobacteria) collection, were isolated from faeces of human and animals and from sewage. Bacteria were maintained as frozen stocks at −80°C in the presence of skim milk as cryoprotective agent. Working cultures were prepared in TPY medium, grown anaerobically at 37°C and harvested at logarithmic phase.
In silico analysis
An in silico analysis was performed for the evaluation of a suitable restriction enzyme. Available hsp60 sequences had been retrieved from cpnDB database and GeneBank, thanks to the work of Jian et al.. In silico digestion analysis was carried out on fragments amplified by universal primers H60F-H60R using two on-line free software: webcutter 2.0 (http://rna.lundberg.gu.se/cutter2) and http://insilico.ehu.es/restriction softwares. Blunt end, frequent cutter enzymes that recognize not degenerated sequences have been considered in order to find a suitable enzyme for all the species (e.g. RsaI, HaeIII, AluI, AccII). However in silico analysis had been performed also on sticky end enzymes (e.g. AatII, Sau3AI, PvuI).
DNA extraction from pure cultures
10 ml of culture were harvested and washed twice with TE buffer (10 mM Tris–HCl, 1 mM EDTA, pH 7.6), resuspended in 1 ml TE containing 15 mg lysozyme and incubated at 37°C overnight. Cells were lysed with 3 ml of lysis buffer (100 mM Tris–HCl, 400 mM NaCl, 2 mM EDTA, pH 8.2), 220 μl SDS (10% w/v) and 150 μl proteinase K (>600 mAU/ml, solution) and incubated for 2 hours in water bath at 60°C. One ml of saturated NaCl solution was added and the suspension was gently inverted twice. Pellets were harvested through centrifugation (5000 × g) at room temperature for 15 minutes. After the transfer of clean supernatants in new tubes, DNA was precipitated with 2.5 volumes of cold ethanol (95%) and resuspended in 300 μl of TE buffer.
Amplification of gene hsp60 and restriction with HaeIII
Universal primers were used to amplify approximately 600 bp of the hsp60 gene in the Bifidobacterium spp. investigated. These primers H60F (5‘-GG(ATGC)GA(CT)GG(ATGC)AC(ATGC)AC(ATGC)AC(ATGC)GC(ATGC)AC(ATGC)GT-3’) and H60R (5’-TC(ATGC)CC(AG)AA(ATGC)CC(ATGC)GG(ATGC)GC(CT)TT(ATGC)AC(ATGC)GC-3’) were designed by Rusanganwa et al. on the basis of the conserved protein sequences GDGTTATV and AVKAPGFGD in HSP60. Amplifications were performed in 20 μl volumes with 1.5 μM of each primer (Eurofins MWG Operon, Ebersberg, Germany), 10 μl 2X HotStarTaq Plus Master Mix (Qiagen, Italy) (1,5 mM MgCl2, 1 U Taq, 0.2 mM dNTP, final concentration) and 150 ng/μl DNA. The PCR cycle consisted of an initial denaturation of 5 min at 95°C followed by 35 cycles of denaturation (30s at 94°C), annealing (30s at 61°C) and extension (45 s at 72°C). The PCR was completed with a final elongation of 10 min at 72°C. The PCR amplification was performed with a PCR Verity 96-well thermal cycler (Applied Biosystems, Milan, Italy). After amplification, the product was visualized via agarose gel (1.3% w/v) in 1X TBE buffer and visualized with ethidium bromide under UV light. A 100 bp DNA ladder (Sigma-Aldrich) was used as a DNA molecular weight marker. Bands were excised from agarose gel (Additional file1: Figure S1) and DNA was eluted with NucleoSpin® Gel and PCR Clean-up (Macherey-Nagel GmbH & Co. KG, Germany) in order to avoid possible non-specific amplifications. 3 μl of the eluted DNA was re-amplified in a 30 μl PCR reaction (see above). BSA was added to the reaction (5% v/v, Fermentas). The PCR products (2 μl) were checked for non-specific amplification on agarose gel. 20 μl (~6 μg) of PCR amplicons were digested with HaeIII enzyme. Restriction digestion was carried out for 2 h at 37°C in 30 μl reaction mixture with 1X SM Restriction Buffer (Sigma-Aldrich), 1.5 μl HaeIII (10 U/μl, Sigma-Aldrich) and water. Digestion products were stained with ethidium bromide and visualized under UV-light (GelDoc™, BioRad), after agarose gel electrophoresis (3.0% agarose (w/v), TBE 1X) at 210 V (3 h). A 20 bp DNA ladder (Sigma-Aldrich) was used. The obtained pictures were elaborated with a free software GNU Image Manipulation Program (Gimp 2.8) only to invert colors and increase contrast.
Precast gradient polyacrylamide gels (4-20%) (Lonza Group Ltd, Switzerland) were also used to obtain RFLP profiles, in order to have a comparison with agarose gels. The vertical electrophoresis apparatus used was P8DS™ Emperor Penguin (Owl, Thermo Scientific) with an adaptor for Lonza precast gels. The run was performed at 100 V in TBE 1X.
A dichotomous key was developed comparing in silico digestion results and the evaluation of visible bands with the use of ImageLab™ 2.0 software (Bio-Rad Laboratories, Inc.).
Results and discussion
In silico analysis
The analysis and comparison of restriction profiles obtained with in silico digestion of bifidobacterial hsp60 sequences allowed the identification of a set of appropriate frequent-cutter endonucleases that recognize non degenerated sequences. The restriction enzyme HaeIII was found to give the clearest and most discriminatory profiles in theoretical PCR-RFLP patterns, discriminating the majority of Bifidobacterium type-strains tested (Table 3). Furthermore, the profiles of other strains, belonging to the investigated species, have been analyzed to confirm the conservation of RFLP profiles within species.
Amplification and restriction analysis of Bifidobacterium spp
Theoretical restriction profiles have been confirmed in vitro on agarose gel. The obtained fragments ranged from 16 bp to 339 bp (Table 3). Fragments lower than 25 bp were not considered as they did not help in species discrimination and in addition they co-migrate with primers. Time course analysis of restricted samples showed the formation of a band of ~200 bp in several species due to an over-digestion (data not shown) and this invalidated the RFLP profiles. For this reason the protocol has been optimized at 2 hours restriction time. Fragments greater than 360 bp were also not considered due to a possible incomplete digestion of such long fragments.
The obtained gels (Figures 1,2,3,4 and5) show species-specific profiles for all type-strains other than B. longum and B. thermacidophilum subspecies. This technique does not allow the identification of the subspecies belonging to these species, which displayed identical RFLP profiles. Matsuki et al.[14, 17] proposed specific primers to differentiate the subspecies of the species B. longum, while B. thermacidophilum subsp. porcinum and B. thermacidophilum subsp. thermacidophilum can be differentiated according to Zhu et al.. The proposed restriction analysis is efficient in discriminating very closely related species and subspecies as B. catenulatum/B. pseudocatenulatum, B. pseudolongum subsp. pseudolongum/B. pseudolongum subsp. globosum and B. animalis subsp. animalis/B. animalis. subsp. lactis.
The same method has been applied with the use of precast gradient polyacrylamide gels. The resolution was greater than that obtained on agarose gels, loading only 4 μl of the restriction reaction instead of the 30 μl used in horizontal electrophoresis. This may allow to reduce the volume of amplification reactions with a consequent reduction of costs.
The comparison between in silico digestion and the obtained gel profiles allowed to develop a dichotomous key (Figure 6) for a faster interpretation of the restriction profiles.
Validation of PCR-RFLP analysis on bifidobacterial isolates
39 strains belonging to 12 different species/subspecies (Table 2) have been investigated to validate the PCR-RFLP technique. Most of the strains tested were previously identified using biochemical tests and in some cases also molecular techniques (species-specific PCR, 16S rDNA sequencing). The obtained data confirmed a conservation of the profiles concerning the species and subspecies tested. Two figures are available as Additional files (Additional file 2: Figure S2: strains belonging to B. animalis subsp. lactis and B. animalis subsp. animalis. Additional file 3: Figure S3: strains belonging to B. longum subsp. longum, B. longum subsp. infantis, B. longum subsp. suis). About 95% of the strains confirmed the taxonomic identification previously assigned. Two strains, B1955 and Su864, previously classified as B. catenulatum and B. longum subsp. suis respectively, gave different profiles from those expected. The RFLP profiles of B1955 turned out to be the same of B. adolescentis ATCC 15703 (T), the dichotomous key confirmed the assignment to the B. adolescentis species. In addition, Su864 was identified as a B. breve strain. These results were also verified through a species-specific PCR .
In this work a PCR-RFLP based method to identify Bifidobacterium spp. was developed and tested on strains belonging to different species. The technique could efficiently differentiate all the 25 species of Bifidobacterium genus and the subspecies belonging to B. pseudolongum and B. animalis, with the support of an easy-to-handle dichotomous key. The technique turned out to be fast and easy, and presented a potential value for a rapid preliminary identification of bifidobacterial isolates.
Polymerase chain reaction
Restriction fragment length polymorphism
Heat-shock protein 60
Amplified ribosomal DNA restriction analysis
Denaturing gradient gel electrophoresis
Enterobacterial repetitive intergenic consensus-PCR
Random amplified polymorphic DNA
- TPY medium:
Tryptone phytone, yeast medium
(Bologna University Scardovi Collection of Bifidobacteria).
Biavati B, Mattarelli P: Genus Bifidobacterium. Bergey’s Manual of systematic bacteriology. Edited by: Goodfellow M, Kampfer P, Busse H-J, Suzuki K-I, Ludwig W, Whitman WB. 2012, New York: Springer, 171-206. Volume 5 2
Gaggìa F, Mattarelli P, Biavati B: Probiotics and prebiotics in animal feeding for safe food production. Int J Food Microbiol. 2010, 141: S15-S28.
Turroni F, Ribbera A, Foroni E, van Sinderen D, Ventura M: Human gut microbiota and bifidobacteria: from composition to functionality. Antonie Van Leeuwenhoek. 2008, 94: 35-50. 10.1007/s10482-008-9232-4.
Endo A, Futagawa-Endo Y, Schumann P, Pukall R, Dicks LM: Bifidobacterium reuteri sp. nov., Bifidobacterium callitrichos sp. nov., Bifidobacterium saguini sp. nov., Bifidobacterium stellenboschense sp. nov. and Bifidobacterium biavatii sp. nov. isolated from faeces of common marmoset (Callithrix jacchus) and red-handed tamarin (Saguinus midas). Syst Appl Microbiol. 2012, 35: 92-97. 10.1016/j.syapm.2011.11.006.
Kim MS, Roh SW, Bae JW: Bifidobacterium stercoris sp. nov., isolated from human faeces. Int J Syst Evol Microbiol. 2010, 60: 2823-2827. 10.1099/ijs.0.019943-0.
Morita H, Nakano A, Onoda H, Toh H, Oshima K, Takami H, Murakami M, Fukuda S, Takizawa T, Kuwahara T, Ohno H, Tanabe S, Hattori M: Bifidobacterium kashiwanohense sp. nov., isolated from healthy infant faeces. Int J Syst Evol Microbiol. 2011, 61: 2610-2615. 10.1099/ijs.0.024521-0.
Aloisio I, Santini C, Biavati B, Dinelli G, Cencič A, Chingwaru W, Mogna L, Di Gioia D: Characterization of Bifidobacterium spp. strains for the treatment of enteric disorders in newborns. App Microbiol Biotechnol. 2012, 96 (6): 1561-1576. 10.1007/s00253-012-4138-5.
Baffoni L, Gaggìa F, Di Gioia D, Santini C, Mogna L, Biavati B: A Bifidobacterium-based synbiotic product to reduce the transmission of C. jejuni along the poultry food chain. Int J Food Microbiol. 2012, 157 (2): 156-161. 10.1016/j.ijfoodmicro.2012.04.024.
Gaggìa F, Di Gioia D, Baffoni L, Biavati B: The role of protective and probiotic cultures in food and feed and their impact in food safety. Trends Foods Sci Tech. 2011, 22: 58-66.
Biavati B, Castagnoli P, Trovatelli LD: Species of the genus Bifidobacterium in the feces of human adults. Microbiologica. 1986, 9: 39-45.
Krízová J, Spanová A, Rittich B: Evaluation of amplified ribosomal DNA restriction analysis (ARDRA) and species-specific PCR for identification of Bifidobacterium species. Syst Appl Microbiol. 2006, 29: 36-44. 10.1016/j.syapm.2005.07.003.
Ventura M, Elli M, Reniero R, Zink R: Molecular microbial analysis of Bifidobacterium isolates from different environments by the species-specific amplified ribosomal DNA restriction analysis (ARDRA). FEMS Microbiol Ecol. 2001, 36: 113-121. 10.1111/j.1574-6941.2001.tb00831.x.
Temmerman R, Masco L, Vanhoutte T, Huys G, Swings J: Development and validation of a nested-PCR-denaturing gradient gel electrophoresis method for taxonomic characterization of bifidobacterial communities. Appl Environ Microbiol. 2003, 69: 6380-6385. 10.1128/AEM.69.11.6380-6385.2003.
Matsuki T, Watanabe K, Tanaka R, Oyaizu H: Rapid identification of human intestinal bifidobacteria by 16S rRNA-targeted species- and group-specific primers. FEMS Microbiol Lett. 1998, 167: 113-121. 10.1111/j.1574-6968.1998.tb13216.x.
Matsuki T, Watanabe K, Tanaka R, Fukuda M, Oyaizu H: Distribution of bifidobacterial species in human intestinal microflora examined with 16S rRNA-gene-targeted species-specific primers. Appl Environ Microbiol. 1999, 65: 4506-4512.
Matsuki T, Watanabe K, Fujimoto J, Kado Y, Takada T, Matsumoto K, Tanaka R: Quantitative PCR with 16S rRNA-gene-targeted species-specific primers for analysis of human intestinal bifidobacteria. Appl Environ Microbiol. 2004, 70: 167-173. 10.1128/AEM.70.1.167-173.2004.
Matsuki T, Watanabe K, Tanaka R: Genus- and species-specific PCR primers for the detection and identification of bifidobacteria. Curr Issues Intest Microbiol. 2003, 4: 61-69.
Ventura M, Canchaya C, Del Casale A, Dellaglio F, Neviani E, Fitzgerald GF, van Sinderen D: Analysis of bifidobacterial evolution using a multilocus approach. Int J Syst Evol Microbiol. 2006, 56: 2783-2792. 10.1099/ijs.0.64233-0.
Shuhaimi M, Ali AM, Saleh NM, Yazid AM: Utilisation of enterobacterial repetitive intergenic consensus (ERIC) sequence-based PCR to fingerprint the genomes of Bifidobacterium isolates and other probiotic bacteria. Biotech Lett. 2001, 23: 731-736. 10.1023/A:1010355003674.
Ventura M, Meylan V, Zink R: Identification and tracing of Bifidobacterium species by use of enterobacterial repetitive intergenic consensus sequences. Appl Environ Microbiol. 2003, 69: 4296-4301. 10.1128/AEM.69.7.4296-4301.2003.
Gómez Zavaglia A, de Urraza P, De Antoni G: Characterization of Bifdobacterium strains using Box primers. Anaerobe. 2000, 6: 169-177. 10.1006/anae.1999.0316.
Masco L, Huys G, Gevers D, Verbrugghen L, Swings J: Identification of Bifidobacterium species using rep-PCR fingerprinting. Syst Appl Microbiol. 2003, 26: 557-563. 10.1078/072320203770865864.
Vincent D, Roy D, Mondou F, Déry C: Characterization of bifidobacteria by random DNA amplification. Int J Food Microbiol. 1998, 43: 185-193. 10.1016/S0168-1605(98)00109-3.
Sidarenka AV, Novik GI, Akimov VN: Application of molecular methods to classification and identification of bacteria of the genus Bifidobacterium. Microbiology. 2008, 77: 251-260. 10.1134/S0026261708030016.
Jian W, Zhu L, Dong X: New approach to phylogenetic analysis of the genus Bifidobacterium based on partial HSP60 gene sequences. Int J Syst Evol Microbiol. 2001, 51: 1633-1638. 10.1099/00207713-51-5-1633.
Blaiotta G, Fusco V, Ercolini D, Aponte M, Pepe O, Villani F: Lactobacillus strain diversity based on partial hsp60 gene sequences and design of PCR-Restriction Fragment Length Polymorphism assays for species identification and differentiation. Appl Environ Microbiol. 2008, 74: 208-215. 10.1128/AEM.01711-07.
Goh SH, Potter S, Wood JO, Hemmingsen SM, Reynolds RP, Chow AW: HSP60 gene sequences as universal targets for microbial species identification: studies with coagulase-negative staphylococci. J Clin Microbiol. 1996, 34: 818-823.
Wong RS, Chow AW: Identification of enteric pathogens by heat shock protein 60kDa (HSP60) gene sequences. FEMS Microbiol Lett. 2002, 206: 107-113. 10.1111/j.1574-6968.2002.tb10994.x.
Hill JE, Penny SL, Crowell KG, Goh SH, Hemmingsen SM: cpnDB: a chaperonin sequence database. Genome Res. 2004, 14: 1669-1675. 10.1101/gr.2649204.
Rusanganwa E, Singh B, Gupta RS: Cloning of HSP60 (GroEL) operon from Clostridium perfringens using a polymerase chain reaction based approach. Biochim Biophys Acta. 1992, 1130: 90-94. 10.1016/0167-4781(92)90468-F.
Bikandi J, San Millán R, Rementeria A, Garaizar J: In silico analysis of complete bacterial genomes: PCR, AFLP-PCR, and endonuclease restriction. Bioinformatics. 2004, 20: 798-799. 10.1093/bioinformatics/btg491.
Rossi M, Altomare L, Rodriguez AG, Brigidi P, Matteuzzi D: Nucleotide sequence, expression and transcriptional analysis of the Bifidobacterium longum MB219 lacZ gene. Arch Microbiol. 2000, 174: 74-80. 10.1007/s002030000178.
Zhu L, Li W, Dong X: Species identification of genus Bifidobacterium based on partial HSP60 gene sequences and proposal of Bifidobacterium thermacidophilum subsp porcinum subsp nov. Int J Syst Evol Microbiol. 2003, 53: 1619-1623. 10.1099/ijs.0.02617-0.
The authors declare that they have no competing interests.
LB conceived the study. LB, VS and ES carried out all the bioinformatics, RFLP analyses, DNA extractions and culture handling. VS conceived the dichotomous key. MM and PM provided some of the strains tested together with the extracted DNA. DDG and FG supervised the work. LB, VS, DDG and FG contributed to paper writing. All authors read and approved the final manuscript. BB supported the research.
Electronic supplementary material
Additional file 1: Figure S1: Example of agarose gel electrophoresis of hsp60 amplicons from different bifidobacterial strains. (JPEG 121 KB)
Additional file 2: Figure S2: Agarose gel electrophoresis of digested hsp60 DNA fragments with HaeIII (negative image). Lane1, ladder 20 bp (Sigma-Aldrich); Lane 2–6, B. animalis subsp.lactis strains Ra20, Ra18, F439, P23, P32; Lane 7–8, B. animalis subsp. animalis strains T169, T6/1; Lane 9, ladder 20 bp (Sigma-Aldrich). (JPEG 467 KB)
Additional file 3: Figure S3: Agarose gel electrophoresis of digested hsp60 DNA fragments with HaeIII (negative image). Lane1, ladder 20 bp (Sigma-Aldrich); Lane 2–4, B. longum subsp. suis strains Su864, Su908, Su932; Lane 5–6, B. longum subsp. longum strains PCB133, ATCC 15707 (T); Lane 7–9, B. longum subsp. infantis strains ATCC 15697 (T), B7740, B7710; Lane 9, ladder 20 bp (Sigma-Aldrich). (JPEG 557 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
About this article
Cite this article
Baffoni, L., Stenico, V., Strahsburger, E. et al. Identification of species belonging to the Bifidobacterium genus by PCR-RFLP analysis of a hsp60 gene fragment. BMC Microbiol 13, 149 (2013). https://doi.org/10.1186/1471-2180-13-149
- Bifidobacterium spp