Holistic integration of omics data reveals the drivers that shape the ecology of microbial meat spoilage scenarios

Background

The use of omics data for monitoring the microbial flow of fresh meat products along a production line and the development of spoilage prediction tools from these data is a promising but challenging task. In this context, we produced a large multivariate dataset (over 600 samples) obtained on the production lines of two similar types of fresh meat products (poultry and raw pork sausages). We describe a full analysis of this dataset in order to decipher how the spoilage microbial ecology of these two similar products may be shaped differently depending on production parameter characteristics.

Methods

Our strategy involved a holistic approach to integrate unsupervised and supervised statistical methods on multivariate data (OTU-based microbial diversity; metabolomic data of volatile organic compounds; sensory measurements; growth parameters), and a specific selection of potential uncontrolled (initial microbiota composition) or controlled (packaging type; lactate concentration) drivers.

Results

Our results demonstrate that the initial microbiota, which is shown to be very different between poultry and pork sausages, has a major impact on the spoilage scenarios and on the effect that a downstream parameter such as packaging type has on the overall evolution of the microbial community. Depending on the process, we also show that specific actions on the pork meat (such as deboning and defatting) elicit specific food spoilers such as Dellaglioa algida, which becomes dominant during storage. Finally, ecological network reconstruction allowed us to map six different metabolic pathways involved in the production of volatile organic compounds involved in spoilage. We were able connect them to the different bacterial actors and to the influence of packaging type in an overall view. For instance, our results demonstrate a new role of Vibrionaceae in isopropanol production, and of Latilactobacillus fuchuensis and Lactococcus piscium in methanethiol/disylphide production. We also highlight a possible commensal behavior between Leuconostoc carnosum and Latilactobacillus curvatus around 2,3-butanediol metabolism.

Conclusion

We conclude that our holistic approach combined with large-scale multi-omic data was a powerful strategy to prioritize the role of production parameters, already known in the literature, that shape the evolution and/or the implementation of different meat spoilage scenarios.

Evolution of fungal community associated with ready-to-eat pineapple during storage under different temperature conditions

The international market of fresh-cut products has witnessed dramatic growth in recent years, stimulated by consumer's demand for healthy, nutritious and convenient foods. One of the main challenging issues for the quality and safety of these products is the potential microbial spoilage that can significantly reduce their shelf-life. The complete identification of fresh-cut product microbiota together with the evaluation of environmental factors impact on microbial composition is of primary importance. We therefore assessed the fungal communities associated with the spoilage of ready-to-eat (RTE) pineapple using a metagenetic amplicon sequencing approach, based on the ITS2 region. Our results revealed a significant variability on fungal species composition between the different batches of RTE pineapple. The initial microbiota composition was the main influencing factor and determined the progress of spoilage. Temperature and storage time were the secondary factors influencing spoilage and their impact was depending on the initial prevalent fungal species, which showed different responses to the various modifications. Our results strongly suggest that further large-scale sampling of RTE pineapple production should be conducted in order to assess the full biodiversity range of fungal community involved in the spoilage process and for unravelling the impact of important environmental factors shaping the initial microbiota.

Heme Uptake in Lactobacillus sakei Evidenced by a New Energy Coupling Factor (ECF)-Like Transport System

