Next Article in Journal
Biogeographic Pattern and Network of Rhizosphere Fungal and Bacterial Communities in Panicum miliaceum Fields: Roles of Abundant and Rare Taxa
Previous Article in Journal
Genomic Diversity of the Rarely Observed Genotype of the Mycobacterium tuberculosis Central Asian (CAS) Lineage 3 from North Brazil
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Distinct Microbiotas Are Associated with Different Production Lines in the Cutting Room of a Swine Slaughterhouse

1
NSERC Industrial Research Chair in Meat Safety (CRSV), Faculté de Médecine Vétérinaire, Université de Montréal, Saint-Hyacinthe, QC J2S 2M2, Canada
2
MESuRS Laboratory (Modelling, Epidemiology and Surveillance of Health Risks), Conservatoire National des Arts et Métiers (Cnam), 75003 Paris, France
3
F. Ménard, Division d’Olymel s.e.c., Ange-Gardien, QC J0E 1E0, Canada
4
CRIPA Swine and Poultry Infectious Diseases Research Center, Faculté de Médecine Vétérinaire, Université de Montréal, Saint-Hyacinthe, QC J2S 2M2, Canada
5
Le Conservatoire National des Arts et Métiers (Cnam), 75003 Paris, France
*
Author to whom correspondence should be addressed.
Microorganisms 2023, 11(1), 133; https://doi.org/10.3390/microorganisms11010133
Submission received: 12 December 2022 / Revised: 31 December 2022 / Accepted: 1 January 2023 / Published: 4 January 2023
(This article belongs to the Special Issue Foodborne Microorganisms and Biofilms)

Abstract

:
The microorganisms found on fresh, raw meat cuts at a slaughterhouse can influence the meat’s safety and spoilage patterns along further stages of processing. However, little is known about the general microbial ecology of the production environment of slaughterhouses. We used 16s rRNA sequencing and diversity analysis to characterize the microbiota heterogeneity on conveyor belt surfaces in the cutting room of a swine slaughterhouse from different production lines (each associated with a particular piece/cut of meat). Variation of the microbiota over a period of time (six visits) was also evaluated. Significant differences of alpha and beta diversity were found between the different visits and between the different production lines. Bacterial genera indicative of each visit and production line were also identified. We then created random forest models that, based on the microbiota of each sample, allowed us to predict with 94% accuracy to which visit a sample belonged and to predict with 88% accuracy from which production line it was taken. Our results suggest a possible influence of meat cut on processing surface microbiotas, which could lead to better prevention, surveillance, and control of microbial contamination of meat during processing.

1. Introduction

Providing safe and organoleptically satisfactory food for the consumer is a major concern for the food industry [1]. In Canada, 11 million tons of food produced each year are preventable waste and four million—or one in eight people—become sick every year from eating contaminated food [2,3]. Spoilage of food by microorganisms is probably responsible for a quarter of the food waste around the world, leading to significant economic losses and underscoring food processing environments as a major concern [4,5]. Moreover, foodborne diseases are the cause of several hundred thousand deaths worldwide every year [6,7]. Reducing food microbial contamination is a continuing battle for the primary food-processing sector, especially with meat products since they provide high nutrient availability for microorganisms and are known vectors of foodborne diseases [4,8].
The quality and safety attributes of final meat products are highly dependent on the raw materials and on their processing, storage conditions, and production environments [9,10]. Indeed, it has been shown that the production environment microbiota—in addition to being potentially involved in the spread of antimicrobial resistance genes [11] and in the alteration of the integrity of processing surfaces—may pose a threat to the attributes of meat [12,13]. Stellato et al. proposed the hypothesis that an equilibrium between the food and the production environment is established since the microbiota found on the surfaces and tools used during processing are often found on the meat [14,15]. Along this vein, current food production sanitary measures mainly rely on low contamination of the production environment while little attention has been paid, as of yet, to the structure of the bacteria communities and the transmission routes of microorganisms [4,16]. Recently, Stellato et al. proposed that since the microbiota of a site is under the influence of specific characteristics such as the type of surface involved, the environmental conditions, and the processes to which the food is subjected, it can be expected that each site will present a specific microbiota [17]. A few studies using diversity analysis found this not to be the case, with environmental processing surface-samples’ clustering according to sampling time [17] and not according to a specific sampling zone [15,17]. These conflicting results highlight the need for studies that address the presence of microenvironments in the food production environment by eliminating, as much as possible, the confounding factors that can influence results and make comparisons between studies difficult.
The recent introduction of high throughput sequencing technology has allowed the efficient characterization of complex and diverse microbial communities [18]. To date, only a limited number of studies have explored the general contamination flows and the microbial organization patterns of food processing environments using 16S rRNA amplicon sequencing [4,15]. The pairing of such sequencing techniques with computational methods has made it possible to infer the sources of origin of microorganisms of interest in a given environment [4]. Several bacterial source-tracking methods have been proposed to detect the origins of bacterial contamination such as Bayesian approaches and random forest algorithms (RF) [19,20].
Slaughterhouses are increasingly identified in the meat processing sector as a key site where unwanted microorganisms are found [4]. At the slaughterhouse, carcasses can be contaminated by the animal’s endogenous microbiota during processing. It has been shown that microorganisms can remain on the surface of a carcass and can survive and detach when they come into contact with surfaces and production equipment [21]. Thus, the contaminated contact surfaces can be sites of cross-contamination for the meat products circulating on them [18]. Indeed, contact surfaces such as conveyor belts have been shown to act as reservoirs for spoilage and pathogen bacteria [21]. Biofilms can settle in irregularities of conveyor belts and act as sources of contamination [22]. In natural environments, biofilms are known to be composed of multiple species [23,24]. The identity of the microorganisms constituting them as well as their interactions have been identified as contributors to the creation of local microbial ecosystems in the production environment [24,25]. Although the nature and quantity of microorganisms found on fresh raw meat at the slaughterhouse represent an initial point that can influence the safety and the spoilage patterns during further transport and storage, and can spread contamination to ready-to-eat facility environments [26,27], very little is known about the general microbial ecology of slaughterhouse production environments, especially with regard to the contamination during further processing stages such as cutting and deboning [16]. During these stages, which take place in the cutting room, the meat carcasses undergo different modifications in order to obtain parts intended to become standard meat cuts. These meat pieces will undergo their final transformation on specific production lines. Thus, a specific production line supports the passage of a specific meat cut. Recently, de Filippis et al. shed light on the significant association between beef microbiota and specific beef cuts [27]. They showed that different cuts from the same carcass can influence the microbial contamination of beef. The influence of the type of meat cut on the microbiota of contact surfaces (production lines) is yet to be established.
Using 16S rRNA sequencing and diversity analysis, the objectives of this study were: (i) to characterize the heterogeneity of the microbiota of conveyor belt contact surfaces present in the cutting room of a swine slaughterhouse in the province of Quebec, Canada, over a period a time (six visits), and to characterize the impact of the production line on the microbiota diversity of those surfaces; (ii) to identify microbial determinants associated with the visit and production line with a focus on bacterial determinants that can include known spoilage species and known foodborne pathogens; finally, (iii) to ascertain if, by using the microbiota of a given sample, it is possible to determine which visit and which production line the sample belongs to using random forest models as a predictive tool. We believe that understanding the composition of the microbiota on food processing surfaces is crucial since these distinct bacterial community dynamics are important in determining the microbial succession patterns that will occur and will influence the quality and safety of the final meat products. This knowledge could lead to more targeted cleaning and disinfecting procedures as well as to a more accurate management of the different meat cuts for further processing. Our study proposes an original approach based on diversity analysis and random forest models to evaluate, for the first time to our knowledge, the microbiota heterogeneity of the potential cross-contamination routes driven by the different production lines and visits in a cutting room of a swine slaughterhouse.

2. Materials and Methods

2.1. Facility Structure, Sampling, Total Microbiota Harvesting, DNA Extraction, Amplicon Library Preparation, and 16S rRNA Sequencing

