Scoary2: rapid association of phenotypic multi-omics data with microbial pan-genomes

Unraveling bacterial gene function drives progress in various areas, such as food production, pharmacology, and ecology. While omics technologies capture high-dimensional phenotypic data, linking them to genomic data is challenging, leaving 40-60% of bacterial genes undescribed. To address this bottleneck, we introduce Scoary2, an ultra-fast microbial genome-wide association studies (mGWAS) software. With its data exploration app and improved performance, Scoary2 is the first tool to enable the study of large phenotypic datasets using mGWAS. As proof of concept, we explore the metabolome of yogurts, each produced with a different Propionibacterium reichii strain and discover two genes affecting carnitine metabolism.

Paenibacillus melissococcoides sp. nov., isolated from a honey bee colony affected by European foulbrood disease

A novel, facultatively anaerobic, Gram-stain-positive, motile, endospore-forming bacterium of the genus Paenibacillus , designated strain 2.1T, was isolated from a colony of Apis mellifera affected by European foulbrood disease in Switzerland. The rod-shaped cells of strain 2.1T were 2.2–6.5 µm long and 0.7–1.1 µm wide. Colonies of strain 2.1T were orange-pigmented under oxic growth conditions on solid basal medium at 35–37 °C. Strain 2.1T showed catalase and cytochrome c oxidase activity. Its polar lipid profile consisted of diphosphatidylglycerol, phosphatidylglycerol, phosphatidylethanolamine, aminophospholipid and phospholipid. The only respiratory quinone was menaquinone 7, and the major cellular fatty acids were anteiso-C15 : 0, anteiso-C17 : 0, iso-C15 : 0, iso-C17 : 0 and palmitic acid (C16 : 0), which is consistent with other members of the genus Paenibacillus . The G+C content of the genomic DNA of strain 2.1T was 53.3 mol%. Phylogenetic analyses based on the 16S rRNA gene sequence similarity showed that strain 2.1T was closely related to Paenibacillus dendritiformis LMG 21716T (99.7 % similarity) and Paenibacillus thiaminolyticus DSM 7262T (98.8 %). The whole-genome average nucleotide identity between strain 2.1T and the type strains of P. dendritiformis and P. thiaminolyticus was 92 and 91 %, respectively, and thus lower than the 95 % threshold value for delineation of genomic prokaryotic species. Based on the results of phylogenetic, genomic, phenotypic and chemotaxonomic analyses we propose the name Paenibacillus melissococcoides sp. nov. for this novel Paenibacillus species. The type strain is 2.1T (=CCOS 2000T=DSM 113619T=LMG 32539T).

Deciphering the global roles of Cold shock proteins in Listeria monocytogenes nutrient metabolism and stress tolerance

Listeria monocytogenes (Lm) accounts for serious public health and food safety problems owing to its stress resilience and pathogenicity. Based on their regulatory involvement in global gene expression events, cold-shock domain family proteins (Csps) are crucial in expression of various stress fitness and virulence phenotypes in bacteria. Lm possesses three Csps (CspA, CspB, and CspD) whose regulatory roles in the context of the genetic diversity of this bacterium are not yet fully understood. We examined the impacts of Csps deficiency on Lm nutrient metabolism and stress tolerance using a set of csp deletion mutants generated in different genetic backgrounds. Phenotype microarrays (PM) analysis showed that the absence of Csps in ∆cspABD reduces carbon (C-) source utilization capacity and increases Lm sensitivity to osmotic, pH, various chemical, and antimicrobial stress conditions. Single and double csp deletion mutants in different Lm genetic backgrounds were used to further dissect the roles of individual Csps in these phenotypes. Selected PM-based observations were further corroborated through targeted phenotypic assays, confirming that Csps are crucial in Lm for optimal utilization of various C-sources including rhamnose and glucose as well as tolerance against NaCl, β-phenyethylamine (PEA), and food relevant detergent stress conditions. Strain and genetic lineage background-based differences, division of labour, epistasis, and functional redundancies among the Csps were uncovered with respect to their roles in various processes including C-source utilization, cold, and PEA stress resistance. Finally, targeted transcriptome analysis was performed, revealing the activation of csp gene expression under defined stress conditions and the impact of Csps on expression regulation of selected rhamnose utilization genes. Overall, our study shows that Csps play important roles in nutrient utilization and stress responses in Lm strains, contributing to traits that are central to the public health and food safety impacts of this pathogen.