Lactobacillus sakei is a nonpathogenic lactic acid bacterium and a natural inhabitant of meat ecosystems. Although red meat is a heme-rich environment, L. sakei does not need iron or heme for growth, although it possesses a heme-dependent catalase. Iron incorporation into L. sakei from myoglobin and hemoglobin was previously shown by microscopy and the L. sakei genome reveals the complete equipment for iron and heme transport. Here, we report the characterization of a five-gene cluster (from lsa1836 to lsa1840 [lsa1836-1840]) encoding a putative metal iron ABC transporter. Interestingly, this cluster, together with a heme-dependent catalase gene, is also conserved in other species from the meat ecosystem. Our bioinformatic analyses revealed that the locus might correspond to a complete machinery of an energy coupling factor (ECF) transport system. We quantified in vitro the intracellular heme in the wild type (WT) and in our Δlsa1836-1840 deletion mutant using an intracellular heme sensor and inductively coupled plasma mass spectrometry for quantifying incorporated 57Fe heme. We showed that in the WT L. sakei, heme accumulation occurs rapidly and massively in the presence of hemin, while the deletion mutant was impaired in heme uptake; this ability was restored by in trans complementation. Our results establish the main role of the L. sakei Lsa1836-1840 ECF-like system in heme uptake. Therefore, this research outcome sheds new light on other possible functions of ECF-like systems.IMPORTANCELactobacillus sakei is a nonpathogenic bacterial species exhibiting high fitness in heme-rich environments such as meat products, although it does not need iron or heme for growth. Heme capture and utilization capacities are often associated with pathogenic species and are considered virulence-associated factors in the infected hosts. For these reasons, iron acquisition systems have been deeply studied in such species, while for nonpathogenic bacteria the information is scarce. Genomic data revealed that several putative iron transporters are present in the genome of the lactic acid bacterium L. sakei In this study, we demonstrate that one of them is an ECF-like ABC transporter with a functional role in heme transport. Such evidence has not yet been brought for an ECF; therefore, our study reveals a new class of heme transport system.

Large microbiota survey reveals how the microbial ecology of cooked ham is shaped by different processing steps

Cooked ham production involves numerous steps shaping the microbial communities of the final product, with consequences on spoilage metabolites production. To identify the main factors driving the ecology of ham and its spoilage, we designed a study encompassing five variables related to ham production: type of storage during meat transportation, churning speed, drain-off time, slicing line and O₂ packaging permeability. About 200 samples from the same facility were obtained and characterized with respect to i) their microbiota based on gyrB amplicon sequencing ii) their production of spoilage-related metabolites based on E-Nose analysis and enzymatic assays. The slicing was the most critical step, shaping two general types of microbiota according to the slicing line: one dominated by Carnobacterium divergens and another one dominated by Leuconostoc carnosum and Serratia proteamaculans. Regarding metabolites production, L. carnosum was associated to d-lactic acid, ethanol and acetic acid production, whereas Serratia proteamaculans was associated to acetic acid production. This last species prevailed with highly O₂-permeable packaging. Within a given slicing line, campaign-based variations were observed, with Lactobacillus sakei, Leuconostoc mesenteroides and Carnobacterium maltaromaticum prevalent in summer. L. sakei was associated with l-lactic acid production and C. maltaromaticum with formic and acetic acid productions.

Large-scale multivariate dataset on the characterization of microbiota diversity, microbial growth dynamics, metabolic spoilage volatilome and sensorial profiles of two industrially produced meat products subjected to changes in lactate concentration and packaging atmosphere

Data in this article provide detailed information on the diversity of bacterial communities present on 576 samples of raw pork or poultry sausages produced industrially in 2017. Bacterial growth dynamics and diversity were monitored throughout the refrigerated storage period to estimate the impact of packaging atmosphere and the use of potassium lactate as chemical preservative. The data include several types of analysis aiming at providing a comprehensive microbial ecology of spoilage during storage and how the process parameters do influence this phenomenon. The analysis includes: the gas content in packaging, pH, chromametric measurements, plate counts (total mesophilic aerobic flora and lactic acid bacteria), sensorial properties of the products, meta-metabolomic quantification of volatile organic compounds and bacterial community metagenetic analysis. Bacterial diversity was monitored using two types of amplicon sequencing (16S rRNA and GyrB encoding genes) at different time points for the different conditions (576 samples for gyrB and 436 samples for 16S rDNA). Sequencing data were generated by using Illumina MiSeq. The sequencing data have been deposited in the bioproject PRJNA522361. Samples accession numbers vary from SAMN10964863 to SAMN10965438 for gyrB amplicon and from SAMN10970131 to SAMN10970566 for 16S.

Detection of an amplification bias associated to Leuconostocaceae family with a universal primer routinely used for monitoring microbial community structures within food products

Objectives

Sequencing of 16S rDNA V3–V4 region is widely applied for food community profiling. However, two different universal forward primers (named here MUYZER-primer1 and KLINDWORTH-primer2) targeting an identical conservative sequence upstream of the V3 region of 16S rRNA gene, and only distinguished by a single mismatch are both used. This study was carried out to compare whether the accuracy of food microbiota analysis would depend on the choice of one of these two primers.

Results

