Next Article in Journal
Rheological Characteristics of Model Gluten-Free Dough with Plantago Seeds and Husk Incorporation
Previous Article in Journal
High-Performance Thin-Layer Chromatography-Immunostaining as a Technique for the Characterization of Whey Protein Enrichment in Edam Cheese
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Foods
Volume 11
Issue 4
10.3390/foods11040535
foods-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
Annamaria Ricciardi
,
Livia Vanessa Storti
,
Marilisa Giavalisco
,
Eugenio Parente
* and
Teresa Zotta
Scuola di Scienze Agrarie, Alimentari, Forestali ed Ambientali (SAFE), Università degli Studi della Basilicata, 85100 Potenza, Italy
*
Author to whom correspondence should be addressed.
Foods 2022, 11(4), 535; https://doi.org/10.3390/foods11040535
Submission received: 22 December 2021 / Revised: 28 January 2022 / Accepted: 10 February 2022 / Published: 13 February 2022
(This article belongs to the Section Food Microbiology)
Browse Figures
Versions Notes

Abstract

:
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 H2O2 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.

1. Introduction

The genus Leuconostoc includes 14 species of great importance in food technology [1]. Some strains of Leuconostoc citreum, Leuc. mesenteroides, and Leuc. lactis have been exploited as starter and/or adjunct cultures to improve the quality and shelf-life of several fermented foods (i.e., kimchi, table olives, sourdough and bakery products, fermented milks, butter, cream, and cheeses) [2,3,4,5]. Leuconostoc strains, in fact, may contribute to the texture and aromatic profile of fermented foods by producing exopolysaccharides [6] and diacetyl, acetoin, acetate, ethanol and 2,3-butylenglycol [7]. However, some species (Leuc. gelidum subsp. gasicomitatum, Leuc. mesenteroides, Leuc. carnosum) are recognized as spoilage agents of meat, fishery, vegetable products, and ready-to-eat meals, as they impair organoleptic features through off-flavor formation, green discoloration, and the production of slime and CO2 [8,9,10,11].
Leuconostoc species are oxygen-tolerant anaerobes with heterofermentative metabolisms. The catabolism of carbohydrates through the phosphoketolase pathway, with the formation of acetyl phosphate as a key intermediate, leads to the production of CO2, lactic acid, acetic acid, or ethanol, depending on the NADH/NAD+ ratio. Heterofermentative species grow poorly when glucose is the sole carbon source, and other sugars, such as sucrose, fructose, and/or maltose, are metabolized for energy production and cell functionality [12]; fructose may be also used as an alternative electron acceptor for the regeneration of reduced NADH. Some members of Leuc. mesenteroides, additionally, are able to convert citrate into oxaloacetate and acetate via citrate lyase-oxaloacetate decarboxylase, supplying the intracellular pools of pyruvate and increasing the production of aroma compounds [13,14]. Like other lactic acid bacteria (LAB), some strains of Leuc. mesenteroides may tolerate oxygen and grow in aerobic conditions. Plihon et al. [15,16] demonstrated that aerobic cultivation of Leuc. mesenteroides improved biomass yield and affected carbon metabolism by shifting ethanol production to acetic acid accumulation.
Several authors [17,18,19] have shown that in some homofermentative LAB, the supplementation with hemin and menaquinone confers further physiological advantages (energy gain, robustness to oxidative and freeze-drying stresses, and synthesis of antioxidant enzymes) compared to unsupplemented aerobiosis.
To date, the respiratory metabolism has been investigated in a few heterofermentative strains of Levilactobacillus brevis [20], Limosilactobacillus reuteri, and Levilactobacillus spicheri [21], confirming the beneficial effects on growth fitness and stress survival. In Leuc. gasicomitatum [22] (actually reclassified as Leuc. gelidum subsp. gasicomitatum) [23], the heme-supplemented aerobiosis had a positive effect on growth rates, biomass accumulation, and the production of aroma compounds. More recently, Zotta et al. [24] verified that potential respiratory phenotypes are distributed in other heterofermentative species and strains, including several members of Leuc. mesenteroides. Among them, Leuc. mesenteroides E30 showed an improved growth and oxidative stress tolerance, as well as manganese-dependent catalase activity, when cultivated in aerobic and respiratory conditions.
In this study, the effect of the atmosphere of incubation on the fitness of Leuc. mesenteroides subsp. cremoris E30 has been further investigated. Specifically, the effect of anaerobic and respiratory cultivations, as well as of citrate metabolism and different pH values, has been evaluated on the growth performances (kinetics and biomass yield), metabolites production, and activities of several enzymes involved in the oxygen utilization and oxidative stress tolerance of Leuc. mesenteroides subsp. cremoris E30 in controlled batch cultivations.

2. Materials and Methods

2.1. Strains and Culture Conditions

Leuconostoc mesenteroides subsp. mesenteroides E30 (isolated from raw milk for the production of Canestrato cheese) was maintained as freeze-dried stock in reconstituted 11% (w/v) skim milk containing 0.1% (w/v) ascorbic acid, in the culture collection of the Laboratory of Industrial Microbiology, Università degli Studi della Basilicata, and was routinely propagated in WMB pH 6.8 [25] for 16 h at 30 °C.

2.2. Fermentation Conditions

