Journal Search Engine
Search Advanced Search Adode Reader(link)
Download PDF Export Citaion korean bibliography PMC previewer
ISSN : 1226-4768(Print)
ISSN : 2288-1247(Online)
Food Engineering Progress Vol.21 No.2 pp.174-179
DOI : https://doi.org/10.13050/foodengprog.2017.21.2.174

Molecular Typing of Leuconostoc citreum Strains Isolated from Korean Fermented Foods Using a Random Amplified Polymorphic DNA Marker

Jasmine Kaur, Sulhee Lee, Anshul Sharma, Young-Seo Park*
Department of Food Science and Biotechnology, Gachon University
Corresponding author: Young-Seo Park, Department of Food Science and Biotechnology, Gachon University, 1342 Seongnam-daero, Seongnam-si, Gyeonggi-do 13120, Korea +82-31-750-5378; +82-31-750-5389ypark@gachon.ac.kr
April 26, 2017 May 13, 2017 May 16, 2017

Abstract

For preliminary molecular typing, PCR-based fingerprinting using random amplified polymorphic DNA (RAPD) is the method of choice. In this study, 14 bacterial strains were isolated from different Korean food sources, identified using 16S rRNA gene sequencing, and characterized through RAPD-PCR. Two PCR primers (239 and KAY3) generated a total of 130 RAPD bands, 14 distinct PCR profiles, 10 polymorphic bands, one monomorphic band, and four unique bands. Dendrogram-based analysis with primer 239 showed that all 14 strains could be divided into seven clades out of which clade VII had the maximum of seven. In contrast, dendrogram analysis with the primer KAY3 divided the 14 L. citreum strains into four clades out of which clade IV consisted of a maximum of 10 strains out of 14. This research identified and characterized bacterial populations associated with different Korean foods. The proposed RAPD-PCR method, based on sequence amplification, could easily identify and discriminate the lactic acid bacteria species at the strain-specific level and could be used as a highly reliable genomic fingerprinting tool.


초록


    Korea Institute of Planning and Evalution for Technology of Food, Agriculture
    314073-03-2-HD040

    Introduction

    Currently, approximately 400 species of lactic acid bacteria (LAB) have been documented, and are typically classified into seven genera including Leuconostoc, Oenococcus, Pediococcus, Lactobacillus, Enterococcus, Streptococcus, and Lactococcus (Carr et al., 2002; Salminen et al., 2004; Zhang et al., 2011). Among these genera, Leuconostoc comprises 24 different species (Kot et al., 2014). This genus contains non-motile, gram-positive, and catalase-negative facultative anaerobes, with DNA exhibiting a relatively low GC content (37-45 mol%). Leuconostoc is associated with meat, poultry, fish, raw milk, dairy products, and various materials of plant-origin (Kim & Chun, 2005; Björkroth & Holzapfel, 2006). It has been determined that different species of Leuconostoc play an important role in fermenting foods such as cheese, kimchi, and sauerkraut. Among LAB, Leuconostoc citreum is one of the most established bacterial strains used for manufacturing kimchi (a Korean traditional dish). It is listed in the generally recognized as safe (GRAS) category; therefore, it plays an important role in various applications in pharmaceutical and biotechnology industries. One important and promising application is the use of these bacteria as a host for the delivery of therapeutic proteins into humans or animals (Son et al., 2016). The accurate identification of LAB is required for both basic research and food-based industrial applications (Tanigawa & Watanabe, 2011). For the rapid identification and analysis of microbial communities, classical biochemical methods can be used, which involve assessing morphological, physiological, and biochemical properties for preliminary identification of bacteria from different sources. However, these techniques are not adequately efficient and sometimes, results can be inconclusive (Tafvizi & Ebrahimi, 2015). Therefore, correct identification and classification of LAB is difficult without the support of genotypic techniques (Kandler & Weiss, 1986; Gevers et al., 2001). Currently, various molecular typing methods such as random amplified polymorphic DNA (RAPD), restriction fragment length polymorphisms (RFLP), pulse-field gel electrophoresis (PFGE), protein fingerprinting, and repetitive element palindromic PCR (rep-PCR) are available for the identification and characterization of LAB species (Villani et al., 1997; Cibik et al., 2000; Pérez et al., 2002; Sánchez et al., 2005; Björkroth & Holzapfel, 2006; Vihavainen & Björkroth, 2009; Nieto-Arribas et al., 2010; Alegría et al., 2013). Genetic fingerprinting techniques can be used to characterize single bacterial isolates or bacterial communities. It was previously published that LAB genotypes can be identified by partial 16S rRNA gene sequencing followed by characterization using species specific PCR (Holzapfel et al., 2001). Among PCR-based methods, RAPD has been used for the differentiation and classification of LAB (Drake et al., 1996; Cocconcelli et al., 1997; Quiberoni et al., 1998; Tailliez et al., 1998; Cibik et al., 2000; Cardamone et al., 2011). The objective of this study was to examine natural LAB populations associated with different food sources in Korea and to report the molecular identification and characterization of 14 bacterial strains.

    Materials and Methods

    Strains and medium

    Fourteen bacterial strains were isolated from seven different food sources, specifically, kimchi (6), salted small octopus (3), pear (1), young radish kimchi (1), salted oyster (1), salted fish (1), and sweet pumpkin (1) from Korea. Strains were grown overnight on MRS (Difco Laboratories, Detroit, MI, USA) plates at 30°C. For long-term preservation, isolated cell suspensions were stored in broth cultures supplemented with 1.5% (w/v) glycerol at -80°C.

    DNA isolation

    For DNA isolation, bacterial cells from 3 mL overnight cultures were pelleted by centrifugation at 16,000×g for 1 min and the pellet was washed with 1 mL of TE buffer (0.1 M Tris- HCl, 0.01 M EDTA, 1 M NaCl). Subsequently, genomic DNA was isolated using an AccuPrep Genomic DNA Extraction kit (Bioneer, Daejeon, Korea) following the manufacturer’s procedure. The quality of DNA was validated using a 0.8% (w/v) agarose gel and 1× TAE buffer.

    Identification of isolated strains

    All isolated strains were identified using a set of two universal primers, specifically, 27F (5'-AGA GTT TGA TCM TGG CTC AG-3') and 1492R (5'-TAC GGY TAC CTT GTT ACG ACT T-3') for amplifying the 16S rRNA gene. The PCR mixtures were subjected to the following thermal cycling conditions: 94°C for 30 s, annealing at 50°C for 1 min, extension at 72°C for 1 min, and a final extension at 72°C for 7 min, in a 20 μL reaction volume, using the Bioneer Premix kit. Amplified products were analyzed using 0.8% agarose gel electrophoresis and purified using a Promega Kit (Gel and PCR cleanup system, Promega, WI, USA) following the manufacturer’s instructions. The purified DNA was sequenced using an ABI 3730 xl DNA Analyzer (Applied Biosystems, Foster City, CA, USA) and service from Macrogen sequencing (Macrogen Inc., Seoul, Korea). The 16S rRNA gene sequence homology was compared using the NCBI database with BLASTN.

    RAPD-PCR

    RAPD primers for this study were commercially synthesized by Macrogen (Macrogen Inc.). Two universal primers, 239 (5'- CTG AAG CGG A-3') and KAY3 (5'-CTG GCG ACT G-3'), which are usually used for the taxonomic discrimination of lactic acid bacteria, were used at a concentration of 10 pmol/ μL for RAPD-PCR typing of the 14 isolated bacterial strains. For PCR amplification, 1 μL of purified genomic DNA, 2 μL of each primer (for two separate reactions), and 17 μL of distilled water were added to the premix kit (Bioneer - 10 mM Tris-HCl (pH 9.0), 40 mM KCl, 1.5 mM MgCl2, 250 μM dNTP mix, and 1 U of Taq polymerase) to a final volume of 20 μL. In an automatic thermocycler (My Cycler, Bio-Rad, CA, USA), the reaction mixtures were subjected to the following temperature cycles: 94°C for 4 min, 2 cycles of denaturation at 94°C for 2 min each, annealing at 39°C for 2 min, extension at 72°C for 2 min, 35 cycles of denaturation at 95°C for 15 s each, annealing at 39°C for 15 s, and extension at 72°C for 20 s, followed by a final elongation at 72°C for 7 min. For electrophoresis, 5 μL of each amplified PCR product was loaded onto a 1.2% (w/v) agarose gel (Lonza, Rockland, NY, USA) in 1× TAE buffer. For DNA staining, 0.5 μL of the ethidium bromide solution (10 mg/mL) was directly added to each 10 mL of agarose solution. Two molecular ladders of 100 bp (Bioneer) and 1 kb (Takara Bio Inc., Shiga, Japan) were used as size markers. After electrophoresis, the gel was photographed using a gel documentation system (BioRad Gel Doc XR+) to verify the PCR results. Finally, the band pattern was analyzed using the Quantity One software and phylogenic trees were constructed by the unweighted pair group method with arithmetic mean (UPGMA) method.

    Results and Discussion

    The objective of this study was to develop a simple, costeffective, easy to understand, RAPD-based molecular tool for typing bacterial strains isolated from different food sources from Korea. For the identification of isolated bacteria, 16S rRNA gene sequencing was performed. Currently, 16S rRNA gene sequencing is the most commonly used method to study bacterial phylogeny and taxonomy. This is possibly because the 16S rRNA gene is present in all bacteria (Woese, 1987). BLAST analysis confirmed that all 14 isolated strains were Leuconostoc citreum. The gene sequences have been submitted to GenBank with the accession numbers (KX286339- KX286352). It has been demonstrated that the use of traditional diagnostic testing (i.e. sugar fermentation patterns) alone, commonly leads to misidentifications among Leuconostoc and other closely related genera such as Lactobacillus, Weissella, and Pediococcus (Kulwichit et al., 2007).

    DNA-based molecular techniques have been confirmed as valuable tools to assess the identity and diversity of different strains (Aznar & Alarcon, 2002). The PCR-based RAPD typing technique uses 10-15 bp primers that bind to multiple genomic sites to detect sequence diversity. The main advantage of this method is that it is fast, easy to perform, and inexpensive, and it has high discriminatory power (Pozo- Bayon et al., 2009). In a previous study (Cibik et al., 2000), a combination of RAPD-PCR and 16S rDNA sequencing with specific primers was applied to analyze the molecular diversity of a large collection of Leuconostoc strains mainly isolated from traditional French cheeses.

    In the present work, a preliminary taxonomic study was performed using 14 isolates from different food sources from Korea. Genomic DNA of 14 L. citreum strains was successfully amplified, and the results obtained using two oligonucleotide primers (239 and KAY3) and RAPD analysis are illustrated in Fig. 1. In total, 130 bands were obtained and seven different RAPD profiles were produced using both primers. The RAPD profile with primer 239 consisted of eight major bands of 2000, 1600, 1200, 1000, 600, 480, 350, and 200 bp (Fig. 1(A)), whereas with primer KAY3, seven major bands of 4000, 2500, 2000, 1600, 1000, 600, and 500 bp were reported (Fig. 1(B)). Many of the strains shared common DNA band sizes with both primers (for example, strains 11037, 11082, 11382, 11391, 11413, 11414, and 11422 with primer 239 and strains 11037, 11076, 11082, 11093, 11382, 11414, and 11447 with primer KAY3) as expected, because all strains were initially identified as L. citreum. Furthermore, with primer 239, the 2000 bp band was only present in strain 11001, and the 500 bp and 600 bp bands were unique to strains 11093 and 11388, respectively. The 350 bp band was prominent in only one sample (11432); similarly, the 200 bp band was identified in two strains (11076 and 11447) of the 14 selected for this study. In contrast to primer 239, a unique band of 500 bp was observed with strain 11388 using primer KAY3 (Fig. 1(B)). The maximum band size observed here was 4000 bp in sample 11432, but this band was very faint. The RAPD typing method has been used for the identification and discrimination of several LAB strains of Pediococcus, Lactobacillus, Enterococcus, and Oenococcus among others (Corroler et al., 1998; Nigatu et al., 1998; Morea et al., 1999; Mora et al., 2000; Suzzi et al., 2000; De Angelis et al., 2001; Moschetti et al., 2001). There has been an increasing need to select microbial strains with commercial applications and to improve quality- and safetyrelated aspects of traditional fermented foods (Fortina et al., 2003).

    Strains with similar RAPD-PCR patterns should not be regarded as identical, even when they come from the same sources, because they might possess dissimilar industrial and/ or biological attributes (Cardamone et al., 2011). To understand the phylogenetic relatedness between isolates, a dendrogram was constructed using the Quantity One software. Dendrogrambased analysis with primer 239 showed that all 14 strains could be divided into seven clades out of which clade VII had the maximum of seven (11382, 11077, 11037, 11082, 11414, 11422, and 11413) strains of 14. Clade V had only two strains (11388 and 11076). Clades I, II, III, IV, and VI were composed of one strain each (11001, 11447, 11391, 11432, and 11093, respectively), as shown in Fig. 2(A). In addition, a single representative sample from each clade revealed a unique banding pattern (Fig. 1(A)). In contrast, dendrogram analysis with the primer KAY3 divided the 14 L. citreum strains into four clades out of which clade IV consisted of a maximum of 10 (11447, 11382, 11413, 11037, 11093, 11082, 11391, 11077, 11414, and 11076) strains of 14. Clade I was represented by two strains (11432 and 11001) (Fig. 2(B)). Clades II and III were composed of strains 11388 and 11422, respectively. Samples in lanes 2, 14, 9, and 13 were representative of clades I, II, and III, which were separated from the major group and showed different RAPD gel band patterns in comparison to the other samples, using primer KAY3 (Fig. 1(B)). It was observed that some isolates grouped together with one primer, but changed groups with another primer. It has been suggested that the problem of RAPD-PCR reproducibility can be overcome by using rigid laboratory protocols. Ramos et al. (2008) developed a procedure to optimize repeatability and avoid bias in sampling loci for genetic analyses based on RAPD data.

    In addition, a single typing method is not sufficient for the identification of bacteria from different food sources. It can be concluded that a multidimensional approach using biochemical characterization, 16S rRNA gene sequencing, and different molecular typing methods could lead to accurate species identification.

    The current study presents the phylogenetic relationships and genetic variations among L. citreum strains isolated from different food products in Korea. RAPD profiles are useful for preliminary screening and typing of bacterial isolates at a lower cost and higher speed. Isolation of novel strains for industrial applications is of utmost importance. For identification of such novel strains, other molecular typing methods are also required. The use of such novel strains in the food industry could offer great potential for augmenting existing processes or for developing new applications.

    Acknowledgements

    This work was supported by the Korea Institute of Planning and Evaluation for Technology in Food, Agriculture, Forestry and Fisheries (IPET) through the High Value-added Food Technology Development Program, funded by the Ministry of Agriculture, Food and Rural Affairs (MAFRA) (grant number: 314073-03-2-HD040).

    Figure

    FOODENGPROG-21-174_F1.gif
    RAPD profiling of 14 strains of Leuconostoc citreum using primers 239 and KAY3.

    (A) RAPD profiling using primer 239. Lane M: marker 1000 bp and 100 bp; lane 2: 11001; lane 3: 11037; lane 4: 11076; lane 5: 11077; lane 6: 11082; lane 7: 11093; lane 8: 11382; lane 9: 11388; lane 10: 11391; lane 11: 11413; lane 12: 11414; lane 13: 11422; lane 14: 11432; lane 15: 11447; lane M: marker 1000 bp and 100 bp. (B) RAPD profiling using primer KAY3. The order of samples is same as using primer 239.

    FOODENGPROG-21-174_F2.gif
    Dendrograms generated after cluster analysis of the Leuconostoc citreum strains using the UPGMA method and the Quantity One software.

    (A) Dendrogram using primer 239. (B) Dendrogram using primer KAY3.

    Table

    Reference

    1. Alegría , Delgado S , Flórez AB , Mayo B (2013) Identification, typing, and functional characterization of Leuconostoc spp. strains from traditional, starter-free cheeses , Dairy Sci. Technol, Vol.93 ; pp.657-673
    2. Aznar R , Alarcon B (2002) On the specificity of PCR detection of Listeria monocytogenes in food a comparison of published primers , Syst. Appl. Microbiol, Vol.25 ; pp.109-119
    3. Björkroth J , Holzapfel W Dworkin M , Falkow S , Rosenberg E , Schleifer KH , Stackebrandt E (2006) Genera Leuconostoc, Oenococcus and Weissella , The Prokaryotes: A Handbook on the Biology of Bacteria: Firmicutes, Cyanobacteria, Springer, ; pp.267-319
    4. Cardamone L , Quiberoni A , Mercanti DJ , Fornasari ME , Reinheimer JA , Guglielmotti DM (2011) Adventitious dairy Leuconostoc strains with interesting technological and biological properties useful for adjunct starters , Dairy Sci. Technol, Vol.91 ; pp.457-470
    5. Carr FJ , Chill D , Maida N (2002) The lactic acid bacteria A literature survey , Crit. Rev. Microbiol, Vol.28 ; pp.281-370
    6. Cibik R , Lepage E , Tailliez P (2000) Molecular diversity of Leuconostoc mesenteroides and Leuconostoc citreum isolated from traditional French cheeses as revealed by RAPD fingerprinting, 16S rDNA sequencing and 16S rDNA fragmentamplification , Syst. Appl. Microbiol, Vol.23 ; pp.267-278
    7. Cocconcelli PS , Parisi MG , Senini L , Bottazzi V (1997) Use of RAPD and 16S rDNA sequencing for the study of Lactobacillus population dynamics in natural whey culture , Lett. Appl. Microbiol, Vol.25 ; pp.8-12
    8. Corroler D , Mangin I , Desmasures N , Gueguen M (1998) An ecological study of lactococci isolated from raw milk in the Camembert cheese registered designation of origin area , ppl. Environ. Microbiol, Vol.64 ; pp.4729-4735
    9. De Angelis M , Corsetti A , Tosti N , Rossi J , Corbo MR , Gobbetti M (2001) Characterization of non-starter lactic acidbacteria from Italian ewe cheeses based on phenotypic, genotypic, and cell wall protein analyses , Appl. Environ. Microbiol, Vol.67 ; pp.2011-2020
    10. Drake M , Small CL , Spence KD , Swanson BG (1996) Rapid detection and identification of Lactobacillus spp. in dairy products by using the polymerase chain reaction , J. Food Prot, Vol.59 ; pp.1031-1036
    11. Fortina MG , Ricci G , Acquti A , Zeppa G , Gandini A , Manachini PL (2003) Genetic characterization of some lactic acid bacteria occurring in an artisanal protected denomination origin (PDO) Italian Cheese, the Toma Piemontese , Food Microbiol, Vol.20 ; pp.397-404
    12. Gevers D , Huys G , Swings J (2001) Applicability of rep-PCR fingerprinting for identification of Lactobacillus species , EMS Microbiol. Lett, Vol.205 ; pp.31-36
    13. Holzapfel WH , Haberer P , Geisen R , Bjorkroth J , Schillinger U (2001) Taxonomy and important features of probiotic microorganisms in food and nutrition , Am. J. Clin. Nutr, Vol.73 ; pp.365S-373S
    14. Kandler O , Weiss N Sneath PHA , Mair NS , Sharpe ME , Holt JG (1986) Regular, non-sporing gram-positive rods , Bergey's Manual of Systematic Bacteriology, Williams and Wilkins Co, ; pp.1209-1234
    15. Kim M , Chun J (2005) Bacterial community structure in kimchi, a Korean fermented vegetable food, as revealed by 16S rRNA gene analysis , Int. J. Food Microbiol, Vol.103 ; pp.91-96
    16. Kot W , Neve H , Heller KJ , Vogensen FK (2014) Bacteriophages of Leuconostoc, Oenococcus, and Weissella , Front Microbiol, Vol.5 ; pp.186
    17. Kulwichit W , Nilgate S , Chatsuwan T , Krajiw S , Unhasuta C , Chongthaleong A (2007) Accuracies of Leuconostoc phenotypic identification a comparison of API systems and conventional phenotypic assays , BMC Infect. Dis, Vol.7 ; pp.69
    18. Mora D , Fortina MG , Parini C , Daffonchio D , Manachini PL (2000) Genomic subpopulations within the species Pediococcus acidilaclici detected by multi locus typing analysis Relationships between pediocin AcH/PA_1 producing and nonproducing strains , Microbiology, Vol.146 ; pp.2027-2038
    19. Morea M , Baruzzi F , Cocconcelli PS (1999) Molecular and physiological characterization of dominant populations in traditional Mozzarella cheese processing , J. Appl. Microbiol, Vol.87 ; pp.574-582
    20. Moschetti G , Blaiotta G , Villani F , Coppola S , Parente E (2001) Comparison of statistical methods for identification of Streptococcus thermophilus, Enterococcus faecalis, and Entercoccus faecium from randomly amplified polymorphic DNA patterns , Appl. Environ. Microbiol, Vol.67 ; pp.2156-2166
    21. Nieto-Arribas P , Seseña S , Poveda JM , Palop L , Cabezas L (2010) Genotypic and technological characterization of Leuconostoc isolates to be used as adjunct starters in Manchego cheese manufacture , Food Microbiol, Vol.27 ; pp.85-93
    22. Nigatu A , Ahrné S , Molin G (1998) Randomly amplified polymorphic DNA for discrimination of Pediococcus pentosaceus and Ped. acidilactici and rapid grouping of Pediococcus isolates , Lett. Appl. Microbiol, Vol.26 ; pp.412-416
    23. Pérez G , Cardell E , Zárate V (2002) Random amplified polymorphic DNA analysis for differentiation of Leuconostoc mesenteroides subspecies isolated from Tenerife cheese , Lett. Appl. Microbiol, Vol.34 ; pp.82-85
    24. Pozo-Bayón M , Pardo I , Ferrer S , Moreno-Arribas MV (2009) Molecular approaches for the identification and characterization of oenological lactic acid bacteria , Afric. J. Biotechnol, Vol.8 ; pp.3995-4001
    25. Quiberoni USA , Tailliez A , Quénée P , Suarez V , Reinheimer J (1998) Genetic (RAPD-PCR) and technological diversities among wild Lactobacillus helveticus strains , J. Appl. Microbiol, Vol.85 ; pp.591-596
    26. Ramos Jr , Telles MPC , Diniz-Filho JAF , Soares TN , Melo DB , Oliveira G (2008) Optimizing reproducibility evaluation for random amplified polymorphic DNA markers , Genet. Mol. Res, Vol.7 ; pp.1384-1391
    27. Salminen MK , Rautelin H , Tynkkynen S , Poussa T , Saxelin M , Valtonen V , Ja¨rvinen A (2004) Lactobacillus bacteremia, clinical significance, and patient outcome, with special focus on probiotic L. Rhamnosus GG , Clin. Infect. Dis, Vol.38 ; pp.62-69
    28. Sánchez JI , Martínez B , Rodríguez A (2005) Rational selection of Leuconostoc strains for mixed starters based on the physiological biodiversity found in raw milk fermentations , Int. J. Food Microbiol, Vol.105 ; pp.377-387
    29. Son YJ , Ryu AJ , Ling L , Han NS , Jeong KJ (2016) Develop ment of a high-copy plasmid for enhanced production of recombinant proteins in Leuconostoc citreum , Microb. Cell Fact, Vol.15 ; pp.12
    30. Suzzi G , Caruso M , Gardini F , Lombardi A , Vannini L , Guerzoni ME , Andrighetto C , Lanorte MT (2000) A survey of the Enterococci isolated from an artisanal Italian goat's cheese (semicotto caprino) , J. Appl. Microbiol, Vol.89 ; pp.267-274
    31. Tafvizi1 F , Ebrahimi MT (2015) Application of repetitive extragenic palindromic elements based on PCR in detectionof genetic relationship of lactic acid bacteria species isolated from traditional fermented food products , J. Agr. Sci. Tech, Vol.17 ; pp.87-98
    32. Tailliez P , Tremblay J , Ehrlich SD , Chopin A (1998) Molecular diversity and culture development for improving the flavor of Proosdij-type cheese , Int. Dairy. J, Vol.13 ; pp.159-168
    33. Tanigawa K , Watanabe K (2011) Multilocus sequence typing reveals a novel sub speciation of Lactobacillus delbrueckii , Microbiology, Vol.157 ; pp.727-738
    34. Vihavainen EJ , Björkroth KJ (2009) Diversity of Leuconostoc gasicomitatum associated with meat spoilage , Int. J. Food Microbiol, Vol.136 ; pp.32-36
    35. Villani F , Moschetti G , Blaiotta G , Coppola S (1997) Characterization of strains of Leuconostoc mesenteroides by analysis of soluble whole-cell protein pattern, DNA fingerprinting and restriction of ribosomal DNA , J. Appl. Microbiol, Vol.82 ; pp.578-588
    36. Woese CR (1987) Bacterial evolution , Microbiol. Rev, Vol.51 ; pp.221-271
    37. Zhang ZG , Ye ZQ , Yu L , Shi P (2011) Phylogenomic reconstruction of lactic acid bacteria an update , BMC Evol. Biol, Vol.11 ; pp.1-12