Alignment of both primers with common food-borne bacteria 16S sequences revealed that the mismatch between both primers might specifically affect the amplification of Leuconostoc, Oenococcus and Fructobacillus species but not Weissella species. Food products containing either Leuconostoc and/or Weissella were selected for a detection test. As expected from our in silico analysis, our study showed that this mismatch induced a strong biased amplification specifically associated to the OTUs belonging to the genus Leuconostoc but not to the genus Weissella. In presence of Muyzer-primer1, none of the sequences expected for Leuconostoc genus was detected whereas those sequences were correctly amplified with Klindworth-primer2. Since Leuconostoc is an important genus in food, agro-environments and in digestive tract of animals, we recommend that Muyzer-primer1 should thus be abandoned for the bacterial characterization of their associated microbiota.

Phylogenomic Analysis of Lactobacillus curvatus Reveals Two Lineages Distinguished by Genes for Fermenting Plant-Derived Carbohydrates

Lactobacillus curvatus is a lactic acid bacterium encountered in many different types of fermented food (meat, seafood, vegetables, and cereals). Although this species plays an important role in the preservation of these foods, few attempts have been made to assess its genomic diversity. This study uses comparative analyses of 13 published genomes (complete or draft) to better understand the evolutionary processes acting on the genome of this species. Phylogenomic analysis, based on a coalescent model of evolution, revealed that the 6,742 sites of single nucleotide polymorphism within the L. curvatus core genome delineate two major groups, with lineage 1 represented by the newly sequenced strain FLEC03, and lineage 2 represented by the type-strain DSM20019. The two lineages could also be distinguished by the content of their accessory genome, which sheds light on a long-term evolutionary process of lineage-dependent genetic acquisition and the possibility of population structure. Interestingly, one clade from lineage 2 shared more accessory genes with strains of lineage 1 than with other strains of lineage 2, indicating recent convergence in carbohydrate catabolism. Both lineages had a wide repertoire of accessory genes involved in the fermentation of plant-derived carbohydrates that are released from polymers of α/β-glucans, α/β-fructans, and N-acetylglucosan. Other gene clusters were distributed among strains according to the type of food from which the strains were isolated. These results give new insight into the ecological niches in which L. curvatus may naturally thrive (such as silage or compost heaps) in addition to fermented food.

Deciphering intra-species bacterial diversity of meat and seafood spoilage microbiota using gyrB amplicon sequencing: A comparative analysis with 16S rDNA V3-V4 amplicon sequencing

Meat and seafood spoilage ecosystems harbor extensive bacterial genomic diversity that is mainly found within a small number of species but within a large number of strains with different spoilage metabolic potential. To decipher the intraspecies diversity of such microbiota, traditional metagenetic analysis using the 16S rRNA gene is inadequate. We therefore assessed the potential benefit of an alternative genetic marker, gyrB, which encodes the subunit B of DNA gyrase, a type II DNA topoisomerase. A comparison between 16S rDNA-based (V3-V4) amplicon sequencing and gyrB-based amplicon sequencing was carried out in five types of meat and seafood products, with five mock communities serving as quality controls. Our results revealed that bacterial richness in these mock communities and food samples was estimated with higher accuracy using gyrB than using16S rDNA. However, for Firmicutes species, 35% of putative gyrB reads were actually identified as sequences of a gyrB paralog, parE, which encodes subunit B of topoisomerase IV; we therefore constructed a reference database of published sequences of both gyrB and pare for use in all subsequent analyses. Despite this co-amplification, the deviation between relative sequencing quantification and absolute qPCR quantification was comparable to that observed for 16S rDNA for all the tested species. This confirms that gyrB can be used successfully alongside 16S rDNA to determine the species composition (richness and evenness) of food microbiota. The major benefit of gyrB sequencing is its potential for improving taxonomic assignment and for further investigating OTU richness at the subspecies level, thus allowing more accurate discrimination of samples. Indeed, 80% of the reads of the 16S rDNA dataset were represented by thirteen 16S rDNA-based OTUs that could not be assigned at the species-level. Instead, these same clades corresponded to 44 gyrB-based OTUs, which differentiated various lineages down to the subspecies level. The increased ability of gyrB-based analyses to track and trace phylogenetically different groups of strains will generate improved resolution and more reliable results for studies of the strains implicated in food processes.