The growth of Leuc. mesenteroides E30 was evaluated for 24 h at 30 °C in modified WMB (mWMB) with 20 g/L of glucose [26], with or without 5 g/L of sodium citrate, in batch cultivations carried out under anaerobic (AN; nitrogen flow at 0.02 L/min) and respiratory (RS; 15% dissolved oxygen concentration, dO2%, supplemented with 1.25 μg/mL of hemin and 1 μg/mL of menaquinone) conditions at pH 5.5 or 6.5, according to a complete 23 factorial design (Table 1).
Bioreactors (3 L working volume for AN and RS conditions, respectively; Applikon, Schiedam, The Netherlands) were inoculated (2% v/v) with an overnight (16 h, 30 °C) WMB anaerobic pre-culture, washed twice with 20 mM potassium phosphate buffer pH 7 (PB7), and standardized to a final absorbance at 650 nm (A650) of 3.0 (Bio-Rad Smart Spec™Plus, Bio-Rad Laboratories Inc., Milan, Italy).
dO2% was measured using a polarographic electrode (Applisens, Applikon Biotechnology, Delft, The Netherlands) and was automatically controlled (ezControl controller, Applikon; set point 15%) by varying the stirrer speed (impeller speed from 300 to 800 rpm; two Rushton turbines, 45 mm diameter) and the opening (from 0% to 100%) of the air flow valve (1 v/v/min maximum air flow). pH was controlled (pH setpoint 5.5 or 6.5) by the automatic addition of a sterile 3:1 NaOH/NaHCO3 solution 4 N, while foaming was controlled by the automatic addition of a sterile 10% (v/v) Antifoam A solution (Fluka, Sigma-Aldrich, St. Louis, MO, USA).
Two independent cultivations were carried out for each growth condition. Samples were aseptically withdrawn at 1 h intervals for the first 9 h and then at 24 h for the measurement of A650. At the end of the stationary phase (24 h), a standard curve was generated to correlate A650 values to cell dry weight (CDW; washed biomass was dried at 105 °C for 24 h). The parameters of the growth curves were estimated with the dynamic model of Baranyi and Roberts [27] using DMFit v 3.5 for Excel [28].

2.3. Oxygen Uptake

Oxygen uptake was measured in standardized (A650 = 1) exponential and stationary phase cell suspensions by using a resazurin reduction assay [29] in anaerobiosis and by monitoring (polarographic electrodes; BioXpert 2 software, Applikon) the decrease in oxygen concentration (dO2%) every 10 s for 5 min in exponential (A650 = 1) and late exponential (9 h) phases in respiration. dO2% values were transformed into µmol of O2 using Henry’s law [30] and calculating the specific oxygen uptake rate (µmol O2/min/g of biomass).

2.4. Biochemical Analyses and Enzymatic Activities

Enzymatic kits (R-Biopharm AG, Darmstadt, Germany) were used to quantify the consumption of glucose and citric acid and the production of ethanol and lactic and acetic acids in culture supernatants collected in exponential and stationary phases.
The amounts of H2O2 in the supernatants, and the activities of pyruvate oxidase (POX), NADH-dependent oxidase (NOX), NADH-dependent peroxidase (NPR), and catalase (CAT) in cell-free extracts (mechanical lysis with FastPrep-24 Instrument, MP Biomedicals, Irvine, CA, USA; 5 cycles of 60 s at speed 6.0) were measured according to Zotta et al. [31], in both exponential and stationary phases. For each biochemical and enzymatic assay, two technical replicates were carried out for each biological condition.

2.5. Statistical Analysis

