# Lactic Starter Dose Shapes S. aureus and STEC O26:H11 Growth, and Bacterial Community Patterns in Raw Milk Uncooked Pressed Cheeses

Justine Piqueras, Christophe Chassard, Cécile Callon, Etienne Rifa, Sébastien

Theil, Annick Lebecque, Céline Delbès

► To cite this version:

Justine Piqueras, Christophe Chassard, Cécile Callon, Etienne Rifa, Sébastien Theil, et al.. Lactic Starter Dose Shapes S. aureus and STEC O26:H11 Growth, and Bacterial Community Patterns in Raw Milk Uncooked Pressed Cheeses. Microorganisms, 2021, 9 (5), pp.1081. 10.3390/microorganisms9051081. hal-03536629

HAL Id: hal-03536629

https://hal.inrae.fr/hal-03536629

Submitted on 20 Jan 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

# Lactic Starter Dose Shapes *S. aureus* and STEC O26:H11 Growth, and Bacterial Community Patterns in Raw Milk Uncooked Pressed Cheeses

Justine Piqueras<fnref n="1" />, Christophe Chassard<fnref n="1" />, Cécile Callon<fnref n="1" />, Etienne Rifa<fnref n="2" />, Sébastien Theil<fnref n="1" /> ID, Annick Lebecque<fnref n="1" /> and Céline Delbès <fnref n="1" />,*, ID

<figure><img src="image_3.png" /></figure>

check for
updates

**Citation:** Piqueras, J.; Chassard, C.; Callon, C.; Rifa, E.; Theil, S.; Lebecque, A.; Delbès, C. Lactic Starter Dose Shapes *S. aureus* and STEC O26:H11 Growth, and Bacterial Community Patterns in Raw Milk Uncooked Pressed Cheeses.
*Microorganisms* **2021**, *9*, 1081.
https://doi.org/10.3390/microorganisms9051081

Academic Editor: Salam A. Ibrahim

Received: 13 April 2021

Accepted: 14 May 2021

Published: 18 May 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_4.png" /></figure>

**Copyright:** © 2021 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Abstract:** Adding massive amounts of lactic starters to raw milk to manage the sanitary risk in the cheese-making process could be detrimental to microbial diversity. Adjusting the amount of the lactic starter used could be a key to manage these adverse impacts. In uncooked pressed cheeses, we investigated the impacts of varying the doses of a lactic starter (the recommended one, 1×, a 0.1× lower and a 2× higher) on acidification, growth of *Staphylococcus aureus* SA15 and Shiga-toxin-producing *Escherichia coli* (STEC) O26:H11 F43368, as well as on the bacterial community patterns. We observed a delayed acidification and an increase in the levels of pathogens with the 0.1× dose. This dose was associated with increased richness and evenness of cheese bacterial community and higher relative abundance of potential opportunistic bacteria or desirable species involved in cheese production. No effect of the increased lactic starter dose was observed. Given that sanitary criteria were paramount to our study, the increase in the pathogen levels observed at the 0.1× dose justified proscribing such a reduction in the tested cheese-making process. Despite this, the effects of adjusting the lactic starter dose on the balance of microbial populations of potential interest for cheese production deserve an in-depth evaluation.

**Keywords:** cheese; raw milk; lactic starter; STEC; CoPS; bacterial diversity; metagenomics

## 1. Introduction

An increasing part of consumers are looking for traditional products, with distinctive sensorial characteristics [1,2], such as raw milk cheeses. The sensorial characteristics of raw milk cheeses are strongly related to the richness of microbiota [3–7], itself dependent on the various initial sources of enrichment [1,5,8,9].

Apart from the variation of sensorial characteristics [1], using raw milk for cheese manufacture can also lead to an increased risk of contamination by food spoilage or pathogenic bacteria [10]. Unwanted bacteria can be responsible for the development of unpleasant flavors or visual defects like slits [11], causing an economic loss for the dairy industry. They can also be involved in sanitary risks. Because of their role in Collective Food Toxi-Infection (CFTI), major pathogens, such as *Salmonella* spp., *Listeria monocytogenes*, Shiga-toxin-producing *Escherichia coli* (STEC) and Coagulase Positive *Staphylococcus* (CoPS), are targeted by controls in food including raw milk cheeses [12].

Worldwide, *Staphylococcus aureus* is the most prevalent pathogen responsible for foodborne intoxication. In France, in 2018, it was the most frequently detected or suspected agent responsible for CFTI and represented 3% of CFTI related to the consumption of dairy products [13], with STEC representing 2.5% of the sources of CFTI [13]. STEC can cause severe hemolytic uremic syndrome, with O26:H11 being the most commonly found

<footnote n="1">¹ UMR 0545 Fromage, Université Clermont Auvergne, INRAE, VetAgro Sup, 20 Côte de Reyne, F-15000 Aurillac, France; justine.piqueras@inrae.fr (J.P.); christophe.chassard@inrae.fr (C.C.); cecile.callon@inrae.fr (C.C.); sebastien.theil@inrae.fr (S.T.); annick.lebecque@orange.fr (A.L.)</footnote>

<footnote n="2">² UMR INSA/INRA 792, Toulouse Biotechnology Institute, INSA/CNRS 5504, 135 Avenue de Rangueil, F-31077 Toulouse, France; etienne.rifa@insa-toulouse.fr</footnote>

<footnote>* Correspondence: celine.delbes@inrae.fr; Tel.: +33-471-456-419</footnote>
serotype in dairy products, including raw milk cheese [14]. Because of the high prevalence of *S. aureus* in dairy products, and the severity of the symptoms associated with STEC infections, these two pathogens are subjected to high vigilance.

Sanitary risk levels also vary according to the cheese type. Thus, due to their short ripening duration, high activity of water and also their slow or moderate acidification, uncooked pressed cheeses present an increased risk for the two previously mentioned pathogens. In these types of cheese, managing acidification by reaching targeted values of pH is a key to controlling sanitary risk [12,15,16] in that it allows for by allowing for an improvement in the protection against the growth of *S. aureus* and STEC.

Acidification is primarily driven by lactic acid bacteria (LAB) that produce lactic acid by lactose degradation. Thus, to ensure early acidification in the cheese-making processes, the milk can be inoculated with starter cultures composed of LAB [10].

Besides traditional inoculation practices, like back-slopping [5,17], multiple commercial solutions have been developed to offer starters adapted to each cheese type [7,10]. These lactic starters are composed of LAB that are generally recognized as safe (GRAS). They are selected for their cheese origin, their synergistic activities with other LAB, as well as their abilities to produce metabolites of interest, from antimicrobial compounds to the metabolites involved in sensorial characteristics. The addition of lactic starters allows for a more constant quality of cheese products [1,10,18], thus reducing the loss of products due to cheese defects. The use of commercial lactic starters has therefore become more and more widespread.

Microbial culture is relevant to evaluate the levels of the main cultivable microbial groups in cheeses and remains the basis of the standard methods for the determination of pathogens in foods (NF EN ISO 6888-2 [19], and ISO/TS 13136:2012 [20]). For more than twenty years, the development of culture-independent methods, based on the sequencing of universally conserved and hypervariable regions of 16S ribosomal DNA, has brought new knowledge about cheese microbial ecosystems, allowing for the rapid identification of the whole bacterial DNA from cheese environments [21]. Thus, such methods allow us to compare the microbial communities between different cheeses [22] or specify their spatial or temporal distribution [23]. Overall, these methods allow us to better understand the importance of the microbial diversity and dynamics in cheeses and to identify their main environmental sources and technological drivers [5]. Potentially negative impacts of the use of lactic starters on the microbial diversity in cheese products have been revealed. Indeed, the use of a lactic starter inoculated at high levels impacts the levels of autochthonous species [24,25]. During the cheese-making process, LAB species from a lactic starter become dominant [26]. When the starter used was produced by a natural method like back-slopping, the dominance of the LAB species identified in the natural lactic starter was also observed [27]. Reduced bacterial diversity, as well as lower richness and evenness, were observed in cheeses manufactured with a lactic starter, compared with cheeses manufactured without the addition of a lactic starter [28].

To our knowledge, only a few studies have looked into the impact of changing the inoculation level of the lactic starter. As an example, Fernandez-García et al. [29] and Khosrowshahi et al. [30] observed that, in Hispanico cheeses (semi-hard) and Iranian White cheeses (close-textured brined cheeses), respectively, the higher the dose of a mesophilic starter, the lower the pH values reached within the first 24 h of cheese-making. On the contrary, Özer et Kesenkaş, [31], did not observe any differences in the pH values in conjunction with the lactic starter dose on the first day in Mihalič cheese (a Turkish hard cheese stored in brine). Among these authors, only Fernandez-García et al. [29] looked at the impacts on the microbial characteristics from the variations in the dose of the lactic starter culture, but they did not observe any modification with regards to the lactic starter dose on mesophilic, Gram-negative bacteria or on lactobacilli counts.

