Lactobacillus casei, Lactobacillus rhamnosus, and Lactobacillus zeae isolates identified by sequence signature and immunoblot phenotype

The diversity and function of sourdough starter microbiomes

Humans have relied on sourdough starter microbial communities to make leavened bread for thousands of years, but only a small fraction of global sourdough biodiversity has been characterized. Working with a community-scientist network of bread bakers, we determined the microbial diversity of 500 sourdough starters from four continents. In sharp contrast with widespread assumptions, we found little evidence for biogeographic patterns in starter communities. Strong co-occurrence patterns observed in situ and recreated in vitro demonstrate that microbial interactions shape sourdough community structure. Variation in dough rise rates and aromas were largely explained by acetic acid bacteria, a mostly overlooked group of sourdough microbes. Our study reveals the extent of microbial diversity in an ancient fermented food across diverse cultural and geographic backgrounds.

groEL Gene-Based Phylogenetic Analysis of Lactobacillus Species by High-Throughput Sequencing

Lactobacillus is a fairly diverse genus of bacteria with more than 260 species and subspecies. Many profiling methods have been developed to carry out phylogenetic analysis of this complex and diverse genus, but limitations remain since there is still a lack of comprehensive and accurate analytical method to profile this genus at species level. To overcome these limitations, a Lactobacillus-specific primer set was developed targeting a hypervariable region in the groEL gene-a single-copy gene that has undergone rapid mutation and evolution. The results showed that this methodology could accurately perform taxonomic identification of Lactobacillus down to the species level. Its detection limit was as low as 10⁴ colony-forming units (cfu)/mL for Lactobacillus species. The assessment of detection specificity using the Lactobacillus groEL profiling method found that Lactobacillus, Pediococcus, Weissella, and Leuconostoc genus could be distinguished, but non-Lactobacillus Genus Complex could not be detected. The groEL gene sequencing and Miseq high-throughput approach were adopted to estimate the richness and diversity of Lactobacillus species in different ecosystems. The method was tested using kurut (fermented yak milk) samples and fecal samples of human, rat, and mouse. The results indicated that Lactobacillus mucosae was the predominant gut Lactobacillus species among Chinese, and L. johnsonii accounted for the majority of lactobacilli in rat and mouse gut. Meanwhile, L. delbrueckii subsp. bulgaricus had the highest relative abundance of Lactobacillus in kurut. Thus, this groEL gene profiling method is expected to promote the application of Lactobacillus for industrial production and therapeutic purpose.

Multilocus sequence typing of Lactobacillus casei isolates from naturally fermented foods in China and Mongolia