It should be noted that this study is part of a larger project which was developed in two phases. The first phase aimed to evaluate spatial and temporal Listeria monocytogenes contamination on the six production lines in a swine slaughterhouse cutting room in Quebec, Canada and to identify microbial determinants of the presence or absence of the bacteria on those surfaces [28]. The current study is the second phase of that project and aims to evaluate the microbiota’s homogeneity on contact surfaces according to processing line (each line is dedicated to a specific cut of meat) and according to different visits using the samples harvested during the first phase. The sampling procedures as well as the processing of the microbiota samples have been described in detail in the work of Shedleur-Bourguignon et al., 2021 [28].
Briefly, the conveyor belt surfaces of six production lines present in the cutting room of a swine slaughterhouse with a capacity of 600 pigs per hour were sampled during midday production on six occasions between September 2017 and March 2018. Surfaces of 900 cm2 were sampled, first by using small brushes in order to mobilize the surface microbiota and second by collecting the biomass with sterile wipes (Innovation Diagnostics, Saint-Eustache, QC, Canada). Six production lines were present in the cutting room, and each production line supported the passage of a specific part of the carcass intended to become a standard meat piece (according to North American standards): half-carcasses, bellies, loins, Bostons, picnics, or hams. Fifty-four surface samples were taken on the main production line (CP), 48 on the belly production line (FL), 48 on the loin production line (LO), 48 on the Boston production line (BO), 48 on the picnic production line (PI), and 48 on the ham production line (FE). One experimental control sample per visit was also taken, which consisted of wipes and brushes that were taken out of their packaging in the cutting room and not allowed to come into contact with any surfaces. Specifically, upon arrival at the laboratory facilities, the heads of the brushes were cut and added to Falcon tubes containing 25 mL of a solution composed of Tris-HCL (10 mM), EDTA (10 mM), and NaCl (0.85%). The used wipes were subsequently cut in half under sterile conditions at the lab and each half wipe was added to the solution and the samples were stomached (Biomérieux Canada, QC, Canada) (the other half wipe was tested in a previous study [28]). Twenty mL was then transferred into two conical tubes of 15 mL (Sarstedt Inc, Saint-Laurent, QC, Canada). The Falcon tubes were then centrifuged (VWR, Saint-Laurent, QC, Canada), the supernatants were discarded, and the bacterial pellets were individually stored at −80 °C until DNA extraction and purification.
The two bacterial pellets belonging to the same sample were pooled. The total DNA of the combined bacterial pellets was then extracted and purified using a modified version of a phenol–chloroform protocol as previously described [28]. Final DNA concentration was measured using the Qubit 3.0 High Sensitivity range assay (Fisher Scientific, Ottawa, ON, Canada). The purity of the DNA was evaluated by using a Nanodrop ND-1000 Spectrophotometer (Thermo Scientific, Wilmington, DE, USA) and by gel electrophoresis (1% of agarose). A 291 bp segment of the V4 hypervariable region of the 16S rRNA gene was amplified using universal primers (515F_Ill and 806R_Ill; Invitrogen, Thermo Fisher Scientific, Waltham, MA, USA) [29,30]. A 30 µL PCR reaction was carried out using the Platinum SuperFi PCR Master Mix (Invitrogen, Burlington, ON, Canada) as described by Shedleur-Bourguignon et al. 2021 [28]. An artificial community (ZymoBIOMICS Microbial Community DNA Standard) (Zymo Research, Irvine, CA, USA) served as a positive control and as an indicator of sequencing quality. The theoretical composition of the artificial community based on 16S sequencing is: 18.4% Lactobacillus; 17.4% Bacillus; 15.5% Staphylococcus; 14.1% Listeria; 10.4% Salmonella; 10.1% Escherichia, 4.2% Pseudomonas, and 9.9% Enterococcus. Six experimentation controls (wipes) and five negative sequencing controls (DNA extraction and PCR steps) were integrated into the plates. The amplicons were then sent to McGill University Genome Quebec Innovation Center (Montreal, QC, Canada) for purification, barcoding, and sequencing by Illumina MiSeq 250 paired-end sequencing.

2.2. Sequencing Data Processing, Diversity, and Statistical Analysis

Diversity analyses were performed in order to characterize the composition and the structure of bacterial populations present on conveyor belt contact surfaces. As previously described by Shedleur-Bourguignon et al., the cleaning and the analyzing of the sequences were executed using Mothur 1.39.5 [28,31]. Alpha and beta diversity analysis were performed using RStudio 3.6.1. The lowest number of sequences in the samples was used as subsampling. For the alpha diversity, the coverage of the subsampling was measured as well as the number of operational taxonomic units (OTUs) in each sample, and their evenness was measured using the inverse Simpson and the Shannon indices. The indices were calculated by using the estimate_richness function of the phyloseq package version 1.24.2 [32]. Comparison statistics were performed between the samples from the different visits and between the different production line samples using Student t-tests and the Kruskal–Wallis test, both with a significance level of 0.05. For the beta diversity, the Jaccard index based on the presence/absence of OTUs and the Bray–Curtis index based on the relative abundance of OTUs were used on the same subsampling. The distance matrices were plotted using 2D nonmetric multidimensional scaling graphs (NMDS) using the ordinate function of the phyloseq package and the ggplot2 package [33], and statistical variations between the different groups (visits and production lines) were assessed using permutational analysis of molecular variance (PERMANOVA) using the ADONIS function of the vegan package [34]. The Bonferroni correction was used in order to correct the significance level in the presence of multiple comparisons. The impact of production lines on the beta diversity of the visits was taken into consideration by adding the production line factor in strata during statistical tests. The same procedure was used to remove the influence of the visits on the beta diversity of the production lines. The Multivariate Association with Linear Models method (MaAsLin version1) was used to identify OTUs that were significantly, positively, or negatively associated with each visit and with each production line in terms of relative abundance using the MaAsLin package [35]. The two variables (visit and production line) were used as a fixed effect and were put into the boost formula to consider the influence of the variables on each other. An association was considered significant at a p-value ≤ 0.05 and a q-value ≤ 0.25.

2.3. Evaluation of the Predictability of a Visit and a Production Line Using Random Forest Models

The models were built in order to classify the samples according to the visit during which they were collected or according to the production line from which they came based on their microbiota. The high number of covariates (OTUs) implicated in the analysis as well as the potential interactions between them guided the choice of using a nonparametric random forest approach as a predictive tool [36]. Two predictive models were built with RStudio 3.6.1 using the RandomForest and the caret packages [37,38]. For each of the two models, eight randomly chosen samples per visit and eight randomly chosen samples per production line were used as a test database to evaluate the performance of the models in order to avoid overfitting. The remaining samples were used to fit the random forest models. As each run of random forest could produce different results, one hundred random forests were performed per model to assess the robustness of the results. A median accuracy (rate of correct predictions) and an empirical 95% confidence interval were estimated to evaluate the performance for each model on the 100 iterations. Confusion matrices of the predictions were also computed on the first iteration of random forest for each model using the yardstick package [39]. The importance of each OTU in the prediction of a sample belonging to a particular production line or visit was determined to identify the most relevant OTUs. The importance was computed using the median of the decrease in accuracy when an OTU is removed from the model among the 100 replications of random forest (Liaw and al., 2002). In order to identify these OTUs, a random noise variable was introduced in the model. OTUs with a higher importance than the random noise in 95% of the replications of the random forest were considered as important in the prediction.

3. Results

3.1. 16S rRNA Gene Amplicon Sequencing

The V4 region of the 16S rRNA gene was sequenced from 310 samples: 294 conveyor belt surface samples were collected on six production lines over six visits, six experimentation controls (one per visit), five negative sequencing controls (PCR and sequencing negative controls), and five DNA positive controls (ZymoBIOMICS Microbial Community DNA Standard). A total of 7,198,363 sequences with an average of 23,437 sequences per sample were obtained after read processing and were grouped into 33 phyla, 84 classes, 152 orders, 313 families, and 952 operational taxonomic units (OTUs). The experimentation controls and the negative sequencing controls contained an average of 16,536 sequences and 7692 sequences, respectively, and no band was visible on gel electrophoresis for these two types of controls prior to sequencing. The ZymoBIOMICS Microbial Community DNA Standard positive controls were composed of the eight bacterial genera expected and showed an average of 17,733 sequences. The average composition of the positive samples based on 16S sequencing was: 16.2% Salmonella; 10.0% Escherichia/Shigella; 13.8% Bacillus; 13.2% Staphylococcus; 17.7% Lactobacillus; 11.5% Listeria; 9.5% Pseudomonas; and 4.4% Enterococcus.

3.2. Production Line Surface Microbiota Description

In order to avoid the description of the sporadic presence of bacterial genera and to instead focus the analyses on results of higher biological value, only the bacterial genera (OTUs) with a relative abundance equal to or greater than 5% in at least five samples from a visit or a production line were compiled in Table 1 and Table 2. Therefore, the 15 most abundant OTUs for all visits combined were from the most abundant to least abundant: Fusobacterium (13.5%); Trueperella (12.9%); Peptoniphilus (12.7%); Acinetobacter (11.9%); Pseudomonas (11.1%); Bacteroides_unclassified (10.6%); Rothia (9.9%); Aerococcaceae_unclassified (9.7%); Porphyromonas (9.7%); Clostridium_sensu_stricto (9.6%); Enhydrobacter (9.5%); Staphylococcus (9.3%); Macrococcus (8.7%); Psychrobacter (8.4%); and Halomonas (8.3%). The 15 most abundant OTUs for all the production lines combined were: Fusobacterium (18.8%); Rothia (15.6%); Trueperella (14.7%); Peptoniphilus (14.6%); Acinetobacter (13.8%); Staphylococcus (12.6%); Enhydrobacter (12.1%); Aerococcaceae_unclassified (11.3%); Peptostreptococcus (10.7%); Bacteroides_unclassified (10.5%); Macrococcus (10.3%); Pseudomonas (10.1%); Porphyromonas (10.1%); Clostridium_sensu_stricto (9.6%), and Psychrobacter (7.6%).

3.3. Microbiota Diversity of Conveyor Belt Surfaces Based on Visit and on Production Line

3.3.1. Alpha Diversity

Alpha diversity analyses were performed in order to describe the bacterial richness and distribution within samples depending on visit and production line using the Observed index, the Shannon Evenness index, and the Inverse Simpson index. Alpha diversity measures of visit samples and production line samples are shown in Table 3 and Table 4. Student’s t-test revealed differences in alpha diversity between several pairs of visits and between multiple pairs of production lines using the three indices. Kruskal–Wallis tests showed significant differences between the visits with the Shannon Evenness and Inverse Simpson indices as well as differences between the production lines with the three indices.

3.3.2. Beta Diversity

