A genome-scale metabolic network of the aroma bacterium Leuconostoc mesenteroides subsp. cremoris

Pangenome reconstruction of Lactobacillaceae metabolism predicts species-specific metabolic traits

Strains across the Lactobacillaceae family form the basis for a trillion-dollar industry. Our understanding of the genomic basis for their key traits is fragmented, however, including the metabolism that is foundational to their industrial uses. Pangenome analysis of publicly available Lactobacillaceae genomes allowed us to generate genome-scale metabolic network reconstructions for 26 species of industrial importance. Their manual curation led to more than 75,000 gene-protein-reaction associations that were deployed to generate 2,446 genome-scale metabolic models. Cross-referencing genomes and known metabolic traits allowed for manual metabolic network curation and validation of the metabolic models. As a result, we provide the first pangenomic basis for metabolism in the Lactobacillaceae family and a collection of predictive computational metabolic models that enable a variety of practical uses.IMPORTANCELactobacillaceae, a bacterial family foundational to a trillion-dollar industry, is increasingly relevant to biosustainability initiatives. Our study, leveraging approximately 2,400 genome sequences, provides a pangenomic analysis of Lactobacillaceae metabolism, creating over 2,400 curated and validated genome-scale models (GEMs). These GEMs successfully predict (i) unique, species-specific metabolic reactions; (ii) niche-enriched reactions that increase organism fitness; (iii) essential media components, offering insights into the global amino acid essentiality of Lactobacillaceae; and (iv) fermentation capabilities across the family, shedding light on the metabolic basis of Lactobacillaceae-based commercial products. This quantitative understanding of Lactobacillaceae metabolic properties and their genomic basis will have profound implications for the food industry and biosustainability, offering new insights and tools for strain selection and manipulation.

Unraveling the carbohydrate metabolic characteristics of Leuconostoc mesenteroides J18 through metabolite and transcriptome analyses

The metabolic characteristics of Leuconostoc mesenteroides subsp. mesenteroides J18, which is mainly responsible for kimchi fermentation, on various carbon sources were investigated through carbon utilization, metabolite, and transcriptome analyses at different culture conditions (10 and 30 °C with/without 2.5% NaCl). The metabolic features of strain J18 were relatively similar across the four culture conditions. However, the metabolic characteristics of strain J18 showed significant variations depending on the carbon source. These distinct metabolic traits of strain J18 on various carbon sources were validated through transcriptomic analyses and the reconstruction of metabolic pathways. The transcriptional expression of the metabolic pathways in response to each carbon source consistently correlated with the production profiles of metabolites, including ethanol, acetoin, diacetyl, and riboflavin, in each carbon source. Our findings suggests that the abundance of Leu. mesenteroides during fermentation and the taste and flavor of fermented food products can be controlled by altering the carbon sources.

Utilization of plant derived lactic acid bacteria for efficient bioconversion of brewers' spent grain into acetoin

Brewers' spent grain (BSG) is a major side-stream from the beer industry, with an annual estimated production of 39 million tons worldwide. Due to its high nutritional value, high abundance and low price, it has been proposed as an ingredient in human food. Here we investigated the ability of different lactic acid bacteria to produce the flavor molecule acetoin in liquid BSG extract, in order to broaden the possibilities of utilization of BSG in human food. All the investigated lactic acid bacteria (LAB) covering the Leuconostoc, Lactobacillus and Lactoccocus species were able to convert the fermentable sugars in liquid BSG into acetoin. Production levels varied significantly between the different LAB species, with Leuconostoc pseudomesenteroides species reaching the highest titers of acetoin with only acetate as the main byproduct, while also being the fastest consumer of the fermentable sugars present in liquid BSG. Surprisingly, the currently best investigated LAB for acetoin production, L. lactis, was unable to consume the maltose fraction of liquid BSG and was therefore deemed unfit for full conversion of the sugars in BSG into acetoin. The production of acetoin in Leu. pseudomesenteroides was pH dependent as previously observed in other LAB, and the conversion of BSG into acetoin was scalable from shake flasks to 1 L bioreactors. While all investigated LAB species produced acetoin under aerobic conditions, Leu. pseudomesenteroides was found to be an efficient and scalable organism for bioconversion of liquid BSG into a safe acetoin rich food additive.

Potential of Cheese-Associated Lactic Acid Bacteria to Metabolize Citrate and Produce Organic Acids and Acetoin

Lactic acid bacteria (LAB) are pivotal in shaping the technological, sensory, and safety aspects of dairy products. The evaluation of proteolytic activity, citrate utilization, milk pH reduction, and the production of organic compounds, acetoin, and diacetyl by cheese associated LAB strains was carried out, followed by Principal Component Analysis (PCA). Citrate utilization was observed in all Leuconostoc (Le.) mesenteroides, Le. citreum, Lactococcus (Lc.) lactis, Lc. garvieae, and Limosilactobacillus (Lm.) fermentum strains, and in some Lacticaseibacillus (Lact.) casei strains. Most strains exhibited proteolytic activity, reduced pH, and generated organic compounds. Multivariate PCA revealed Le. mesenteroides as a prolific producer of acetic, lactic, formic, and pyruvic acids and acetoin at 30 °C. Enterococcus sp. was distinguished from Lact. casei based on acetic, formic, and pyruvic acid production, while Lact. casei primarily produced lactic acid at 37 °C. At 42 °C, Lactobacillus (L.) helveticus and some L. delbrueckii subsp. bulgaricus strains excelled in acetoin production, whereas L. delbrueckii subsp. bulgaricus and Streptococcus (S.) thermophilus strains primarily produced lactic acid. Lm. fermentum stood out with its production of acetic, formic, and pyruvic acids. Overall, cheese-associated LAB strains exhibited diverse metabolic capabilities which contribute to desirable aroma, flavor, and safety of dairy products.

Leuconostoc performance in soy-based fermentations - Survival, acidification, sugar metabolism, and flavor comparisons

Leuconostoc spp. is often regarded as the flavor producer, responsible for the production of acetoin and diacetyl in dairy cheese. In this study, we investigate seven plant-derived Leuconostoc strains, covering four species, in their potential as a lyophilized starter culture for flavor production in fermented soy-based cheese alternatives. We show that the process of lyophilization of Leuconostoc can be feasible using a soy-based lyoprotectant, with survivability up to 63% during long term storage. Furthermore, the storage in this media improves the subsequent growth in a soy-based substrate in a strain specific manner. The utilization of individual raffinose family oligosaccharides was strain dependent, with Leuconostoc pseudomesenteroides NFICC99 being the best consumer. Furthermore, we show that all investigated strains were able to produce a range of volatile flavor compounds found in dairy cheese products, as well as remove certain dairy off-flavors from the soy-based substrate like hexanal and 2-pentylfuran. Also here, NFICC99 was strain producing most cheese-related volatile flavor compounds, followed by Leuconostoc mesenteroides NFICC319. These findings provide initial insights into the development of Leuconostoc as a potential starter culture for plant-based dairy alternatives, as well as a promising approach for generation of stable, lyophilized cultures.

Differences in Aroma Metabolite Profile, Microstructure, and Rheological Properties of Fermented Milk Using Different Cultures

Texture and flavour are the key attributes determining sensory quality and are highly affected by starter cultures. A selection of phenotypic strains is needed to create diverse texture and flavour to meet consumers' preferences. In this study, the use of five lactic acid bacteria strains in the production of fermented milk, along with the metabolite profiles, microstructure, and rheological properties of the fermented milk samples, was investigated. Our results showed that Lactobacillus helveticus (LH) and Streptococcus thermophilus (ST) had a stronger acidification during fermentation but resulted in products with a coarser protein network compared to Lactococcus lactis (BL1) and Leuconostoc mesenteroides (CL3). Milk fermented by LH had the highest viscosity and exopolysaccharide concentration, while milk fermented by ST had the highest concentration of diacetyl. Although Leuconostoc pseudomesenteroides (CL3ST) had a minimal acidification capability, it produced high levels of ethyl-derived compounds associated with sweet, fruity, and floral fragrances. The results demonstrated that LH and ST could be used as starter cultures targeting fermented milks with different viscosities, while BL1, CL3, and CL3ST are suitable as adjunct cultures to impact different acidic sharpness and flavour notes.