All statistical and graphic analyses were performed using the free software R (https://www.r-project.org/, accessed date on 1 December 2021), version 3.4.2 [32], while package tidyverse [33] was used for tidying data and results and for generating graphs.

3. Results

3.1. Growth Kinetics and Oxygen Uptake

The kinetics of growth of Leuc. mesenteroides E30 was evaluated at different growth conditions (AN vs. RS), carbon source (glucose vs. citrate), and pH values (5.5 vs. 6.5) (Figure 1). The dynamic model of Baranyi and Roberts [27], used to estimate the growth parameters (Table 2), provided a good fit for all cultivations (R2 from 0.998 to 0.972).
The maximum specific growth rate (µmax) was affected by each examined factor (see Supplementary Table S1 for significance levels). The highest µmax values were found in the respiratory cultures compared to anaerobic ones, regardless of pH setpoint and citrate supplementation. Respiration also affected the duration of the lag phase, which was generally higher than anaerobic cultivations. Citrate boosted the growth rate only in anaerobiosis; on the contrary, its supplementation decreased the µmax in respiratory conditions. As expected, pH 5.5 impaired the growth rate compared to pH 6.5 (optimal values for cultivation).
The production of biomass (Figure 2) was affected by the type of incubation (p < 0.001) in both growth phases; interaction between the type of incubation and citrate (p < 0.05) influenced the biomass value in the exponential phase, whereas the interaction between the type of incubation and pH (p < 0.01) affected the biomass production in the stationary phase (see Supplementary Table S2 for significance levels). Respiratory conditions increased the biomass production compared with anaerobic cultivations, and the most abundant production was measured at pH 6.5 and in the absence of citrate. Cultivation at pH 5.5 had a different effect depending on the type of incubation; in fact, the suboptimal pH promoted the biomass production under anaerobic incubation but impaired it in respiration.
As expected, oxygen consumption was observed only in respiratory cultures; oxygen uptake, however, was not affected by citrate, pH, or growth phase, and no difference in the specific oxygen consumption rate was observed in cells grown under different respiratory conditions.

3.2. Consumption of Substrates and Production of Metabolites

Results related to the consumption of substrates and production of metabolites are shown in Table 3. Respiration increased the biomass yield of Leuc. mesenteroides E30 compared with anaerobiosis (see Supplementary Table S2 for significance levels). In the exponential growth phase, the strain consumed about 10–40% of glucose initially present in the culture medium, exhausting it completely only in the stationary phase. Citrate was totally consumed already in the exponential phase.
Respiration reduced lactate and ethanol formation and concurrently increased acetate production, compared to anaerobic incubation (see Supplementary Table S3 for significance levels). In all conditions, citrate slightly increased the lactate concentration, whose maximum yield was observed during anaerobic cultivation at pH 6.5. The production of acetate was dependent on the interaction of different factors (i.e., citrate supplementation, pH values, and growth phase). Ethanol, which was an abundant metabolite in anaerobic growing cells, reached the highest concentration when the pH was controlled at pH 6.5 and in the absence of citrate.

3.3. Activities of Enzymes Related to Oxygen Metabolism and Oxidative Stress

The activity of pyruvate oxidase (POX), the main enzyme of aerobic metabolism, was lower than the sensitivity limit of the assay (0.005 µkatal/mg protein) for all growth conditions. The activities of NADH-dependent oxidase (NOX), NADH-dependent peroxidase (NPR), and catalase (CAT) in exponential (A650 = 1) and stationary (24 h of incubation) phases were reported in Figure 3A–C.
Both in exponential and stationary phases, the enzymatic levels of NOX in respiratory cultures were higher (see Supplementary Table S4 for significance levels) than those measured in anaerobic ones. In respiratory conditions, the highest levels of NOX were detected in the exponential growth phase; in the stationary phase, the enzymatic activity was slightly impaired at pH 6.5 regardless of citrate supplementation. In the exponential phase, the activity of NPR was affected (p < 0.01) by citrate and by interactions between the type of incubation and citrate and between the type of incubation and pH. In the stationary phase, instead, pH (p < 0.01) and interaction between the type of incubation and pH (p < 0.05) showed a significant effect on NPR activity (Supplementary Table S4). The catalase activity of Leuc. mesenteroides E30 was mainly affected by the type of incubation (p < 0.01) in both exponential and stationary growth phases (Supplementary Table S4), and the anaerobically growing cells showed the highest activities. H2O2 was not detected in any of the tested conditions (data not shown), probably due to its degradation by catalase and NPR activities.

4. Discussion

In this study, for the first time, the effect of respiratory cultivation (i.e., oxygen and supplementation with hemin and menaquinone), pH (5.5 and 6.5), and citrate (as an alternative carbon source) was evaluated on the growth, metabolite, and enzymatic activities (POX, NOX, NPR, CAT) of Leuc. mesenteroides under controlled batch conditions.
Respiration provided several physiological and metabolic advantages for some LAB (mainly belonging to the homofermentative species Lactococcus lactis, Lactiplantibacillus plantarum, Lacticaseibacillus casei [17,18,19,20]), compared to anaerobic cultivation. In this study, we proved that respiratory cultivation improved some physiological properties also for Leuc. mesenteroides E30.
Our data, in fact, demonstrated that the strain was able to perform a respiratory metabolism exhibiting the typical traits of respiratory phenotypes (i.e., increase in growth rates and biomass accumulation and redirection of the pyruvate metabolism towards acetate production). In other heterofermentative LAB, the respiratory growth resulted in physiological benefits (improved biomass, long-term survival, oxygen uptake, changes in carbon metabolism, and the accumulation of acetate and aroma compounds).
Specifically, in Levilactobacillus spicheri and Limosilactobacillus reuteri cultivated under uncontrolled conditions, the heme-supplemented aerobiosis increased the maximum specific growth rates and the final cell density compared to anaerobic incubation [21]. Similarly, in Leuc. gelidum subsp. gasicomitatum (former Leuc. gasicomitatum), the oxygen and heme, naturally present in meat products, positively affected the growth and biomass yield with important consequences also on the metabolic profile [22]. Oxygen, in fact, represses the expression of alcohol dehydrogenase, thus shifting ethanol to acetate production. However, lactate remains the major metabolite.
Aerobic conversion of pyruvate to acetate could be mediated either by pyruvate dehydrogenase (PDH) or by POX. The latter is the key enzyme for aerobic accumulation of acetate in homofermentative LAB and some heterofermentative LAB, while the effect of PDH on aerobic condition remains unclear [19,21,34]. In Leuc. mesenteroides E30, the activity of POX was not found; this could be probably attributed to the presence of the POX gene in the strains of Leuc. mesenteroides subsp. mesenteroides and subsp. dextranicum (IMG database; NCBI database), but not in those belonging to Leuc. mesenteroides subsp. cremoris. Leuc. mesenteroides E30 was identified at the subspecies level as subsp. cremoris ([1], see label M148); therefore, the lack of POX gene is expected.
In our strain, the capability of using oxygen, detoxifying reactive oxygen species (ROS), and coping with oxidative conditions may be due to the presence of other enzymes, such as NOX and NPR. In LAB, NOX and NPR activities contribute to oxidative stress resistance and to intracellular redox balance because these enzymes use NADH to regenerate NAD+. Specifically, NOX reacts with oxygen to produce either water or H2O2. The latter compound is reduced to water by NPR activity. In Leuc. mesenteroides E30, NOX and NPR were found in both anaerobic and respiratory cultures; however, the respiratory growth of Leuc. mesenteroides E30 remarkably increased NOX activity.
Similarly, an oxygen-induced function of NOX was detected in Limosilactobacillus panis. Furthermore, oxygen availability did not affect the gene expression but increased the enzymatic activity with remarkable effects on energy yield and metabolism [35]. In Fructilactobacillus sanfranciscensis, the oxygen availability did not affect NOX activity; however, a NOX-negative strain showed a fructose-dependent growth response under aerobic conditions [36].
The fact that NPR activity of Leuc. mesenteroides E30 was not affected by the respiratory condition could suggest the presence of only water-forming NOX, as previously described in Leuc. mesenteroides subsp. mesenteroides by Sakamoto et al. [37]. The absence of H2O2 in all conditions could prove the effectiveness of the NOX/NPR coupled system in the oxidative stress protection of Leuc. mesenteroides E30.
The degradation of H2O2 in Leuc. mesenteroides E30 could be due to catalase activity, even if its production was not induced only by heme supplementation, suggesting possible activity of a manganese-dependent isoform. Heme- and Mn-dependent catalase activities have been previously described in some heterofermentative LAB [21,24,38]. Moreover, some annotated genomes of Leuc. mesenteroides exclusively harbor gene encoding for Mn catalase (IMG database).
In this study, the effect of glucose and citrate, as carbon sources, was also evaluated. Specifically, citrate, used in co-metabolism to glucose, affected the biomass production and metabolic production of Leuc. mesenteroides E30, depending on the type of incubation. Citrate had a positive impact on the growth performances of anaerobically grown cells, but it impaired those of respiratory cells. In some LAB, citrate may be converted into succinate via the reductive tricarboxylic acid cycle or supply the pool of pyruvate by citrate lyase-oxaloacetate decarboxylase activities [19]. In anaerobic conditions, pyruvate is reduced into lactate by lactate dehydrogenase activity; in aerobic conditions, the accumulation of pyruvate resulting also from citrate metabolism, could interfere with the growth of respiratory cultures. Citrate, instead, had no effect on acetic acid and ethanol yields.
Heterofermentative LAB may use several carbon sources, i.e., maltose, fructose, and sucrose [12], including citrate, with a direct effect on their physiology [36]. Özcan et al. [39], through an in silico study, observed a stimulation of Leuc. mesenteroides growth during the co-utilization of citrate and glucose, in agreement with other experimental data [40,41]. The presence of citrate, contributing to the pyruvate pool through oxaloacetate decarboxylase, caused an increased production of lactate and, consequently, an increased oxidation of NADH. Therefore, because the requirement of ethanol production for the re-oxidation of NADH decreases, citrate contributes to the acetate pool by acetate kinase activity with additional ATP yield.
To date, the potential of respiratory phenotypes has been industrially exploited only by the CHR Hansen Company (Denmark) to improve the production of a lactococcal starter culture (direct-to-vat F-DVS pHageControl™ R-604). The benefits of the respiratory metabolism, –, are different. Furthermore, in previous studies, Reale et al. [42,43] demonstrated that some strains of Lacticaseibacillus casei cultivated in respiratory conditions affected the biochemical and organoleptic features of some foods (i.e., reduction in oxidative processes, increase in secondary proteolysis and changes in the aroma compound in Cheddar-type cheeses and in sourdoughs).
Leuc. mesenteroides plays a positive role in the production of several fermented foods. The capability to produce gas, exopolysaccharides, and several aroma compounds, in fact, may contribute to the development of a desired texture and flavour in different products (e.g., fermented vegetables, bakery products, and some cheese varieties). Therefore, the understanding of metabolic pathways and the exploitation of strains with improved phenotypic features and stress resistance, also for this species, may be of practical relevance for different food-related applications.

5. Conclusions

This study demonstrated that respiratory cultivation improved the growth rate and biomass production, and affected the metabolic pattern, energetic state, and activity of antioxidant enzymes of Leuc. mesenteroides E30, compared to anaerobic conditions.
Therefore, the respiratory phenotype of Leuc. mesenteroides E30 could be exploited as a starter and/or an adjunct culture to improve the organoleptic and nutritional properties of several foods, especially vegetable- and cereal-based products.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/foods11040535/s1, Table S1. Analysis of variance on the effect of atmosphere of incubation, pH, citrate and their interactions on specific growth rate (μmax), Table S2: Analysis of variance on the effect of atmosphere of incubation, pH, citrate and their interactions; Table S3. Analysis of variance on the effect of atmosphere of incubation, pH, citrate and their interactions; Table S4. Analysis of variance on the effect of atmosphere of incubation, pH, citrate and their interactions.

Author Contributions

Conceptualization, E.P., A.R. and T.Z.; formal analysis, L.V.S. and M.G.; data curation, E.P.; investigation, L.V.S. and M.G.; writing—original draft preparation, E.P.; writing—review and editing, A.R. and T.Z.; supervision, E.P., A.R. and T.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ricciardi, A.; Storti, L.V.; Zotta, T.; Felis, G.E.; Parente, E. Analysis of rpoB polymorphism and PCR-based approaches for the identification of Leuconostoc mesenteroides at the species and subspecies level. Int. J. Food Microbiol. 2020, 318, 108474. [Google Scholar] [CrossRef]
  2. Botta, C.; Langerholc, T.; Cencič, A.; Cocolin, L. In vitro selection and characterization of new probiotic candidates from table olive microbiota. PLoS ONE 2014, 9, e94457. [Google Scholar] [CrossRef] [PubMed]
  3. Corona, O.; Alfonzo, A.; Ventimiglia, G.; Nasca, A.; Francesca, N.; Martorana, A.; Moschetti, G.; Settanni, L. Industrial application of selected lactic acid bacteria isolated from local semolinas for typical sourdough bread production. Food Microbiol. 2016, 59, 43–56. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Moon, S.H.; Kim, C.R.; Chang, H.C. Heterofermentative lactic acid bacteria as a starter culture to control kimchi fermentation. LWT-Food Sci. Technol. 2018, 88, 181–188. [Google Scholar] [CrossRef]
  5. Pogačić, T.; Maillard, M.B.; Leclerc, A.; Hervé, C.; Chuat, V.; Valence, F.; Thierry, A. Lactobacillus and Leuconostoc volatilomes in cheese conditions. Appl. Microbiol. Biotechnol. 2016, 100, 2335–2346. [Google Scholar] [CrossRef]
  6. Xu, Y.; Coda, R.; Shi, Q.; Tuomainen, P.; Katina, K.; Tenkanen, M. Exopolysaccharides production during the fermentation of soybean and fava bean flours by Leuconostoc mesenteroides DSM 20343. J. Agric. Food Chem. 2017, 65, 2805–2815. [Google Scholar] [CrossRef]
  7. D’Angelo, L.; Cicotello, J.; Zago, M.; Guglielmotti, D.; Quiberoni, A.; Suárez, V. Leuconostoc strains isolated from dairy products: Response against food stress conditions. Food Microbiol. 2017, 66, 28–39. [Google Scholar] [CrossRef]
  8. Comi, G.; Iacumin, L. Identification and process origin of bacteria responsible for cavities and volatile off-flavour compounds in artisan cooked ham. J. Food Sci. Technol. 2012, 47, 114–121. [Google Scholar] [CrossRef]
  9. Comi, G.; Andyanto, D.; Manzano, M.; Iacumin, L. Lactococcus lactis and Lactobacillus sakei as bio-protective culture to eliminate Leuconostoc mesenteroides spoilage and improve the shelf life and sensorial characteristics of commercial cooked bacon. Food Microbiol. 2016, 58, 16–22. [Google Scholar] [CrossRef]
  10. Pothakos, V.; Nyambi, C.; Zhang, B.Y.; Papastergiadis, A.; De Meulenaer, B.; Devlieghere, F. Spoilage potential of psychrotrophic lactic acid bacteria (LAB) species: Leuconostoc gelidum subsp. gasicomitatum and Lactococcus piscium on sweet bell pepper (SBP) simulation medium under different gas compositions. Int. J. Food Microbiol. 2014, 178, 120–129. [Google Scholar] [CrossRef]
  11. Pothakos, V.; Aulia, Y.A.; van der Linden, I.; Uyttendaele, M.; Devlieghere, F. Exploring the strain-specific attachment of Leuconostoc gelidum subsp. gasicomitatum on food contact surfaces. Int. J. Food Microbiol. 2015, 199, 41–46. [Google Scholar] [CrossRef] [PubMed]
  12. Gänzle, M.G. Lactic metabolism revisited: Metabolism of lactic acid bacteria in food fermentations and food spoilage. Curr. Opin. Food Sci. 2015, 2, 106–117. [Google Scholar] [CrossRef]
  13. Smid, E.J.; Kleerebezem, M. Production of aroma compounds in lactic fermentations. Annu. Rev. Food Sci. Technol. 2014, 5, 313–326. [Google Scholar] [CrossRef] [PubMed]
  14. Zaunmüller, T.; Eichert, M.; Richter, H.; Unden, G. Variations in the energy metabolism of biotechnologically relevant heterofermentative lactic acid bacteria during growth on sugars and organic acids. Appl. Microbiol. Biotechnol. 2006, 72, 421–429. [Google Scholar] [CrossRef]
  15. Plihon, F.; Taillandier, P.; Strehaiano, P. Oxygen effect on batch cultures of Leuconostoc mesenteroides: Relationship between oxygen uptake, growth and end-products. Appl. Microbiol. Biotechnol. 1995, 43, 117–122. [Google Scholar] [CrossRef]
  16. Plihon, F.; Taillandier, P.; Strehaiano, P. Oxygen effect on lactose catabolism by a Leuconostoc mesenteroides strain: Modeling of general O2-dependent stoichiometry. Biotechnol. Bioeng. 1996, 49, 63–69. [Google Scholar] [CrossRef]
  17. Lechardeur, D.; Cesselin, B.; Fernandez, A.; Lamberet, G.; Garrigues, C.; Pedersen, M.; Gaudu, P.; Gruss, A. Using heme as an energy boost for lactic acid bacteria. Curr. Opin. Biotech. 2011, 22, 143–149. [Google Scholar] [CrossRef]
  18. Pedersen, M.B.; Gaudu, P.; Lechardeur, D.; Petit, M.A.; Gruss, A. Aerobic respiration metabolism in lactic acid bacteria and uses in biotechnology. Annu. Rev. Food Sci. Technol. 2012, 3, 37–58. [Google Scholar] [CrossRef]
  19. Zotta, T.; Ricciardi, A.; Parente, E. Aerobic metabolism in the genus Lactobacillus: Impact on stress response and potential applications in the food industry. J. Appl. Microbiol. 2017, 122, 857–869. [Google Scholar] [CrossRef] [Green Version]
  20. Brooijmans, R.J.W.; Smit, B.; Santos, F.; van Riel, J.; de Vos, W.M.; Hugenholtz, J. Heme and menaquinone induced electron transport in lactic acid bacteria. Microb. Cell Fact. 2009, 8, 28. [Google Scholar] [CrossRef] [Green Version]
  21. Ianniello, R.G.; Zheng, J.; Zotta, T.; Ricciardi, A.; Gänzle, M.G. Biochemical analysis of respiratory metabolism in the heterofermentative Lactobacillus spicheri and Lactobacillus reuteri. J. Appl. Microbiol. 2015, 119, 763–775. [Google Scholar] [CrossRef] [PubMed]
  22. Jääskeläinen, E.; Johansson, P.; Kostiainen, O.; Nieminen, T.; Schmidt, G.; Somervuo, P.; Mohsina, M.; Vanninen, P.; Auvinen, P.; Björkroth, J. Significance of heme-based respiration in meat spoilage caused by Leuconostoc gasicomitatum. Appl. Environ. Microbiol. 2013, 79, 1078–1085. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Rahkila, R.; De Bruyne, K.; Johansson, P.; Vandamme, P.; Björkroth, J. Reclassification of Leuconostoc gasicomitatum as Leuconostoc gelidum subsp. gasicomitatum comb. nov. description of Leuconostoc gelidum subsp. aenigmaticum subsp. nov. designation of Leuconostoc gelidum subsp. gelidum subsp. nov. and emended description of Leuconostoc gelidum. Int. J. Syst. Evol. Microbiol. 2014, 64, 1290–1295. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Zotta, T.; Ricciardi, A.; Ianniello, R.G.; Storti, L.V.; Glibota, N.A.; Parente, E. Aerobic and respirative growth of heterofermentative lactic acid bacteria: A screening study. Food Microbiol. 2018, 76, 117–127. [Google Scholar] [CrossRef] [PubMed]
  25. Zotta, T.; Ricciardi, A.; Guidone, A.; Sacco, M.; Muscariello, L.; Mazzeo, M.F.; Cacace, G.; Parente, E. Inactivation of ccpA and aeration affect growth, metabolite production and stress tolerance in Lactobacillus plantarum WCFS1. Int. J. Food Microbiol. 2012, 155, 51–59. [Google Scholar] [CrossRef] [PubMed]
  26. Ianniello, R.G.; Zotta, T.; Matera, A.; Genovese, F.; Parente, E.; Ricciardi, A. Investigation of factors affecting aerobic and respiratory growth in the oxygen-tolerant strain Lactobacillus casei N87. PLoS ONE 2016, 11, e0164065. [Google Scholar] [CrossRef]
  27. Baranyi, J.; Roberts, T.A. A dynamic approach to predicting bacterial growth in food. Int. J. Food Microbiol. 1994, 23, 277–294. [Google Scholar] [CrossRef]
  28. Baranyi, J. DMFit Manual Version 3.5 Institute of Food Research; Norwich Research Park: Norwich, UK, 2015. [Google Scholar]
  29. Ricciardi, A.; Ianniello, R.G.; Tramutola, A.; Parente, E.; Zotta, T. Rapid detection assay for oxygen consumption in the Lactobacillus casei group. Ann. Microbiol. 2014, 64, 1861–1864. [Google Scholar] [CrossRef]
  30. Sander, R. Compilation of Henry’s Law constants (Version 4.0) for water as solvent. Atmos. Chem. Phys. 2015, 15, 4399–4981. [Google Scholar] [CrossRef] [Green Version]
  31. Zotta, T.; Ricciardi, A.; Ianniello, R.G.; Parente, E.; Reale, A.; Rossi, F.; Iacumin, L.; Comi, G.; Coppola, R. Assessment of aerobic and respiratory growth in the Lactobacillus casei group. PLoS ONE 2014, 9, e99189. [Google Scholar] [CrossRef] [Green Version]
  32. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2021; Available online: https://www.R-project.org/ (accessed on 1 December 2021).
  33. Wickham, H. Tidyverse: Easily Install and Load ‘Tidyverse’ Packages. R Package Version 1.1.1. 2017. Available online: https://CRAN.R-project.org/package=tidyverse (accessed on 1 December 2021).
  34. Guo, T.; Zhang, L.; Xin, Y.; Xu, Z.; He, H.; Kong, J. Oxygen-inducible conversion of lactate to acetate in heterofermentative Lactobacillus brevis ATCC 367. Appl. Environ. Microbiol. 2017, 83, 21. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Kang, T.S.; Korber, D.R.; Tanaka, T. Influence of oxygen on NADH recycling and oxidative stress resistance systems in Lactobacillus panis PM1. AMB Express 2013, 3, 10–19. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Jänsch, A.; Freiding, S.; Behr, J.; Vogel, R.F. Contribution of the NADH-oxidase (Nox) to the aerobic life of Lactobacillus sanfranciscensis DSM20451T. Food Microbiol. 2011, 28, 29–37. [Google Scholar] [CrossRef] [PubMed]
  37. Sakamoto, M.; Uchimura, T.; Komagata, K. Comparison of H2O-forming NADH oxidase from Leuconostoc mesenteroides subsp. mesenteroides NRIC 1541T and H2O2-forming NADH oxidase from Sporolactobacillus inulinus NRIC 1133T. J. Biosci. Bioeng. 1996, 82, 531–537. [Google Scholar] [CrossRef]
  38. Whittenbury, R. Hydrogen peroxide formation and catalase activity in the lactic acid bacteria. J. Gen. Microbiol. 1964, 35, 13–26. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Özcan, E.; Selvi, S.S.; Nikerel, E.; Teusink, B.; Öner, E.T.; Çakır, T. A genome-scale metabolic network of the aroma bacterium Leuconostoc mesenteroides subsp. cremoris. Appl. Microbiol. Biotechnol. 2019, 103, 3153–3165. [Google Scholar] [CrossRef]
  40. Marty-Teysset, C.; Lolkema, J.S.; Schmitt, P.; Diviès, C.; Konings, W.N. The citrate metabolic pathway in Leuconostoc mesenteroides: Expression, amino acid synthesis, and a-ketocarboxylate transport. J. Bacteriol. 1996, 178, 6209–6215. [Google Scholar] [CrossRef] [Green Version]
  41. Schmitt, P.; Divies, C.; Cardona, R. Origin of end-products from the co-metabolism of glucose and citrate by Leuconostoc mesenteroides subsp. cremoris. Appl. Microbiol. Biotechnol. 1992, 36, 679–683. [Google Scholar] [CrossRef]
  42. Reale, A.; Ianniello, R.G.; Ciocia, F.; Di Renzo, T.; Boscaino, B.; Ricciardi, A.; Coppola, R.; Parente, E.; Zotta, T.; McSweeney, P.L.H. Effect of respirative and catalase-positive Lactobacillus casei adjuncts on the production and quality of Cheddar-type cheese. Int. Dairy J. 2016, 63, 78–87. [Google Scholar] [CrossRef]
  43. Reale, A.; Di Renzo, T.; Zotta, T.; Preziuso, M.; Boscaino, F.; Ianniello, R.; Storti, L.V.; Tremonte, P.; Coppola, R. Effect of respirative cultures of Lactobacillus casei on model sourdough fermentation. LWT-Food Sci. Technol. 2016, 73, 622–629. [Google Scholar] [CrossRef]
Foods 11 00535 g001 550
Figure 1. Growth kinetics of Leuconostoc mesenteroides E30. The conditions of each fermentation (F1 to F8) are detailed in Table 1. X-axis, time (hours) of incubation; Y-axis, growth (lnX/X0), dissolved oxygen (dO2%), and NaOH:NaHCO3 (meq/L) consumption; black symbols, biomass production; blue line, kinetics of growth (ln X/X0) estimated with the model of Baranyi and Roberts [27] where X and X0 are, respectively, the biomass measured at each withdrawal time and at start of fermentation; green line, dO2% measured through polarographic electrode; red line, milliequivalents/L (meq/L) of NaOH:NaHCO3 consumed during fermentation. Code of fermentation name: AN, anaerobiosis; RS, respiration; 0, absence of sodium citrate in mWMB; 1, supplementation of mWMB with 5 g/L of sodium citrate; 5.5, pH controlled at 5.5; 6.5, pH controlled at 6.5.
Figure 1. Growth kinetics of Leuconostoc mesenteroides E30. The conditions of each fermentation (F1 to F8) are detailed in Table 1. X-axis, time (hours) of incubation; Y-axis, growth (lnX/X0), dissolved oxygen (dO2%), and NaOH:NaHCO3 (meq/L) consumption; black symbols, biomass production; blue line, kinetics of growth (ln X/X0) estimated with the model of Baranyi and Roberts [27] where X and X0 are, respectively, the biomass measured at each withdrawal time and at start of fermentation; green line, dO2% measured through polarographic electrode; red line, milliequivalents/L (meq/L) of NaOH:NaHCO3 consumed during fermentation. Code of fermentation name: AN, anaerobiosis; RS, respiration; 0, absence of sodium citrate in mWMB; 1, supplementation of mWMB with 5 g/L of sodium citrate; 5.5, pH controlled at 5.5; 6.5, pH controlled at 6.5.
Foods 11 00535 g001
Foods 11 00535 g002 550
Figure 2. RS, respiratory growth (with 15% dissolved oxygen (dO2%), supplementation with 1.25 μg/mL of hemin and 1 μg/mL of menaquinone); AN, anaerobic growth (with nitrogen flow at 0.02 L/min). Conditions of each fermentation (F1 to F8) are detailed in Table 1. Code of fermentation name: AN, anaerobiosis; RS, respiration; 0, absence of sodium citrate in mWMB; 1, supplementation of sodium citrate (5 g/L) in mWMB; 5.5, controlled pH at 5.5; 6.5, pH controlled at 6.5. Growth phase: E, exponential phase; S, stationary phase. Mean values and standard errors are reported.
Figure 2. RS, respiratory growth (with 15% dissolved oxygen (dO2%), supplementation with 1.25 μg/mL of hemin and 1 μg/mL of menaquinone); AN, anaerobic growth (with nitrogen flow at 0.02 L/min). Conditions of each fermentation (F1 to F8) are detailed in Table 1. Code of fermentation name: AN, anaerobiosis; RS, respiration; 0, absence of sodium citrate in mWMB; 1, supplementation of sodium citrate (5 g/L) in mWMB; 5.5, controlled pH at 5.5; 6.5, pH controlled at 6.5. Growth phase: E, exponential phase; S, stationary phase. Mean values and standard errors are reported.
Foods 11 00535 g002
Foods 11 00535 g003 550
Figure 3. Enzymatic activity of NADH oxidase (NOX, (A)), NADH peroxidase (NPR, (B)), and catalase (CAT, (C)), respectively, expressed in µkatal/mg protein. Growth conditions: RS, respiratory growth (with 15% dissolved oxygen (dO2%), supplementation with 1.25 μg/mL of hemin and 1 μg/mL of menaquinone); AN, anaerobic growth (with nitrogen flow at 0.02 L/min). Conditions of each fermentation (F1 to F8) are detailed in Table 1. Code of fermentation name: AN, anaerobiosis; RS, respiration; 0, absence of sodium citrate in mWMB; 1, supplementation of sodium citrate (5 g/L) in mWMB; 5.5, controlled pH at 5.5; 6.5, pH controlled at 6.5. Growth phase: E, exponential phase; S, stationary phase. Mean values and standard errors are reported.
Figure 3. Enzymatic activity of NADH oxidase (NOX, (A)), NADH peroxidase (NPR, (B)), and catalase (CAT, (C)), respectively, expressed in µkatal/mg protein. Growth conditions: RS, respiratory growth (with 15% dissolved oxygen (dO2%), supplementation with 1.25 μg/mL of hemin and 1 μg/mL of menaquinone); AN, anaerobic growth (with nitrogen flow at 0.02 L/min). Conditions of each fermentation (F1 to F8) are detailed in Table 1. Code of fermentation name: AN, anaerobiosis; RS, respiration; 0, absence of sodium citrate in mWMB; 1, supplementation of sodium citrate (5 g/L) in mWMB; 5.5, controlled pH at 5.5; 6.5, pH controlled at 6.5. Growth phase: E, exponential phase; S, stationary phase. Mean values and standard errors are reported.
Foods 11 00535 g003
Table
Table 1. Experimental design used to evaluate the effect of the type of incubation, carbon source, and pH on growth, biomass, metabolites production, and enzymatic activities.
Table 1. Experimental design used to evaluate the effect of the type of incubation, carbon source, and pH on growth, biomass, metabolites production, and enzymatic activities.
FermentationIncubation aCitrate bpH
F1AN06.5
F2AN16.5
F3AN05.5
F4AN15.5
F5RS06.5
F6RS16.5
F7RS05.5
F8RS15.5
a AN, anaerobic growth with nitrogen flow at 0.02 L/min; RS, respiratory growth with 15% dissolved oxygen (dO2%), supplementation with 1.25 μg/mL of hemin and 1 μg/mL of menaquinone. b 0, absence of sodium citrate in mWMB; 1, supplementation of mWMB with 5 g/L of sodium citrate.
Table
Table 2. Parameters of growth kinetics of Leuconostoc mesenteroides E30.
Table 2. Parameters of growth kinetics of Leuconostoc mesenteroides E30.
Growth Conditionµmax (h−1)lag (h)R2 of Fit
F1_AN_0_6.5_a0.46 ± 0.020.00 ± 0.000.972
F1_AN_0_6.5_b0.46 ± 0.011.38 ± 0.210.996
F2_AN_1_6.5_a0.64 ± 0.011.05 ± 0.110.998
F2_AN_1_6.5_b0.57 ± 0.011.13 ± 0.140.997
F3_AN_0_5.5_a0.40 ± 0.021.92 ± 0.250.995
F3_AN_0_5.5_b0.42 ± 0.023.09 ± 0.220.997
F4_AN_1_5.5_a0.50 ± 0.010.00 ± 0.000.997
F4_AN_1_5.5_b0.47 ± 0.011.54 ± 0.190.997
F5_RS_0_6.5_a0.77 ± 0.032.17 ± 0.190.997
F5_RS_0_6.5_b0.82 ± 0.060.36 ± 3.570.992
F6_RS_1_6.5_a0.72 ± 0.074.58 ± 0.320.988
F6_RS_1_6.5_b0.71 ± 0.032.58 ± 0.210.995
F7_RS_0_5.5_a0.63 ± 0.032.11 ± 0.180.996
F7_RS_0_5.5_b0.67 ± 0.032.41 ± 0.200.996
F8_RS_1_5.5_a0.55 ± 0.063.24 ± 0.500.974
F8_RS_1_5.5_b0.49 ± 0.033.00 ± 0.350.988
Growth condition: RS, respiratory growth (with 15% dissolved oxygen (dO2%), supplementation with 1.25 μg/mL of hemin and 1 μg/mL of menaquinone); AN, anaerobic growth (with nitrogen flow at 0.02 L/min). Conditions of each fermentation (F1 to F8) are detailed in Table 1. Code of fermentation name: AN, anaerobiosis; RS, respiration; 0, absence of sodium citrate in mWMB; 1, supplementation of mWMB with 5 g/L of sodium citrate; 5.5, controlled pH at 5.5; 6.5, pH controlled at 6.5; µmax, maximum specific growth rate; lag (h), duration of lag phase; R2 of fit quantified the goodness of modeling.
Table
Table 3. The consumption of glucose and citrate and the production of metabolites of Leuconostoc mesenteroides E30.
Table 3. The consumption of glucose and citrate and the production of metabolites of Leuconostoc mesenteroides E30.
Growth Condition aGrowth Phase b Biomass Yield cGlucose
Consumed d		Lactic Acid Yield eAcetic ACID Yield fEthanol Yield g
F1_AN_0_6.5E0.074 ± 0.03234.8 ± 3.3560.595 ± 0.2480.038 ± 0.0100.450 ± 0.0186
S0.072 ± 0.009131.93 ± 4.5120.979 ± 0.0430.019 ± 0.0210.682 ± 0.046
F2_AN_1_6.5E0.102 ± 0.07027.16 ± 21.2531.707 ± 1.3391.205 ± 0.9780.158 ± 0.065
S0.076 ± 0.006116.81 ± 5.7611.209 ± 0.0700.196 ± 0.0120.666 ± 0.108
F3_AN_0_5.5E0.074 ± 0.00532.43 ± 0.7460.763 ± 0.0030.036 ± 0.0240.606 ± 0.074
S0.082 ± 0.009130.78 ± 1.3610.843 ± 0.0010.010 ± 0.0030.653 ± 0.087
F4_AN_1_5.5E0.140 ± 0.00617.01 ± 2.4241.706 ± 0.2781.530 ± 0.1760.247 ± 0.185
S0.082 ± 0.006122.02 ± 2.8341.186 ± 0.2870.177 ± 0.0380.590 ± 0.148
F5_RS_0_6.5E0.232 ± 0.03613.71 ± 4.8470.802 ± 0.1860.726 ± 0.1810.090 ± 0.098
S0.135 ± 0.051126.43 ± 1.5100.613 ± 0.0110.634 ± 0.0080.130 ± 0.146
F6_RS_1_6.5E0.145 ± 0.02224.78 ± 1.4910.924 ± 0.1640.000 ± 0.0000.000 ± 0.000
S0.128 ± 0.004124.88 ± 1.9200.783 ± 0.0690.611 ± 0.0280.105 ± 0.059
F7_RS_0_5.5E0.077 ± 0.02743.77 ± 14.1690.244 ± 0.0880.147 ± 0.0730.000 ± 0.000
S0.113 ± 0.005126.91 ± 4.7540.468 ± 0.0700.499 ± 0.0220.089 ± 0.097
F8_RS_1_5.5E0.455 ± 0.17411.21 ± 0.5591.877 ± 0.0531.649 ± 0.0940.000 ± 0.000
S0.102 ± 0.033122.41 ± 4.6610.940 ± 0.0090.469 ± 0.0250.100 ± 0.008
a Growth conditions: RS, respiratory growth (with 15% dissolved oxygen (dO2%), supplementation with 1.25 μg/mL of hemin and 1 μg/mL of menaquinone); AN, anaerobic growth (with nitrogen flow at 0.02 L/min). The conditions of each fermentation (F1 to F8) are detailed in Table 1. Code of fermentation name: AN, anaerobiosis; RS, respiration; 0, absence of sodium citrate in mWMB; 1, supplementation of mWMB with 5 g/L of sodium citrate; 5.5, controlled pH at 5.5; 6.5, pH controlled at 6.5. b Growth phase: E, exponential phase; S, stationary phase. Mean values of two biological replicates ± standards errors are reported. c Biomass yield: calculated as (X-X0)/(S0-S), i.e., biomass production (X-X0), g/L, relative to consumed glucose (S0-S), g/L. d Glucose consumed: consumed glucose (S0-S), mM. e Lactate yield: production of DL-lactic acid (P-P0), g/L, relative to consumed glucose (S0-S), g/L. f Acetate yield: production of acetic acid (A-A0), g/L, relative to consumed glucose (S0-S), g/L. g Ethanol yield: production of acetic acid (E-E0), g/L, relative to consumed glucose (S0-S), g/L.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ricciardi, A.; Storti, L.V.; Giavalisco, M.; Parente, E.; Zotta, T. The Effect of Respiration, pH, and Citrate Co-Metabolism on the Growth, Metabolite Production and Enzymatic Activities of Leuconostoc mesenteroides subsp. cremoris E30. Foods 2022, 11, 535. https://doi.org/10.3390/foods11040535

AMA Style

Ricciardi A, Storti LV, Giavalisco M, Parente E, Zotta T. The Effect of Respiration, pH, and Citrate Co-Metabolism on the Growth, Metabolite Production and Enzymatic Activities of Leuconostoc mesenteroides subsp. cremoris E30. Foods. 2022; 11(4):535. https://doi.org/10.3390/foods11040535

Chicago/Turabian Style

Ricciardi, Annamaria, Livia Vanessa Storti, Marilisa Giavalisco, Eugenio Parente, and Teresa Zotta. 2022. "The Effect of Respiration, pH, and Citrate Co-Metabolism on the Growth, Metabolite Production and Enzymatic Activities of Leuconostoc mesenteroides subsp. cremoris E30" Foods 11, no. 4: 535. https://doi.org/10.3390/foods11040535

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop