Comparative genomics analysis of genus Leuconostoc resolves its taxonomy and elucidates its biotechnological importance

The Metabolism of Leuconostoc Genus Decoded by Comparative Genomics

Leuconostoc encompasses a number of species that frequently appear in foods where they play different roles, ranging from ripening to spoiling. The number of available Leuconostoc genomes has recently increased and enabled the precise taxonomic and phylogenetic delineation of species. Nonetheless, a thorough investigation of the functions and the metabolic potential of Leuconostoc species has never been accomplished. In this study, all the currently available 553 Leuconostoc genomes were downloaded from NCBI GenBank and annotated utilizing specific tools in order to reconstruct the metabolic potential of the genus in terms of carbohydrate hydrolysis and fermentative pathways, transporters, and anabolic potential. The analysis revealed that species cluster based on their metabolic potential, showing unique adaptation and ecological roles. Pentose phosphate and phosphoketolase pathways were highlighted as the main ones of central metabolism. The various identified PTS and ABC transporters showed adaptability to different sugars. The metabolic diversity described in this study not only supports the role of Leuconostoc spp. in natural ecosystems but also highlights their potential in industrial applications, particularly in the fermentation industry where their ability to metabolize a wide range of substrates can be harnessed for the production of various fermented foods and bioproducts.

Synergistic effect of the coculture of Leuconostoc pseudomesenteroides and Lactococcus lactis, isolated from honeybees, on the generation of plant-based dairy alternatives based on soy, pea, oat, and potato drinks

The production of plant-based dairy alternatives has been majorly focused on the improvement of sensorial, technological and nutritional properties, to be able to mimic and replace milk-based fermented products. The presence of off-flavours and antinutrients, the lack of production of dairy-like flavours or the metabolic inaccessibility of plant proteins are some of the challenges to overcome to generate plant-based dairy alternatives. However, in the present study, it is demonstrated how the synergistic effect of two LAB strains, when cocultured, can simultaneously solve those challenges when fermenting in four different plant-based raw materials: soy, pea, oat, and potato drinks (SPOP). The fermentation was performed through the mono- and co-culture of the two LAB strains isolated from Apis mellifera (honeybee): Leuconostoc pseudomesenteroides NFICC 2004 and Lactococcus lactis NFICC 2005. Firstly, the coculture of both strains demonstrated to increase the acidification rate of the four plant matrices. Moreover, L. pseudomesenteroides (LP) demonstrated to in situ produce high concentrations of mannitol when fructose was present as C-source. Furthermore, L. pseudomesenteroides, which encoded for PII-proteinase, demonstrated to break down SPOP proteins, releasing free amino acids that were used by L.lactis (LL) for growth and metabolism. Lastly, the analysis of their co-metabolic volatile performance showed the principal ability of removal of the main off-flavours found in SPOP, such as hexanal, 1-octen-3-ol, 2-pentylfuran, pentanal, octanal, heptanal, and nonanal, mainly led by L. pseudomesenteroides, as well as the production of dairy-like flavours, such as diacetyl and 3-methyl-1-butanol, triggered by L. lactis metabolism. Overall, these findings endorsed the use of honeybee isolated strains as starter cultures, demonstrated the potential of coupling genotypes and phenotypes of multiple strains to improve the organoleptic properties suggesting a potential of combining plant-based matrices for the generation of future high-quality plant-based dairy alternatives.

Metagenomic Insights into the Taxonomic and Functional Features of Traditional Fermented Milk Products from Russia

Fermented milk products (FMPs) contain probiotics that are live bacteria considered to be beneficial to human health due to the production of various bioactive molecules. In this study, nine artisanal FMPs (kefir, ayran, khurunga, shubat, two cottage cheeses, bryndza, khuruud and suluguni-like cheese) from different regions of Russia were characterized using metagenomics. A metagenomic sequencing of ayran, khurunga, shubat, khuruud and suluguni-like cheese was performed for the first time. The taxonomic profiling of metagenomic reads revealed that Lactococcus species, such as Lc. lactis and Lc. cremoris prevailed in khuruud, bryndza, one sample of cottage cheese and khurunga. The latter one together with suluguni-like cheese microbiome was dominated by bacteria, affiliated to Lactobacillus helveticus (32-35%). In addition, a high proportion of sequences belonging to the genera Lactobacillus, Lactococcus and Streptococcus but not classified at the species level were found in the suluguni-like cheese. Lactobacillus delbrueckii, as well as Streptococcus thermophilus constituted the majority in another cottage cheese, kefir and ayran metagenomes. The microbiome of shubat, produced from camel's milk, was significantly distinctive, and Lentilactobacillus kefiri, Lactobacillus kefiranofaciens and Bifidobacterium mongoliense represented the dominant components (42, 7.4 and 5.6%, respectively). In total, 78 metagenome-assembled genomes with a completeness ≥ 50.2% and a contamination ≤ 8.5% were recovered: 61 genomes were assigned to the Enterococcaceae, Lactobacillaceae and Streptococcaceae families (the Lactobacillales order within Firmicutes), 4 to Bifidobacteriaceae (the Actinobacteriota phylum) and 2 to Acetobacteraceae (the Proteobacteria phylum). A metagenomic analysis revealed numerous genes, from 161 to 1301 in different products, encoding glycoside hydrolases and glycosyltransferases predicted to participate in lactose, alpha-glucans and peptidoglycan hydrolysis as well as exopolysaccharides synthesis. A large number of secondary metabolite biosynthetic gene clusters, such as lanthipeptides, unclassified bacteriocins, nonribosomal peptides and polyketide synthases were also detected. Finally, the genes involved in the synthesis of bioactive compounds like β-lactones, terpenes and furans, nontypical for fermented milk products, were also found. The metagenomes of kefir, ayran and shubat was shown to contain either no or a very low count of antibiotic resistance genes. Altogether, our results show that traditional indigenous fermented products are a promising source of novel probiotic bacteria with beneficial properties for medical and food industries.