Consequently, adjusting the lactic starter to the necessary and sufficient dose could be a key to improving cheese production, to managing sanitary risks while promoting microbial diversity.
In this regard, the principal aim of our study was to investigate first whether a reduced dose of lactic starter, compared with the recommended dose, would satisfy sanitary criteria while minimizing the effects on the microbial diversity. We also investigated whether the protection against pathogens could be improved using a higher lactic starter dose, and at what cost for the bacterial diversity. Taking the dose recommended by the lactic starter manufacturer as a reference, we also simultaneously tested a 10-fold lower (0.1×) and a 2-fold higher (2×) dose. These were added into raw milk originating from two different farms. The effects of these three different doses on the growth of *S. aureus* SA15 and STEC O26:H11 F43368, as well as on the microbial community composition, were assessed in uncooked pressed cheese models.

## 2. Materials and Methods

### 2.1. Strains and Culture Conditions

A unique batch of the lyophilized commercial culture MY800 (*Streptococcus thermophilus*, *Lactobacillus delbrueckii* spp. *bulgaricus*, *Lb. delbrueckii* subsp. *lactis*, Danisco, Copenhagen, Denmark) was used as a lactic starter culture for all cheese manufactures. The lyophilized lactic starter was stored at -20 °C before use. The recommendation for use is 5 Danisco Culture Units (DCU) for each 100 L of milk. The lyophilized lactic starter was revived in sterile reconstituted skimmed milk (LACTALIS Ingredients, Bourgbarré, France) for 2 h 30 min at 42 °C, and then was stored at 4 °C until use within 3 days.