Extensive diversity and rapid turnover of phage defense repertoires in cheese-associated bacterial communities

BACKGROUND

Phages are key drivers of genomic diversity in bacterial populations as they impose strong selective pressure on the evolution of bacterial defense mechanisms across closely related strains. The pan-immunity model suggests that such diversity is maintained because the effective immune system of a bacterial species is the one distributed across all strains present in the community. However, only few studies have analyzed the distribution of bacterial defense systems at the community-level, mostly focusing on CRISPR and comparing samples from complex environments. Here, we studied 2778 bacterial genomes and 188 metagenomes from cheese-associated communities, which are dominated by a few bacterial taxa and occur in relatively stable environments.

RESULTS

We corroborate previous laboratory findings that in cheese-associated communities nearly identical strains contain diverse and highly variable arsenals of innate and adaptive (i.e., CRISPR-Cas) immunity systems suggesting rapid turnover. CRISPR spacer abundance correlated with the abundance of matching target sequences across the metagenomes providing evidence that the identified defense repertoires are functional and under selection. While these characteristics align with the pan-immunity model, the detected CRISPR spacers only covered a subset of the phages previously identified in cheese, providing evidence that CRISPR does not enable complete immunity against all phages, and that the innate immune mechanisms may have complementary roles.

CONCLUSIONS

Our findings show that the evolution of bacterial defense mechanisms is a highly dynamic process and highlight that experimentally tractable, low complexity communities such as those found in cheese, can help to understand ecological and molecular processes underlying phage-defense system relationships. These findings can have implications for the design of robust synthetic communities used in biotechnology and the food industry. Video Abstract.

Postprandial Responses on Serum Metabolome to Milk and Yogurt Intake in Young and Older Men

The identification and validation of biomarkers of food intake (BFIs) is a promising approach to develop more objective and complementary tools to the traditional dietary assessment methods. Concerning dairy, their evaluation in terms of intake is not simple, given the variety of existing foods, making it difficult to establish the association between specific dairy products consumption and the effects on human health, which is also dependent on the study population. Here, we aimed at identifying BFI of both milk (M) and yogurt (Y) in 14 healthy young (20–35 years) and 14 older (65–80 years). After a 3-week run-in period of dairy exclusion from the diet, the subjects acutely consumed 600 ml of M or Y. Metabolomics analyses were conducted on serum samples during the following 6 h (LC-MS and GC-MS). Several metabolites showing increased iAUC after milk or yogurt intake were considered as potential BFI, including lactose (M > Y, 2-fold), galactitol (M > Y, 1.5-fold), galactonate (M > Y, 1.2-fold), sphingosine-1-phosphate (M > Y from 2.1-fold), as well as an annotated disaccharide (Y > M, 3.6-fold). Delayed serum kinetics were also observed after Y compared to M intake lysine (+22 min), phenylalanine (+45 min), tyrosine (+30min), threonine (+38 min) 3-phenyllactic acid (+30 min), lactose (+30 min), galactitol (+45min) and galactonate (+30 min). The statistical significance of certain discriminant metabolites, such as sphingosine-1-phosphate and several free fatty acids, was not maintained in the older group. This could be related to the physiological modifications induced by aging, like dysregulated lipid metabolism, including delayed appearance of dodecanoic acid (+60 min) or altered postprandial appearance of myristic acid (+70% Cmax), 3-dehydroxycarnitine (−26% Cmin), decanoylcarnitine (−51% Cmin) and dodecanoylcarnitine (−40% Cmin). In conclusion, candidate BFI of milk or yogurt could be identified based on the modified postprandial response resulting from the fermentation of milk to yogurt. Moreover, population specificities (e.g., aging) should also be considered in future studies to obtain more accurate and specific BFI.

Functional strain redundancy and persistent phage infection in Swiss hard cheese starter cultures

Undefined starter cultures are poorly characterized bacterial communities from environmental origin used in cheese making. They are phenotypically stable and have evolved through domestication by repeated propagation in closed and highly controlled environments over centuries. This makes them interesting for understanding eco-evolutionary dynamics governing microbial communities. While cheese starter cultures are known to be dominated by a few bacterial species, little is known about the composition, functional relevance, and temporal dynamics of strain-level diversity. Here, we applied shotgun metagenomics to an important Swiss cheese starter culture and analyzed historical and experimental samples reflecting 82 years of starter culture propagation. We found that the bacterial community is highly stable and dominated by only a few coexisting strains of Streptococcus thermophilus and Lactobacillus delbrueckii subsp. lactis. Genome sequencing, metabolomics analysis, and co-culturing experiments of 43 isolates show that these strains are functionally redundant, but differ tremendously in their phage resistance potential. Moreover, we identified two highly abundant Streptococcus phages that seem to stably coexist in the community without any negative impact on bacterial growth or strain persistence, and despite the presence of a large and diverse repertoire of matching CRISPR spacers. Our findings show that functionally equivalent strains can coexist in domesticated microbial communities and highlight an important role of bacteria-phage interactions that are different from kill-the-winner dynamics.