Synergistic Immunostimulatory Activities of Probiotic Strains, Leuconostoc lactis and Weissella cibaria, and the Prebiotic Oligosaccharides They Produce

Synbiotics contain health-beneficial bacteria, i.e., probiotics and prebiotics selectively utilized by the probiotics. Herein, three probiotic strains, Leuconostoc lactis CCK940, L. lactis SBC001, and Weissella cibaria YRK005, and the oligosaccharides produced by these strains (CCK, SBC, and YRK, respectively) were used to prepare nine synbiotic combinations. Macrophages (RAW 264.7) were treated with these synbiotic combinations and the corresponding lactic acid bacteria and oligosaccharides alone to evaluate the treatments' immunostimulatory activities. The level of nitric oxide (NO) production was significantly higher in the macrophages treated with the synbiotics than in those treated with the corresponding probiotic strains and the oligosaccharide alone. The immunostimulatory activities of the synbiotics increased regardless of the probiotic strain and the type of oligosaccharide used. The expressions of tissue necrosis factor-α, interleukin-1β, cyclooxygenase-2, inducible NO synthase genes, and extracellular-signal-regulated and c-Jun N-terminal kinases were significantly higher in the macrophages treated with the three synbiotics than in those treated with the corresponding strains or with the oligosaccharides alone. These results indicate that the synergistic immunostimulatory activities of probiotics and the prebiotics they produced in the studied synbiotic preparations resulted from the activation of the mitogen-activated protein-kinase-signaling pathway. This study suggests the combined use of these probiotics and prebiotics in the development of synbiotic preparations as health supplements.

Metataxonomic insights in the distribution of Lactobacillaceae in foods and food environments

Members of the family Lactobacillaceae, which now includes species formerly belonging to the genera Lactobacillus and Pediococcus, but also Leuconostocaceae, are of foremost importance in food fermentations and spoilage, but also as components of animal and human microbiota and as potentially pathogenic microorganisms. Knowledge of the ecological distribution of a given species and genus is important, among other things, for the inclusion in lists of microorganisms with a Qualified Presumption of Safety or with beneficial use. The objective of this work is to use the data in FoodMicrobionet database to obtain quantitative insights (in terms of both abundance and prevalence) on the distribution of these bacteria in foods and food environments. We first explored the reliability of taxonomic assignments using the SILVA v138.1 reference database with full length and partial sequences of the 16S rRNA gene for type strain sequences. Full length 16S rRNA gene sequences allow a reasonably good classification at the genus and species level in phylogenetic trees but shorter sequences (V1-V3, V3-V4, V4) perform much worse, with type strains of many species sharing identical V4 and V3-V4 sequences. Taxonomic assignment at the genus level of 16S rRNA genes sequences and the SILVA v138.1 reference database can be done for almost all genera of the family Lactobacillaceae with a high degree of confidence for full length sequences, and with a satisfactory level of accuracy for the V1-V3 regions. Results for the V3-V4 and V4 region are still acceptable but significantly worse. Taxonomic assignment at the species level for sequences for the V1-V3, V3-V4, V4 regions of the 16S rRNA gene of members of the family Lactobacillaceae is hardly possible and, even for full length sequences, and only 49.9 % of the type strain sequences can be unambiguously assigned to species. We then used the FoodMicrobionet database to evaluate the prevalence and abundance of Lactobacillaceae in food samples and in food related environments. Generalist and specialist genera were clearly evident. The ecological distribution of several genera was confirmed and insights on the distribution and potential origin of rare genera (Dellaglioa, Holzapfelia, Schleiferilactobacillus) were obtained. We also found that combining Amplicon Sequence Variants from different studies is indeed possible, but provides little additional information, even when strict criteria are used for the filtering of sequences.

Identification of Leuconostoc species based on novel marker genes identified using real-time PCR via computational pangenome analysis

Leuconostoc species are important microorganisms in food fermentation but also cause food spoilage. Although these species are commercially important, their taxonomy is still based on inaccurate identification methods. Here, we used computational pangenome analysis to develop a real-time PCR-based method for identifying and differentiating the 12 major Leuconostoc species found in food. Analysis of pan and core-genome phylogenies showed clustering of strains into 12 distinct groups according to the species. Pangenome analysis of 130 Leuconostoc genomes from these 12 species enabled the identification of each species-specific gene. In silico testing of the species-specific genes against 143 publicly available Leuconostoc and 100 other lactic acid bacterial genomes showed that all the assays had 100% inclusivity/exclusivity. We also verified the specificity for each primer pair targeting each specific gene using 23 target and 124 non-target strains and found high specificity (100%). The sensitivity of the real-time PCR method was 10² colony forming units (CFUs)/ml in pure culture and spiked food samples. All standard curves showed good linear correlations, with an R² value of ≥0.996, suggesting that screened targets have good specificity and strong anti-interference ability from food sample matrices and non-target strains. The real-time PCR method can be potentially used to determine the taxonomic status and identify the Leuconostoc species in foods.

Comparative genome analysis of four Leuconostoc strains with a focus on carbohydrate-active enzymes and oligosaccharide utilization pathways

•

Comparative genomic analysis of four Leuconostoc strains was performed.

•

Leuconostoc spp. shared genomic similarity, but their genetic content differed.

•

Leuconostoc spp. showed different genes encoding CAZymes.

•

Oligosaccharide’s utilization and folate biosynthesis pathways were investigated.

Leuconostoc is mostly found in food, plants, and dairy products. Due to their innate genomic features, such as the presence of carbohydrate-active enzymes, bacteriocins, and plasmids, Leuconostoc spp. have great biotechnological potential. In this study, four strains were isolated and identified as Leuconostoc mesenteroides SG315 (LA), L. citreum SG255 (LB), L. lactis CCK940 (LC), and L. lactis SBC001 (LD). Comparative analysis was performed using their draft genome sequences. Differences among the four strains were analyzed using the average nucleotide identity, dot plot, and multiple alignments of conserved genomic sequences. Functional profiling revealed 2134, 1917, 1751, and 1816 open reading frames; 2023, 1823, 1655, and 1699 protein-coding genes; 60, 57, 83, and 82 RNA-coding genes; and GC content of 37.5 %, 38.8 %, 43.3 %, and 43.2 %, in LA, LB, LC, and LD, respectively. The total number of genes encoding carbohydrate-active enzymes was 76 (LA), 73 (LB), 57 (LC), and 67 (LD). These results indicate that the four strains shared a large number of genes, but their gene content is different. Furthermore, most genes with unknown functions were observed in the prophage regions of the genome. This study also elucidated the oligosaccharide utilization and folate biosynthesis pathways in Leuconostoc spp. Taken together, our findings provide useful information on the genomic diversity of CAZymes in the four Leuconostoc strains and suggest that these species could be used for potent exploitation.

Phylogenomic analysis of the genus Leuconostoc

Leuconostoc is a genus of saccharolytic heterofermentative lactic acid bacteria that inhabit plant-derived matrices and a variety of fermented foods (dairy products, dough, milk, vegetables, and meats), contributing to desired fermentation processes or playing a role in food spoilage. At present, the genus encompasses 17 recognized species. In total, 216 deposited genome sequences of Leuconostoc were analyzed, to check the delineation of species and to infer their evolutive genealogy utilizing a minimum evolution tree of Average Nucleotide Identity (ANI) and the core genome alignment. Phylogenomic relationships were compared to those obtained from the analysis of 16S rRNA, pheS, and rpoA genes. All the phylograms were subjected to split decomposition analysis and their topologies were compared to check the ambiguities in the inferred phylogenesis. The minimum evolution ANI tree exhibited the most similar topology with the core genome tree, while single gene trees were less adherent and provided a weaker phylogenetic signal. In particular, the 16S rRNA gene failed to resolve several bifurcations and Leuconostoc species. Based on an ANI threshold of 95%, the organization of the genus Leuconostoc could be amended, redefining the boundaries of the species L. inhae, L. falkenbergense, L. gelidum, L. lactis, L. mesenteroides, and L. pseudomesenteroides. Two strains currently recognized as L. mesenteroides were split into a separate lineage representing a putative species (G16), phylogenetically related to both L. mesenteroides (G18) and L. suionicum (G17). Differences among the four subspecies of L. mesenteroides were not pinpointed by ANI or by the conserved genes. The strains of L. pseudomesenteroides were ascribed to two putative species, G13 and G14, the former including also all the strains presently belonging to L. falkenbergense. L. lactis was split into two phylogenetically related lineages, G9 and G10, putatively corresponding to separate species and both including subgroups that may correspond to subspecies. The species L. gelidum and L. gasicomitatum were closely related but separated into different species, the latter including also L. inhae strains. These results, integrating information of ANI, core genome, and housekeeping genes, complemented the taxonomic delineation with solid information on the phylogenetic lineages evolved within the genus Leuconostoc.