Lactobacillus casei is a lactic acid bacterium used in manufacturing of many fermented food products. To investigate the genetic diversity and population biology of this food-related bacterium, 224 Lb. casei isolates and 5 reference isolates were examined by multilocus sequence typing (MLST). Among them, 224 Lb. casei isolates were isolated from homemade fermented foods, including naturally fermented dairy products, acidic gruel, and Sichuan pickles from 38 different regions in China and Mongolia. The MLST scheme was developed based on the analysis of 10 selected housekeeping genes (carB, clpX, dnaA, groEL, murE, pyrG, pheS, recA, rpoC, and uvrC). All 229 isolates could be allocated to 171 unique sequence types, including 25 clonal complexes and 71 singletons. The high index of association value (1.3524) and standardized index of association value (0.1503) indicate the formation of an underlying clonal population by all the isolates. However, split-decomposition, relative frequency of occurrence of recombination and mutation, and relative effect of recombination and mutation in the diversification values confirm that recombination may have occurred, and were more frequent than mutation during the evolution of Lb. casei. Results from Structure analyses (version 2.3; http://pritch.bsd.uchicago.edu/structure.html) demonstrated that there were 5 lineages in the Lb. casei isolates, and the overall relatedness built by minimum spanning tree showed no clear relationship between the clonal complexes with either the isolation sources or sampling locations of the isolates. Our newly developed MLST scheme of Lb. casei was an easy and valuable tool that, together with the construction of an MLST database, will contribute to further detailed studies on the evolution and population genetics of Lb. casei from various niches.

Culture-Independent Metagenomic Surveillance of Commercially Available Probiotics with High-Throughput Next-Generation Sequencing

Millions of people consume dietary supplements either following a doctor's recommendation or at their own discretion to improve their overall health and well-being. This is a rapidly growing trend, with an associated and expanding manufacturing industry to meet the demand for new health-related products. In this study, we examined the contents and microbial viability of several popular probiotic products on the United States market. Culture-independent methods are proving ideal for fast and efficient analysis of foodborne pathogens and their associated microbial communities but may also be relevant for analyzing probiotics containing mixed microbial constituents. These products were subjected to next-generation whole-genome sequencing and analyzed by a custom in-house-developed k-mer counting method to validate manufacturer label information. In addition, the batch variability of respective products was examined to determine if any changes in their formulations and/or the manufacturing process occurred. Overall, the products we tested adhered to the ingredient claims and lot-to-lot differences were minimal. However, there were a few discrepancies in the naming of closely related Lactobacillus and Bifidobacterium species, whereas one product contained an apparent Enterococcus contaminant in two of its three lots. With the microbial contents of the products identified, we used traditional PCR and colony counting methods to comparatively assess our results and verify the viability of the microbes in these products with regard to the labeling claims. Of all the supplements examined, only one was found to be inaccurate in viability. Our use of next-generation sequencing as an analytical tool clearly demonstrated its utility for quickly analyzing commercially available products containing multiple microbes to ensure consumer safety. IMPORTANCE The rapidly growing supplement industry operates without a formal premarket approval process. Consumers rely on product labels to be accurate and true. Those products containing live microbials report both identity and viability on most product labels. This study used next-generation sequencing technology as an analytical tool in conjunction with classic culture methods to examine the validity of the labels on supplement products containing live microbials found in the United States marketplace. Our results show the importance of testing these products for identity, viability, and potential contaminants, as well as introduce a new culture-independent diagnostic approach for testing these products. Podcast: A podcast concerning this article is available.

Stress responses in probiotic Lactobacillus casei

Survival in harsh environments is critical to both the industrial performance of lactic acid bacteria (LAB) and their competitiveness in complex microbial ecologies. Among the LAB, members of the Lactobacillus casei group have industrial applications as acid-producing starter cultures for milk fermentations and as specialty cultures for the intensification and acceleration of flavor development in certain bacterial-ripened cheese varieties. They are amongst the most common organisms in the gastrointestinal (GI) tract of humans and other animals, and have the potential to function as probiotics. Whether used in industrial or probiotic applications, environmental stresses will affect the physiological status and properties of cells, including altering their functionality and biochemistry. Understanding the mechanisms of how LAB cope with different environments is of great biotechnological importance, from both a fundamental and applied perspective: hence, interaction between these strains and their environment has gained increased interest in recent years. This paper presents an overview of the important features of stress responses in Lb. casei, and related proteomic or gene expression patterns that may improve their use as starter cultures and probiotics.

High resolution melting analysis (HRM) as a new tool for the identification of species belonging to the Lactobacillus casei group and comparison with species-specific PCRs and multiplex PCR

The correct identification and characterisation of bacteria is essential for several reasons: the classification of lactic acid bacteria (LAB) has changed significantly over the years, and it is important to distinguish and define them correctly, according to the current nomenclature, avoiding problems in the interpretation of literature, as well as mislabelling when probiotic are used in food products. In this study, species-specific PCR and HRM (high-resolution melting) analysis were developed to identify strains belonging to the Lactobacillus casei group and to classify them into L. casei, Lactobacillus paracasei and Lactobacillus rhamnosus. HRM analysis confirmed to be a potent, simple, fast and economic tool for microbial identification. In particular, 201 strains, collected from International collections and attributed to the L. casei group, were examined using these techniques and the results were compared with consolidated molecular methods, already published. Seven of the tested strains don't belong to the L. casei group. Among the remaining 194 strains, 6 showed inconsistent results, leaving identification undetermined. All the applied techniques were congruent for the identification of the vast majority of the tested strains (188). Notably, for 46 of the strains, the identification differed from the previous attribution.

Identification of Lactobacillus strains of goose origin using MALDI-TOF mass spectrometry and 16S-23S rDNA intergenic spacer PCR analysis

The objective of our study was to identify Lactobacillus sp. strains of goose origin using MALDI-TOF mass spectrometry, ITS-PCR and ITS-PCR/RFLP. All three techniques proved to be valuable tools for identification of avian lactobacilli and produced comparable classification results. Lactobacillus strains were isolated from 100% of geese aged 3 weeks to 4 years, but from only 25% of chicks aged 1-10 days. Among the 104 strains isolated, we distinguished 14 Lactobacillus species. The dominant species was Lactobacillus salivarius (35.6%), followed by Lactobacillus johnsonii (18.3%), Lactobacillus ingluviei (11.5%) and Lactobacillus agilis (7.7%). The intact-cell MALDI-TOF mass spectrometry enabled rapid species identification of the lactobacilli with minimal pretreatment. However, it produced more than one identification result for 11.5% examined strains (mainly of the species L. johnsonii). ITS-PCR distinguished 12 genotypes among the isolates, but was not able to differentiate closely related strains, i.e. between Lactobacillus amylovorus and Lactobacillus kitasatonis and between Lactobacillus paracasei, Lactobacillus rhamnosus and Lactobacillus zeae. These species were differentiated by ITS-PCR/RFLP using the restriction enzymes TaqI and MseI. The results obtained indicate that ITS-PCR and ITS-PCR/RFLP assays could be used not only for interspecific, but also for intraspecific, typing.

16S rRNA PCR-Denaturing Gradient Gel Electrophoresis of Oral Lactobacillus casei Group and Their Phenotypic Appearances

This study aimed to develop a 16S rRNA PCR-denaturing gradient gel electrophoresis (DGGE) to identify the species level of Lactobacillus casei group and to investigate their characteristics of acid production and inhibitory effect. PCR-DGGE has been developed based on the 16S rRNA gene, and a set of HDA-1-GC and HDA-2, designed at V2-V3 region, and another set of CARP-1-GC and CARP-2, designed at V1 region, have been used. The bacterial strains included L. casei ATCC 393, L. paracasei CCUG 32212, L. rhamnosus ATCC 7469, L. zeae CCUG 35515, and 46 clinical strains of L. casei/paracasei/rhamnosus. Inhibitory effect against Streptococcus mutans and acid production were examined. Results revealed that each type species strain and identified clinical isolate showed its own unique DGGE pattern using CARP1-GC and CARP2 primers. HDA1-GC and HDA2 primers could distinguish the strains of L. paracasei from L. casei. It was found that inhibitory effect of L. paracasei was stronger than L. casei and L. rhamnosus. The acid production of L. paracasei was lower than L. casei and L. rhamnosus. In conclusion, the technique has been proven to be able to differentiate between closely related species in L. casei group and thus provide reliable information of their phenotypic appearances.

Comparative genome analysis of Lactobacillus casei strains isolated from Actimel and Yakult products reveals marked similarities and points to a common origin

The members of the Lactobacillus genus are widely used in the food and feed industry and show a remarkable ecological adaptability. Several Lactobacillus strains have been marketed as probiotics as they possess health-promoting properties for the host. In the present study, we used two complementary next-generation sequencing technologies to deduce the genome sequences of two Lactobacillus casei strains LcA and LcY, which were isolated from the products Actimel and Yakult, commercialized as probiotics. The LcA and LcY draft genomes have, respectively, an estimated size of 3067 and 3082 Mb and a G+C content of 46.3%. Both strains are close to identical to each other and differ by no more than minor chromosomal re-arrangements, substitutions, insertions and deletions, as evident from the verified presence of one insertion-deletion (InDel) and only 29 single-nucleotide polymorphisms (SNPs). In terms of coding capacity, LcA and LcY are predicted to encode a comparable exoproteome, indicating that LcA and LcY are likely to establish similar interactions with human intestinal cells. Moreover, both L. casei LcA and LcY harboured a 59.6 kb plasmid that shared high similarities with plasmids found in other L. casei strains, such as W56 and BD-II. Further analysis revealed that the L. casei plasmids constitute a good evolution marker within the L. casei species. The plasmids of the LcA and LcY strains are almost identical, as testified by the presence of only three verified SNPs, and share a 3.5 kb region encoding a remnant of a lactose PTS system that is absent from the plasmids of W56 and BD-II but conserved in another smaller L. casei plasmid (pLC2W). Our observations imply that the results obtained in animal and human experiments performed with the Actimel and Yakult strains can be compared with each other as these strains share a very recent common ancestor.

Funding Information The present work was supported by the Center of Excellence in Microbial Food Safety Research (Academy of Finland, Grant 141140), Grant ERC 250172 – Microbes Inside from the European Research Council and Grants 137389 and 141130 from the Academy of Finland. F.P.D. was funded by a postdoctoral research fellowship (Academy of Finland, Grant 252123).

A novel method for simultaneous Enterococcus species identification/typing and van genotyping by high resolution melt analysis

In order to develop a typing and identification method for van gene containing Enterococcus faecium, two multiplex PCR reactions were developed for use in HRM-PCR (High Resolution Melt-PCR): (i) vanA, vanB, vanC, vanC23 to detect van genes from different Enterococcus species; (ii) ISR (intergenic spacer region between the 16S and 23S rRNA genes) to detect all Enterococcus species and obtain species and isolate specific HRM curves. To test and validate the method three groups of isolates were tested: (i) 1672 Enterococcus species isolates from January 2009 to December 2009; (ii) 71 isolates previously identified and typed by PFGE (pulsed-field gel electrophoresis) and MLST (multi-locus sequence typing); and (iii) 18 of the isolates from (i) for which ISR sequencing was done. As well as successfully identifying 2 common genotypes by HRM from the Austin Hospital clinical isolates, this study analysed the sequences of all the vanB genes deposited in GenBank and developed a numerical classification scheme for the standardised naming of these vanB genotypes. The identification of Enterococcus faecalis from E. faecium was reliable and stable using ISR PCR. The typing of E. faecium by ISR PCR: (i) detected two variable peaks corresponding to different copy numbers of insertion sequences I and II corresponding to peak I and II respectively; (ii) produced 7 melt profiles for E. faecium with variable copy numbers of sequences I and II; (iii) demonstrated stability and instability of peak heights with equal frequency within the patient sample (36.4±4.5 days and 38.6±5.8 days respectively for 192 patients); (iv) detected ISR-HRM types with as much discrimination as PFGE and more than MLST; and (v) detected ISR-HRM types that differentiated some isolates that were identical by PFGE and MLST. In conjunction with the rapid and accurate van genotyping method described here, this ISR-HRM typing and identification method can be used as a stable identification and typing method with predictable instability based on recombination and concerted evolution of the rrn operon that will complement existing typing methods.

Discrimination of dairy industry isolates of the Lactobacillus casei group

Lactobacilli are a major part of the microflora of the gut and of many fermented dairy products, and are found in a variety of environments. Lactobacillus casei, Lactobacillus paracasei, Lactobacillus rhamnosus, and Lactobacillus zeae form a closely related taxonomic group within the facultatively heterofermentative lactobacilli. The classification and nomenclature of these bacteria are controversial. In this study, relationships between these species were investigated using type strains and dairy industry isolates examined with DNA-based techniques and conventional carbohydrate use tests. Carbohydrate use patterns gave poor discrimination of some species, but DNA PCR using specific primers targeted to sequences of the 16S rRNA gene discriminated 4 types consistent with the currently recognized species. Pulsed-field agarose gel electrophoresis of chromosomal NotI restriction fragments identified 18 different band patterns from 21 independent Lactobacillus isolates and confirmed the identity of L. casei strains from 2 culture collections (CSCC 5203 and ASCC 290), both representing the type strain of L. casei. Some isolates were reclassified as L. rhamnosus, suggesting that the prevalence of L. rhamnosus as a natural component of the microflora of dairy foods and dairy environments has previously been underestimated. These methods can provide a practical basis for discrimination of the species and identification of individual industrial strains.

Taxonomic and strain-specific identification of the probiotic strain Lactobacillus rhamnosus 35 within the Lactobacillus casei group

Lactobacilli are lactic acid bacteria that are widespread in the environment, including the human diet and gastrointestinal tract. Some Lactobacillus strains are regarded as probiotics because they exhibit beneficial health effects on their host. In this study, the long-used probiotic strain Lactobacillus rhamnosus 35 was characterized at a molecular level and compared with seven reference strains from the Lactobacillus casei group. Analysis of rrn operon sequences confirmed that L. rhamnosus 35 indeed belongs to the L. rhamnosus species, and both temporal temperature gradient gel electrophoresis and ribotyping showed that it is closer to the probiotic strain L. rhamnosus ATCC 53103 (also known as L. rhamnosus GG) than to the species type strain. In addition, L. casei ATCC 334 gathered in a coherent cluster with L. paracasei type strains, unlike L. casei ATCC 393, which was closer to L. zeae; this is evidence of the lack of relatedness between the two L. casei strains. Further characterization of the eight strains by pulsed-field gel electrophoresis repetitive DNA element-based PCR identified distinct patterns for each strain, whereas two isolates of L. rhamnosus 35 sampled 40 years apart could not be distinguished. By subtractive hybridization using the L. rhamnosus GG genome as a driver, we were able to isolate five L. rhamnosus 35-specific sequences, including two phage-related ones. The primer pairs designed to amplify these five regions allowed us to develop rapid and highly specific PCR-based identification methods for the probiotic strain L. rhamnosus 35.

Lactobacillus strain diversity based on partial hsp60 gene sequences and design of PCR-restriction fragment length polymorphism assays for species identification and differentiation

A phylogenetic tree showing diversities among 116 partial (499-bp) Lactobacillus hsp60 (groEL, encoding a 60-kDa heat shock protein) nucleotide sequences was obtained and compared to those previously described for 16S rRNA and tuf gene sequences. The topology of the tree produced in this study showed a Lactobacillus species distribution similar, but not identical, to those previously reported. However, according to the most recent systematic studies, a clear differentiation of 43 single-species clusters was detected/identified among the sequences analyzed. The slightly higher variability of the hsp60 nucleotide sequences than of the 16S rRNA sequences offers better opportunities to design or develop molecular assays allowing identification and differentiation of either distant or very closely related Lactobacillus species. Therefore, our results suggest that hsp60 can be considered an excellent molecular marker for inferring the taxonomy and phylogeny of members of the genus Lactobacillus and that the chosen primers can be used in a simple PCR procedure allowing the direct sequencing of the hsp60 fragments. Moreover, in this study we performed a computer-aided restriction endonuclease analysis of all 499-bp hsp60 partial sequences and we showed that the PCR-restriction fragment length polymorphism (RFLP) patterns obtainable by using both endonucleases AluI and TacI (in separate reactions) can allow identification and differentiation of all 43 Lactobacillus species considered, with the exception of the pair L. plantarum/L. pentosus. However, the latter species can be differentiated by further analysis with Sau3AI or MseI. The hsp60 PCR-RFLP approach was efficiently applied to identify and to differentiate a total of 110 wild Lactobacillus strains (including closely related species, such as L. casei and L. rhamnosus or L. plantarum and L. pentosus) isolated from cheese and dry-fermented sausages.

Multilocus sequence typing of Lactobacillus casei reveals a clonal population structure with low levels of homologous recombination

Robust genotyping methods for Lactobacillus casei are needed for strain tracking and collection management, as well as for population biology research. A collection of 52 strains initially labeled L. casei or Lactobacillus paracasei was first subjected to rplB gene sequencing together with reference strains of Lactobacillus zeae, Lactobacillus rhamnosus, and other species. Phylogenetic analysis showed that all 52 strains belonged to a single compact L. casei-L. paracasei sequence cluster, together with strain CIP107868 (= ATCC 334) but clearly distinct from L. rhamnosus and from a cluster with L. zeae and CIP103137(T) (= ATCC 393(T)). The strains were genotyped using amplified fragment length polymorphism, multilocus sequence typing based on internal portions of the seven housekeeping genes fusA, ileS, lepA, leuS, pyrG, recA, and recG, and tandem repeat variation (multilocus variable-number tandem repeats analysis [MLVA] using nine loci). Very high concordance was found between the three methods. Although amounts of nucleotide variation were low for the seven genes (pi ranging from 0.0038 to 0.0109), 3 to 12 alleles were distinguished, resulting in 31 sequence types. One sequence type (ST1) was frequent (17 strains), but most others were represented by a single strain. Attempts to subtype ST1 strains by MLVA, ribotyping, clustered regularly interspaced short palindromic repeat characterization, and single nucleotide repeat variation were unsuccessful. We found clear evidence for homologous recombination during the diversification of L. casei clones, including a putative intragenic import of DNA into one strain. Nucleotides were estimated to change four times more frequently by recombination than by mutation. However, statistical congruence between individual gene trees was retained, indicating that recombination is not frequent enough to disrupt the phylogenetic signal. The developed multilocus sequence typing scheme should be useful for future studies of L. casei strain diversity and evolution.