Improving acetoin production through construction of a genome-scale metabolic model

Acetoin was widely used in food, medicine, and other industries, because of its unique fragrance. Bacillus amyloliquefaciens was recognized as a safe strain and a promising acetoin producer in fermentation. However, due to the complexity of its metabolic network, it had not been fully utilized. Therefore, a genome-scale metabolic network model (iJYQ746) of B. amyloliquefaciens was constructed in this study, containing 746 genes, 1736 reactions, and 1611 metabolites. The results showed that Mg²⁺, Mn²⁺, and Fe²⁺ have inhibitory effects on acetoin. When the stirring speed was 400 rpm, the maximum titer was 49.8 g L^-1. Minimization of metabolic adjustments (MOMA) was used to identify potential metabolic modification targets 2-oxoglutarate aminotransferase (serC, EC 2.6.1.52) and glucose-6-phosphate isomerase (pgi, EC 5.3.1.9). These targets could effectively accumulate acetoin by increasing pyruvate content, and the acetoin synthesis rate was increased by 610% and 10%, respectively. This provides a theoretical basis for metabolic engineering to reasonably transform B. amyloliquefaciens and produce acetoin.

Microbial communities and metabolic pathways involved in reductive decolorization of an azo dye in a two-stage AD system

Multiple stage anaerobic system was found to be an effective strategy for reductive decolorization of azo dyes in the presence of sulfate. Bulk color removal (56-90%) was achieved concomitant with acidogenic activity in the 1st-stage reactor (R1), while organic matter removal (≤100%) and sulfate reduction (≤100%) occurred predominantly in the 2nd-stage reactor (R2). However, azo dye reduction mechanism and metabolic routes involved remain unclear. The involved microbial communities and conditions affecting the azo dye removal in a two-stage anaerobic digestion (AD) system were elucidated using amplicon sequencing (16S rRNA, fhs, dsrB and mcrA) and correlation analysis. Reductive decolorization was found to be co-metabolic and mainly associated with hydrogen-producing pathways. We also found evidence of the involvement of an azoreductase from Lactococcus lactis. Bacterial community in R1 was sensitive and shifted in the presence of the azo dye, while microorganisms in R2 were more protected. Higher diversity of syntrophic-acetate oxidizers, sulfate reducers and methanogens in R2 highlights the role of the 2nd-stage in organic matter and sulfate removals, and these communities might be involved in further transformations of the azo dye reduction products. The results improve our understanding on the role of different microbial communities in anaerobic treatment of azo dyes and can help in the design of better solutions for the treatment of textile effluents.

The Effect of Respiration, pH, and Citrate Co-Metabolism on the Growth, Metabolite Production and Enzymatic Activities of Leuconostoc mesenteroides subsp. cremoris E30

Leuconostoc mesenteroides includes strains used as starter and/or adjunct cultures for the production of several fermented foods. In this study, the effect of anaerobic and respiratory cultivations, as well as of citrate supplementation and different pH values, was evaluated on growth, biomass, metabolite, and enzymatic activities (pyruvate oxidase, POX; NADH-dependent oxidase, NOX; NADH-dependent peroxidase, NPR) of Leuconostoc mesenteroides subsp. cremoris E30. We compared the respiration-increased growth rate and biomass production of Leuc. mesenteroides E30 to anaerobic cultivation. A supplementation of citrate impaired the growth rate of the respiratory cells. As expected, anaerobic cultures did not consume oxygen, and a similar trend in oxygen uptake was observed in respiratory cultures. The aerobic incubation caused changes in the metabolic pattern, reducing the production of ethanol in favour of acetic acid. Citrate was already exhausted in the exponential phase and did not affect the yields in acetic acid and ethanol. NOX activity increased in the presence of oxygen, while catalase was also detected in the absence of hemin. The absence of H₂O₂ suggested its degradation by NPR and catalase. Respiratory cultivation provided benefits (increase in growth rate, biomass, and activity in antioxidant enzymes) for Leuc. mesenteroides E30. Therefore, the exploitation of respiratory phenotypes may be useful for the formulation of competitive starter or adjunct cultures.

Developments in effective use of volatile organic compound analysis to assess flavour formation during cheese ripening

In the burgeoning demand for optimization of cheese production, ascertaining cheese flavour formation during the cheese making process has been the focal point of determining cheese quality. In this research reflection, we have highlighted how valuable volatile organic compound (VOC) analysis has been in assessing contingent cheese flavour compounds arising from non-starter lactic acid bacteria (NSLAB) along with starter lactic acid bacteria (SLAB), and whether VOC analysis associated with other high-throughput data might help provide a better understanding the cheese flavour formation during cheese process. It is widely known that there is a keen interest to merge all omics data to find specific biomarkers and/or to assess aroma formation of cheese. Towards that end, results of VOC analysis have provided valuable insights into the cheese flavour profile. In this review, we are pinpointing the effective use of flavour compound analysis to perceive flavour-forming ability of microbial strains that are convenient for dairy production, intertwining microbiome and metabolome to unveil potential biomarkers that occur during cheese ripening. In doing so, we summarised the functionality and integration of aromatic compound analysis in cheese making and gave reflections on reconsidering what the role of flavour-based analysis might have in the future.

Phenotypic Traits and Immunomodulatory Properties of Leuconostoc carnosum Isolated From Meat Products

Twelve strains of Leuconostoc carnosum from meat products were investigated in terms of biochemical, physiological, and functional properties. The spectrum of sugars fermented by L. carnosum strains was limited to few mono- and disaccharides, consistently with the natural habitats of the species, including meat and fermented vegetables. The strains were able to grow from 4 to 37°C with an optimum of approximately 32.5°C. The ability to grow at temperatures compatible with refrigeration and in presence of up to 60 g/L NaCl explains the high loads of L. carnosum frequently described in many meat-based products. Six strains produced exopolysaccharides, causing a ropy phenotype of colonies, according to the potential involvement on L. carnosum in the appearance of slime in packed meat products. On the other side, the study provides evidence of a potential protective role of L. carnosum WC0321 and L. carnosum WC0323 against Listeria monocytogenes, consistently with the presence in these strains of the genes encoding leucocin B. Some meat-based products intended to be consumed without cooking may harbor up to 10⁸ CFU/g of L. carnosum; therefore, we investigated the potential impact of this load on health. No strains survived the treatment with simulated gastric juice. Three selected strains were challenged for the capability to colonize a mouse model and their immunomodulatory properties were investigated. The strains did not colonize the intestine of mice during 10 days of daily dietary administration. Intriguingly, despite the loss of viability during the gastrointestinal transit, the strains exhibited different immunomodulatory effect on the maturation of dendritic cells in vivo, the extent of which correlated to the production of exopolysaccharides. The ability to stimulate the mucosal associated immune system in such probiotic-like manner, the general absence of antibiotic resistance genes, and the lack of the biosynthetic pathways for biogenic amines should reassure on the safety of this species, with potential for exploitation of selected starters.

Genetic Diversity of Leuconostoc mesenteroides Isolates from Traditional Montenegrin Brine Cheese

In many dairy products, Leuconostoc spp. is a natural part of non-starter lactic acid bacteria (NSLAB) accounting for flavor development. However, data on the genomic diversity of Leuconostoc spp. isolates obtained from cheese are still scarce. The focus of this study was the genomic characterization of Leuconostoc spp. obtained from different traditional Montenegrin brine cheeses with the aim to explore their diversity and provide genetic information as a basis for the selection of strains for future cheese production. In 2019, sixteen Leuconostoc spp. isolates were obtained from white brine cheeses from nine different producers located in three municipalities in the northern region of Montenegro. All isolates were identified as Ln. mesenteroides. Classical multilocus sequence tying (MLST) and core genome (cg) MLST revealed a high diversity of the Montenegrin Ln. mesenteroides cheese isolates. All isolates carried genes of the bacteriocin biosynthetic gene clusters, eight out of 16 strains carried the citCDEFG operon, 14 carried butA, and all 16 isolates carried alsS and ilv, genes involved in forming important aromas and flavor compounds. Safety evaluation indicated that isolates carried no pathogenic factors and no virulence factors. In conclusion, Ln. mesenteroides isolates from Montenegrin traditional cheeses displayed a high genetic diversity and were unrelated to strains deposited in GenBank.

How water-soluble saccharides drive the metabolism of lactic acid bacteria during fermentation of brewers' spent grain

We proposed a novel phenomic approach to track the effect of short-term exposures of Lactiplantibacillus plantarum and Leuconostoc pseudomesenteroides to environmental pressure induced by brewers' spent grain (BSG)-derived saccharides. Water-soluble BSG-based medium (WS-BSG) was chosen as model system. The environmental pressure exerted by WS-BSG shifted the phenotypes of bacteria in species- and strains-dependent way. The metabolic drift was growth phase-dependent and likely underlay the diauxic profile of organic acids production by bacteria in response to the low availability of energy sources. Among pentosans, metabolism of arabinose was preferred by L. plantarum and xylose by Leuc. pseudomesenteroides as confirmed by the overexpression of related genes. Bayesian variance analysis showed that phenotype switching towards galactose metabolism suffered the greatest fluctuation in L. plantarum. All lactic acid bacteria strains utilized more intensively sucrose and its plant-derived isomers. Sucrose-6-phosphate activity in Leuc. pseudomesenteroides likely mediated the increased consumption of raffinose. The increased levels of some phenolic compounds suggested the involvement of 6-phospho-β-glucosidases in β-glucosides degradation. Expression of genes encoding β-glucoside/cellobiose-specific EII complexes and phenotyping highlighted an increased metabolism for cellobiose. Our reconstructed metabolic network will improve the understanding of how lactic acid bacteria may transform BSG into suitable food ingredients.

Role of Lactic Acid Bacteria Phospho-β-Glucosidases during the Fermentation of Cereal by-Products

Bioprocessing using lactic acid bacteria (LAB) is a powerful means to exploit plant-derived by-products as a food ingredient. LAB have the capability to metabolize a large variety of carbohydrates, but such metabolism only relies on few metabolic routes, conferring on them a high fermentation potential. One example of these pathways is that involving phospho-β-glucosidase genes, which are present in high redundancy within LAB genomes. This enzymatic activity undertakes an ambivalent role during fermentation of plant-based foods related to the release of a wide range of phenolic compounds, from their β-D-glycosylated precursors and the degradation of β-glucopyranosyl derived carbohydrates. We proposed a novel phenomic approach to characterize the metabolism drift of Lactiplantibacillus plantarum and Leuconostoc pseudomesenteroides caused by a lignocellulosic by-product, such as the brewers’ spent grain (BSG), in contrast to Rich De Man, Rogosa and Sharpe (MRS) broth. We observed an increased metabolic activity for gentiobiose, cellobiose and β-glucoside conjugates of phenolic compounds during BSG fermentation. Gene expression analysis confirmed the importance of cellobiose metabolism while a release of lignin-derived aglycones was found during BSG fermentation. We provided a comprehensive view of the important role exerted by LAB 6-phospho-β-glucosidases as well the major metabolic routes undertaken during plant-based fermentations. Further challenges will consider a controlled characterization of pbg gene expression correlated to the metabolism of β-glucosides with different aglycone moieties.

Dynamic co-culture metabolic models reveal the fermentation dynamics, metabolic capacities and interplays of cheese starter cultures

In this study, we have investigated the cheese starter culture as a microbial community through a question: can the metabolic behaviour of a co-culture be explained by the characterized individual organism that constituted the co-culture? To address this question, the dairy-origin lactic acid bacteria Lactococcus lactis subsp. cremoris, Lactococcus lactis subsp. lactis, Streptococcus thermophilus and Leuconostoc mesenteroides, commonly used in cheese starter cultures, were grown in pure and four different co-cultures. We used a dynamic metabolic modelling approach based on the integration of the genome-scale metabolic networks of the involved organisms to simulate the co-cultures. The strain-specific kinetic parameters of dynamic models were estimated using the pure culture experiments and they were subsequently applied to co-culture models. Biomass, carbon source, lactic acid and most of the amino acid concentration profiles simulated by the co-culture models fit closely to the experimental results and the co-culture models explained the mechanisms behind the dynamic microbial abundance. We then applied the co-culture models to estimate further information on the co-cultures that could not be obtained by the experimental method used. This includes estimation of the profile of various metabolites in the co-culture medium such as flavour compounds produced and the individual organism level metabolic exchange flux profiles, which revealed the potential metabolic interactions between organisms in the co-cultures.

Genome scale metabolic models and analysis for evaluating probiotic potentials

Probiotics are live beneficial microorganisms that can be consumed in the form of dairy and food products as well as dietary supplements to promote a healthy balance of gut bacteria in humans. Practically, the main challenge is to identify and select promising strains and formulate multi-strain probiotic blends with consistent efficacy which is highly dependent on individual dietary regimes, gut environments, and health conditions. Limitations of current in vivo and in vitro methods for testing probiotic strains can be overcome by in silico model guided systems biology approaches where genome scale metabolic models (GEMs) can be used to describe their cellular behaviors and metabolic states of probiotic strains under various gut environments. Here, we summarize currently available GEMs of microbial strains with probiotic potentials and propose a knowledge-based framework to evaluate metabolic capabilities on the basis of six probiotic criteria. They include metabolic characteristics, stability, safety, colonization, postbiotics, and interaction with the gut microbiome which can be assessed by in silico approaches. As such, the most suitable strains can be identified to design personalized multi-strain probiotics in the future.

Current Status and Applications of Adaptive Laboratory Evolution in Industrial Microorganisms

Adaptive laboratory evolution (ALE) is an evolutionary engineering approach in artificial conditions that improves organisms through the imitation of natural evolution. Due to the development of multi-level omics technologies in recent decades, ALE can be performed for various purposes at the laboratory level. This review delineates the basics of the experimental design of ALE based on several ALE studies of industrial microbial strains and updates current strategies combined with progressed metabolic engineering, in silico modeling and automation to maximize the evolution efficiency. Moreover, the review sheds light on the applicability of ALE as a strain development approach that complies with non-recombinant preferences in various food industries. Overall, recent progress in the utilization of ALE for strain development leading to successful industrialization is discussed.

A metabolic reconstruction of Lactobacillus reuteri JCM 1112 and analysis of its potential as a cell factory

Background

Lactobacillus reuteri is a heterofermentative Lactic Acid Bacterium (LAB) that is commonly used for food fermentations and probiotic purposes. Due to its robust properties, it is also increasingly considered for use as a cell factory. It produces several industrially important compounds such as 1,3-propanediol and reuterin natively, but for cell factory purposes, developing improved strategies for engineering and fermentation optimization is crucial. Genome-scale metabolic models can be highly beneficial in guiding rational metabolic engineering. Reconstructing a reliable and a quantitatively accurate metabolic model requires extensive manual curation and incorporation of experimental data.

Results

A genome-scale metabolic model of L. reuteri JCM 1112^T was reconstructed and the resulting model, Lreuteri_530, was validated and tested with experimental data. Several knowledge gaps in the metabolism were identified and resolved during this process, including presence/absence of glycolytic genes. Flux distribution between the two glycolytic pathways, the phosphoketolase and Embden–Meyerhof–Parnas pathways, varies considerably between LAB species and strains. As these pathways result in different energy yields, it is important to include strain-specific utilization of these pathways in the model. We determined experimentally that the Embden–Meyerhof–Parnas pathway carried at most 7% of the total glycolytic flux. Predicted growth rates from Lreuteri_530 were in good agreement with experimentally determined values. To further validate the prediction accuracy of Lreuteri_530, the predicted effects of glycerol addition and adhE gene knock-out, which results in impaired ethanol production, were compared to in vivo data. Examination of both growth rates and uptake- and secretion rates of the main metabolites in central metabolism demonstrated that the model was able to accurately predict the experimentally observed effects. Lastly, the potential of L. reuteri as a cell factory was investigated, resulting in a number of general metabolic engineering strategies.

Conclusion

We have constructed a manually curated genome-scale metabolic model of L. reuteri JCM 1112^T that has been experimentally parameterized and validated and can accurately predict metabolic behavior of this important platform cell factory.