Functional strain redundancy and persistent phage infection in Swiss hard cheese starter cultures

Undefined starter cultures are poorly characterized bacterial communities from environmental origin used in cheese making. They are phenotypically stable and have evolved through domestication by repeated propagation in closed and highly controlled environments over centuries. This makes them interesting for understanding eco-evolutionary dynamics governing microbial communities. While cheese starter cultures are known to be dominated by a few bacterial species, little is known about the composition, functional relevance, and temporal dynamics of strain-level diversity. Here, we applied shotgun metagenomics to an important Swiss cheese starter culture and analyzed historical and experimental samples reflecting 82 years of starter culture propagation. We found that the bacterial community is highly stable and dominated by only a few coexisting strains of Streptococcus thermophilus and Lactobacillus delbrueckii subsp. lactis. Genome sequencing, metabolomics analysis, and co-culturing experiments of 43 isolates show that these strains are functionally redundant, but differ tremendously in their phage resistance potential. Moreover, we identified two highly abundant Streptococcus phages that seem to stably coexist in the community without any negative impact on bacterial growth or strain persistence, and despite the presence of a large and diverse repertoire of matching CRISPR spacers. Our findings show that functionally equivalent strains can coexist in domesticated microbial communities and highlight an important role of bacteria-phage interactions that are different from kill-the-winner dynamics.

Formation of 3-Methylbutanal and 3-Methylbutan-1-ol Recognized as Malty during Fermentation in Swiss Raclette-Type Cheese, Reconstituted Milk, and de Man, Rogosa, and Sharpe Broth

This work aimed to determine the formation over time of 3-methylbutanal and 3-methylbutan-1-ol recognized as malty during the manufacture of Raclette-type cheese and the fermention of reconstituted skim milk, and filter-sterilized MRS broth. Using dynamic headspace-vacuum transfer in trap extraction followed by gas chromatography coupled with mass spectrometry-olfactometry (DHS-VTT-GC-MS-O) as a screening method for the malty compounds, five compounds (2-methylpropanal, 2- and 3-methylbutanal, and 2- and 3-methylbutan-1-ol) were identified as potential compounds causing the malty aroma in starter culture development and Raclette-type cheeses. Focus on compounds having a predominant sensorial effect (3-methylbutanal and 3-methylbutan-1-ol), spikings of leucine, ¹³C-labeled leucine, α-ketoisocaproic acid, and α-ketoglutaric acid provided a better understanding of their formation pathway. This study highlighted the discrepancies in the formation of 3-methylbutanal and 3-methylbutan-1-ol between the growth media; namely, despite the presence of free leucine available in MRS and the addition of an excess, no increase of the target compounds was observed. The concentration of these compounds in MRS increased only when α-ketoglutaric acid or α-ketoisocaproic acid was added, and a preference for the pathway to α-hydroxyisocaproic acid instead of 3-methylbutanal was shown. In addition, a formation of 3-methylbutanal when the bacteria were not yet active was observed when spiking α-ketoisocaproic acid, which potentially indicates that this part of the metabolism could take place extracellularly. These results could potentially unveil other, not-yet-identified reactants, directly influencing the production of compounds responsible for the malty aroma in Raclette cheese.

Evolution of Listeria monocytogenes During a Persistent Human Prosthetic Hip Joint Infection

Listeria monocytogenes associated prosthetic joint infections (PJI) are a rare but increasing clinical problem of listeriosis. We characterized two isolates of the same L. monocytogenes strain isolated within five years of each other from a recurrent human prosthetic joint infection. The two isolates although clonally identical were phenotypically distinct confirming that the original infection strain had evolved within the human host PJI environment giving rise to a phenotypically distinct variant. The recurrent PJI isolate displayed various phenotypic differences compared to the parental original PJI isolate including diminished growth and carbon source metabolism, as well as altered morphology and increased stress sensitivity. The PJI isolates were both diminished in virulence due to an identical truncation mutation in the major virulence regulator PrfA. Genome wide sequence comparison provided conclusive evidence that the two isolates were identical clonal descendants of the same L. monocytogenes strain that had evolved through acquisition of various single nucleotide polymorphisms (SNPs) as well as insertion and deletion events (InDels) during a persistent human PJI. Acquired genetic changes included a specific mutation causing premature stop codon (PMSC) and truncation of RNAse J1 protein. Based on analysis of this naturally truncated as well as other complete RNAse J1 deletion mutants we show that the long-term survival of this specific L. monocytogenes strain within the prosthetic joint might in part be explained by the rnjA PMSC mutation that diminishes virulence and activation of the host immune system in a zebrafish embryo localized infection model. Overall our analysis of this special natural case provides insights into random mutation events and molecular mechanisms that might be associated with the adaptation and short-term evolution of this specific L. monocytogenes strain within a persistent human PJI environment.

Variable Carbon Source Utilization, Stress Resistance, and Virulence Profiles Among Listeria monocytogenes Strains Responsible for Listeriosis Outbreaks in Switzerland

A combination of phenotype microarrays, targeted stress resistance and virulence assays and comparative genome analysis was used to compare a set of Listeria monocytogenes strains including those involved in previous Swiss foodborne listeriosis outbreaks. Despite being highly syntenic in gene content these strains showed significant phenotypic variation in utilization of different carbon (C)-sources as well as in resistance of osmotic and pH stress conditions that are relevant to host and food associated environments. An outbreak strain from the 2005 Swiss Tomme cheese listeriosis outbreak (Lm3163) showed the highest versatility in C-sources utilized whereas the strain responsible for the 1983 to 1987 Vacherin Montd'or cheese listeriosis outbreak (LL195) showed the highest tolerance to both osmotic and pH stress conditions among the examined strains. Inclusion of L-norvaline led to enhanced resistance of acidic stress in all the examined strains and there were strain-strain-specific differences observed in the ability of other amino acids and urea to enhance acid stress resistance in L. monocytogenes. A strain dependent inhibition pattern was also observed upon inclusion of β-phenylethylamine under alkaline stress conditions. In targeted phenotypic analysis the strain-specific differences in salt stress tolerance uncovered in phenotypic microarrays were corroborated and variations in host cell invasion and virulence among the examined strains were also revealed. Outbreak associated strains representing lineage I serotype 4b showed superior pathogenicity in a zebrafish infection model whilst Lm3163 a lineage II serotype 1/2a outbreak strain demonstrated the highest cellular invasion capacity amongst the tested strains. A genome wide sequence comparison of the strains only revealed few genetic differences between the strains suggesting that variations in gene regulation and expression are largely responsible for the phenotypic differences revealed among the examined strains. Our results have generated data that provides a potential basis for the future design of improved Listeria specific media to enhance routine detection and isolation of this pathogen as well as provide knowledge for developing novel methods for its control in food.

Metabolic Footprinting of Fermented Milk Consumption in Serum of Healthy Men

Background

Fermentation is a widely used method of natural food preservation that has consequences on the nutritional value of the transformed food. Fermented dairy products are increasingly investigated in view of their ability to exert health benefits beyond their nutritional qualities.

Objective

To explore the mechanisms underpinning the health benefits of fermented dairy intake, the present study followed the effects of milk fermentation, from changes in the product metabolome to consequences on the human serum metabolome after its ingestion.

Methods

A randomized crossover study design was conducted in 14 healthy men [mean age: 24.6 y; mean body mass index (in kg/m2): 21.8]. At the beginning of each test phase, serum samples were taken 6 h postprandially after the ingestion of 800 g of a nonfermented milk or a probiotic yogurt. During the 2-wk test phases, subjects consumed 400 g of the assigned test product daily (200 g, 2 times/d). Serum samples were taken from fasting participants at the end of each test phase. The serum metabolome was assessed through the use of LC-MS-based untargeted metabolomics.

Results

Postprandial serum metabolomes after milk or yogurt intake could be differentiated [orthogonal projections to latent structures discriminant analysis (OPLS-DA) Q2 = 0.74]. Yogurt intake was characterized by higher concentrations of 7 free amino acids (including proline, P = 0.03), reduced concentrations of 5 bile acids (including glycocholic acid, P = 0.04), and modulation of 4 indole derivative compounds (including indole lactic acid, P = 0.01). Fasting serum samples after 2 wk of daily intake of milk or yogurt could also be differentiated based on their metabolic profiles (OPLS-DA Q2 = 0.56) and were discussed in light of the postprandial results.

Conclusion

Metabolic pathways related to amino acids, indole derivatives, and bile acids were modulated in healthy men by the intake of yogurt. Further investigation to explore novel health effects of fermented dairy products is warranted.This trial was registered at clinicaltrials.gov as NCT02230345.

Isothermal micro calorimetry--a new method for MIC determinations: results for 12 antibiotics and reference strains of E. coli and S. aureus

Background

Antimicrobial susceptibility testing of microorganisms is performed by either disc diffusion or broth dilution tests. In clinical use, the tests are often still performed manually although automated systems exist. Most systems, however, are based on turbidometric methods which have well-known drawbacks.

Results

In this study we evaluated isothermal micro calorimetry (IMC) for the determination of minimal inhibitory concentrations (MICs) of 12 antibiotics for five micro-organisms. Here we present the data for the 12 antibiotics and two representative microorganisms E. coli (a Gram^-) and S. aureus (a Gram⁺). IMC was able to determine the MICs correctly according to CLSI values. Since MICs require 24 hours, time was not reduced. However, IMC provided new additional data – a continuous record of heat-producing bacterial activity (e.g. growth) in calorimetry ampoules at subinhibitory antibiotic concentrations. Key features of the heatflow (P) and aggregate heat (Q) vs. time curves were identified (t_delay and ΔQ/Δt respectively). Antibiotics with similar modes of action proved to have similar effects on t_delay and/or ΔQ/Δt.

Conclusion

IMC can be a powerful tool for determining the effects of antibiotics on microorganisms in vitro. It easily provides accurate MICs – plus a potential means for analyzing and comparing the modes of action of antibiotics at subinhibitory concentrations. Also IMC is completely passive, so after evaluation, ampoule contents (media, bacteria, etc.) can be analyzed by any other method desired.

Evaluation of the passage of Lactobacillus gasseri K7 and bifidobacteria from the stomach to intestines using a single reactor model

Background

Probiotic bacteria are thought to play an important role in the digestive system and therefore have to survive the passage from stomach to intestines. Recently, a novel approach to simulate the passage from stomach to intestines in a single bioreactor was developed. The advantage of this automated one reactor system was the ability to test the influence of acid, bile salts and pancreatin.

Lactobacillus gasseri K7 is a strain isolated from infant faeces with properties making the strain interesting for cheese production. In this study, a single reactor system was used to evaluate the survival of L. gasseri K7 and selected bifidobacteria from our collection through the stomach-intestine passage.

Results

Initial screening for acid resistance in acidified culture media showed a low tolerance of Bifidobacterium dentium for this condition indicating low survival in the passage. Similar results were achieved with B. longum subsp. infantis whereas B. animalis subsp. lactis had a high survival.

These initial results were confirmed in the bioreactor model of the stomach-intestine passage. B. animalis subsp. lactis had the highest survival rate (10%) attaining approximately 5 × 10⁶ cfu ml^-1 compared to the other tested bifidobacteria strains which were reduced by a factor of up to 10⁶. Lactobacillus gasseri K7 was less resistant than B. animalis subsp. lactis but survived at cell concentrations approximately 1000 times higher than other bifidobacteria.

Conclusion

In this study, we were able to show that L. gasseri K7 had a high survival rate in the stomach-intestine passage. By comparing the results with a previous study in piglets we could confirm the reliability of our simulation. Of the tested bifidobacteria strains, only B. animalis subsp. lactis showed acceptable survival for a successful passage in the simulation system.

Detection of the pediocin gene pedA in strains from human faeces by real-time PCR and characterization of Pediococcus acidilactici UVA1

Background

Bacteriocin-producing lactic acid bacteria are commonly used as natural protective cultures. Among them, strains of the genus Pediococcus are particularly interesting for their ability to produce pediocin, a broad spectrum antimicrobial peptide with a strong antagonistic activity against the food-borne pathogen Listeria monocytogenes. Furthermore, there is increasing interest in isolating new bacteriocin-producing strains of human intestinal origin that could be developed for probiotic effects and inhibition of pathogenic bacteria in the gut. In this work, we typed a new strain, co-isolated from baby faeces together with a Bifidobacterium thermophilum strain, and characterized its proteinaceous compound with strong antilisterial activity.

Results

The newly isolated strain UVA1 was identified as a Pediococcus acidilactici by carbohydrate fermentation profile, growth at 50°C and 16S rDNA sequencing. The partially purified bacteriocin was heat resistant up to 100°C, active over a wide range of pH (2 to 9) and susceptible to proteolytic enzymes. The molecular weight, estimated by SDS-PAGE, was similar to that of pediocin AcH/PA-1 (4.5 kDa). P. acidilactici UVA1 harboured a 9.5-kb plasmid that could be cured easily, which resulted in the loss of the antimicrobial activity. Southern hybridization using the DIG-labelled pedA-probe established that the bacteriocin gene was plasmid-borne as for all pediocin described so far. Nucleotide sequence of the whole operon (3.5 kb) showed almost 100 % similarity to the pediocin AcH/PA-1 operon. The mRNA transcript for pedA could be detected in P. acidilactici UVA1 but not in the cured derivative, confirming the expression of the pedA-gene in UVA1. Using a new real-time PCR assay, eleven out of seventeen human faecal samples tested were found to contain pedA-DNA.

Conclusion

We identified and characterised the first pediocin produced by a human intestinal Pediococcus acidilactici isolate and successfully developed a new real-time PCR assay to show the large distribution of pedA-containing strains in baby faecal samples.

Classification of a moderately oxygen-tolerant isolate from baby faeces as Bifidobacterium thermophilum

Background

Bifidobacteria are found at varying prevalence in human microbiota and seem to play an important role in the human gastrointestinal tract (GIT). Bifidobacteria are highly adapted to the human GIT which is reflected in the genome sequence of a Bifidobacterim longum isolate. The competitiveness against other bacteria is not fully understood yet but may be related to the production of antimicrobial compounds such as bacteriocins. In a previous study, 34 Bifidobacterium isolates have been isolated from baby faeces among which six showed proteinaceous antilisterial activity against Listeria monocytogenes. In this study, one of these isolates, RBL67, was further identified and characterized.

Results

Bifidobacterium isolate RBL67 was classified and characterized using a polyphasic approach. RBL67 was classified as Bifidobacterium thermophilum based on phenotypic and DNA-DNA hybridization characteristics, although 16S rDNA analyses and partial groEL sequences showed higher homology with B. thermacidophilum subsp. porcinum and B. thermacidophilum subsp. thermacidophilum, respectively. RBL67 was moderately oxygen-tolerant and was able to grow at pH 4 and at a temperature of 47°C.

Conclusion

In order to assign RBL67 to a species, a polyphasic approach was used. This resulted in the classification of RBL67 as a Bifidobacterium thermophilum strain. To our knowledge, this is the first report about B. thermophilum isolated from baby faeces since the B. thermophilum strains were related to ruminants and swine faeces before. B. thermophilum was previously only isolated from animal sources and was therefore suggested to be used as differential species between animal and human contamination. Our findings may disapprove this suggestion and further studies are now conducted to determine whether B. thermophilum is distributed broader in human faeces. Furthermore, the postulated differentiation between human and animal strains by growth above 45°C is no longer valid since B. thermophilum is able to grow at 47°C. In our study, 16S rDNA and partial groEL sequence analysis were not able to clearly assign RBL67 to a species and were contradictory. Our study suggests that partial groEL sequences may not be reliable as a single tool for species differentiation.

Detection of antifungal properties in Lactobacillus paracasei subsp. paracasei SM20, SM29, and SM63 and molecular typing of the strains

Lactobacilli isolated from different food and feed samples such as raw milk, cheese, yoghurt, olives, sour dough, as well as corn and grass silage, were screened for their antifungal activities. Out of 1,424 isolates tested, 82 were shown to be inhibitory to different yeasts (Candida spp. and Zygosaccharomyces bailii) and a Penicillium sp., which were previously isolated from spoiled yoghurt and fruits. Carbohydrate fermentation patterns suggested that a substantial portion, 25%, belonged to the Lactobacillus casei group, including L. casei, L. paracasei, and L. rhamnosus. The isolates SM20 (DSM14514), SM29 (DSM14515), and SM63 (DSM14516) were classified by PCR using species-specific primers to target the corresponding type strains (L. casei, L. paracasei, and L. rhamnosus) as controls. Further molecular typing methods such as randomly amplified polymorphic DNA, pulsed-field gel electrophoresis, and sequencing analysis of the 16S rRNA gene allowed classifying strains SM20, SM29, and SM63 as L. paracasei subsp. paracasei in accordance with the new reclassification of the L. casei group proposed by Collins et al.