The similarity of the microbiota’s structure was compared at the OTU level between the samples according to the different visits using the Bray–Curtis index and the Jaccard index. The similarity of the microbiota’s structure was also compared between the samples according to the different production lines using the same indices. Statistically significant differences were found between all pairs of visits (p ˂ 0.001) with both indices. Figure 1—which illustrates the structure of the microbiota between samples belonging to the different visits using the Bray–Curtis index—does not show apparent clustering between samples belonging to the same visit (see Figure 1). However, when each visit was plotted in two-by-two NMDS plots, more important differences emerged between some pairs of visits: V1 and V3, V1 and V4, V1 and V5, V3 and V4, and V4 and V5 (see Figure 2). Statistically significant differences were also found between all pairs of production lines (p ˂ 0.001) with the Jaccard and the Bray–Curtis indices. These differences in bacterial community composition of the different production lines can be clearly visualized in an NMDS graph using the Bray–Curtis index (see Figure 3). Figure 3 shows a clustering between samples from the same production line. This clustering is even more apparent when the different production lines were plotted two-by-two (see Figure 4). Similar distance matrix NMDS plots were obtained using the Jaccard index (see Figures S1–S4 (Supplementary Materials)).

3.3.3. Multivariate Association with Linear Model Analysis

Multivariate association with linear model analysis (MaAsLin 1) was assessed to identify OTUs positively or negatively associated with the different visits and the different production lines. The number of total OTUs positively and negatively associated with the different visits in terms of relative abundance as well as a summary of associated OTUs of interest (bacterial genera, which can include known spoilage species and known foodborne pathogens species) are detailed in Table 5. Production line associations from the MaAsLin test are presented in Table 6. The information is presented in the same way as the MaAsLin visit data. The complete list of the taxonomic attribution of OTUs associated with the various visits and production lines is available in the Supplementary Materials (see Tables S1–S4).

3.4. Random Forest

The two models based on random forests were used to classify the samples according to the visit during which they were collected or according to the production line from which they came based on their microbiota. The confusion matrices obtained for the visit model and for the production line model are presented in Figure 5 and Figure 6, respectively. The visit model had an average prediction rate of 94.3% (89.6–100.0%). Visit one showed the best prediction rate (99.0%) and visit three showed the lowest rate of prediction (90.0%). Visits two, four, five, and six showed prediction rates of 95.8%, 91.5%, 92.9%, and 96.9%, respectively. The production line model has an average production rate of 88.1% (81.3–93.9%). The ham production line showed the best prediction rate (96.8%), while the loin production line showed the lowest prediction rate (77.5%). The main production line, the belly production line, the Boston production line, and the picnic production line, respectively, showed prediction rates of 87.4%, 92.5%, 96.3%, and 78.0%. The importance of variables (OTUs) in the prediction of a visit or a production line was evaluated. Among the 952 OTUs determined in the study, only a small proportion of them—143 and 182—were considered important to predict the visits and the production lines, respectively. The complete lists of these OTUs can be seen in Tables S5 and S6.

4. Discussion

In our study, a total of 294 samples of production line surfaces associated with the circulation of different meat cuts were collected in the cutting room of a swine slaughterhouse. The microbiota of the samples was analyzed by 16S rRNA amplicon sequencing in order to evaluate the impact of visit and production line on the diversity and structure of bacterial communities found on contact surfaces. The majority of the bacterial sequences collected in our study belonged to the phyla Proteobacteria (36.4%), Firmicutes (35.9%), Fusobacteria (31.4%), Actinobacteria (30.0%), and Bacteroidetes (27.2%). The dominance of these phyla in pig slaughterhouse environments has been reported by previous studies, with the exception of the Fusobacteria phylum [40,41]. The Fusobacteria phylum is known to be a normal constituent of the oropharyngeal and gastrointestinal microbiota of pigs [42]. A recent study by Wylensek et al. revealed pig-specific species within the Fusobacterium genus [43]. Thus, bacteria belonging to this phylum could have been transferred from the gastrointestinal tract of pigs to the carcasses during evisceration and then transferred to surfaces in contact with the meat [21,44]. The dominance of the phylum Fusobacteria can potentially be explained by the specific context of our study since our study was interested in the microbiota of conveyer belt surfaces in a cutting room during production while most of the available studies focused on multiple slaughterhouse compartments/rooms [12]. Our results showed that the Fusobacteria phylum had a significantly higher relative abundance during the fourth and the sixth visits and that the phylum was important for the identification of the visit and production line in the random forest models.
Several genera found among the 15 most abundant bacteria across all visits and production lines combined have been characterized in the literature as having phenotypes allowing them to counteract the stress conditions encountered in the production environment. Indeed, the dominant presence of psychotropic genera such as Pseudomonas, Acinetobacter, and Psychrobacter was expected due to the refrigeration temperatures maintained in the cutting room [12]. A study by Botta et al. correlated the greater abundance and persistence of these three genera in secondary processing rooms with the lowest temperatures of these environments [26]. In addition, several of the most abundant phyla in our study, such as Staphylococcus, are known to be able to form biofilms. Biofilms can confer a greater resistance to disinfection products to bacteria but also to shear forces produced by moving parts in production [45]. A recurrent introduction of microorganisms by the carcasses following the contact of meat products with processing surfaces in the cutting room can also not be excluded. Several bacterial genera with potentially pathogenic species and/or potential spoilage species were among the OTUs with a relative abundance equal to or greater than five percent in at least five samples from a visit or production line. The presence of these OTUs on fresh pork meat have been previously described [10]. However, the presence of the bacterial genus Clostridium among the 15 most abundant OTUs during visits and on production lines and its strong presence in the cutting room was not expected. Even though the Clostridium genus is not known to be a problem in the environmental processing environment of swine slaughterhouses regarding its potential pathogenicity for humans, the high representation of the genus Clostridium in our study suggests that the presence of this potential foodborne pathogen may be underestimated.
Variability in the mean relative abundance of frequent bacterial genera (OTUs) was identified. Indeed, several genera showed a higher relative abundance during certain visits or on certain production lines. For instance, the genus Staphylococcus showed a higher relative abundance during the third visit. The genera Peptostreptococcus (V2) and Macrococcus (V4) are also examples of bacteria overrepresented during a single visit. In addition, our study showed statistical differences in terms of alpha diversity between several pairs of visits with the three indices. Moreover, statistically significant differences of beta diversity were found between all pairs of visits, indicating that the structures of major and minor bacterial populations are affected by the day of production. Put together, these results, as expected, suggest that the number of different bacterial genera and the uniformity of their distribution on surfaces are affected by the day of production. Indeed, in a different context, Stellato et al. investigated the bacterial biogeographical patterns on surfaces and tools in a hospital cooking center and showed that, based on the composition of the microbiota, samples from food production environments grouped according to sampling time (two months separated the two samplings) [17]. As the slaughter process is thought to present very little variation, it is possible to hypothesize that the differences in the mean relative abundance as well as in alpha and beta diversity are mainly supported by the arrival of different inputs (carcasses) in the cutting room, by the staff working during these different days, and by differences in the application of washing and disinfection measures. However, a study by Braley et al., conducted in the same facility as our study, revealed no significant difference of microbiota between pig carcass batches [46]. The fact that the study conducted by Braley et al. was realized during a short period of time (one day) and the sampling in our study took place at monthly intervals could explain these variations in the impact of time. Another study by Lim et al. showed that depending on the season, changes in the abundance of certain bacterial genera could be observed on food-contact and non-food contact surfaces in a foodservice facility [18]. The impact of season on the microbiota composition of fresh beef meat was also highlighted in a study by Hwang et al. [47]. As our sampling took place over six months (one sampling each month) and covered three seasons (summer, autumn, and winter, which are notably distinct in Canada), it is possible that conditions specific to the seasons have influenced the quantity and diversity of the microorganisms harvested. It should be noted that the impact of the visits on the diversity of the processing surfaces microbiotas associated with production line was taken into consideration in the present study.
Variability in the mean relative abundance of bacterial genera (OTUs) was also identified in relation to the production line. Indeed, the genera Trueperella (FL), Pseudomonas (FE), Acinetobacter (LO), Rothia (LO), and Staphylococcus (LO) presented a higher mean relative abundance on certain production lines. In the same perspective, a study by Biasino et al. mapped the distribution of microbiological contamination of pig carcasses [44]. They showed significant differences, using a culture-dependent methodology, in the contamination level of the different carcass areas regarding the total aerobic bacteria Enterobacteriaceae and Salmonella. De Filippis et al., for their part, have highlighted a significant association between beef microbiota and specific beef cuts using 16S rRNA amplicon sequencing [27]. The results of these two studies that highlight the influence of different meat cuts on the microbial contamination of meat are in line with our observations but leave unanswered the link between the microbiota of the production surfaces and the pieces of meat they support.
Our study showed statistical differences in terms of alpha diversity between several pairs of production lines. The highest alpha diversity measures have been associated with the CP and the FE production lines. The BO production line showed the lowest alpha diversity values while the FL, LO, and PI production lines showed variable levels of diversity depending on the index used. In addition, statistically significant differences of beta diversity were found between all pairs of production lines. This is the first time to our knowledge that a study has revealed the presence of bacterial microbiotas distinct in number, uniformity, and structure between production lines in a cutting room of a slaughterhouse. As the only difference that exists between the conveyors associated with the different production lines is the type of meat cut circulating on them, our study represents a first step in the exploration of the influence of meat cut on surface contact microbiota. However, these results were obtained in a single slaughterhouse, therefore it would be interesting and relevant to study whether these observations are present in other cutting rooms of pig slaughterhouses. These results concur with several studies which have highlighted the presence of different microbiotas in the food production environment. Zwirzitz et al. applied the software SourceTracker on full-length 16S rRNA gene sequencing data from a swine slaughterhouse and revealed that contact surfaces such as the polishing tunnel and a railing used for the classification of the carcasses strongly contributed to the composition of the microbiota of meat samples [4]. The study showed that many bacterial genera were unique to specific sites in a swine slaughterhouse, indicating that specific microorganisms occupy environmental niches in the facility. The authors suggested that the spatial distribution of the microorganisms is the result of specific transmission routes within the facility. Botta et al. conducted a study in three red-meat slaughterhouses and revealed significant differences in the meat processing-plant surface microbiota between the environments (the three different plants), the temporal phases (before or after cleaning and sanitizing), and between the type of room (deboning rooms and processing room). Although these studies have explored the general contamination flows and the microbial organization patterns [1,4] of food processing environments, none of them focused on a single compartment/room in a food plant.
Multivariate association with linear model analysis (MaAsLin) was assessed to identify OTUs positively or negatively associated with the different visits and the different production lines. Among the determinants positively or negatively associated with the different visits, 12 have been previously identified as possibly containing known spoilage species and known foodborne pathogen species: Acinetobacter, Flavobacterium, Pseudomonas, Psychrobacter, Lactobacillus, Lactococcus, Streptococcus, Brochothrix, Carnobacterium, Clostridium, Enterococcus, and Aerococcus. The identification of visits during which these microbial genera showed significantly higher relative abundance on contact surfaces raises the question as to whether there are differences in the relative abundance of these bacterial genera on the carcasses. Indeed, it has been shown that microorganisms can remain on the surface of the carcass, survive, and detach when they come in contact with surfaces and production equipment [21]. Six of those bacterial genera currently recognized as spoilage determinants were lactic acid bacteria genera (LAB): Lactobacillus, Lactococcus, Carnobacterium, Aerococcus, Streptococcus, and Enterococcus. LAB are recognized as important competitors of other spoilage microorganisms and have been reported to be associated with the souring, formation of slime, bone taint, and greening of refrigerated fresh raw meat [10,48,49]. In addition to the lactic acid bacteria, the genus Brochothrix (Otu00245), another Gram-positive bacterium, emerged as a determinant in the MaAsLin tests. The species B. thermosphacta is known as an important spoiler of various food matrixes, including the formation of slime on pork meat [10,48]. A study by Hultman and al. suggested that refrigerated food-processing environments provide conditions for the persistence of Brochothrix spp and LAB and that meat-contact surfaces could act as reservoirs for those bacteria [15]. The Acinetobacter and Pseudomonas genera, for their part, are considered as important sources of meat spoilage products in aerobic conditions [50], while several species of the genus Clostridium have been identified as causative agents for defects of meat in anaerobic conditions [49].
Regarding the production lines, 12 have been previously identified as common potential meat spoilage bacterial genera. In our study, several LAB (Lactococcus, Lactobacillus, Carnobacterium, and Enterococcus) showed a positive association in terms of relative abundance with the CP, LO, and FE production lines, suggesting predicted locations of LAB reservoirs on production lines. The genus Clostridium was the only genus identified as a determinant of five production lines (CP, FL, LO, PI, and FE). However, it is interesting to note that the distribution of the OTUs (OTUs 00016, 00059, 00136, 00203, 00307, and 00441) associated with this bacterial genus is differentially organized. For example, the OTU 00203 and the OTUs 00307 and 00441 were found only on the PI and CP production lines, respectively. This suggests that particular Clostridium species would be specifically associated with certain production lines. Thus, since only a fraction of the Clostridium species is identified as a causative agent for defects in vacuum-packed meat [49], precisely identifying the species associated with the different production lines could help identify the cuts of meat best-suited for this type of packaging. Except for the genus Clostridium, no OTU identified as potential spoilage bacteria showed a positive association with the FL production line. On the contrary, four bacterial genera associated with meat spoilage showed a negative association with the FL production line: Psychrobacter, Moraxella, Acinetobacter, and Pseudomonas [14,51]. The PI production line also showed a few positive associations (Clostridium and Psychrobacter genera) but also three negative associations (Acinetobacter, Pseudomonas, and Salmonella genera). While Psychrobacter and Moraxella possess a low spoilage potential due to their lack of several important biochemical attributes such as proteolysis, Pseudomonas, the species P. fluorescens, P. putida, and P. fragi, in particular contribute to the spoilage of raw meat to a large extent [52]. No association was found with the BO production line. These results suggest, in the context of our study, that the BO, the FL, and the PI production lines are the least-concerned production lines in terms of relative abundance regarding bacteria capable of deteriorating meat quality.
Four bacterial genera potentially pathogenic to humans were found among the bacterial determinants positively or negatively associated with the different visits. The potential pathogenic genus Campylobacter, of which species C. coli frequently contaminates pigs and consequently food of porcine origin, was identified as a microbial determinant of visits five and six [53]. The genus Listeria, of which the species L. monocytogenes has been linked to several outbreaks associated with the consumption of pork products, was associated with the second visit [54,55]. Based on these results, the introduction of these two potential pathogens seems to have been sporadic. It is therefore possible to hypothesize that the inputs (carcasses) introduced into the cutting room on these two dates had a higher relative abundance of Campylobacter or Listeria than the inputs from the other visits and could therefore have transferred these pathogens in higher concentrations to the contact surfaces. This observation is reinforced by the fact that the OTU associated with the genus Listeria was unique to visit two.
The Campylobacter genus also showed positive associations with the CP and the FE production lines in the same way as the genus Escherichia/Shigella, an indicator of fecal contamination, with the FL and the FE production lines [56]. The Salmonella genus, a major concern of pork meat-product contamination, showed a negative association with the PI production line [57,58]. The results of our study are partially in accordance with a study by Biasino et al. [44], which showed that the ventral and anterior areas (foreleg, head, sternum, and throat) were the most contaminated parts of pig carcasses by Enterobacteriaceae and Salmonella. These results are consistent with our study, which showed a lower relative abundance of the genus Salmonella on the PI production line. However, for E. coli, our study showed a higher relative abundance of the Escherichia/Shigella genera on the FE production line [44]. This divergence in results can be explained by differences between the production facilities sampled, such as the evisceration techniques and the cleaning and disinfection procedures, by the use of different detection methods, and also by the fact that the Biasino et al. study was conducted on carcass parts while our study was conducted on surfaces in contact with these carcass parts [44]. Unfortunately, few studies have examined the transmission of Campylobacter in the pork industry, which does not allow the comparison of our results. The Staphylococcus genus showed a unique positive association with the LO production line. This bacterial genus is a potential opportunistic pathogen of which pigs and humans are potential hosts [59]. As no association was found with any visit, the hypothesis of an introduction via workers on this particular production line is probable. The presence on several conveyors of the genus Clostridium, within which the C. perfingens species represents a risk to humans, has been discussed above. The identification (using the MaAsLin test) of microbial determinants positively or negatively associated with the different production lines shows for the first time to our knowledge the existence of significant differences in terms of relative abundance in the organization of microbial communities on the production lines associated with specific meat cuts at the slaughterhouse.
Having identified bacterial determinants associated with the different visits and the different production lines, two predictive models were built based on random forests as proof of concept in order to establish, through the predictive character of these models, the existence of distinct microbiotas in a cutting room of a slaughterhouse. Based on the microbiota of a sample, our models were able to predict its visit number at a rate of 94% and which production line it came from at a rate of 88%. The proof of concept provided by the creation of our predictive models suggests their usage on a larger scale and their potential integration into larger surveillance systems based on blockchain technology. Using these models could be helpful in predicting the origin of a problematic sample quickly and efficiently. The predictive value of surface provenance based on a sample of the microbiota of a piece of meat remains to be explored.
We believe that the results of our study, which identified for the first time to our knowledge the presence of distinct microbiotas on contact surfaces related to the circulation of specific meat cuts, represent a crucial step in the understanding of the microbial ecology of a slaughterhouse. The distinct bacterial community dynamics which take place on the production lines in a cutting room are important in determining the microbial succession pattern that will occur and ultimately influence the quality and safety of the final meat products. Therefore, this study is a step forward in the prevention, surveillance, and control of the microbial contamination of meat products.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/microorganisms11010133/s1, Table S1: OTUs positively associated with visit according to the MaAsLin test; Table S2: OTUs negatively associated with visit according to the MaAsLin test; Table S3: OTUs positively associated with production line according to the MaAsLin test; Table S4: OTUs negatively associated with production line according to the MaAsLin test; Table S5: OTUs hierarchized according to their mean importance for the identification of a visit in the random forest visit model; Table S6: OTUs hierarchized according to their mean importance for the identification of a production line in the random forest production line model.

Author Contributions

Conceptualization, P.F. and F.S.-B.; methodology, P.F., F.S.-B., T.D., M.N.H. and W.P.T.; software, F.S.-B., T.D., M.N.H. and A.T.; validation, F.S.-B., P.F., T.D., M.N.H. and A.T.; formal analysis, F.S.-B. and T.D.; resources, W.P.T. and J.L.; writing—original draft preparation, F.S.-B. and P.F.; writing—review and editing, F.S.-B., P.F., T.D., W.P.T., J.L., A.T. and M.N.H.; visualization, F.S.-B. and T.D.; supervision, P.F. and A.T.; project administration, P.F.; funding acquisition, P.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Natural Sciences and Engineering Research Council of Canada (NSERC, No.: RDCPJ 520873-17), F. Ménard, division d’Olymel s.e.c., financial partner in the Industrial Research Chair in Meat Safety and the Swine and Poultry Infectious Diseases Research Center (CRIPA).

Data Availability Statement

The data presented in this study are openly available in the NCBI Sequence Read Archive under accession number PRJNA758928 (https://www.ncbi.nlm.nih.gov/sra/PRJNA758928 (accessed on 1 September 2021).

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study, in the collection, analyses, or interpretation of data, in the final writing of the manuscript, or in the decision to publish the results.

References

  1. Bagge-Ravn, D. The microbial ecology of processing equipment in different fish industries—Analysis of the microflora during processing and following cleaning and disinfection. Int. J. Food Microbiol. 2003, 87, 239–250. [Google Scholar] [CrossRef]
  2. Public Health Agency of Canada. Food-Related Illnesses, Hospitalizations and Deaths in Canada; Government of Canada: Gatineau, QC, Canada, 2016.
  3. Environmental and Climate Change Canada. Taking Stock: Reducing Food Loss and Waste in Canada; Government of Canada: Gatineau, QC, Canada, 2019.
  4. Zwirzitz, B.; Wetzels, S.U.; Dixon, E.D.; Stessl, B.; Zaiser, A.; Rabanser, I.; Thalguter, S.; Pinior, B.; Roch, F.-F.; Strachan, C.; et al. The sources and transmission routes of microbial populations throughout a meat processing facility. NPJ Biofilms Microbiomes 2020, 6, 1–12. [Google Scholar] [CrossRef]
  5. Huis in’t Veld, J.H.J. Microbial and biochemical spoilage of foods: An overview. Int. J. Food Microbiol. 1996, 33, 1–18. [Google Scholar] [CrossRef]
  6. Havelaar, A.H.; Kirk, M.D.; Torgerson, P.R.; Gibb, H.J.; Hald, T.; Lake, R.J.; Praet, N.; Bellinger, D.C.; de Silva, N.R.; Gargouri, N.; et al. World Health Organization Global Estimates and Regional Comparisons of the Burden of Foodborne Disease in 2010. PLoS Med. 2015, 12, e1001923. [Google Scholar] [CrossRef] [Green Version]
  7. Zwirzitz, B.; Wetzels, S.U.; Dixon, E.D.; Fleischmann, S.; Selberherr, E.; Thalguter, S.; Quijada, N.M.; Dzieciol, M.; Wagner, M.; Stessl, B. Co-Occurrence of Listeria spp. and Spoilage Associated Microbiota During Meat Processing Due to Cross-Contamination Events. Front. Microbiol. 2021, 12, 632935. [Google Scholar] [CrossRef]
  8. Authority, E.F.S.J.E.J. The European Union summary report on trends and sources of zoonoses, zoonotic agents and food-borne outbreaks in 2016. EFSA J. 2017, 15, e05077. [Google Scholar]
  9. Pothakos, V.; Devlieghere, F.; Villani, F.; Björkroth, J.; Ercolini, D. Lactic acid bacteria and their controversial role in fresh meat spoilage. Meat Sci. 2015, 109, 66–74. [Google Scholar] [CrossRef]
  10. Nychas, G.-J.E.; Skandamis, P.N.; Tassou, C.C.; Koutsoumanis, K.P. Meat spoilage during distribution. Meat Sci. 2008, 78, 77–89. [Google Scholar] [CrossRef]
  11. Verraes, C.; Van Boxstael, S.; Van Meervenne, E.; Van Coillie, E.; Butaye, P.; Catry, B.; De Schaetzen, M.A.; Van Huffel, X.; Imberechts, H.; Dierick, K.; et al. Antimicrobial resistance in the food chain: A review. Int. J. Environ. Res. Public Health 2013, 10, 2643–2669. [Google Scholar] [CrossRef] [Green Version]
  12. Møretrø, T.; Langsrud, S. Residential Bacteria on Surfaces in the Food Industry and Their Implications for Food Safety and Quality. Compr. Rev. Food Sci. Food Saf. 2017, 16, 1022–1041. [Google Scholar] [CrossRef] [Green Version]
  13. Giaouris, E.; Heir, E.; Desvaux, M.; Hébraud, M.; Møretrø, T.; Langsrud, S.; Doulgeraki, A.; Nychas, G.-J.; Kačániová, M.; Czaczyk, K.; et al. Intra- and inter-species interactions within biofilms of important foodborne bacterial pathogens. Front. Microbiol. 2015, 6, 841. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Stellato, G.; La Storia, A.; De Filippis, F.; Borriello, G.; Villani, F.; Ercolini, D. Overlap of Spoilage-Associated Microbiota between Meat and the Meat Processing Environment in Small-Scale and Large-Scale Retail Distributions. Appl. Environ. Microbiol. 2016, 82, 4045–4054. [Google Scholar] [CrossRef] [Green Version]
  15. Hultman, J.; Rahkila, R.; Ali, J.; Rousu, J.; Bjorkroth, J. Meat Processing Plant Microbiome and Contamination Patterns of Cold-Tolerant Bacteria Causing Food Safety and Spoilage Risks in the Manufacture of Vacuum-Packaged Cooked Sausages. Appl. Environ. Microbiol. 2015, 81, 7088–7097. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Choi, Y.M.; Park, H.J.; Jang, H.I.; Kim, S.A.; Imm, J.Y.; Hwang, I.G.; Rhee, M.S. Changes in microbial contamination levels of porcine carcasses and fresh pork in slaughterhouses, processing lines, retail outlets, and local markets by commercial distribution. Res. Vet. Sci. 2013, 94, 413–418. [Google Scholar] [CrossRef]
  17. Stellato, G.; La Storia, A.; Cirillo, T.; Ercolini, D. Bacterial biogeographical patterns in a cooking center for hospital foodservice. Int. J. Food Microbiol. 2015, 193, 99–108. [Google Scholar] [CrossRef]
  18. Lim, E.S.; Kim, J.J.; Sul, W.J.; Kim, J.-S.; Kim, B.; Kim, H.; Koo, O.K. Metagenomic Analysis of Microbial Composition Revealed Cross-Contamination Pathway of Bacteria at a Foodservice Facility. Front. Microbiol. 2021, 12, 636329. [Google Scholar] [CrossRef]
  19. Smith, A.; Sterba-Boatwright, B.; Mott, J. Novel application of a statistical technique, Random Forests, in a bacterial source tracking study. Water Res. 2010, 44, 4067–4076. [Google Scholar] [CrossRef]
  20. Bokulich, N.A.; Bergsveinson, J.; Ziola, B.; Mills, D.A. Mapping microbial ecosystems and spoilage-gene flow in breweries highlights patterns of contamina-tion and resistance. eLife 2015, 4, e04634. [Google Scholar] [CrossRef]
  21. Midelet, G.; Carpentier, B. Transfer of microorganisms, including Listeria monocytogenes, from various materials to beef. Appl. Environ. Microbiol. 2002, 68, 4015–4024. [Google Scholar] [CrossRef] [Green Version]
  22. Fagerlund, A.; Møretrø, T.; Heir, E.; Briandet, R.; Langsrud, S. Cleaning and Disinfection of Biofilms Composed of Listeria monocytogenes and Background Microbiota from Meat Processing Surfaces. Appl. Environ. Microbiol. 2017, 83, e01046-17. [Google Scholar] [CrossRef] [Green Version]
  23. Bridier, A.; Sanchez-Vizuete, P.; Guilbaud, M.; Piard, J.-C.; Naïtali, M.; Briandet, R. Biofilm-associated persistence of food-borne pathogens. Food Microbiol. 2015, 45, 167–178. [Google Scholar] [CrossRef] [PubMed]
  24. Giaouris, E.; Heir, E.; Hébraud, M.; Chorianopoulos, N.; Langsrud, S.; Møretrø, T.; Habimana, O.; Desvaux, M.; Renier, S.; Nychas, G.-J. Attachment and biofilm formation by foodborne bacteria in meat processing environments: Causes, implications, role of bacterial interactions and control by alternative novel methods. Meat Sci. 2014, 97, 298–309. [Google Scholar] [CrossRef] [PubMed]
  25. Tolker-Nielsen, T.; Molin, S. Spatial Organization of Microbial Biofilm Communities. Microb. Ecol. 2000, 40, 75–84. [Google Scholar] [CrossRef] [PubMed]
  26. Botta, C.; Ferrocino, I.; Pessione, A.; Cocolin, L.; Rantsiou, K. Spatiotemporal Distribution of the Environmental Microbiota in Food Processing Plants as Impacted by Cleaning and Sanitizing Procedures: The Case of Slaughterhouses and Gaseous Ozone. Appl. Environ. Microbiol. 2020, 86, e01861-20. [Google Scholar] [CrossRef]
  27. De Filippis, F.; La Storia, A.; Villani, F.; Ercolini, D. Exploring the Sources of Bacterial Spoilers in Beefsteaks by Culture-Independent High-Throughput Se-quencing. PLoS ONE 2013, 8, e70222. [Google Scholar] [CrossRef] [Green Version]
  28. Shedleur-Bourguignon, F.; Thériault, W.P.; Longpré, J.; Thibodeau, A.; Fravalo, P. Use of an Ecosystem-Based Approach to Shed Light on the Heterogeneity of the Contamina-tion Pattern of Listeria monocytogenes on Conveyor Belt Surfaces in a Swine Slaughterhouse in the Province of Quebec, Canada. Pathogens 2021, 10, 1368. [Google Scholar] [CrossRef] [PubMed]
  29. Apprill, A.; McNally, S.; Parsons, R.; Weber, L. Minor revision to V4 region SSU rRNA 806R gene primer greatly increases detection of SAR11 bacterio-plankton. Aquat. Microb. Ecol. 2015, 75, 129–137. [Google Scholar] [CrossRef] [Green Version]
  30. Parada, A.E.; Needham, D.M.; Fuhrman, J.A. Every base matters: Assessing small subunit rRNA primers for marine mi-crobiomes with mock communities, time series and global field samples. Environ. Microbiol. 2016, 18, 1403–1414. [Google Scholar] [CrossRef]
  31. Larivière-Gauthier, G.; Thibodeau, A.; Letellier, A.; Yergeau, A.; Fravalo, P. Reduction of Salmonella Shedding by Sows during Gestation in Relation to Its Fecal Microbiome. Front. Microbiol. 2017, 8, 2219. [Google Scholar] [CrossRef] [Green Version]
  32. McMurdie, P.J.; Holmes, S. phyloseq: An R package for reproducible interactive analysis and graphics of microbiome census data. PLoS ONE 2013, 8, e61217. [Google Scholar] [CrossRef] [Green Version]
  33. Wickham, H. ggplot2: Elegant Graphics for Data Analysis; Springer: New York, NY, USA, 2016. [Google Scholar]
  34. Oksanen, J.; Blanchet, F.G.; Kindt, R.; Legendre, P.; Minchin, P.R.; O’hara, R.B.; Simpson, G.L.; Solymos, P.; Stevens, M.H.; Wagner, H. Vegan: Community Ecology Package; R Package Version 2.5-6. 2018. Available online: https://CRAN.R-project.org/package=vegan (accessed on 28 December 2022).
  35. Morgan, X.C.; Tickle, T.L.; Sokol, H.; Gevers, D.; Devaney, K.L.; Ward, D.V.; Reyes, J.A.; Shah, S.A.; LeLeiko, N.; Snapper, S.B.; et al. Dysfunction of the intestinal microbiome in inflammatory bowel disease and treatment. Genome Biol. 2012, 13, R79. [Google Scholar] [CrossRef] [PubMed]
  36. Duchemin, T.; Bar-Hen, A.; Lounissi, R.; Dab, W.; Hocine, M.N. Hierarchizing Determinants of Sick Leave: Insights From a Survey on Health and Well-being at the Workplace. J. Occup. Environ. Med. 2019, 61, e340–e347. [Google Scholar] [CrossRef] [PubMed]
  37. Breiman, L. Random Forests. Mach. Learn. 2001, 45, 5–32. [Google Scholar] [CrossRef] [Green Version]
  38. Kuhn, M. Building Predictive Models in R Using the caret Package. J. Stat. Softw. 2008, 28, 1–26. [Google Scholar] [CrossRef] [Green Version]
  39. Kuhn, M.; Vaughan, D. Yardstick: Tidy Characterizations of Model Performance. R package Version 0.0.8. 2021. Available online: https://CRAN.R-project.org/package=yardstick (accessed on 28 December 2022).
  40. Calero, G.C.; Gómez, N.C.; Lerma, L.L.; Benomar, N.; Knapp, C.W.; Abriouel, H. In silico mapping of microbial communities and stress responses in a porcine slaughterhouse and pork products through its production chain, and the efficacy of HLE disinfectant. Food Res. Int. 2020, 136, 109486. [Google Scholar] [CrossRef] [PubMed]
  41. Bridier, A.; Le Grandois, P.; Moreau, M.-H.; Prénom, C.; Le Roux, A.; Feurer, C.; Soumet, C. Impact of cleaning and disinfection procedures on microbial ecology and Salmonella antimicrobial re-sistance in a pig slaughterhouse. Sci. Rep. 2019, 9, 12947. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. De Witte, C.; Flahou, B.; Ducatelle, R.; Smet, A.; De Bruyne, E.; Cnockaert, M.; Taminiau, B.; Daube, G.; Vandamme, P.; Haesebrouck, F. Detection, isolation and characterization of Fusobacterium gastrosuis sp. nov. colonizing the stomach of pigs. Syst. Appl. Microbiol. 2016, 40, 42–50. [Google Scholar] [CrossRef] [PubMed]
  43. Wylensek, D.; Hitch, T.C.A.; Riedel, T.; Afrizal, A.; Kumar, N.; Wortmann, E.; Liu, T.; Devendran, S.; Lesker, T.R.; Hernández, S.B.; et al. A collection of bacterial isolates from the pig intestine reveals functional and taxonomic diversity. Nat. Commun. 2020, 11, 1–26. [Google Scholar] [CrossRef]
  44. Biasino, W.; De Zutter, L.; Mattheus, W.; Bertrand, S.; Uyttendaele, M.; Van Damme, I. Correlation between slaughter practices and the distribution of Salmonella and hygiene indicator bacteria on pig carcasses during slaughter. Food Microbiol. 2018, 70, 192–199. [Google Scholar] [CrossRef] [PubMed]
  45. Aggarwal, S.; Stewart, P.S.; Hozalski, R.M. Biofilm Cohesive Strength as a Basis for Biofilm Recalcitrance: Are Bacterial Biofilms Overdesigned? Microbiol. Insights 2015, 8, S31444. [Google Scholar] [CrossRef]
  46. Braley, C.; Fravalo, P.; Gaucher, M.-L.; Larivière-Gauthier, G.; Shedleur-Bourguignon, F.; Longpré, J.; Thibodeau, A. Similar Carcass Surface Microbiota Observed Following Primary Processing of Different Pig Batches. Front. Microbiol. 2022, 13, 849883. [Google Scholar] [CrossRef]
  47. Hwang, B.K.; Choi, H.; Choi, S.H.; Kim, B.-S. Analysis of Microbiota Structure and Potential Functions Influencing Spoilage of Fresh Beef Meat. Front. Microbiol. 2020, 11, 1657. [Google Scholar] [CrossRef] [PubMed]
  48. Nychas, G.-J.; Marshall, D.L.; Sofos, J. Meat, Poultry, and Seafood. In Food Microbiology: Fundamentals and Frontiers; ASM Press: Washington, DC, USA, 2007; pp. 105–140. [Google Scholar]
  49. Remenant, B.; Jaffrès, E.; Dousset, X.; Pilet, M.-F.; Zagorec, M. Bacterial spoilers of food: Behavior, fitness and functional properties. Food Microbiol. 2015, 45, 45–53. [Google Scholar] [CrossRef] [PubMed]
  50. Iulietto, M.F.; Sechi, P.; Borgogni, E.; Cenci-Goga, B.T. Meat Spoilage: A Critical Review of a Neglected Alteration Due to Ropy Slime Producing Bacteria. Ital. J. Anim. Sci. 2015, 14, 4011. [Google Scholar] [CrossRef]
  51. Doulgeraki, A.I.; Ercolini, D.; Villani, F.; Nychas, G.-J.E. Spoilage microbiota associated to the storage of raw meat in different conditions. Int. J. Food Microbiol. 2012, 157, 130–141. [Google Scholar] [CrossRef] [PubMed]
  52. Betts, G. 23—Other Spoilage Bacteria, in Food Spoilage Microorganisms; Blackburn, C.d.W., Ed.; Woodhead Publishing: Sawston, UK, 2006; pp. 668–693. [Google Scholar]
  53. Moore, J.E.; Madden, R.H. Survival of Campylobacter coli in porcine liver. Food Microbiol. 2001, 18, 1–10. [Google Scholar] [CrossRef]
  54. Norton, D.M.; McCamey, M.A.; Gall, K.L.; Scarlett, J.M.; Boor, K.J.; Wiedmann, M. Molecular Studies on the Ecology of Listeria monocytogenes in the Smoked Fish Processing Industry. Appl. Environ. Microbiol. 2001, 67, 198–205. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Laksanalamai, P.; Joseph, L.A.; Silk, B.J.; Burall, L.S.; LTarr, C.; Gerner-Smidt, P.; Datta, A.R. Genomic characterization of Listeria monocytogenes strains involved in a multistate listeriosis out-break associated with cantaloupe in US. PLoS ONE 2012, 7, e42448. [Google Scholar] [CrossRef]
  56. Barco, L.; Belluco, S.; Roccato, A.; Ricci, A. A systematic review of studies on Escherichia coli and Enterobacteriaceae on beef carcasses at the slaughter-house. Int. J. Food Microbiol. 2015, 207, 30–39. [Google Scholar] [CrossRef]
  57. Martínez-Avilés, M.; Garrido-Estepa, M.; Álvarez, J.; de la Torre, A. Salmonella Surveillance Systems in Swine and Humans in Spain: A Review. Veter Sci. 2019, 6, 20. [Google Scholar] [CrossRef] [Green Version]
  58. Self, J.L.; Luna-Gierke, R.E.; Fothergill, A.; Holt, K.G.; Vieira, A.R. Outbreaks attributed to pork in the United States, 1998–2015. Epidemiol. Infect. 2017, 145, 2980–2990. [Google Scholar] [CrossRef] [PubMed]
  59. Savariraj, W.R.; Ravindran, N.B.; Kannan, P.; Paramasivam, R.; Senthilkumar, T.M.A.; Kumarasamy, P.; Rao, V.A. Prevalence, antimicrobial susceptibility and virulence genes of Staphylococcus aureus isolated from pork meat in retail outlets in India. J. Food Saf. 2019, 39, e12589. [Google Scholar] [CrossRef]
Figure 1. Non-Metric Multidimensional Scaling Graph (NMDS) of the microbiota structures of the six visit samples using Bray–Curtis index. Each point represents a sample, and each visit is represented by the combination of a symbol and a color. Blue square: Visit one; Orange circle: Visit two; Pink triangle: Visit three; Turquoise star: Visit four; Yellow x: Visit five; Purple diamond: Visit six.
Figure 1. Non-Metric Multidimensional Scaling Graph (NMDS) of the microbiota structures of the six visit samples using Bray–Curtis index. Each point represents a sample, and each visit is represented by the combination of a symbol and a color. Blue square: Visit one; Orange circle: Visit two; Pink triangle: Visit three; Turquoise star: Visit four; Yellow x: Visit five; Purple diamond: Visit six.
Microorganisms 11 00133 g001
Figure 2. Non-Metric Multidimensional Scaling Graphs (NMDS) using Bray–Curtis index illustrating sample microbiome beta-diversity according to visit and paired two by two. Blue square: Visit one; Orange circle: Visit two; Pink triangle: Visit three; Turquoise star: Visit four; Yellow x: Visit five; Purple diamond: Visit six.
Figure 2. Non-Metric Multidimensional Scaling Graphs (NMDS) using Bray–Curtis index illustrating sample microbiome beta-diversity according to visit and paired two by two. Blue square: Visit one; Orange circle: Visit two; Pink triangle: Visit three; Turquoise star: Visit four; Yellow x: Visit five; Purple diamond: Visit six.
Microorganisms 11 00133 g002
Figure 3. Non-Metric Multidimensional Scaling Graph (NMDS) of the microbiota structures of the six production line samples using Bray–Curtis index. Each point represents a sample, and each production line is represented by the combination of a symbol and a color. Turquoise star: Boston production line (BO); Blue square: Main production line (CP); Purple diamond: Ham production line (FE); Orange circle: Belly production line (FL); Pink triangle: Loin production line (LO); Yellow x: Picnic production line (PI).
Figure 3. Non-Metric Multidimensional Scaling Graph (NMDS) of the microbiota structures of the six production line samples using Bray–Curtis index. Each point represents a sample, and each production line is represented by the combination of a symbol and a color. Turquoise star: Boston production line (BO); Blue square: Main production line (CP); Purple diamond: Ham production line (FE); Orange circle: Belly production line (FL); Pink triangle: Loin production line (LO); Yellow x: Picnic production line (PI).
Microorganisms 11 00133 g003
Figure 4. Non-Metric Multidimensional Scaling Graphs (NMDS) using Bray–Curtis index illustrating sample microbiome beta-diversity according to production line and paired two by two. Turquoise star: Boston production line (BO); Blue square: Main production line (CP); Purple diamond: Ham production line (FE); Orange circle: Belly production line (FL); Pink triangle: Loin production line (LO); Yellow x: Picnic production line (PI).
Figure 4. Non-Metric Multidimensional Scaling Graphs (NMDS) using Bray–Curtis index illustrating sample microbiome beta-diversity according to production line and paired two by two. Turquoise star: Boston production line (BO); Blue square: Main production line (CP); Purple diamond: Ham production line (FE); Orange circle: Belly production line (FL); Pink triangle: Loin production line (LO); Yellow x: Picnic production line (PI).
Microorganisms 11 00133 g004
Figure 5. Confusion matrix of the visit model. The x-axis line indicates the identity of the six visits and the y-axis line shows the prediction of a sample visit based on its microbiota. 1: Visit one; 2: Visit two; 3: Visit three; 4: Visit four; 5: Visit five; 6: Visit six.
Figure 5. Confusion matrix of the visit model. The x-axis line indicates the identity of the six visits and the y-axis line shows the prediction of a sample visit based on its microbiota. 1: Visit one; 2: Visit two; 3: Visit three; 4: Visit four; 5: Visit five; 6: Visit six.
Microorganisms 11 00133 g005
Figure 6. Confusion matrix of the production line model. The x-axis line indicates the identity of the six production lines and the y-axis line shows the prediction of a sample production line based on its microbiota. BO: Boston production line; CP: main production line; FE: Ham production line; FL: Belly production line; LO: Loin production line; PI: Picnic production line.
Figure 6. Confusion matrix of the production line model. The x-axis line indicates the identity of the six production lines and the y-axis line shows the prediction of a sample production line based on its microbiota. BO: Boston production line; CP: main production line; FE: Ham production line; FL: Belly production line; LO: Loin production line; PI: Picnic production line.
Microorganisms 11 00133 g006
Table 1. Bacterial genera with an average relative abundance equal to or greater than 5% according to the six visits independent of production line.
Table 1. Bacterial genera with an average relative abundance equal to or greater than 5% according to the six visits independent of production line.
GenusV1V2V3V4V5V6
%Nbr%Nbr %Nbr %Nbr %Nbr %Nbr
Fusobacterium16.664210.472414.393912.722112.973413.9434
Trueperella11.272313.172013.622611.981713.502713.8231
Pseudomonas10.311712.452311.12113.49328.882410.1520
Acinetobacter14.322211.562012.512212.22348.802012.1920
Peptoniphilus13.321515.29913.01910.36513.501110.9331
Rothia11.191312.81710.1257.4358.5669.365
Enhydrobacter10.96810.48915.927
Staphylococcus 23.06910.347
Bacteroides_unclassified12.12249.88117.541112.53511.45269.7912
Aerococcaceae_unclassified 6.77612.451212.581613.0012
Psychrobacter 7.0789.12510.4010
Porphyromonas11.35119.35810.84197.22712.45156.785
Parvimonas 5.95 7.5717
Clostridium_sensu_stricto 9.345 9.6469.588
Peptostreptococcus 16.696
Macrococcus8.988 18.589
Chryseobacterium 6.746 7.259
Epilithonimonas11.67
Prevotella 7.4310
Clostridiaceae_1_unclassified 6.645
Only the bacterial genera (OTUs) with a relative abundance equal to or greater than 5% in at least five samples from a visit were compiled. V1: Visit one; V2: Visit two; V3: Visit three; V4: Visit four; V5: Visit five; V6: Visit six. %: Average relative abundance of OTU in samples present at 5% or more. Nbr: Number of samples in which OTU is present at 5% or more. OTUs with an increase in mean relative abundance of at least 5% compared to other visits are identified in green.
Table 2. Bacterial genera with an average relative abundance equal to or greater than 5% according to the six production lines and independent of visit.
Table 2. Bacterial genera with an average relative abundance equal to or greater than 5% according to the six production lines and independent of visit.
GenusCPFLLOBOPIFE
%Nbr%Nbr%Nbr%Nbr%Nbr%Nbr
Fusobacterium14.632824.463317.242116.324423.483516.3633
Trueperella15.391924.283611.35612.174114.512710.4915
Pseudomonas9.17327.68108.03229.24229.341317.0438
Acinetobacter10.573210.341321.392214.752412.31713.530
Peptoniphilus10.19516.582010.46815.241519.991514.9717
Rothia10.2728 34.357 10.98116.925
Enhydrobacter 14.52915.74514510.825
Staphylococcus 16.3415
Bacteroides_unclassified8.131210.26209.89512.232015.53276.655
Aerococcaceae_unclassified10.381617.511514.8710 8.078
Psychrobacter10.2316 5.738
Porphyromonas8.5814 14.9269.37148.731010.1217
Parvimonas 7.066 7.936 6.4311
Clostridium_sensu_stricto10.4618 9.265
Peptostreptococcus 13.828
Macrococcus 11.9124
Streptococcus8.025
Chryseobacterium7.879
Epilithonimonas 11.438
Rhodopseudomonas5.989
Prevotella 7.946
Clostridiaceae_1_unclassified 6.625
Only the bacterial genera (OTUs) with a relative abundance equal to or greater than 5% in at least five samples from a production line were compiled. CP: Main production line; FL: Belly production line; LO: Loin production line; BO: Boston production line; PI: Picnic production line; FE: Ham production line. %: Average relative abundance of OTU in samples present at 5% or more. Nbr: Number of samples in which OTU is present at 5% or more. OTUs with an increase in mean relative abundance of at least 5% on one production line over any other lines are shown in green.
Table 3. Comparison of the alpha-diversity indices of the microbiota samples from the six visits.
Table 3. Comparison of the alpha-diversity indices of the microbiota samples from the six visits.
ObservedShannonInv. Simpson
V1193 *V3, V4, V5, V62.97 *V2, V3, V4, V59.85 *V2, V4, V5
V22183.28 *V1, V613.64 *V1, V6
V3221 *V13.16 *V1, V4, V511.17 *V4, V5
V4239 *V13.39 *V1, V3, V614.51 *V1, V3, V6
V5245 *V13.47 *V1, V3, V616.5 *V1, V3, V6
V6226 *V13.07 *V2, V4, V510.35 *V2, V4, V5
All visits p-values0.08380.0008 **0.0003 **
Diversity indices were calculated with 1000 iterations based on 10,087 sequences. Each line in the graph corresponds to the alpha diversity measures associated with a visit. Pairs of visits showing a significant difference of alpha diversity according to Student’s t-tests are indicated with superscript after the symbol * at the appropriate line. For example, for the Observed Index, visit one is significantly different from visits three, four, five, and six. V1: Visit one; V2: Visit two; V3: Visit three; V4: Visit four; V5: Visit five; V6: Visit 6. ** Significant differences with the Kruskal-Wallis test.
Table 4. Comparison of alpha-diversity indices of the microbiota samples from the six production lines.
Table 4. Comparison of alpha-diversity indices of the microbiota samples from the six production lines.
ObservedShannonInv. Simpson
CP372 *LO, FL, BO, PI, FE3.73 *LO, FL, BO, PI, FE18.5 *LO, FL, BO, PI, FE
LO169 *CP, PI, FE3.22 *CP, FL, BO11.87 *CP, FL
FL180 *CP, FE2.95 *CP, LO, FE9.18 *CP, LO, FE
BO167 *CP, PI, FE3.03 *CP, LO, FE10.46 *CP, FE
PI196 *CP, LO, BO, FE3.06 *CP, FE11.51 *CP
FE237 *CP, LO, FL, BO, PI3.29 *CP, FL, BO, PI13.52 *CP, FL, BO
All production lines p-values2.2 × 10−16 **4.64 × 10−10 **1.32 × 10−8 **
Diversity indices were calculated with 1000 iterations based on 10,087 sequences. Each line in the graph corresponds to the alpha diversity measures associated with a production line. Pairs of production lines showing a significant difference of alpha diversity according to Student’s t-tests are indicated with superscript after the symbol * at the appropriate line. For example, for the observed index, the FL production line is significantly different from the CP and FE production lines. CP: Main production line; FL: Belly production line; LO: Loin production line; BO: Boston production line; PI: Picnic production line; FE: Ham production line. ** Significant differences with the Kruskal-Wallis test.
Table 5. Relevant bacteria (i.e., commonly known bacterial genera that can include potential spoilage species or potential foodborne pathogen species) extracted from the list of OTUs identified using the MaAsLin test and associated with the different visits.
Table 5. Relevant bacteria (i.e., commonly known bacterial genera that can include potential spoilage species or potential foodborne pathogen species) extracted from the list of OTUs identified using the MaAsLin test and associated with the different visits.
Positive AssociationsNegative Bacterial Associations
TotalUniqueRelevant Bacterial GeneraTotalUniqueRelevant Bacterial Genera
V100 00
V2144Listeria_80__Otu00067201Escherichia_Shigella_Otu00046
Brochothrix_Otu00245
Streptococcus_Otu00029
V33910Acinetobacter_Otu00084202Lactococcus_Otu00132
Flavobacterium_Otu00354Streptococcus_Otu00029
Pseudomonas_Otu00173
Psychrobacter_Otu00047
V43911Acinetobacter_Otu00084231Escherichia_Shigella_Otu00046
Enterococcus_77__Otu00071Aerococcus_Otu00053
Flavobacterium_Otu00138Lactococcus_Otu00132
Lactobacillus_Otu00092
Psychrobacter_Otu00047
Streptococcus_Otu00024
V54912Campylobacter_Otu00197221Escherichia_Shigella_Otu00046
Brochothrix_Otu00245Aerococcus_Otu00053
Enterococcus_77__Otu00071Lactococcus_Otu00132
Flavobacterium_Otu00354Streptococcus_Otu00029
Flavobacterium_Otu00138
Flavobacterium_Otu00160
Lactobacillus_Otu00092
Pseudomonas_Otu00173
Psychrobacter_Otu00047
V6338Campylobacter_Otu00197244Escherichia_Shigella_Otu00046
Campylobacter_Otu00233Brochothrix_Otu00245
Campylobacter_Otu00024Carnobacterium_58__Otu00055
Campylobacter_Otu00233Pseudomonas_Otu00037
Flavobacterium_Otu00160Streptococcus_Otu00029
Flavobacterium_Otu00354
Flavobacterium_Otu00138
Lactobacillus_Otu00092
Psychrobacter_Otu00047
Column: 1—Visit ID; 2—Total number of OTUs positively associated with one or more visits; 3—Number of OTUs positively associated with only one particular visit; 4—Potential spoilage or pathogen OTUs positively associated with a visit (potential pathogen bacteria are identified in bold); 5—Total number of OTUs negatively associated with one or more visits; 6—Number of OTUs negatively associated with only one particular visit; 7—Potential spoilage or pathogen OTUs negatively associated with one or more visits (potential pathogen bacteria are identified in bold). V1: Visit one; V2: Visit two; V3: Visit three; V4: Visit four; V5: Visit five; V6: Visit 6.
Table 6. Relevant bacteria (i.e., commonly known bacterial genera that can include potential spoilage species or potential foodborne pathogen species) extracted from the list of OTUs identified using the MaAsLin test and associated with the different production lines.
Table 6. Relevant bacteria (i.e., commonly known bacterial genera that can include potential spoilage species or potential foodborne pathogen species) extracted from the list of OTUs identified using the MaAsLin test and associated with the different production lines.
Positive AssociationsNegative Associations
TotalUniqueRelevant Bacterial GeneraTotalUniqueRelevant Bacterial Genera
CP16997Campylobacter_Otu00493132Acinetobacter_Otu00020
Clostridium_sensu_stricto_Otu00016
Clostridium_sensu_stricto_Otu00136
Clostridium_sensu_stricto_Otu00307
Clostridium_sensu_stricto_Otu00441
Clostridium_sensu_stricto_Otu00059
Clostridium_XlVa_75__Otu00385
Enterococcus_77__Otu00071
Flavobacterium_Otu00480
Flavobacterium_Otu00333
Flavobacterium_Otu00160
Lactobacillus_Otu00092, Otu00099
Lactococcus_Otu00262, Otu00132
Moraxella_Otu00259
Pseudomonas_Otu00179, Otu00003
Psychrobacter_Otu00047
Psychrobacter_Otu00012
Psychrobacter_Otu00093
FL4214Clostridium_sensu_stricto_Otu00016422Acinetobacter_Otu00011, Otu0020
Pseudomonas_Otu00003, Otu00037
Psychrobacter_Otu00012, Otu00047
Moraxella_Otu00045
LO3611Clostridium_sensu_stricto_Otu00016404Acinetobacter_Otu00020, Otu00011
Enterococcus_77__Otu00071Moraxella_Otu00045
Flavobacterium_Otu00138
Staphylococcus_99__Otu00008
BO00 00
PI295Clostridium_sensu_stricto_97__Otu00203165Salmonella_66__Otu00026
Psychrobacter_Otu00047Acinetobacter_Otu00020, Otu00011
Pseudomonas_Otu00023
FE458Clostridium_sensu_stricto_Otu00059210Moraxella_Otu00045
Clostridium_sensu_stricto_Otu00016Morganella_91__Otu00337
Clostridium_sensu_stricto_Otu00136
Campylobacter_Otu00493
Escherichia_Shigella_Otu00046
Carnobacterium_58__Otu00055
Lactococcus_Otu00262
Psychrobacter_Otu00012
Column: 1—Production lines ID; 2—Total number of OTUs positively associated with one or more production lines; 3—Numbers of OTUs positively associated with only one particular production line; 4—Potential spoilage or pathogen OTUs positively associated with a production line (potential pathogen bacteria are identified in bold); 5—Total number of OTUs negatively associated with one or more production lines; 6—Number of OTUs negatively associated with only one particular production line; 7—Potential spoilage or pathogen OTUs negatively associated with a production line (potential pathogen bacteria are identified in bold); CP: Main production line; FL: Belly production line; LO: Loin production line; BO: Boston production line; PI: Picnic production line; FE: Ham production line.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Shedleur-Bourguignon, F.; Duchemin, T.; P. Thériault, W.; Longpré, J.; Thibodeau, A.; Hocine, M.N.; Fravalo, P. Distinct Microbiotas Are Associated with Different Production Lines in the Cutting Room of a Swine Slaughterhouse. Microorganisms 2023, 11, 133. https://doi.org/10.3390/microorganisms11010133

AMA Style

Shedleur-Bourguignon F, Duchemin T, P. Thériault W, Longpré J, Thibodeau A, Hocine MN, Fravalo P. Distinct Microbiotas Are Associated with Different Production Lines in the Cutting Room of a Swine Slaughterhouse. Microorganisms. 2023; 11(1):133. https://doi.org/10.3390/microorganisms11010133

Chicago/Turabian Style

Shedleur-Bourguignon, Fanie, Tom Duchemin, William P. Thériault, Jessie Longpré, Alexandre Thibodeau, Mounia N. Hocine, and Philippe Fravalo. 2023. "Distinct Microbiotas Are Associated with Different Production Lines in the Cutting Room of a Swine Slaughterhouse" Microorganisms 11, no. 1: 133. https://doi.org/10.3390/microorganisms11010133

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