Intragenomic 16S rRNA gene heterogeneity in Lactococcus lactis subsp. cremoris

High-throughput qPCR and 16S rRNA gene amplicon sequencing as complementary methods for the investigation of the cheese microbiota

Background

Next-generation sequencing (NGS) methods and especially 16S rRNA gene amplicon sequencing have become indispensable tools in microbial ecology. While they have opened up new possibilities for studying microbial communities, they also have one drawback, namely providing only relative abundances and thus compositional data. Quantitative PCR (qPCR) has been used for years for the quantification of bacteria. However, this method requires the development of specific primers and has a low throughput. The constraint of low throughput has recently been overcome by the development of high-throughput qPCR (HT-qPCR), which allows for the simultaneous detection of the most prevalent bacteria in moderately complex systems, such as cheese and other fermented dairy foods. In the present study, the performance of the two approaches, NGS and HT-qPCR, was compared by analyzing the same DNA samples from 21 Raclette du Valais protected designation of origin (PDO) cheeses. Based on the results obtained, the differences, accuracy, and usefulness of the two approaches were studied in detail.

Results

The results obtained using NGS (non-targeted) and HT-qPCR (targeted) show considerable agreement in determining the microbial composition of the cheese DNA samples studied, albeit the fundamentally different nature of these two approaches. A few inconsistencies in species detection were observed, particularly for less abundant ones. The detailed comparison of the results for 15 bacterial species/groups measured by both methods revealed a considerable bias for certain bacterial species in the measurements of the amplicon sequencing approach. We identified as probable origin to this PCR bias due to primer mismatches, variations in the number of copies for the 16S rRNA gene, and bias introduced in the bioinformatics analysis.

Conclusion

As the normalized microbial composition results of NGS and HT-qPCR agreed for most of the 21 cheese samples analyzed, both methods can be considered as complementary and reliable for studying the microbial composition of cheese. Their combined application proved to be very helpful in identifying potential biases and overcoming methodological limitations in the quantitative analysis of the cheese microbiota.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12866-022-02451-y.

Lactococcus lactis Diversity Revealed by Targeted Amplicon Sequencing of purR Gene, Metabolic Comparisons and Antimicrobial Properties in an Undefined Mixed Starter Culture Used for Soft-Cheese Manufacture

The undefined mixed starter culture (UMSC) is used in the manufacture of cheeses. Deciphering UMSC microbial diversity is important to optimize industrial processes. The UMSC was studied using culture-dependent and culture-independent based methods. MALDI-TOF MS enabled identification of species primarily from the Lactococcus genus. Comparisons of carbohydrate metabolism profiles allowed to discriminate five phenotypes of Lactococcus (n = 26/1616). The 16S sequences analysis (V1–V3, V3–V4 regions) clustered the UMSC microbial diversity into two Lactococcus operational taxonomic units (OTUs). These clustering results were improved with the DADA2 algorithm on the housekeeping purR sequences. Five L. lactis variants were detected among the UMSC. The whole-genome sequencing of six isolates allowed for the identification of the lactis subspecies using Illumina^® (n = 5) and Pacbio^® (n = 1) technologies. Kegg analysis confirmed the L. lactis species-specific niche adaptations and highlighted a progressive gene pseudogenization. Then, agar spot tests and agar well diffusion assays were used to assess UMSC antimicrobial activities. Of note, isolate supernatants (n = 34/1616) were shown to inhibit the growth of Salmonella ser. Typhimurium CIP 104115, Lactobacillus sakei CIP 104494, Staphylococcus aureus DSMZ 13661, Enterococcus faecalis CIP103015 and Listeria innocua CIP 80.11. Collectively, these results provide insightful information about UMSC L. lactis diversity and revealed a potential application as a bio-protective starter culture.

Identification of single nucleotide polymorphisms (SNPs) in the 16S rRNA gene of foodborne Bacillus spp

The main goal of this work was the identification of single nucleotide polymorphisms (SNPs) in the 16S rRNA gene of foodborne Bacillus spp. that may be useful for typing purposes. These species include, among others, Bacillus cereus, an important pathogenic species involved in food poisoning, and Bacillus licheniformis, Bacillus subtilis and Bacillus pumilus, which are causative agents of food spoilage described as responsible for foodborne disease outbreaks. With this purpose in mind, 52 Bacillus strains isolated from culture collections and fresh and processed food were considered. SNP type "Y" at sites 212 and 476 appeared in the majority of B. licheniformis studied strains. SNP type "R" at site 278 was detected in many strains of the B. subtilis/Bacillus amyloliquefaciens group, while polymorphism "Y" at site 173 was characteristic of the majority of strains of B. cereus/Bacillus thuringiensis group. The analysis of SNPs provided more intra-specific information than phylogenetic analysis in the cases of B. cereus and B. subtilis. Moreover, this study describes novel SNPs that should be considered when designing 16S rRNA-based primers and probes for multiplex-PCR, Real-Time PCR and microarray systems for foodborne Bacillus spp.

Development of a Lactobacillus specific T-RFLP method to determine lactobacilli diversity in complex samples

Terminal restriction fragment length polymorphism (T-RFLP) analysis has been widely used for studying microbial communities. However, most T-RFLP assays use 16S rDNA as the target and are unable to accurately characterize a microbial subpopulation. In this study, we developed a novel T-RFLP protocol based on Lactobacillus hsp60 to rapidly characterize and compare lactobacilli composition. The theoretical terminal restriction fragment (TRF) profiles were calculated from 769 Lactobacillus hsp60 sequences from online databases. In silico digestion with restriction endonucleases AluI and TacI on hsp60 amplicons generated 83 distinct TRF patterns, of which, 70 were species specific. To validate the assay, five previously sequenced lactobacilli were cultured independently, mixed at known concentrations and subjected to analysis by T-RFLP. All five strains generated the predicted TRFs and a qualitative consistent relationship was revealed. We performed the T-RFLP protocol on fecal samples from mice fed 6 different diets (n=4). Principal component analysis and agglomerative hierarchical clustering revealed that the lactobacilli community was strongly connected to dietary supplementation. Our study demonstrates the potential for using Lactobacillus specific T-RFLP to characterize lactobacilli communities in complex samples.