Reducing Salt in Raw Pork Sausages Increases Spoilage and Correlates with Reduced Bacterial Diversity

Raw sausages are perishable foodstuffs; reducing their salt content raises questions about a possible increased spoilage of these products. In this study, we evaluated the influence of salt reduction (from 2.0% to 1.5% [wt/wt]), in combination with two types of packaging (modified atmosphere [50% mix of CO2-N2] and vacuum packaging), on the onset of spoilage and on the diversity of spoilage-associated bacteria. After 21 days of storage at 8°C, spoilage was easily observed, characterized by noticeable graying of the products and the production of gas and off-odors defined as rancid, sulfurous, or sour. At least one of these types of spoilage occurred in each sample, and the global spoilage intensity was more pronounced in samples stored under modified atmosphere than under vacuum packaging and in samples with the lower salt content. Metagenetic 16S rRNA pyrosequencing revealed that vacuum-packaged samples contained a higher total bacterial richness (n = 69 operational taxonomic units [OTUs]) than samples under the other packaging condition (n = 46 OTUs). The core community was composed of 6 OTUs (Lactobacillus sakei, Lactococcus piscium, Carnobacterium divergens, Carnobacterium maltaromaticum, Serratia proteamaculans, and Brochothrix thermosphacta), whereas 13 OTUs taxonomically assigned to the Enterobacteriaceae, Enterococcaceae, and Leuconostocaceae families comprised a less-abundant subpopulation. This subdominant community was significantly more abundant when 2.0% salt and vacuum packaging were used, and this correlated with a lower degree of spoilage. Our results demonstrate that salt reduction, particularly when it is combined with CO2-enriched packaging, promotes faster spoilage of raw sausages by lowering the overall bacterial diversity (both richness and evenness).

IMPORTANCE

Our study takes place in the context of raw meat product manufacturing and is linked to a requirement for salt reduction. Health guidelines are calling for a reduction in dietary salt intake. However, salt has been used for a very long time as a hurdle technology, and salt reduction in meat products raises the question of spoilage and waste of food. The study was conceived to assess the role of sodium chloride reduction in meat products, both at the level of spoilage development and at the level of bacterial diversity, using 16S rRNA amplicon sequencing and raw pork sausage as a meat model.

Selection of bioprotective cultures for preventing cold-smoked salmon spoilage

Biopreservation is a natural technology of food preservation, which consists of inoculating food with microorganisms selected for their antibacterial properties. The objective of this study was to select lactic acid bacteria (LAB) to improve the quality of cold-smoked salmon (CSS). In this work, different strains representative of the 4 dominant species, identified in a previous study by pyrosequencing the 16S rRNA gene, were isolated and their spoiling potential in CSS blocks, sterilized by ionization, was assessed by twelve trained panelists along the vacuum storage at 8°C. Photobacterium phosphoreum, Brochothrix thermosphacta and Serratia proteamaculans released strong off-odors whereas the spoiling potential of Carnobacterium divergens was weaker. The spoiling capacity of Lactococcus piscium EU2241, Leuconostoc gelidum EU2247, Lactobacillus sakei EU2885, Staphylococcus equorum S030674 and 4 commercial starters was tested by the same method and 2 strains were eliminated due to off-odor production. The effect of the 6 selected LAB against the 4 specific spoiling organisms (SSOs) selected was tested by challenge tests in sterile CSS blocks. The protective effect of the LAB differed from one SSO to another and no correlation could be established between the sensory improvement, SSO inhibition, and the implantation or acidification of protective cultures (PCs). All the PCs except L. piscium reduced the off-odors released by P. phosphoreum although some of them had no effect on its growth. S. equorum, which did not grow in CSS, favored the implantation of P. phosphoreum but prevented its off-odor formation. L. piscium was the only strain that prevented the spoilage of B. thermosphacta and S. proteamaculans although it did not grow very well and did not acidify the product. L. gelidum EU2247 inhibited the growth of these 2 SSOs and lowered the pH but had no effect on the sensory quality. Finally, L. piscium was tested in 2 naturally contaminated products, with a positive effect on 1 batch. This effect was not correlated with the microbial ecosystem as determined by acultural and cultural techniques. Based on these results, the selection strategy is discussed.