STEC O26:H11 strain F43368 and coagulase (+) *S. aureus* strain SA15, isolated from "tome de pays" in 2009, and raw milk, respectively, were provided by UMRF (Unité Mixte de Recherche sur le Fromage, INRAE, Aurillac, France). Both pathogenic strains were revived in Brain Heart Infusion (BHI) broth (BIOKAR Diagnostics, Pantin, France), for 24 h at 37 °C. The strain of STEC O26:H11 was enumerated on BHI agar medium. The *S. aureus* strain was enumerated on Baird-Parker medium, supplemented with Rabbit Plasma Fibrinogen (RPF medium) (BIOKAR Diagnostics, Pantin, France) [32]. Aliquots of known concentrations were stored at -80 °C, until use in the cheese-making assays, 2 weeks after storage for the last. Inoculated cheeses samples on D1 were tested for the presence of staphylococcal enterotoxins by the Enzyme Linked Fluorescent Assay based Vidas Set 2 system (AgroLab's, Aurillac, France), which targets staphylococcal enterotoxins A to E with a limit of detection of 0.25 ng toxin per gram of food, according to NF EN ISO 19020:2017 [33].

### 2.2. Experimental Design

In January 2019, raw milk was obtained from two different farms ("F_15" and "F_38"), located no more than 20 km from the cheese-making plant. It was collected early in the morning from a refrigerated milk tank on the farm and used within 1 h to produce Saint-Nectaire-like uncooked pressed cheeses of reduced size (500 g) at the UMRF experimental cheese plant.

On each day of the cheese-making process, the milk from both farms was tested under the same conditions of inoculation, and each pair was repeated three times, over a 2-week period. Two cheeses were made from each batch of milk: a "control" cheese, without pathogen inoculation, and an "assay" cheese, with controlled inoculation of pathogens.

Pathogens were inoculated into the milk prior to the cheese-making process. The strains *S. aureus* SA15 and STEC O26:H11 F43368 were inoculated at the expected levels of 10³ CFU/mL and 10² CFU/mL, respectively, together in the assay milk vats. In control vats, no pathogen was added.

Both the assay and the control vats were inoculated either with the recommended dose of lactic starter MY800 (6 × 10⁶ CFU/g, referred to as "1×"), with a 10-fold lower dose (6 × 10⁵ CFU/g, referred to as "0.1×") or with a 2-fold higher dose (1 × 10⁷ CFU/g, referred to as "2×").
## 2.3. Cheese-Making and Sampling

Vats were filled with 5 L of raw milk, and all the vats were heated at 33 °C. The lactic starter (MY800) was inoculated at the chosen dose in each vat and then these vats were inoculated with pathogens. After being mixed for 15 min, the milk was inoculated with the ripening starter (30 mL per 500 L of milk, composed of *Brevibacterium linens*, *Debaryomyces hansenii*, *Penicillium fuscoglaucum*, and *Fusarium domesticum* (SN MEL N27 500 L, LIP Aurillac)), and rennet (40 mL per 100 L of milk, containing 520 mg/L of active chymosin (Fabre LCP Food Ingredients, Prat, France)). The inoculated milk was processed according to the uncooked pressed cheese technology as previously described [34,35] and salted at 5% with coarse salt. The cheeses were then ripened in UMRF's ripening cellars at 8 °C and 95% relative humidity for 28 days. During the ripening, the cheeses were washed at 6 and 9 days of ripening with the solution of the ripening starter. They were turned over once a week.

The raw milk was sampled before the cheese was made to evaluate the bacterial and chemical composition and again after inoculation to control the dose of the lactic starter and the level of pathogens. Cheese cores were sampled after 1 (D1), and 8 (D8) days. After 28 days (D28), a last sampling was realized during which the rind (with a thickness of 1.5 mm) and the core from the same cheese were separately sampled. At each sampling, a 100 g-piece of cheese was taken and separated into equivalent portions intended for microbiological, metagenomic and biochemical analyses. The sample dedicated to the metagenomic analyses was roughly sliced with a knife before storage at -20 °C.

## 2.4. Physico-Chemical Analyses

Several physico-chemical parameters were determined in the fresh raw milk (Agro-Lab's, Aurillac, France). Urea, fat, protein, lactose and casein levels, as well as fatty acid profiles, were estimated by infra-red-light spectrometry. Somatic cell counts were estimated by epifluorescence and the lipolysis rate by Copper Soap method.

pH measures were carried out with a 926 VTV pH-meter with Ingold electrode 406 MX (Mettler-Toledo S.A., Viroflay, France), in the milk, on the fresh cheese at 6 h (D0.25) and at the time of sampling (D1 and D28). They were taken from three different points on each sampled cheese, at room temperature.

The dry matter content was calculated in the milk, in one-day-old cheese (D1, at the entry in the cellar) and at the end of ripening, according to ISO 5534:2004 [IDF 4:2004] [36].

The rheological method used was uniaxial compression at a constant displacement rate, as described by Frétin et al. [37]. Resistance to penetration measures of the cheese core piece were taken on seven points uniformly distributed over the surface of the control cheese sample. The Young's modulus (stress/relative strain), proportional to cheese firmness and inversely related to cheese elasticity and creaminess [38], was calculated.

## 2.5. Microbiological Analyses

Ten grams of cheese were ground in phosphate buffer solution using a homogenizer mill (BagMixer, Interscience, Saint Nom la Brétèche, France), as described by Frétin et al. [37]. To evaluate the levels of microbial groups, both the raw and inoculated milk, as well as the ground cheese, were plated directly without or with dilution in Ringer's solution. Inoculations were carried out with EasySpiral (Interscience) in exponential mode, with serial dilutions in Ringer's solution.

In the fresh raw milk, the levels of *E. coli* and total microorganisms were evaluated by AgroLab's, by counting on TBX selective media and by flow cytometry (Bactoscan FC, FOSS, Hilleroed, Denmark), respectively, while the level of autochthonous coagulase positive *Staphylococcus* was estimated by counting on an RPF medium, and the absence of autochthonous STEC was verified by enrichment in peptone water followed by plating on ChromID coli at 42 °C (Biomerieux, Marcy l'Etoile, France).

After revival in sterile reconstituted skimmed milk, the lactic starter was enumerated on M17 agar [39] at 42 °C, for 24 h to determine the volume needed to inoculate the milk
at the desired dose during the manufacturing of the cheese. Both batches of inoculated milk were then verified for the levels of *Streptococcus*, STEC O26:H11 and *S. aureus*, on M17 agar at 42 °C, ChromID coli at 42 °C, and RPF medium at 37 °C, respectively. The levels of the pathogens were controlled in the cheese cores, at D1 and D28, by plating on the same media. As well, the level of thermophilic streptococci was also checked at D28 in the cheese cores and in the rinds. The levels of the other microbial groups such as lactic and ripening bacteria, yeasts and molds were evaluated by counting on MRS agar [40] at 30 °C, CRBM agar [41] at 25 °C, and OGA medium [42] at 25 °C, respectively, at D28 in the core and the rind of the cheeses.

## 2.6. Amplicon Sequencing

The milk (100 mL) was centrifuged at 5285 × g for 30 min, at 4 °C (Sigma 6K15, Kana-
gawa, Japan). The fat and the supernatant were discarded, and the pellet was resuspended
with 1mL of PBS. After a 5-min centrifugation at 13,000 × g (Eppendorf 5415D, Hamburg,
Germany), the supernatant was removed, and the pellet was stored at -20 °C until extrac-
tion. A thawed, two-gram sample of the control cheese was ground with a homogenizer
mill (BagMixer, Interscience, Saint Nom la Brétèche, France) with 2 mL of 4 M guanidine
thiocyanate -0.1 M Tris (pH 7.5) and 500 µL of the mixture were transferred into a 2 mL
tube, on which an extraction was performed.

Extractions were performed with a series of 24 sample tubes including one blank extraction tube, in which ground cheese or the milk pellet were replaced by PCR grade water (Merck KGaA, Darmstadt, Germany). 40 µL of 10% N-lauroylsarcosine and 200 mg of zirconium beads were added to each tube and then mixed with a bead mill homogenizer (Precellys Evolution, Bertin Technologies SAS, Ozyme, France) (one 20-s run at a speed of 6500 m/s). The following steps were the same as described by Frétin et al. [34], except the final precipitation which was performed as described by Duthoit et al. [43].

The DNA was quantified using the Qubit Fluorometric Quantitation (Invitrogen, Thermo Fisher Scientific, Waltham, MA, USA) method with a High Sensitivity kit for DNA extracts from milk and a Broad Range kit for DNA extracts from cheeses (rind and core).

The variable region V3–V4 of the 16S rRNA gene (~510 bp) was targeted using the primers MSQ-16SV3F (5’-TACGGRAGGCWGCAG-3’) [44] and PCR1R-460 (5’-TTACCAGG GTATCTAATCCT-3’) and amplified as previously described by Frétin et al. [8]. The amplified products were used for Illumina paired-end library preparation and cluster generation, followed by 250 bp paired-end sequencing on an Illumina MiSeq instrument (INRAE, GeT-PLaGE platform, Toulouse, France).

The raw data, provided by sequencing, were processed with the pipeline rANOMALY [45], which relies on the dada2 R package [46] to produce amplicon sequence variants (ASV) as taxonomic units. A step of decontamination was carried out, based both on the prevalence of contaminant ASVs, as identified in the blank samples, and on the DNA concentration, as described by Theil et Rifa [45]. The filtered ASV count table was used to perform statistical analyses.

## 2.7. Statistical Analyses

Significance was declared at *p* < 0.05. *p*-values between 0.05 and 0.10 were considered a tendency, while *p* > 0.10 was considered not significant. Statistical analyses were performed with R for Windows (version 3.6.3, R Foundation for Statistical Computing, Vienna, Austria). The packages "dunn.test", "ggplot2" were used for pairwise comparison of the lactic starter dose. Pathogen counts and pH values were analyzed with a Dunn's test, allowing multiple comparisons using rank sum, with a Bonferroni's correction. The farm effect was tested by a Wilcoxon's test. The absence of interaction with the farm was confirmed with an ANOVA, where the variables, lactic starter dose, farm and their interactions, were the fixed factors.

To study how modulating the dose of the lactic starter impacts the bacterial commu-
nities in raw milk-cheeses, the rANOMALY interface was used. To compare evenness
and richness (α-diversity) between tested conditions, different indexes were used, such as
Shannon, Simpson, and InvSimpson, which highlight evenness, as opposed to Observed, Chao1, ACE or Fisher indexes, which are more related to richness. Different distance methods were used to estimate the community composition differences (β-diversity) between two conditions, as Bray Curtis, Jaccard and weighted Unifrac. The farm and the lactic starter dose modulation effects on alpha-diversity and beta-diversity were tested by PERMANOVA and Adonis' tests, respectively. The differences of dispersion of bacterial profiles to the centroid between conditions were tested by TukeyHSD test. Differential analyses of families, genus and species abundances in samples, according to the farm or the dose, were performed on taxa with a minimum relative abundance of 0.001 and a minimum prevalence of 4 in the milk and of 3 in other samples. Three methods were used: DeSeq, MetaGenomeSeq and MetaCoder. A taxon was considered as differentially abundant between tested conditions when it was differentially found with a minimum of two methods (*p*-value < 0.05).

### 3. Results

#### 3.1. *Quality of Raw Milks*

The major biochemical characteristics of the raw milk, originating from each of the two farms and used to make the cheese, are summarized in the Supplementary Table S1. None of the measured descriptors (pH, levels of lactose, fat, protein, urea, casein and the number of somatic cells) showed a significant difference according to the farm they originated from.

The level of total microorganisms tends to be higher in the milk from the farm F_15 (1.30 × 10⁴ CFU/mL) than in the milk from the farm F_38 (6.33 × 10³ CFU/mL) (Table S1). The mean levels of autochthonous coagulase positive *Staphylococcus* were 1.04 × 10² and 6.16 × 10² CFU/mL in the raw milk from the farms F_15 and F_38, respectively, with no significant difference between the farms. The levels of total *E. coli* were lower than 10 CFU/mL whichever the farm, and no STECs were detected considering the detection threshold in the enrichment broth (<5 CFU/mL).

#### 3.2. *Levels of Pathogens*

After the inoculation of the pathogens into the milks, mean levels of *S. aureus* were 9.30 × 10² and 6.37 × 10² CFU/mL in the milk from the F_15 and the F_38, respectively. The levels of *S. aureus* in most of the control cheeses were either below the detection threshold, or below those in the corresponding inoculated cheeses whatever the ripening stage. The level of *S. aureus* in the inoculated cheeses increased over the first 24 h of cheese-making (Figure 1), to reach mean levels of 1.95 × 10⁶, 3.09 × 10⁵ and 1.25 × 10⁵ CFU/g for the lactic starter doses 0.1×, 1× and 2×, respectively. At D1, no significant difference in *S. aureus* level was found in relation to the lactic starter dose, whether the farms were treated together or separately (Table S2). In addition, none of the staphylococcal enterotoxin A to E was detected whatever the lactic starter dose (data not shown). At the end of ripening (D28), the levels of *S. aureus* in the inoculated cheeses were significantly lower than at D1, reaching 1.53 × 10⁵, 2.28 × 10⁴, and 2.33 × 10⁴ CFU/g for the lactic starter doses 0.1×, 1× and 2×, respectively (Table S2). The level of *S. aureus* was significantly higher with the 0.1× dose than with both of the other tested doses (Figure 1), with mean differences of 0.8 Log [CFU/g] compared with both the 1× (*p*-value = 0.004) and the 2× doses (*p*-value = 0.024). No significant difference was reported either with respect to the farm of origin, the ripening time or the dose interaction. Whatever the ripening time, the levels of *S. aureus* did not differ between the 1× and 2× lactic starter doses.

After the inoculation of the pathogens into the milk, mean levels of STEC O26:H11 were 1.04 × 10² and 1.13 × 10² CFU/mL in the milk from the F_15 and the F_38, respectively. In the control cheeses, the STEC O26:H11 level was below the detection threshold (Table S2), whatever the lactic starter dose, the farm or the ripening time. In the inoculated cheeses, over the first 24 h of the cheese-making process, the level of STEC O26:H11 increased to reach mean levels of 2.75 × 10⁷, 7.75 × 10⁵ and 1.96 × 10⁵ CFU/g for the lactic starter doses, 0.1×, 1× and 2×, respectively (Table S2). At D1, the STEC O26:H11 level was higher
with the 0.1× dose than with both of the other doses (Figure 2), with mean differences of 1.8 Log [CFU/g] compared with the 1× dose (p-value = 0.060) and of 2.3 Log [CFU/g] compared with the 2× dose (p-value = 0.001).

<figure><img src="image_5.png" /></figure>

**Figure 1.** Level in Log [CFU/g] of *S. aureus* in milk (D0) and in cheeses at D1 and D28, according to the lactic starter dose (n = 6). The differences according to the dose were estimated by a multiple comparison test, a Dunn's test. \* *p* < 0.05; \*\* *p* < 0.01; n.s. non-significant. For milk samples (n = 6), no statistical tests were performed.

<figure><img src="image_6.png" /></figure>

**Figure 2.** Level in Log [CFU/g] of STEC O26:H11 in milk (D0) and in cheeses at D1 and D28, according to the lactic starter dose (n = 6). The differences according to the dose were estimated by a multiple comparison test, a Dunn's test. \*\* *p* < 0.01; \* tendency; n.s. non-significant. For milk samples (n = 6), no statistical tests were performed.

The STEC O26:H11 levels in the cheeses at the end of ripening (D28) remained stable compared with those in the cheeses at D1, with mean values of 3.44 × 10⁷, 4.20 × 10⁵ and 1.10 × 10⁵ CFU/g for the 0.1×, 1× and 2× doses, respectively (Table S2). On D28, the STEC O26:H11 level was higher with the 0.1× dose than with both of the others, with mean difference of 2 Log [CFU/g] compared with the 1× dose (p-value = 0.069) and of 2.5 Log [CFU/g] compared to the 2× dose (p-value = 0.001). Whatever the ripening time, no difference in the levels of STEC O26:H11 was found either between the 1× and 2× lactic starter doses, or between farms.

### 3.3. Physico-Chemical Parameters

The pH values in the pathogen-inoculated cheeses and in the control cheeses were compared and no significant difference was reported, whatever the ripening time or the lactic starter dose (Table S2). The results detailed below (Figure 3) correspond to the inoculated cheeses. Whatever the lactic starter dose, a significant decrease in the pH values
was observed in the milk (with a mean value of 6.66) at D0.25, where mean values reached 6.53, 5.93 and 5.74, for 0.1×, 1× and 2× doses, respectively. Between D0.25 and D1, pH values dropped to 5.47, 5.23 and 5.19 for 0.1×, 1× and 2× doses, respectively. The pH value in the cheeses with the 0.1× dose dropped further to 5.24 at D28. The acidification was delayed with the 0.1× dose compared with both of the other doses at the early stages of the cheese-making process (D0.25 and D1). At D0.25, the pH value in the cheeses with the 0.1× dose exceeded those in the cheeses at 1× and 2× dose by 0.60 (p-value = 0.059) and 0.79 (p-value = 0.001), respectively. At D1, the pH value in the cheeses with the 0.1× dose exceeded those in the 1× and the 2× dose cheeses by 0.24 (p-value = 0.096) and by 0.28 (p-value = 0.006), respectively. At D28, the pH values did not significantly differ with respect to the lactic starter dose. However, the pH values in the cheeses manufactured with the milk from the F_15 were significantly lower than those in the cheeses manufactured with the milk from the F_38, with respective mean values of 5.21 and 5.37 (Table S2). Whatever the ripening time, no difference in the pH values was observed between the 1× and 2× lactic starter doses.

<figure><img src="image_7.png" /><figcaption>**Figure 3.** pH values in milk (D0) and in cheeses at D0.25, D1 and D28, according to the lactic starter dose (n = 6). The differences according to the dose were estimated by a multiple comparison test, a Dunn's test. **p** < 0.01; • tendency; n.s. non-significant. For milk samples (n = 6), no statistical tests were performed.</figcaption></figure>

The dry matter content in the pathogen-inoculated cheeses and in the control cheeses were compared and no significant difference was reported, whatever the ripening time or the lactic starter dose (Table S2). The results detailed below correspond to the inoculated cheeses. At D1, a lower dry matter content was observed with the 0.1× dose compared to the 2× dose, and then it remained stable between D1 and D28 in the 0.1×-dose cheeses, but significantly increased in the cheeses manufactured with 1× and 2× doses. No difference in dry matter content was found between either the 1× and 2× lactic starter doses, or in relation to the farm, whatever the ripening time.

The Young's modulus values (Table S2), obtained from rheological analyses, were significantly lower with the 0.1× dose than with the 2× dose of the lactic starter, which involved the higher the lactic starter dose, the firmer the cheese. Difference relating to the farm was also observed, and the mean Young's modulus value in cheeses manufactured with milk from the F_38 was significantly higher (mean value of 0.401 MPa) than in cheese manufactured with milk from the F_15 (mean value of 0.327 MPa).

### 3.4. *Microbial Counts*

At the end of ripening time (D28), the levels of thermophilic streptococci in the cores of the control cheeses (Table S3) were significantly higher in the cheeses manufactured with the 2× dose of the lactic starter than those with the 0.1× dose (p-value = 0.044). In the cheese cores, no significant difference was observed with respect to the farm. In the
rinds, no significant difference in the level of thermophilic streptococci was observed, with respect to the lactic starter dose, and mean levels reached 2.66 × 10⁸ CFU/g. However, higher levels of thermophilic streptococci were noticed in the cheeses manufactured with the milk from the F_38 than in those manufactured with the milk from the F_15.

The levels of the other microbial groups (Table S3) were similar whatever the lactic starter dose, except that the higher the lactic starter dose the lower the ripening bacteria level in the cheese cores. Nevertheless, several microbial groups showed different levels in the cheeses in relation to the farm: the levels of the yeasts and of ripening bacteria varied both in the cheese cores and in the cheese rinds, while the levels of molds and of lactic bacteria differed only in the cheese cores or in the cheese rinds, respectively.

### 3.5. Bacterial Diversity

The differences in the bacterial community evenness and diversity, according to the farm of origin and to the lactic starter dose, were evaluated. The bacterial community in the milk and in the cheese cores differed in richness depending on the farm they originated from, as indicated by the Observed, Chao1, ACE and Fisher indexes (Table S4). Milk from F_15 showed a higher richness than that from the F_38 but, during the ripening, for all lactic starter doses combined, the difference was reversed, with the richness index higher with the F_38 than with the F_15. Considering the data of the two farms combined, the most significant differences in the bacterial community richness and evenness were linked to the lactic starter dose, especially in the cheese cores ($p$-value ≤ 0.001). The values of the richness and the evenness indexes were increased when the lactic starter dose was decreased, whatever the index (Shannon, Simpson or InvSimpson) mainly in the cheese cores (D8 and D28) and, to a lesser extent, in the cheese rinds. These indexes, whatever the dose or the farm, also increased between D8 and D28 in the cheese cores.

The ordination of the bacterial community profiles through non-metric multidimensional scaling (NMDS) according to the farm or to the lactic starter dose was evaluated. The bacterial profiles in the raw milk (Table S5), and in the cheese cores at D28 (Figure 4D, Table S5), were clustered by farm, as shown by the Bray-Curtis, Jaccard and wUF methods. The variability in the bacterial composition of the milk (Figure S1) was higher in the milk originating from the farm F_38 than in that from the farm F_15, as shown by the difference in profiles dispersion to centroid. The bacterial profiles of the cheeses were clustered principally by lactic starter dose (Figure 4), as shown by the Bray-Curtis, Jaccard and wUF methods. We noted a cluster composed of the cheeses manufactured with the 0.1× dose, significantly distant from both of the other doses, in the cheese cores (at D8 and D28) and, to a lesser extent, in the cheese rinds. The bacterial profiles in the cores of the cheeses (at D8 and D28), manufactured with the 0.1× dose, were also more scattered than those manufactured with both of the other doses, but the dispersion to centroid for 1× and 2× doses increased during the cheese ripening and was higher in the cheese rinds than in the cheese cores. No difference was observed between 1× and 2× dose.

The *Streptococcaceae* family, *Streptococcus* genus and *S. thermophilus* species were the most predominant bacterial taxa in the cheese cores, both at D8 and D28, as well as in the cheese rinds. However, their relative abundance decreased during ripening and was lower in the rinds than in the cores. None of the DNA sequences detected in the raw milk or in the control cheeses was affiliated to *S. aureus*, or to *E. coli*. Differential analyses were performed on the relative abundance of bacterial taxa depending on the tested conditions in the raw milk, the cheese cores at D8 and D28 and the cheese rinds. For each ripening time, the family, genus and species ranks were investigated. Mean relative abundance of taxa are displayed in the Table S6. Only taxa differentially abundant between conditions are reported in the table. Whatever the taxonomic rank, most of the differences in relation to the farm were observed in the milk, and none was observed in the cheese rinds. In the raw milk samples, differentially abundant taxa were generally more abundant in the F_15 than in the F_38. Thus, the *Streptococcaceae* family, in this case mainly represented by the species *Lactococcus lactis*, were more abundant in the raw milk from the F_15, followed
by the *Lactobacillaceae*, *Bifidobacteriaceae* and *Corynebacteriaceae* families. Conversely, the higher mean relative abundance of species and genus belonging to the *Pseudomonadaceae* and *Bacillaceae* families and the *Streptococcus* genus were observed in the raw milk samples from the F_38. At D8, only the *Lactobacillaceae* family, genus and species still showed higher mean relative abundance in the F_15 cheeses. In addition, the *Staphylococcus* genus (identified as *S. xylosus*, a non-CoPS species) and *Enterococcus faecalis* presented higher mean relative abundance in the F_38. Finally, at D28, additional taxa belonging to *Enterococcaceae*, *Enterococcus*, *E. faecalis* and *Serratia* sp. showed higher mean relative abundance in the cores of the F_38 cheeses.

<figure><img src="image_8.png" /><figcaption>**Figure 4.** Beta-diversity of (a) Cheese core samples at D8 (n = 18), (b) Cheese core samples at D28 (n = 18), and (c) Cheese rind samples at D28 (n = 18) clustered according to the lactic starter dose; and (d) Cheese core samples at D28 (n = 18) clustered according to the farm. Non-Metric multidimensional scaling (NMDS) based on the Bray Curtis algorithm of bacterial communities. *p*-value was obtained after permutational MANOVA analysis (Adonis statistical test) and indicates significance between sample groups.</figcaption></figure>

Differences in the relative abundance of some families, genera and species were observed with respect to the lactic starter dose in the cheese cores, whatever the time of ripening, as well as in the cheese rinds. The relative abundance of the *Streptococcaceae* family, *Streptococcus* genus and *S. thermophilus* increased with the lactic starter dose in the cheese cores both at D8 and D28, with a maximum difference of 1.8-fold between the 0.1× and the 2× doses at D28. The same trend was observed in the cheese rinds, but with a higher difference in mean relative abundance between the 0.1× dose and both of the other doses (with multiplying factors of 10× and 16× for the 1× and 2× doses, respectively). The families, genera and species belonging to the *Enterobacteriaceae* and *Enterococcaceae* families and to *L. lactis* were differentially abundant in relation to the lactic starter dose at D8. Their mean relative abundance was significantly higher with the 0.1× dose than with both of the other doses. Then, at D28 in the cheese cores, the mean relative abundance of *Enterobacteriaceae* and *Enterococcaceae* families dropped to an equivalent abundance
regardless of the lactic starter dose. On the contrary, the mean relative abundance of *L. lactis*
increased between D8 and D28 and was still significantly higher with the 0.1× dose than
with both of the other doses. In the cheese rinds, the higher the dose, the lower the mean
relative abundance of *Brevibacterium* genus and species, of the *Enterobacteriaceae* and, to a
lesser extent of the *Staphylococcaceae* families, with significant differences between the 0.1×
and the 2× doses.

4. Discussion

Managing sanitary protection while promoting microbial diversity is a major challenge for raw milk cheese manufacture improvement. Several studies have reported differences in the microbial or physico-chemical characteristics of raw milk cheeses with or without added lactic starter. To our knowledge, few studies have evaluated the impacts on acidification of modulating the lactic starter dose from the recommended one, or on the level of pathogens or the microbial diversity in the cheese. We hypothesized that adjusting the dose of lactic starter used could be a key to managing sanitary protection while promoting bacterial diversity in raw milk cheeses. We compared the impacts of the lactic starter used as recommended by the manufacturer with those of a 10-fold lower and a 2-fold higher dose on the levels of *S. aureus* and STEC O26:H11 and on the microbial diversity in uncooked pressed cheese made from raw milk. We tested a high-risk situation, combining high levels of inoculated pathogens and a cheese process which could allow their growth. The levels of pathogens were defined according to previous studies in literature, allowing for comparisons between growth patterns: initial inoculation levels of 2 Log [CFU/mL] of STEC O26:H11 were used to describe its growth in uncooked pressed model cheeses [47,48]. Medved'ová et al. [49] studied the impact of lactic starter addition on growth of *S. aureus* and *E. coli* at 3 Log [CFU/mL] each, in pasta-filata cheeses.

We observed a delay in acidification with the 0.1× dose compared to both the 1× and the 2× doses, with differences of 0.60 after 6 h (D0.25) and of 0.24 after one day (D1). The 10-fold reduction in the lactic starter dose was associated with a less significant delay in acidification compared to that observed without a lactic starter added, as reported in the literature. Thus, Carafa et al. [24] reported pH values reduced by 0.39 to 0.53 at 24 h in Malga Traditional Mountain cheeses (an Italian semi-hard cheese manufactured with raw milk) prepared with an autochthonous thermophilic starter added at 5 × 10⁶ CFU/mL, compared with cheeses without the addition of the lactic starter. An even larger reduction in the pH was observed after 24 h (0.6 units) by Medved'ová et al. [49], between the cheeses manufactured with a thermophilic starter added at 10⁵ CFU/mL and the control pasta-filata cheeses without a lactic starter added.

In our study, no significant difference in pH values was observed between the cheeses
manufactured with the 1× or the 2× dose. These results are consistent with those reported
in the work of Khosrowshahi et al. [30], which studied the impact of three different doses
(1× as the reference, 2× and 4×) of a mesophilic lactic starter on acidification during the
manufacture of Iranian White Cheese. They observed reductions in pH values at renneting
by 0.08 and 0.28 between the 1× dose and the 2× or 4× doses, respectively.

Moreover, the recommended dose and the 2× dose allowed Saint-Nectaire-like cheeses
to reach targeted pH value at D0.25 (around 5.8, or below, at 6 h), which limits the growth
of CoPS at 24 h [15], conversely to the 0.1× dose. This result highlighted the importance
of the recommended dose as a minimal lactic starter dose to ensure correct acidification,
possibly contributing to the management of sanitary risks.

To our knowledge, as regards both pathogens evaluated in the present study, the
scientific literature has only reported their behavior in cheeses manufactured with or
without the addition of a lactic starter.

Although the initial level of STEC can influence the maximal level reached in cheese, a
similar growth pattern was reported by different authors. In our study, we observed a fast
increase in the STEC level by around 4 Log [CFU/g] during the first 24 h. This is similar
to the results reported by Callon et al. [48], Delbès-Paus et al. [47], or Frétin et al. [34], in
uncooked pressed cheese models with a standard dose of thermophilic starter (around 10⁷ CFU/mL) and initial STEC inoculation levels of either 0.5 or 2 Log [CFU/mL]. We observed a persistence until 28 days of ripening, whatever the tested lactic starter dose. This was also observed by the abovementioned authors. Yet some studies on different cheese types, with longer duration of ripening from 15 weeks to 7 months, showed a reduction in the level or even a non-detection of *E. coli* at the end of ripening [24,28]. At D1 in our study, decreased levels of STEC O26:H11 by 1.8 Log [CFU/g] and 2.3 Log [CFU/g] were observed with 1× and 2× doses compared with the 0.1× dose, respectively. The same trend was observed at D28, with levels of STEC O26:H11 decreased by 2 Log [CFU/g] and 2.5 Log [CFU/g], with the 1× and the 2× doses, respectively, compared with the 0.1× dose. The differences in STEC counts in cheese at D1 fall into the same range as those reported by Medved'ová et al. [49], showing a decrease in the levels of *E. coli* by 2.10 Log [CFU/g] at D1 in cheeses with the addition of a lactic starter in comparison with the control cheeses with no lactic starter added.

The growth patterns of *S. aureus* in our uncooked pressed cheese models were consistent with those previously reported. Delbès et al. [15] described a fast growth of *S. aureus* during the first 6 h of the cheese-making process (+3 Log [CFU/g]) and then a slight reduction during ripening (-0.55 Log [CFU/g] at D30). In addition, we showed that the *S. aureus* levels were increased by 0.8 Log [CFU/g] at D28 in the cheeses manufactured with the 0.1× dose compared with those manufactured with both of the other doses. The threshold of 10⁴ CFU/g for CoPS counts in cheeses products, as defined by the Commission Regulation (EC) No 2073/2005 of 15 November 2005 on microbiological criteria for foodstuffs (last update: 8 March 2020) [50], was consequently overpassed with the 0.1× lactic starter dose, which was not the case with the recommended dose. *S. aureus* SA15 produced staphylococcal enterotoxins in vitro in Brain Heart Infusion incubated at 30 °C for 48 h [51] and was found to possess the enterotoxin-encoding *sec4* and *sel2* genes [52]. Yet, no staphylococcal enterotoxin was detected in the inoculated cheeses on D1, whatever the lactic starter dose. These results suggest that neither the SA15 strain, nor the *S. aureus* populations naturally present in the raw milk used, produced enterotoxins in detectable concentrations in the tested cheeses, despite *S. aureus* levels reached 6 Log [CFU/g] and pH was around 6.5 at D0.25 at the 0.1× dose. These results were consistent with those from previous studies in Saint-Nectaire-model cheese made from cows' raw milk [15], similar to our cheese models, in which no staphylococcal enterotoxin was detected, although coagulase-positive staphylococci counts were around 5 Log CFU/g at D0.25, a propitious condition to SE accumulation. The increase in the *S. aureus* counts with the 10-fold reduced lactic starter dose are similar to those reported by Al-Nabulsi et al. [53] in white brined cheeses manufactured without a lactic starter. Using standard brine (10% of NaCl) and a storage temperature of 25 °C, the authors showed reductions by 1.17 Log [CFU/g] and by 0.56 Log [CFU/g] at D1 and at D28, respectively, between the cheeses manufactured without and with a lactic starter. Medved'ová et al. [49] also observed a reduction by 0.67 Log [CFU/g] of CoPS levels at D1 in pasta-filata cheeses with an added lactic starter compared with those without a lactic starter. At the end of storage (D28), they obtained a reduction of 2.04 Log [CFU/g] between the cheeses with and without a lactic starter.

Within the range of the lactic starter doses tested in the present work, the best protection against pathogens (as STEC O26:H11 or *S. aureus*) was obtained with the lactic starter added as recommended by the manufacturer. While we observed no difference, whatever the pathogen, between the 1× and 2× doses, a 10-fold lower dose did not guarantee protection against sanitary risks. This protecting effect of the lactic starter was associated with a faster acidification and lower pH values at D1, but it cannot be excluded that it could be related to antimicrobial properties of *S. thermophilus*. As shown by several recent studies, *S. thermophilus* strains possess anti-pathogenic properties associated with the production of bacteriocins [54]. However, we used a commercial and complex lactic starter aimed at controlling acidification in a wide range of cheese processes and subject to the manufacturer's culture strain rotation strategy to combat phages in the dairy industry. To the
best of our knowledge, no information is available on the ability of the strains included in
this lactic starter to produce bacteriocins, but we think it unlikely that strains with specific
anti-pathogenic properties (including production of bacteriocins) will be added in rotation
in this type of lactic starter.

Despite the decreasing relative abundance of *S. thermophilus* sequences during ripen-
ing, this species was still found in high relative abundance at the end of ripening, both in
the cheese cores (from 46% of the total 16S rRNA gene sequences with the 0.1× dose to
85% with the 2× dose) and in the cheese rinds (up to 75% with the 2× dose). We noticed
a 10-fold lower relative abundance of *S. thermophilus* in the cheese rinds with the 0.1×
dose compared with both of the other doses. It cannot be excluded that the abundance
of *S. thermophilus* may be overestimated after metabarcoding analyses due to the possi-
ble amplification of DNA from dead cells. However, the counts in viable and cultivable
streptococci in the cheese cores at D28 were still high, increasing with the lactic starter
dose and ranging from 5.7 × 10⁸ to 1.41 × 10⁹ CFU/g. In the cheese rinds, the levels of
thermophilic streptococci were not different in relation to the lactic starter dose and were
lower than those obtained in the cheese core (reached 2.38 × 10⁸ CFU/g). On the whole,
these results are consistent with those reported by Silvetti et al. [55] in Silter cheese (Italian
hard cheese), who noticed an increase from 7.3 to 77% of *Streptococcus* spp. between curd
and end of ripening (200 days) in cheeses with a thermophilic starter added at 10⁶ CFU/mL.
Conversely, several authors showed a reduction or an absence of detection of the added
streptococci in fresh cheese [28] and in semi-hard cheeses matured for up to 90 days [31],
or even 120 days [25].

In our work, at D8, we observed that, in the cheese cores, minority taxa such as *Lacti- caseibacillus*, *Staphylococcus* genera and species, or, at D28, *Enterococcus* and *Serratia* genera and species, were differentially abundant depending on the farm. However, the bacterial community composition of the cheese was mainly driven by the lactic starter dose. Re- duced richness and evenness (α-diversity) or dispersion of bacterial profiles to the centroid (beta-diversity), were observed with the 1× and the 2× lactic starter dose compared with the 0.1× dose. These results were consistent with those obtained by Choi et al. [28] who described, in semi-hard Cheddar cheese, higher α-diversity index (Shannon and Chao1) without an added lactic starter than with the addition of a thermophilic starter (inoculated at 1 × 10⁶ CFU/mL). Minority taxa (<7%) affiliated to the *Enterobacteriaceae* family, which includes potential opportunistic species, showed higher relative abundance in the cheeses with the lowest amount of lactic starter. Carafa et al. [24] also observed a higher relative abundance of *Enterobacteriaceae* in their control cheeses than in the cheeses with an added lactic starter. Furthermore, Silvetti et al. [55] reported lower specific richness in cheese with the addition of an autochthonous lactic starter, including unwanted bacteria as *Pseudomonas* or *Serratia* spp. Unlike other studies, we found that the dominant taxa, corresponding to genus or species with potential interest in cheese, also showed higher relative abundance with the 0.1× dose than with both of the other doses. Notably, the relative abundance of autochthonous *L. lactis* dropped from 17% to 0.2% in the cheese cores at D8, and from 21% to 0.8% at D28, by increasing the thermophilic starter dose. Similarly, in the cheese rinds, the relative abundance of *Brevibacterium* genus and species decreased from around 26% to 4% between the 0.1× dose and the 2× dose. In the early stages, *L. lactis* contributes to acidification and flavor development through the involvement of the lactose and citrate metabolisms [3,56] and can also produce antimicrobial compounds such as nisin [57]. Some species belonging to *Brevibacterium* genus, as *B. linens*, *B. antiquum* or *B. aurantiacum*, are involved in the cheese-ripening process and participate in the development of sensorial typicity (flavor and color) [58,59], as well as in proteolysis [60]. While the counts in cult- ivable LAB were comparable regardless of the lactic starter dose, the counts in ripening bacteria in the cheese cores decreased as the added lactic starter dose increased (Table S3). Given the major role of LAB, and of ripening bacteria as *Brevibacterium* sp. in cheese flavor development, the changes in their relative abundance bring into question the potential impact of the lactic starter dose on the cheese sensory qualities. Considering the slower
acidification and draining of the 0.1 × dose fresh cheeses, that resulted in lower dry matter content at D1 and D28, it is very likely that the lactic acid starter dose affected the ripened cheese texture and flavor. This assumption is supported by the increased Young's modulus values, indicating an increased cheese firmness with increased lactic starter dose, and likely negatively correlated to the proteolysis level [61].

**5. Conclusions**

To sum up, the 10-fold lower lactic starter dose enhanced the bacterial diversity in the cheese as well as the relative abundance of potentially helpful taxa (*Brevibacterium* or *Lactococcus* species) along with that of other potentially opportunistic ones (*Enterobacteriaceae* family). Microbial diversity can be associated with a protecting effect against pathogens in cheese [34,62,63]. In the present study, despite an increased bacterial diversity, the delayed acidification observed with the 0.1 × dose was likely responsible for the higher levels of the inoculated STEC and *S. aureus* in the cheeses.

Given that sanitary criteria were paramount to our study, it can be said that the 10-fold lower dose is to be proscribed. Nevertheless, with lower and more realistic levels of inoculated pathogens, an intermediate reduced dose of lactic starter, between the 0.1 × and the standard dose, may be sufficient to establish an efficient balance between acidification, pathogen growth and microbial diversity.

No difference in the level of pathogens or in bacterial diversity in the cheese was observed between the 1× and 2× doses. Further studies would be needed to determine if this was a result of an insufficient difference between the two tested doses or from a plateau effect.

Despite the reduced dispersion of cheese bacterial patterns with an increase in the dose of the lactic starter cultures, our study revealed differences in the bacterial community composition in relation to the farm of origin, not only in the raw milk, but also in the cheeses. Since the microbiota of raw milk is highly variable from one farm to another, the effects of adjusting the lactic starter dose should be evaluated with a larger range of raw milk.

To conclude, we confirmed that the modulation of the amount of the lactic starter could be a key to managing the qualities and the sanitary risks in cheese, although specific studies would be necessary to finely adjust the dose depending on the type of the lactic starter and on the type of the cheese. Insofar as the balance of the indigenous microbial populations of potential interest for cheese production has been altered in relation to the lactic starter dose, it would be important to assess the associated effects on the sensory qualities of the cheese.

**Supplementary Materials:** The following are available online at https://www.mdpi.com/article/10.3390/microorganisms9051081/s1, **Table S1:** Biochemical and microbial characteristics in the raw milk used (n = 3). **Table S2:** pH values, dry matter content (%), coagulase + *S. aureus* and STEC O26:H11 levels (CFU/g) and Young's modulus values (MPa) in the inoculated and control cheeses at D0.25 (pH only) D1 and D28. **Table S3:** Microbial counts (CFU/g) at D28 in the control cheese cores and rinds. **Table S4:** α-Diversity indexes calculated from the bacterial profiles in the raw milk, the cheese cores and the cheese rinds according to the farm and the lactic starter dose. **Table S5:** Statistical assessment (p-values) of bacterial profiles clustering and of intra-cluster dispersion calculated according to the farm and the dose, in the raw milk, the cheese cores and the cheese rinds. **Figure S1:** Relative abundance of the 10 major genera detected in the raw milk, organized by farms and days of cheese-making. **Table S6:** Mean relative abundance for taxa differentially abundant between tested conditions (farm or lactic starter dose).

**Author Contributions:** Conceptualization, C.D., C.C. (Christophe Chassard), A.L., and J.P.; Resources, C.C. (Cécile Callon), E.R. and S.T.; Visualization, J.P. and C.D.; Project administration, C.D.; Methodology, C.C. (Cécile Callon) and C.D.; Software, E.R. and S.T.; Validation, C.D. and C.C. (Christophe Chassard); Formal Analysis, J.P.; Investigation, J.P.; Data Curation, E.R. and S.T.; Writing—Original Draft Preparation, J.P. and C.D.; Writing—Review & Editing, C.D., C.C. (Christophe
Chassard) and C.C. (Cécile Callon); Supervision, C.D., A.L., C.C. (Christophe Chassard); Funding acquisition, C.D. All authors have read and agreed to the published version of the manuscript.

**Funding:** This work is part of the project CPER Qualimont. It is co-funded by The Auvergne-Rhône-Alpes Region (Project n° 17 011 032 01—22454) and the European Regional Development Fund (ERDF FC5093).

**Institutional Review Board Statement:** Not applicable.

**Data Availability Statement:** Raw sequence data are available at the European Nucleotide Archive (ENA) database under the accession number PRJEB44120.

**Acknowledgments:** The authors thank René Lavigne (INRAE UMR 0545 Fromage) for the making of the cheese, Romain Pigearias for his technical help in DNA extractions, Beatrice Desserre (INRAE UMR 0545 Fromage) for laboratory management, Isabelle Constant (INRAE UMR1213 Herbivores) for her assistance in the rheological measures and Céline Gres (INRAE UMR 0545 Fromage) for administrative support. They also thank Bruno Martin (INRAE UMR1213 Herbivores), Hélène Falentin (INRAE UMR 1253 Science et Technologie du Lait et de l'Œuf) and Anne Saint-Eve (AgroParisTech UMR 782 Génie et microbiologie des procédés alimentaires) for their helpful feedback on the experiment design, Patricia Mabrut for her very helpful English revisions and thank Alain Kondjoyan (INRAE UR0370 Qualité des Produits Animaux) for the management of the project CPER Qualimont.

**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 writing of the manuscript, or in the decision to publish the results.

## References

1. Settanni, L.; Moschetti, G. New trends in technology and identity of traditional dairy and fermented meat production processes: Preservation of typicality and hygiene. *Trends Food Sci. Technol.* **2014**, *37*, 51–58. [CrossRef]

2. Zindy, P.; Halawany-Darson, R.; Hauwuy, A. Consumers’ perception and attitude towards regional “terroir” products: Focus on “terroir” cheeses. *INRA Prod. Anim.* **2017**, *12*, 229–240.

3. Blaya, J.; Barzideh, Z.; LaPointe, G. Symposium review: Interaction of starter cultures and nonstarter lactic acid bacteria in the cheese environment. *J. Dairy Sci.* **2018**, *101*, 3611–3629. [CrossRef] [PubMed]

4. Callon, C.; Berdagué, J.L.; Dufour, E.; Montel, M.C. The effect of raw milk microbial flora on the sensory characteristics of salers-type cheeses. *J. Dairy Sci.* **2005**, *88*, 3840–3850. [CrossRef]

5. Gobbetti, M.; Di Cagno, R.; Calasso, M.; Neviani, E.; Fox, P.F.; De Angelis, M. Drivers that establish and assembly the lactic acid bacteria biota in cheeses. *Trends Food Sci. Technol.* **2018**, *78*, 244–254. [CrossRef]

6. Gobbetti, M.; De Angelis, M.; Di Cagno, R.; Mancini, L.; Fox, P.F. Pros and cons for using Non-Starter Lactic Acid Bacteria (NSLAB) as secondary/adjunct starters for cheese ripening. *Trends Food Sci. Technol.* **2015**, *45*, 167–178. [CrossRef]

7. Kongo, J.M. Lactic acid bacteria as starter-cultures for cheese processing: Past, present and future developments. In *Lactic Acid Bacteria—R & D for Food, Health and Livestock Purposes*; InTech Open Limited: London, UK, 2013. [CrossRef]

8. Frétin, M.; Martin, B.; Rifa, E.; Isabelle, V.-M.; Pomiès, D.; Ferlay, A.; Montel, M.-C.; Delbès, C. Bacterial community assembly from cow teat skin to ripened cheeses is influenced by grazing systems. *Sci. Rep.* **2018**, *8*, [CrossRef] [PubMed]

9. Montel, M.-C.; Buchin, S.; Mallet, A.; Delbes-Paus, C.; Vuitton, D.A.; Desmasures, N.; Berthier, F. Traditional cheeses: Rich and diverse microbiota with associated benefits. *Int. J. Food Microbiol.* **2014**, *177*, 136–154. [CrossRef]

10. Johnson, M.E. A 100-year review: Cheese production and quality. *J. Dairy Sci.* **2017**, *100*, 9952–9965. [CrossRef]

11. O'Sullivan, D.J.; Giblin, L.; McSweeney, P.L.H.; Sheehan, J.J.; Cotter, P.D. Nucleic acid-based approaches to investigate microbial-related cheese quality defects. *Front. Microbiol.* **2013**, *4*. [CrossRef]

12. Donnelly, C. *Review of Controls for Pathogen Risks in Scottish Artisan Cheeses Made from Unpasteurised Milk*; Food Standards Scotland: Aberdeen, UK, 2018.

13. Fournet, N.; Laurent, E.; Jones, G.; Jourdan Da Silva, N.; Tourdjman, M.; Chereau, F.; Nisavanh, A.; De Valk, H. *Surveillance of Collective Food Toxi-Infections. Data from the Mandatory Declaration 2018*; Santé Publique France: Paris, France, 2020.

14. Bonanno, L.; Delubac, B.; Michel, V.; Auvray, F. Influence of stress factors related to cheese-making process and to STEC detection procedure on the induction of stx phages from STEC O26:H11. *Front. Microbiol.* **2017**, *8*. [CrossRef]

15. Delbes, C.; Alomar, J.; Chougui, N.; Martin, J.-F.; Montel, M.-C. *Staphylococcus aureus* Growth and enterotoxin production during the manufacture of uncooked, semihard cheese from cows’ raw milk. *J. Food Prot.* **2006**, *69*, 2161–2167. [CrossRef] [PubMed]

16. Miszczycha, S.D.; Perrin, F.; Ganet, S.; Jamet, E.; Tenenhaus-Aziza, F.; Montel, M.-C.; Thevenot-Sergentet, D. Behavior of different shiga toxin-producing *Escherichia coli* serotypes in various experimentally contaminated raw-milk cheeses. *Appl. Environ. Microbiol.* **2013**, *79*, 150–158. [CrossRef]

17. Parente, E.; Guidone, A.; Matera, A.; De Filippis, F.; Mauriello, G.; Ricciardi, A. Microbial community dynamics in thermophilic undefined milk starter cultures. *Int. J. Food Microbiol.* **2016**, *217*, 59–67. [CrossRef]
18. Lucchini, R.; Cardazzo, B.; Carraro, L.; Negrinotti, M.; Balzan, S.; Novelli, E.; Fasolato, L.; Fasoli, F.; Farina, G. Contribution of natural milk culture to microbiota, safety and hygiene of raw milk cheese produced in alpine malga. *Ital. J. Food Saf.* **2018**, *7*. [CrossRef]

19. International Standard Organisation (ISO). ISO 6888-2:1999(en) Microbiology of Food and Animal Feeding Stuffs—Horizontal Method for the Enumeration of Coagulase-Positive Staphylococci (Staphylococcus aureus and Other Species)—Part 2: Technique Using Rabbit Plasma Fibrinogen Agar Medium; CEN Comite Europeen de Normalisation: Brussels, Belgium, 1999; Volume 7.

20. International Standard Organisation (ISO). ISO/TS 13136:2012(en) Microbiology of Food and Animal Feed—Real-Time Polymerase Chain Reaction (PCR)-Based Method for the Detection of Food-Borne Pathogens—Horizontal Method for the Detection of Shiga Toxin-Producing Escherichia coli (STEC) and the Determination of O157, O111, O26, O103 and O145 Serogroups; CEN Comite Europeen de Normalisation: Brussels, Belgium, 2012; Volume 22.

21. Mitra, S. Multiple data analyses and statistical approaches for analyzing data from metagenomic studies and clinical trials. In *Evolutionary Genomics: Statistical and Computational Methods*; Methods in Molecular Biology; Anisimova, M., Ed.; Springer: New York, NY, USA, 2019; pp. 605–634, ISBN 978-1-4939-9074-0.

22. Quigley, L.; O'Sullivan, O.; Beresford, T.P.; Ross, R.P.; Fitzgerald, G.F.; Cotter, P.D. High-Throughput sequencing for detection of subpopulations of bacteria not previously associated with artisanal cheeses. *Appl. Environ. Microbiol.* **2012**, *78*, 5717–5723. [CrossRef]

23. O'Sullivan, D.J.; Cotter, P.D.; O'Sullivan, O.; Giblin, L.; McSweeney, P.L.H.; Sheehan, J.J. Temporal and spatial differences in microbial composition during the manufacture of a continental-type cheese. *Appl. Environ. Microbiol.* **2015**, *81*, 2525–2533. [CrossRef] [PubMed]

24. Carafa, I.; Stocco, G.; Franceschi, P.; Summer, A.; Tuohy, K.M.; Bittante, G.; Franciosi, E. Evaluation of autochthonous lactic acid bacteria as starter and non-starter cultures for the production of traditional mountain cheese. *Food Res. Int.* **2019**, *115*, 209–218. [CrossRef] [PubMed]

25. Feutry, F.; Oneca, M.; Berthier, F.; Torre, P. Biodiversity and growth dynamics of lactic acid bacteria in artisanal PDO ossau-iraty cheeses made from raw ewe's milk with different starters. *Food Microbiol.* **2012**, *29*, 33–42. [CrossRef]

26. Settanni, L.; Di Grigoli, A.; Tornambé, G.; Bellina, V.; Francesca, N.; Moschetti, G.; Bonanno, A. Persistence of wild *Streptococcus thermophilus* strains on wooden vat and during the manufacture of a traditional caciocavallo type cheese. *Int. J. Food Microbiol.* **2012**, *155*, 73–81. [CrossRef]

27. Kamimura, B.A.; Cabral, L.; Noronha, M.F.; Baptista, R.C.; Nascimento, H.M.; Sant'Ana, A.S. Amplicon sequencing reveals the bacterial diversity in milk, dairy premises and Serra Da Canastra artisanal cheeses produced by three different farms. *Food Microbiol.* **2020**, *89*, 103453. [CrossRef]

28. Choi, J.; In Lee, S.; Rackerby, B.; Frojen, R.; Goddik, L.; Ha, S.-D.; Park, S.H. Assessment of overall microbial community shift during cheddar cheese production from raw milk to aging. *Appl. Microbiol. Biotechnol.* **2020**, *104*, 6249–6260. [CrossRef]

29. Fernandez-García, E.; Tomillo, J.; Nuñez, M. Formation of biogenic amines in raw milk Hispánico cheese manufactured with proteinases and different levels of starter culture. *J. Food Prot.* **2000**, *63*, 1551–1555. [CrossRef] [PubMed]

30. Khosrowshahi, A.; Madadlou, A.; Ebrahim Zadeh Mousavi, M.; Emam-Djomeh, Z. Monitoring the chemical and textural changes during ripening of Iranian white cheese made with different concentrations of starter. *J. Dairy Sci.* **2006**, *89*, 3318–3325. [CrossRef]

31. Özer, E.; Kesenkaş, H. The effect of using different starter culture combinations on ripening parameters, microbiological and sensory properties of Mihaliç cheese. *J. Food Sci. Technol.* **2019**, *56*, 1202–1211. [CrossRef] [PubMed]

32. De Buyser, M.L.; Lombard, B.; Schulten, S.M.; In't Veld, P.H.; Scotter, S.L.; Rollier, P.; Lahellec, C. Validation of EN ISO standard methods 6888 Part 1 and Part 2: 1999—Enumeration of coagulase-positive staphylococci in foods. *Int. J. Food Microbiol.* **2003**, *83*, 185–194. [CrossRef]

33. International Standard Organisation (ISO). ISO 19020:2017(en) Microbiology of the Food Chain—Horizontal Method for the Immunoenzymatic Detection of Staphylococcal Enterotoxins in Foodstuffs; CEN Comite Europeen de Normalisation: Brussels, Belgium, 2017; Volume 21.

34. Frétin, M.; Chassard, C.; Delbès, C.; Lavigne, R.; Rifa, E.; Theil, S.; Fernandez, B.; Laforce, P.; Callon, C. Robustness and efficacy of an inhibitory consortium against *E. coli* O26:H11 in raw milk cheeses. *Food Control* **2020**, *115*, 107282. [CrossRef]

35. Millet, L.; Saubusse, M.; Didienne, R.; Tessier, L.; Montel, M.C. Control of *Listeria monocytogenes* in raw-milk cheeses. *Int. J. Food Microbiol.* **2006**, *108*, 105–114. [CrossRef] [PubMed]

36. International Standard Organisation (ISO). ISO 5534:2004 (en) (IDF 4:2004) Cheese and Processed Cheese—Determination of the Total Solids Content (Reference Method); CEN Comite Europeen de Normalisation: Brussels, Belgium, 2004; Volume 7.

37. Frétin, M.; Ferlay, A.; Verdier-Metz, I.; Fournier, F.; Montel, M.-C.; Farruggia, A.; Delbès, C.; Martin, B. The effects of low-input grazing systems and milk pasteurisation on the chemical composition, microbial communities, and sensory properties of uncooked pressed cheeses. *Int. Dairy J.* **2017**, *64*, 56–67. [CrossRef]

38. Noël, Y.; Ardö, Y.; Pocher, S.; Hunter, A.; Lavanchy, P.; Luginbuhl, W.; Bars, D.L.; Polychroniadou, A.; Pellegrino, L. Characterisation of protected denomination of origin cheeses: Relationships between sensory texture and instrumental data. *Lait* **1998**, *78*, 569–588. [CrossRef]

39. Terzaghi, B.E.; Sandine, W.E. Improved medium for lactic Streptococci and their bacteriophages. *Appl. Microbiol.* **1975**, *29*, 807–813. [CrossRef]

40. Man, J.C.D.; Rogosa, M.; Sharpe, M.E. A medium for the cultivation of lactobacilli. *J. Appl. Bacteriol.* **1960**, *23*, 130–135. [CrossRef]
41. Denis, C.; Gueguen, M.; Henry, E.; Levert, D. New media for the numeration of cheese surface bacteria. *Lait* **2001**, *81*, 365–379.
[CrossRef]

42. Mossel, D.A.; Kleynen-Semmeling, A.M.; Vincentie, H.M.; Beerens, H.; Catsaras, M. Oxytetracycline-Glucose-Yeast extract agar for selective enumeration of moulds and yeasts in foods and clinical material. *J. Appl. Bacteriol.* **1970**, *33*, 454–457. [CrossRef] [PubMed]

43. Duthoit, F.; Godon, J.-J.; Montel, M.-C. Bacterial community dynamics during production of registered designation of origin salers cheese as evaluated by 16S RNA gene single-strand conformation polymorphism analysis. *Appl. Environ. Microbiol.* **2003**, *69*, 3840–3848. [CrossRef] [PubMed]

44. Poirier, S.; Rué, O.; Peguilhan, R.; Coeuret, G.; Zagorec, M.; Champomier-Vergès, M.-C.; Loux, V.; Chaillou, S. 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. *PLoS ONE* **2018**, *13*, e0204629. [CrossRef]

45. Theil, S.; Rifa, E. Amplicon workflow for microbial community analysis. *F1000Research* **2021**, *10*, 7. [CrossRef]

46. Callahan, B.J.; McMurdie, P.J.; Rosen, M.J.; Han, A.W.; Johnson, A.J.A.; Holmes, S.P. DADA2: High-Resolution sample inference from illumina amplicon data. *Nat. Methods* **2016**, *13*, 581–583. [CrossRef]

47. Delbès-Paus, C.; Miszczycha, S.; Ganet, S.; Helinck, S.; Veisseire, P.; Pochet, S.; Thévenot, D.; Montel, M.-C. Behavior of *Escherichia coli* O26:H11 in the Presence of *Hafnia alvei* in a model cheese ecosystem. *Int. J. Food Microbiol.* **2013**, *160*, 212–218. [CrossRef]

48. Callon, C.; Arliguie, C.; Montel, M.-C. Control of shigatoxin-producing *Escherichia coli* in cheese by dairy bacterial strains. *Food Microbiol.* **2016**, *53*, 63–70. [CrossRef]

49. Medved'ová, A.; Koňuchová, M.; Kvočiková, K.; Hatalová, I.; Valík, L. Effect of lactic acid bacteria addition on the microbiological safety of pasta-filata types of cheeses. *Front. Microbiol.* **2020**, *11*. [CrossRef]

50. European Commission. Commission of the European Communities commission regulation (EC) No 2073/2005 of 15 November 2005 on microbiological criteria for foodstuffs. *Off. J. Eur. Union* **2005**, *L338*, 1–26.

51. Alomar, J.; Lebert, A.; Montel, M.C. Effect of temperature and PH on growth of *Staphylococcus aureus* in co-culture with *Lactococcus garvieae*. *Curr. Microbiol.* **2008**, *56*, 408–412. [CrossRef] [PubMed]

52. Delpech, P.; Bornes, S.; Alaterre, E.; Bonnet, M.; Gagne, G.; Montel, M.-C.; Delbès, C. *Staphylococcus aureus* transcriptomic response to inhibition by H₂O₂-producing *Lactococcus garvieae*. *Food Microbiol.* **2015**, *51*, 163–170. [CrossRef] [PubMed]

53. Al-Nabulsi, A.A.; Osaili, T.M.; AbuNaser, R.A.; Olaimat, A.N.; Ayyash, M.; Al-Holy, M.A.; Kadora, K.M.; Holley, R.A. Factors affecting the viability of *Staphylococcus aureus* and production of enterotoxin during processing and storage of white-brined cheese. *J. Dairy Sci.* **2020**, *103*, 6869–6881. [CrossRef]

54. Evivie, S.E.; Ogwu, M.C.; Abdelazez, A.; Bian, X.; Liu, F.; Li, B.; Huo, G. Suppressive effects of *Streptococcus thermophilus* KLDS 3.1003 on some foodborne pathogens revealed through in vitro, in vivo and genomic insights. *Food Funct.* **2020**, *11*, 6573–6587. [CrossRef]

55. Silvetti, T.; Capra, E.; Morandi, S.; Cremonesi, P.; Decimo, M.; Gavazzi, F.; Giannico, R.; De Noni, I.; Brasca, M. Microbial population profile during ripening of Protected Designation of Origin (PDO) siller cheese, produced with and without autochthonous starter culture. *LWT* **2017**, *84*, 821–831. [CrossRef]

56. Marilley, L.; Casey, M.G. Flavours of cheese products: Metabolic pathways, analytical tools and identification of producing strains. *Int. J. Food Microbiol.* **2004**, *90*, 139–159. [CrossRef]

57. Bukvicki, D.; Siroli, L.; D'Alessandro, M.; Cosentino, S.; Fliss, I.; Said, L.B.; Hassan, H.; Lanciotti, R.; Patrignani, F. Unravelling the potential of *Lactococcus lactis* strains to be used in cheesemaking production as biocontrol agents. *Foods* **2020**, *9*, 1815. [CrossRef]

58. Anast, J.M.; Dzieciol, M.; Schultz, D.L.; Wagner, M.; Mann, E.; Schmitz-Esser, S. Brevibacterium from Austrian hard cheese harbor a putative histamine catabolism pathway and a plasmid for adaptation to the cheese environment. *Sci. Rep.* **2019**, *9*, 6164. [CrossRef]

59. Bockelmann, W.; Willems, K.P.; Neve, H.; Heller, K.H. Cultures for the ripening of smear cheeses. *Int. Dairy J.* **2005**, *15*, 719–732. [CrossRef]

60. Boyaval, P.; Desmazaud, M.J. Current state of research on Brevibacterium linens. *Lait* **1983**, *63*, 187–216. [CrossRef]

61. Famelart, M.H.; Graet, Y.L.; Michel, F.; Richoux, R.; Riaublanc, A. Evaluation of the methods of measurement for functional properties of Emmental cheeses from the west of France. *Lait* **2002**, *82*, 225–245. [CrossRef]

62. Saubusse, M.; Millet, L.; Delbès, C.; Callon, C.; Montel, M.C. Application of single strand conformation polymorphism—PCR method for distinguishing cheese bacterial communities that inhibit *Listeria monocytogenes*. *Int. J. Food Microbiol.* **2007**, *116*, 126–135. [CrossRef]

63. Gay, M.; Amgar, A. Factors moderating *Listeria monocytogenes* growth in raw milk and in soft cheese made from raw milk. *Lait* **2005**, *85*, 153–170. [CrossRef]