# Lactococcus lactis Diversity Revealed by Targeted Amplicon Sequencing of purR Gene, Metabolic Comparisons and Antimicrobial Properties in an Undefined Mixed Starter Culture Used for Soft-Cheese Manufacture

Sabrina Saltaji, Olivier Rué, Valérie Sopena, Sophie Sablé, Fatoumata Tambadou, Sandrine Didelot, Romain Chevrot

► To cite this version:

Sabrina Saltaji, Olivier Rué, Valérie Sopena, Sophie Sablé, Fatoumata Tambadou, et al., Lactococcus lactis Diversity Revealed by Targeted Amplicon Sequencing of purR Gene, Metabolic Comparisons and Antimicrobial Properties in an Undefined Mixed Starter Culture Used for Soft-Cheese Manufacture. Foods, 2020, 9 (5), pp.622. 10.3390/foods9050622. hal-03516147

HAL Id: hal-03516147

https://hal.science/hal-03516147

Submitted on 26 Aug 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
Article

*Lactococcus lactis* Diversity Revealed by Targeted
Amplicon Sequencing of *purR* Gene, Metabolic
Comparisons and Antimicrobial Properties in an
Undefined Mixed Starter Culture Used for
Soft-Cheese Manufacture

Sabrina Saltaji <fnref n="1" />,<fnref n="2" />, Olivier Rué <fnref n="3" />,<fnref n="4" />✉, Valérie Sopena <fnref n="1" />, Sophie Sablé <fnref n="1" />, Fatoumata Tambadou <fnref n="2" />, Sandrine Didelot <fnref n="1" />,† and Romain Chevrot <fnref n="1" />,*,†

Received: 10 April 2020; Accepted: 7 May 2020; Published: 13 May 2020

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

check for
updates

**Abstract:** The undefined mixed starter culture (UMSC) is used in the manufacture of cheeses. Deciphering UMSC microbial diversity is important to optimize industrial processes. The UMSC was studied using culture-dependent and culture-independent based methods. MALDI-TOF MS enabled identification of species primarily from the *Lactococcus* genus. Comparisons of carbohydrate metabolism profiles allowed to discriminate five phenotypes of *Lactococcus* (n = 26/1616). The 16S sequences analysis (V1–V3, V3–V4 regions) clustered the UMSC microbial diversity into two *Lactococcus* operational taxonomic units (OTUs). These clustering results were improved with the DADA2 algorithm on the housekeeping *purR* sequences. Five *L. lactis* variants were detected among the UMSC. The whole-genome sequencing of six isolates allowed for the identification of the *lactis* subspecies using Illumina® (n = 5) and Pacbio® (n = 1) technologies. Kegg analysis confirmed the *L. lactis* species-specific niche adaptations and highlighted a progressive gene pseudogenization. Then, agar spot tests and agar well diffusion assays were used to assess UMSC antimicrobial activities. Of note, isolate supernatants (n = 34/1616) were shown to inhibit the growth of *Salmonella* ser. Typhimurium CIP 104115, *Lactobacillus sakei* CIP 104494, *Staphylococcus aureus* DSMZ 13661, *Enterococcus faecalis* CIP103015 and *Listeria innocua* CIP 80.11. Collectively, these results provide insightful information about UMSC *L. lactis* diversity and revealed a potential application as a bio-protective starter culture.

**Keywords:** starter cultures; dairy *Lactococcus*; carbohydrate metabolism; DADA2 algorithm; antimicrobial activities; lactococcin B-like bacteriocin

**1. Introduction**

Artisanal soft cheeses are a class of handmade products following ancient practices and
manufactured with a limited mechanization process. These cheeses include the use of undefined
mixed starter cultures (UMSC) for traditional dairy-based food fermentations [1-3]. UMSC are
composed of undefined multiple-strain cultures of lactic acid bacteria (LAB). They harbor complex
microbiota, which naturally proceed from raw milk, and are specific to the producing region [4].

<footnote n="1">¹ UMR CNRS 7266 LIENSs, Université de La Rochelle, 17000 La Rochelle, France; sabrina.saltaji1@univ-lr.fr (S.S.); vgauthie@univ-lr.fr (V.S.); sophie.sable@univ-lr.fr (S.S.); sandrine.didelot@univ-lr.fr (S.D.)</footnote>

<footnote n="2">² Flores de Terroirs, 17440 Aytré, France; fatoumata.tambadou@flores-de-terroirs.com</footnote>

<footnote n="3">³ INRAE, Université Paris-Saclay, 78350 Jouy-en-Josas, France; olivier.rue@inra.fr</footnote>

<footnote n="4">⁴ INRAE, BioinfOmics, MIGALE Bioinformatics Facility, Université Paris-Saclay, 78350 Jouy-en-Josas, France</footnote>

<footnote>* Correspondence: romain.chevrot01@univ-lr.fr</footnote>

<footnote>† Authors contributed equally to this study.</footnote>
Most important bacteria in dairy UMSC are *Lactococcus lactis*, *Leuconostoc mesenteroides*, *Streptococcus* spp. and *Lactobacillus* spp. [1,3,5,6]. However, limited knowledge is currently available on the genetic diversity of UMSC [7]. Changes caused by intrinsic milk properties, but also by other factors such as the dairy environment and cheese making procedures, can influence the UMSC microbial composition. Therefore its complexity is specific beyond the species/subspecies level [8–10]. From the beginning of dairying, Lactococci have been widely used as a principal component of starter cultures [6,11]. Both subspecies *L. lactis* subsp. *lactis* and *L. lactis* subsp. *cremoris* are the upmost contributors to the acidification process following the rapid fermentation of lactose [7]. The diversity of *Lactococcus* strains provides an important reservoir for starters and may contribute to flavor differences and specific features in a fermented milk-product [12]. Lactococci participate in food preservation, texture formation and flavor by producing aromatic compounds (alcohols, aldehydes, ketones) or by metabolizing substrates such as citrate, amino acid or fats [13]. The growing interest to preserve traditional cheeses has stimulated the use of raw milk UMSC and thus, a better understanding of the diversity of *Lactococcus lactis* strains is required. High-throughput sequencing (HTS) and whole-genome sequencing (WGS) technologies have opened a new window to decipher such microbial diversity. Lately, industrial and technological relevance has motivated extensive genomic studies [14–16] to better predict production performance and advise culture strategies to avoid fermentation failures [17]. Comparative genomics has also provided insightful information into genotype–phenotype associations. Phenotypic data and Lactococci strain-specific information are effective methods for the assignment of functions underlying industrial traits or niche adaptation [18,19].

Starter culture use implies the respect of safety and quality patterns in the dairy industry. Spoilage bacteria and foodborne pathogens are of major concern in food production and processing environments as well as for consumers since they cause economic losses and foodborne illnesses. Some LAB species, especially *L. lactis*, can exert a food preservative effect due to a nutritional competition ability or production of organic acids or antimicrobial compounds. The pH reduction makes the dairy matrix unfavorable for the growth of food-spoilage or pathogenic microorganisms. LAB can also produce metabolic end-products including hydrogen peroxide (H₂O₂), CO₂, diacetyl, acetaldehydes, D-isomers of amino acids, reuterin and bacteriocins [20,21]. These antimicrobial strategies are increasingly used for limiting the growth of unwanted microorganisms [22]. In a raw milk starter culture, the food micro-structure is in a liquid form and favors the accumulation of biomass and metabolites causing chemical changes such as a decline of pH. In a gel food microstructure (such as in cheese), microorganisms are immobilized and constrained to form colonies accompanied by local accumulation of metabolic end-products. The immobilization of microorganisms in a cheese sample might result in interactions between colonies and local pH changes to stimulate an acid tolerance response. Therefore, food preservation depends also on the structure of the product [23].

In this study, genotype–phenotype associations and bacterial screenings were carried out to
characterize the UMSC bacterial diversity to better understand UMSC carbohydrate metabolism and
to evaluate antimicrobial activities.

**2. Materials and Methods**

2.1. *Sample Collection*

The UMSC was obtained following raw milk treatment as described in the Bionatif® patented
process [24]. The raw milk was obtained from a farm located in Normandy (Isigny-sur-mer, France),
and the UMSC was used daily for soft-cheese manufacture. The UMSC was kept frozen at -80 °C
(Flores de terroirs, Aytré, France).
## 2.2. Gene Sequencing of the UMSC Microbiota

### 2.2.1. DNA Extraction

Samples of UMSC were aseptically homogenized and centrifuged at 13,000× g for 10 min at 4 °C. DNA was extracted with the DNeasy PowerFood Microbial Kit (Qiagen, Hilden, Germany), which combines mechanical (beat beating) and chemical lysis of the cell. The purified DNA was stored at -80 °C.

### 2.2.2. Gene Sequencing and Data Analysis of 16S rRNA

The V1–V3 and the V3–V4 regions of the 16S rRNA genes were amplified by PCR. The V1–V3 bacterial universal primers were selected [25,26]. Primers 27F (3’-AGAGTTTGATCATGGCTCAG-5’) and 1492R (3’-TACGGYTACCTTGTTACGACTT-5’) were used for the V3–V4 amplification fragment. For each sample, an amplified DNA quality check was performed. The amplified DNA was then sequenced by the Eurofins Genomics platform (Ebersberg, Germany) using an Illumina® MiSeq Technology.

The FROGS (find rapidly OTUs with Galaxy solution) pipeline (version 3.1) was used to retrieve the operational taxonomic units (OTUs) under the Galaxy environment [27]. This pipeline includes raw FASTQ files demultiplexing and a preprocess step to filter and delete sequences with unexpected lengths, with ambiguous bases (N) and the ones that do not have primer sequences at both 3' and 5' ends. The clustering of sequences was performed by Swarm using an iterative growth process with a local clustering threshold of 3 [28]. Chimera sequences were identified and removed using Vsearch [29]. OTUs accounting for more than 0.5‰ of the total reads were considered. The reference sequences of each OTU were blasted against the SILVA database (version 128) for taxonomic assignation [30].

### 2.2.3. *purR* Gene Sequencing and Data Analysis

The *purR* gene was amplified by PCR with primers purR-324F (5’-YACTCCATCAAATC TTCGTAAAAT-3’), purR-811R (5’-TGTCATTAAATATATTTCCCAATTGAACA-3’) described by Frantzen et al. [7]. Forward (5’-TCGTCGGCAGCGTCAGATGTGTATAAGAGACAG) and reverse (5’-GTCTCGTGGGCTCGGAGATGTGTATAAGAGACAG) Illumina® adapter overhangs were added to the 5’ ends of primers to enable Nextera XT DNA indexing of the PCR products. The libraries were sequenced using an Illumina® MiSeq technology on the Eurofins Genomics platform (Ebersberg, Germany). The attempted amplicon size was 487 pb [7].

The data were analyzed using the DADA2 pipeline (v. 1.12, http://benjjneb.github.io/dada2/tutorial.html), which retrieves biological sequences from reads by modelling the Illumina-sequencing errors [31]. Once read quality profiles were validated, the core-denoising algorithm of DADA2 was performed on the forward and reverse reads separately. Amplicon sequence variants (ASV) were inferred and pairs merged. ASVs table was constructed and chimeras were identified and removed only if the parent sequence was at least twice abundant than a left/right-segment recombination. The taxonomic affiliation was assigned using a local BLAST database produced using the lactococcal genomes available on the NCBI database.

## 2.3. Bacterial Library

Presumptive LAB were enumerated on MRS, BHI or M17 agar (Biokar, Allonne, France) according to the conditions reported in Table 1. Bacterial isolates were isolated (n = 1616), purified twice and maintained frozen at -80 °C in the corresponding growth medium containing 20% (v/v) glycerol (Table 1).
<table><caption>**Table 1.** Growth conditions used for lactic acid bacteria (LAB) and non-LAB isolation.</caption><thead><tr><th colspan="2">Growth Conditions</th></tr></thead><tbody><tr><td rowspan="6">Aerobic atmosphere</td><td rowspan="6">28 °C</td><td>BHI agar</td></tr><tr><td>MRS agar</td></tr><tr><td>M17 agar</td></tr><tr><td>MRS agar + 2% NaCl</td></tr><tr><td>MRS agar + 4% NaCl</td></tr><tr><td>MRS agar + 6.5% NaCl</td></tr><tr><td>BHI agar + 6.5% NaCl</td></tr><tr><td rowspan="3">Microaerophilic atmosphere</td><td rowspan="3">28 °C</td><td>BHI agar</td></tr><tr><td>MRS agar</td></tr><tr><td>M17 agar</td></tr><tr><td rowspan="3">Aerobic atmosphere</td><td>22 °C</td><td>BHI agar</td></tr><tr><td>37 °C</td><td>MRS agar</td></tr><tr><td>44 °C</td><td>MRS agar</td></tr></tbody></table>

Each isolate was seeded twice on BHI agar plates (Biokar Diagnostics, Allonne, France), at 28 °C for 30 h.

## 2.4. Phenotypic Profiles and LAB Isolates Identification

### 2.4.1. Bacterial Carbohydrate Metabolism Profile

Two isolates issued from growth conditions were randomly selected. Their carbohydrate metabolism profile was determined twice using the API 50 CH system. Overall, it consisted of 49 biochemical tests. API 50 CH was used in conjunction with API 50 CHL medium based on the manufacturer's instructions (Biomerieux, Marcy l'Etoile, France).

### 2.4.2. MALDI-TOF MS Bacterial Identification

MALDI-TOF MS (Biotyper®, Bruker, Palaiseau, France) identifications were performed on each bacterial isolate studied for their carbohydrate metabolism profile. MALDI-TOF MS analysis was performed by Eurofins Microbiologie Ouest based on their internal procedures (Nantes, France).

### 2.4.3. DNA Extraction, Routine PCR Procedure and Full-Length 16S rRNA Gene Sequencing

Full-length 16S rRNA gene sequencing was performed by the Eurofins Genomics Europe Applied (Ebersberg, Germany) according to the company's internal test ID BJ00M Bacterial Species Identification. The obtained sequences were aligned using the MEGA7 program package [32] and compared on the NCBI database using the BLAST program (https://blast.ncbi.nlm.nih.gov/Blast.cgi) [33].

## 2.5. Genome Sequencing and Data Analysis

### 2.5.1. Pure Culture DNA Extraction

Six *L. lactis* isolates were chosen following API 50 CH results to perform genotype-phenotype associations. Overnight cultures from these isolates were centrifuged at 13,000× g for 10 min at 4 °C. These Genomic DNA was extracted using the DNeasy blood and tissue kit (Qiagen, Germany) according to the manufacturer's instruction. Genomic DNA was quantified using the NanoDrop ND-1000 (ThermoFisher Scientific, Illkirch-Graffenstaden, France).

### 2.5.2. Whole Genome Sequencing and Gene Annotation

Whole-genome sequence data were generated by PacBio® RS II and Illumina® iSeq sequencing technology in order to produce de novo assemblies. The PacBio® RS II data were obtained and
assembled using the Hierarchical Genome Assembly Process (HGAP) assembly algorithm [34], on the Eurofins Genomics platform (Ebersberg, Germany). The Illumina® iSeq data were produced and assembled using the SPAdes algorithm (Mérieux Nutrisciences, Nantes, France) [35]. Chromosomal organization of the resulting contigs were predicted by aligning them against the *L. lactis* subsp. *lactis* IL1403 reference genome using the MicroScope platform [36]. The individual genome assemblies were deposited in the NCBI database with accession codes listed in Table 2. Overall genomic relatedness between genomes was assessed using the Orthologous Average Nucleotide Identity (OrthoANI) tool [37]. Carbohydrate metabolism-related pathway information was identified from those described in the Kyoto Encyclopedia of Genes and Genomes (KEGG, http://www.KEGG.jp/) [38]. The mining of genes clusters related to the production of antimicrobial peptides was performed using BAGEL4 [39].

<table><caption>**Table 2.** Genomic features of subspecies in the genus *Lactococcus*.</caption><thead><tr><th>Strains</th><th>Sequencing Technology</th><th>Scaffold Number</th><th>Genome Size (Mb)</th><th>GC %</th><th>Predicated Coding Genes</th><th>GenBank Accession Number</th></tr></thead><tbody><tr><td><i>L. lactis</i> subsp. <i>lactis</i> L1</td><td>Illumina® iSeq</td><td>177</td><td>2.75</td><td>34.96</td><td>2945</td><td>JAATTT000000000</td></tr><tr><td><i>L. lactis</i> subsp. <i>lactis</i> L2</td><td>Illumina® iSeq</td><td>193</td><td>2.78</td><td>34.99</td><td>3003</td><td>JAATTS000000000</td></tr><tr><td><i>L. lactis</i> subsp. <i>lactis</i> L3</td><td>PacBio® RS II</td><td>7</td><td>2.84</td><td>35.22</td><td>3024</td><td>JAATTU000000000</td></tr><tr><td><i>L. lactis</i> subsp. <i>lactis</i> L8</td><td>Illumina® iSeq</td><td>242</td><td>2.81</td><td>34.99</td><td>3036</td><td>JAATTR000000000</td></tr><tr><td><i>L. lactis</i> subsp. <i>lactis</i> L9</td><td>Illumina® iSeq</td><td>229</td><td>2.80</td><td>34.98</td><td>3014</td><td>JAATTQ000000000</td></tr><tr><td><i>L. lactis</i> subsp. <i>lactis</i> L14</td><td>Illumina® iSeq</td><td>178</td><td>2.78</td><td>34.98</td><td>2966</td><td>JAATTP000000000</td></tr><tr><td><i>L. lactis</i> subsp. <i>lactis</i> IL1403</td><td>[40]</td><td>[40]</td><td>2.37</td><td>35.33</td><td>2473</td><td>AE005176.1</td></tr></tbody></table>

## 2.6. Accession Number(s)

The 16S and *purR* amplicon data were deposited at DDBJ/ENA/GenBank under BioProject number PRJNA615620. The whole-genome project was deposited at DDBJ/ENA/GenBank under BioProject number PRJNA615395.

## 2.7. Antimicrobial Activity

### 2.7.1. Bacterial Strains and Growth Conditions

Indicator strains were selected from culture collections and maintained at -80 °C in their growth medium supplemented with 20% glycerol (v/v). They were cultured twice in growth conditions specific to each strain (Table 3).

<table><caption>**Table 3.** Growth media and growth conditions of food-spoilage and foodborne bacteria used as indicator strains.</caption><thead><tr><th>Indicator Strains</th><th>Strain</th><th>Culture Collection</th><th>Growth Medium</th><th>Growth Conditions (°C)</th><th>Microaerophilic Condition</th></tr></thead><tbody><tr><td>Bacillus cereus</td><td>CIP 66.24T</td><td>CIP<sup>1</sup></td><td>BHI<sup>5</sup></td><td>30</td><td>No</td></tr><tr><td>Bacillus subtilis</td><td>CIP 52.65T</td><td>CIP</td><td>BHI</td><td>30</td><td>No</td></tr><tr><td>Brochothrix thermosphacta</td><td>DSMZ 20171T</td><td>DSMZ<sup>2</sup></td><td>BHI</td><td>25</td><td>No</td></tr><tr><td>Carnobacterium maltaromaticum</td><td>CIP 103158</td><td>CIP</td><td>MRS<sup>6</sup></td><td>30</td><td>Yes</td></tr><tr><td>Cronobacter sakazakii</td><td>CIP 103183T</td><td>CIP</td><td>BHI</td><td>30</td><td>No</td></tr><tr><td>Escherichia coli</td><td>DSMZ 1103</td><td>DSMZ</td><td>BHI</td><td>37</td><td>No</td></tr><tr><td>Enterococcus faecalis</td><td>CIP 103015</td><td>CIP</td><td>BHI</td><td>37</td><td>No</td></tr><tr><td>Leuconostoc mesenteroides</td><td>ATCC 14935</td><td>ATCC<sup>3</sup></td><td>MRS</td><td>30</td><td>Yes</td></tr><tr><td>Listeria innocua</td><td>CIP 80.11</td><td>CIP</td><td>BHI</td><td>37</td><td>No</td></tr></tbody></table>
<table><caption>**Table 3. Cont.**</caption><thead><tr><th>Indicator Strains</th><th>Strain</th><th>Culture Collection</th><th>Growth Medium</th><th>Growth Conditions (°C)</th><th>Microaerophilic Condition</th></tr></thead><tbody><tr><td>Listeria ivanovii</td><td>CIP 78.42T</td><td>CIP</td><td>BHI</td><td>37</td><td>No</td></tr><tr><td>Listeria monocytogenes</td><td>DSMZ 19094</td><td>DSMZ</td><td>BHI</td><td>37</td><td>No</td></tr><tr><td>Listeria monocytogenes</td><td>ATCC 35152</td><td>ATCC</td><td>BHI</td><td>37</td><td>No</td></tr><tr><td>Mucor circinelloides</td><td>UBOCC-A-112187</td><td>UBOCC<sup>4</sup></td><td>MEA<sup>7</sup></td><td>25</td><td>No</td></tr><tr><td>Mucor racemosus</td><td>UBOCC-A-111130</td><td>UBOCC</td><td>MEA</td><td>25</td><td>No</td></tr><tr><td>Penicillium solitum</td><td>UBOCC-A-111055</td><td>UBOCC</td><td>MEA</td><td>25</td><td>No</td></tr><tr><td>Pseudomonas aeruginosa</td><td>ATCC 27853</td><td>ATCC</td><td>BHI</td><td>30</td><td>No</td></tr><tr><td>Pseudomonas fluorescens</td><td>CIP 69.13</td><td>CIP</td><td>BHI</td><td>30</td><td>No</td></tr><tr><td>Staphylococcus aureus</td><td>DSMZ 13661</td><td>DSMZ</td><td>BHI</td><td>37</td><td>No</td></tr><tr><td>Staphylococcus aureus</td><td>CIP 76.25</td><td>CIP</td><td>BHI</td><td>37</td><td>No</td></tr><tr><td>Salmonella ser. Enteritidis</td><td>CIP 82.97</td><td>CIP</td><td>BHI</td><td>37</td><td>No</td></tr><tr><td>Salmonella ser.</td><td>CIP 104115</td><td>CIP</td><td>BHI</td><td>37</td><td>No</td></tr><tr><td>Typhimurium</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Yarrowia lipolytica</td><td>UBOCC-A-211004</td><td>UBOCC</td><td>YPD<sup>8</sup></td><td>25</td><td>No</td></tr><tr><td>Yersinia enterocolitica</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>subsp. enterocolitica</td><td>DSMZ 4780</td><td>DSMZ</td><td>BHI</td><td>26</td><td>No</td></tr></tbody></table>

## 2.7.2. Determination of Inhibitory Properties of LAB Isolates against Foodborne Pathogenic Bacteria and Spoilage Bacteria

### Agar Spot Test

An agar spot test was used to determine the antimicrobial potential of LAB isolates against the indicator strains (Table 3). Five microliters of an overnight culture (approximately $10^8$ CFU/mL) of each LAB isolate were spotted onto BHI agar plates (15 g/L agar) and left to dry for 20 min at room temperature under sterile conditions. The plates were then overlaid with 10 mL of soft agar (8 g/L) seeded with $10^6$ CFU/mL of the indicator strain (Table 3). Plates were incubated overnight at the optimum growth temperature of each indicator strain and checked for inhibition halos. As positive control, *Lactobacillus sakei* CIP 104,494 was used because of its large antimicrobial spectrum [41].

### Preparations of Cell-Free Supernatant (CFS) and Neutralized CFS (NCFS)

Cell-free supernatants (CFSs) were obtained by centrifugation at 8600 g for 10 min at 4 °C to separate bacterial cells and supernatants. The neutralized CFS (NCFS) were prepared by adjusting the pH of CFS to 7.0 with NaOH 1M to exclude the antimicrobial effect of organic acids. Samples were heated at 100 °C for 10 min to inhibit the enzyme activity. The antimicrobial properties of CFSs and NCFSs were tested on the indicator strains by the agar well diffusion assay.

### Agar Well Diffusion Assay

BHI agar was seeded with $10^6$ CFU/mL of indicator strains (Table 3), homogenized and poured into sterile Petri dishes. After drying under sterile conditions, sterile glass wells were placed on agar plates. Aliquots of non-diluted CFS or NCFS produced by LAB isolates were used to fill in the wells (100 µL). Plates were incubated overnight at the optimum growth temperature of each indicator strain. Inhibition halos were observed and recorded.

### Inhibitory Effect during Co-Cultivation

The ability of LAB isolates to inhibit the growth of *S. aureus* DSMZ 13,661 and *S. ser. Typhimurium* CIP 104,115 was further studied during co-cultivation assays in a double concentrated BHI media (pH 7) [42]. The co-cultivation process consisted of two different culture batches. Batch 1 corresponded to the co-inoculation of *S. aureus* DSMZ 13,661 with the selected LAB isolate and batch 2, the co-inoculation of *S. ser. Typhimurium* CIP 104,115 with the LAB isolate. Growth control was performed using a single

<footnote n="1">¹ Collections of Institut Pasteur (CIP), ² Leibniz Institute Deutsche Sammlung von Mikroorganismen und Zellkulturen,</footnote>

<footnote n="3">³ American Type Culture Collection, ⁴ Cultures de l'Université de Bretagne Occidentale, ⁵ brain heart infusion,</footnote>

<footnote n="6">⁶ de Man Rogosa and Sharpe, ⁷ malt extract agar, ⁸ yeast potato dextrose.</footnote>
culture for each indicator strain or the LAB isolate. An inoculum of 1% (v/v) of LAB isolate overnight culture (10⁹ CFU/mL in BHI broth) and 1% (v/v) of indicator strain overnight culture (10⁹ CFU/mL in BHI broth) were seeded to perform Batch 1 or Batch 2 co-culture assay. All batches were incubated at 28 °C for 48 h. The CFS from each batch culture was collected every 2 h. Antimicrobial activities were detected by the agar well diffusion assay as described above.

# 3. Results and Discussion

## 3.1. UMSC Bacterial Diversity Studies

Cheese manufacture behavior, processing environment and starter cultures are known to harbor an indigenous microbiota that influences organoleptic features on artisanal cheeses. The bacterial composition of starter cultures supports organoleptic and biochemical changes during ripening [43]. Revealing information on microbial consortia are essential for artisanal cheesemakers, considering the willingness of the consumer market to require locally manufactured products with different organoleptic properties [9,43,44]. Part of a daily routine, the studied UMSC is used for soft-cheese manufacture. Its microbial transfer is applied during rennet addition, a few hours after milk enrichment with acidifying starters. This cheese-making practice, combined with a high diversity of microbial activities, is key to developing gustatory characteristics. In this study, the microbial consortium of the UMSC was deciphered using culture-independent and culture-dependent based methods.

First, the UMSC microbial diversity was assessed using three different culture-independent methods: (i) analysis of 16S rRNA gene fragment sequences, (ii) *purR* gene sequence analysis, (iii) whole-genome sequencing analysis. Afterward, two culture-based methods were used on bacterial isolates recovered from the UMSC: MALDI-TOF MS and full-length 16S rRNA gene sequencing analysis.

### 3.1.1. Analysis of 16S rRNA Gene Fragment Sequences

An assessment of the microbial diversity started with targeted amplicon sequencing of 16S rRNA gene fragments. Read counts (total of 65,447 reads for V3–V4 and 61,383 reads for V1–V3) were filtered and relative abundances were calculated for each OTU. One dominant genus, *Lactococcus* sp., encompassing 99.9% of relative abundance was identified for each 16S rRNA gene region. With a clustering threshold of three, two OTUs were obtained among the *Lactococcus* genus.

### 3.1.2. *purR* Gene Sequence Analysis

The housekeeping gene *purR* was then studied to complete the assessment of microbial diversity. It was previously selected after pan-genome analysis for the differentiation of clades beyond the *L. lactis* subspecies level [7]. The *purR* gene encodes for a transcriptional regulator of purine biosynthesis. The use of a single nucleotide inference method highlighted five amplicon sequence variants (ASVs) after raw data processing using the DADA2 algorithm (Table 4). A total of 108,455 reads were obtained with an exact *purR* amplicon size of 435 nucleotides for each clustered ASV [45].

<table><caption>**Table 4.** Amplicon sequence variants (ASVs) matrix count (%) determined with *purR* gene analysis.</caption><thead><tr><td></td><td>Abundance (%)</td></tr></thead><tbody><tr><td>ASV_1</td><td>99.86</td></tr><tr><td>ASV_2</td><td>0.13</td></tr><tr><td>ASV_3</td><td>0.005</td></tr><tr><td>ASV_4</td><td>0.005</td></tr><tr><td>ASV_5</td><td>0.002</td></tr></tbody></table>

A close strain similarity was revealed and each ASV variant was represented by a single nucleotide sequence type. Nucleotide variant analysis showed that the most abundant sequence (ASV_1) was 100% identical to the corresponding reference sequence *L. lactis* subsp. *lactis* UC08 (GenBank: CP015903.1)
or UC11 (GenBank: CP015904.1). A matrix table highlighted nucleotide substitutions between each ASVs, from 1 up to 18 different bases (Figure 1A). The nucleotide translation into amino acid sequences identified three different ASVs (see matrix table, Figure 1B). Amino acids sequences of variants ASV_1, ASV_2 and ASV_3 were identical. Variants ASV_4 and ASV_5 showed three amino acid substitutions compared to other ASV (or the reference sequence). A comparison between ASV_4 and ASV_5 showed one amino acid substitution. Substitutions include amino acid changes in position 7: Glu (ASV_1), Ala (ASV_4) and Ser (ASV_5), position 21: Thr (ASV_1) and Pro (ASV_4 and ASV_5) and position 22: Lys (ASV_1) and Gln (ASV_4 and ASV_5). Dendrograms were constructed using nucleotide alignments (Figure 1A) or protein alignments (Figure 1B).

<figure><img src="image_3.png" /><figcaption>**Figure 1.** Dendrograms showing sequence similarities of *L. lactis* strains based on a cluster analysis of the *purR* amplicon nucleotide variants (A) or amino acid variants (B). *purR* genes were sequenced from the undefined mixed starter culture (UMSC) global genomic DNA extraction. The tree was constructed using Phylogeny.fr [46]. Matrix tables highlight differences in the composition of nucleotide (A) or amino acid (B) sequences.</figcaption></figure>

These five ASVs sequence reads were all affiliated to the species level (100% identity) as *L. lactis* using BLAST searches against a local database including any sequences of lactococcal genomes part of the NCBI database (built on 10 February 2020). Interestingly, the phylogenetic analysis of *purR* did not identify *L. lactis* to the subspecies level, as it did for Frantzen et al. [7]. In their work, primers designed by Frantzen et al. were used for relative quantification of the microbial community in starter cultures. Medium amplicon size was then required. Therefore, the *purR* amplicon size (435 nucleotides) was shorter than the full-length *purR* gene (816 nucleotides), identified in the whole-genome sequences. The targeted amplicon *purR* showed an alignment between position 347 to 782. A BLAST search performed on the full-length *purR* gene affiliated the genome sequences to the subspecies level (100% identity) as *L. lactis* subsp. *lactis*. These results suggested that new primers devised to amplify longer fragments of the *purR* gene would be needed to explore *L. lactis* diversity to the subspecies level in the UMSC. Another approach was attempted using the *eps* operon, both located on plasmids and chromosomes [7]. Unfortunately, irregular amplicon sizes were obtained, and it was impossible to differentiate isolates to the species/subspecies level (data not shown).

The ASV method infers the biological sequences in the sample before the introduction of amplification and sequencing errors. It distinguishes sequence variants differing by as little as one nucleotide [47]. The OTU assignment methods gather clusters of reads that differ by less than a fixed clustering threshold of three [28]. In a previous study, Lactococci strains were confirmed to have multiple copies
of the 16S rRNA gene. These copies harbor intra-genomic heterogeneity, which makes it difficult to differentiate bacteria to the species/subspecies level using 16S rRNA amplicon [45,48]. ASV methods better discriminate ecological patterns capturing all biological variation [31]. Furthermore, DNA sequences of protein-coding genes are more effective than 16S rRNA genes when distinguishing between very closely related bacteria [7,49]. The targeted amplicon comparison highlighted a noteworthy difference between results obtained with fragments of the 16S rRNA gene (2 OTUs) or with the *purR* gene (5 ASVs). The housekeeping *purR* gene seems to be the best candidate since it enabled the differentiation of clades beyond the species level.

### 3.1.3. Whole-Genome Sequencing Analysis

The ANI value is used to estimate the genetic distance between isolates, which reflects the sequence identities of the conserved regions between two genomes [50]. Whole-genome sequencing and assemblies were performed. The pairwise OrthoANI values of six in-house *Lactococcus* and one reference genome (*L. Lactis* subsp. *lactis* IL1403) were calculated to assess their genetic relatedness. One cluster was retrieved (Supplementary Materials Figure S2). Based on the guidelines set out by Richter and Rosselló-Móra, *L. lactis* subsp. *lactis* strains are a different species to subspecies *cremoris* strains possessing an ANI value below 95% [51]. Here, the ANI identity was >99.8% for all of the *L. lactis* genomes, suggesting they all belong to the *lactis* subspecies group.

### 3.1.4. Bacterial Library

Second, a culture-dependent based method was performed to complete the study of the UMSC bacterial diversity. A bacterial library was established from the UMSC with a set of 1616 bacterial isolates, according to 13 growth conditions (Table 5).

<table><caption>**Table 5.** Growth conditions and number of collected isolates used to build the bacterial library.</caption><thead><tr><th></th><th>Growth Conditions</th><th></th><th>Number of Collected Isolates</th></tr></thead><tbody><tr><td rowspan="7">Aerobic atmosphere</td><td rowspan="7">28 °C</td><td>BHI agar</td><td>128</td></tr><tr><td>MRS agar</td><td>192</td></tr><tr><td>M17 agar</td><td>192</td></tr><tr><td>MRS agar + 2% NaCl</td><td>96</td></tr><tr><td>MRS agar + 4% NaCl</td><td>96</td></tr><tr><td>MRS agar + 6.5% NaCl</td><td>96</td></tr><tr><td>BHI agar + 6.5% NaCl</td><td>96</td></tr><tr><td rowspan="3">Microaerophilic atmosphere</td><td rowspan="3">28 °C</td><td>BHI agar</td><td>48</td></tr><tr><td>MRS agar</td><td>192</td></tr><tr><td>M17 agar</td><td>192</td></tr><tr><td rowspan="3">Aerobic atmosphere</td><td>22 °C</td><td>BHI agar</td><td>96</td></tr><tr><td>37 °C</td><td>MRS agar</td><td>96</td></tr><tr><td>44 °C</td><td>MRS agar</td><td>96</td></tr><tr><td colspan="2">Total</td><td></td><td>1616</td></tr></tbody></table>

### 3.1.5. MALDI-TOF MS and Full-Length 16S rRNA Gene Sequencing Analysis

Two representative isolates of each growth condition were randomly selected for MALDI-TOF MS and full-length 16S rRNA gene sequencing characterizations. A set of 26 bacterial isolates (n = 26/1616) was analyzed as a sample that could represent the bacterial diversity of each growth condition. Four species identifications were obtained: *Lactococcus lactis* (L) (n = 18/26), *Lactococcus raffinolactis* (R) (n = 1/26), *Enterococcus faecalis* (E) (n = 5/26) and *Staphylococcus warneri* (S) (n = 2/26). The culture-dependent analysis showed similarities in the UMSC diversity and agreed with culture-independent studies for the *L. lactis* identification. Results also highlighted the presence of other species besides *L. lactis*. These latter could be part of the remaining 0.1% of relative abundance outlined by 16S rRNA targeted amplicon sequencing. In terms of genus richness and effective number
of genera, the differences observed among 16S rRNA analytical method were most likely a result of unclassified reads in the bioinformatic pipeline [52].

## 3.2. UMSC Strains Metabolic and Functional Characteristics

### 3.2.1. UMSC Genotype-Phenotype Associations

API 50 CH systems were used to explore carbohydrates metabolism. Height carbohydrate metabolism profiles were retrieved from the previous 26 examined LAB isolates (Supplementary Materials Figure S1). Three dairy *L. lactis* type strains were also tested (Table 6). Carbohydrate metabolism plays a key role in the growth of *Lactococcus* species. On average, 12.5% of genes are associated with carbohydrate metabolism and involved in glycolysis/gluconeogenesis and the pentose phosphate pathway (PPP) [53]. Four *L. lactis* phenotypic groups were obtained: (A) for isolates L1 to L6, L10 to L13 and L15 to L17 (*n* = 13); (B) for isolates L7 and L8 (*n* = 2); (C) for isolates L14 and L18 (*n* = 2) and (D) for isolate L9 (*n* = 1) (Table 6).

<table><caption>**Table 6.** Phenotypic characteristics of UMSC *Lactococcus* isolates. Phenotypic characteristics of two strain types, ATCC 19,435 and ATCC 19,257; data were extracted from Buron-Moles et al. and Yu et al. [14,53]. Phenotypic data of the type strain IL1403 was extracted from Siezen et al. [54].</caption><thead><tr><th>Carbohydrate</th><th>A Group<br>(n = 13)<sup>a</sup></th><th>B Group<br>(n = 2)<sup>b</sup></th><th>C Group<br>(n = 2)<sup>c</sup></th><th>D Group<br>(n = 1)<sup>d</sup></th><th>ATCC<br>19435<sup>e</sup></th><th>ATCC<br>19257<sup>f</sup></th><th>IL1403<sup>h</sup></th></tr></thead><tbody><tr><td>Strain origin</td><td>Raw milk</td><td>Raw milk</td><td>Raw milk</td><td>Raw milk</td><td>Milk</td><td>Cream</td><td>Dairy</td></tr><tr><td>Control</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>D-ribose</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>-</td><td>+</td></tr><tr><td>L-arabinose</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>D-xylose</td><td>+</td><td>-</td><td>+</td><td>+</td><td>+</td><td>-</td><td>-</td></tr><tr><td>D-galactose</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>-</td><td>+</td></tr><tr><td>D-glucose</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td></tr><tr><td>D-fructose</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td></tr><tr><td>D-mannose</td><td>+</td><td>+</td><td>+</td><td>+</td><td>w<sup>g</sup></td><td>+</td><td>+</td></tr><tr><td>D-mannitol</td><td>-</td><td>+</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td></tr><tr><td>D-sorbitol</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Arbutin</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>-</td><td>+</td></tr><tr><td>Salicin</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>-</td><td>+</td></tr><tr><td>D-cellobiose</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>-</td><td>+</td></tr><tr><td>D-maltose</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>-</td><td>+</td></tr><tr><td>D-lactose</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>w</td></tr><tr><td>D-melibiose</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>D-sucrose</td><td>-</td><td>+</td><td>+</td><td>+</td><td>-</td><td>-</td><td>-</td></tr><tr><td>D-trehalose</td><td>+</td><td>+</td><td>+</td><td>+</td><td>+</td><td>-</td><td>+</td></tr><tr><td>D-raffinose</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Starch</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>D-tagatose</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Percentage of<br>carbohydrates (n = 49)<br>utilized per strain (%)</td><td>33</td><td>45</td><td>35</td><td>35</td><td>29</td><td>10</td><td>29</td></tr></tbody></table>

Lₓ: *Lactococcus lactis* species, <sup>a</sup> A Group: L1 to L6, L10 to L13 and L15 to L17, <sup>b</sup> B Group: L7 and L8, <sup>c</sup> C Group: L14 and L18, <sup>d</sup> D Group: L9, <sup>e</sup> *L. lactis* subsp. *lactis* ATCC 19,435 [14,52], <sup>f</sup> *L. lactis* subsp. *cremoris* ATCC 19,257 [14,52], <sup>g</sup> weakly positive, <sup>h</sup> *L. lactis* subsp. *lactis* IL1403 [54].

A genotype-phenotype association was devised to study possible relations between genes and phenotypes observed in the API 50 CH systems. A representative set of six *L. lactis* isolates were selected for WGS studies: three isolates were part of the A group (L1, L2, L3), one isolate was part of the B group (L8), one isolate was part of the C group (L14) and one isolate was part of the D group (L9). An Illumina iSeq® platform (L1, L2, L8, L9, L14 genomes) and a PacBio® RS II technology (L3 genome) were used for genome assemblies. The role of *L. lactis* as an important industrial starter strain is mainly due to the rapid conversion of lactose to lactic acid. The KEGG pathways (https://www.genome.jp/kegg-bin) of Lactococci cell metabolism summarize two possible catabolic ways for lactose utilization: (i) hydrolysis of lactose to give α-D-glucose and D-galactose
performed by a β-galactosidase (EC 3.2.1.23), and (ii) the conversion into a lactose 6-phosphate (Lac-6P) by a lactose phosphotransferase (EC 2.7.1.207). The Lac-6P is subsequently hydrolyzed by 6-phospho-β-galactosidase (EC 3.2.1.85) to give D-galactose-6P (Gal-6P), and then Gal-6P is catabolized by the tagatose 6-phosphate pathway [55]. Genome sequence analysis highlighted both catabolic pathways with genes encoding these key enzymes. Conversely, phenotypic results indicated that the B group was able to degrade D-sorbitol, although it lacked all gene families found to be significantly associated with this degradation. Buron-Moles et al. reported the same results with dairy *L. lactis* isolates [14].

Regarding D-tagatose degradation, phenotypic results observed in the API 50 CH systems revealed that isolates of the B group could use this carbohydrate source whilst other isolates (or type strains) could not. D-tagatose is hydrolyzed by a tagatose-6-phosphate kinase (EC 2.7.1.144) to give D-tagatose-6P (Tag-6P). Tag-6P is subsequently hydrolyzed by a tagatose-1,6-diphosphate aldolase (EC 4.1.2.40) to give Tag-6P₂, and then it enters in the PPP. These enzymes were found in each studied genome even though the B group was the only one able to grow on D-tagatose. By studying genotype-phenotype associations between gene families and their carbohydrate metabolism, Buron-Moles et al. reported that D-tagatose can only be degraded by a few LAB members of *Lactobacillus* and *Pediococcus* strains [14]. Siezen et al. also reported the inability of *Lactococcus* strains to grow on D-tagatose [54]. Here, the *L. lactis* isolates were highlighted for their uncommon D-tagatose degradation.

D-xylose units are components of the β-1,4-xylan, a hemicellulosic polysaccharide of plant materials [54]. The β-1,4-xylosidase (EC 3.2.1.37) is a key enzyme involved in the β-1,4-xylan degradation to obtain D-xylose oligomers. D-xylose is then converted by a xylose isomerase (EC 5.3.1.5) to form D-xylulose (D-xyl). D-xyl is then catabolized by a xylulokinase (EC 2.7.1.17) to enter in the PPP. Genes encoding these enzymes were found in each analyzed genome. The carbohydrate metabolism study demonstrated that isolates were able to use D-xylose, except for those that were part of the B group, the *cremoris* subspecies type strain and the dairy strain *L. lactis* IL1403. Similarly, phenotypic results highlighted isolates able to grow on D-sucrose except for group A and *L. lactis* type strains, which were sucrose-negative. A previous study also reported the inability of dairy *Lactococcus* strains to utilize D-sucrose [14], unlike strains from plant origin that could grow on the D-sucrose [18,54]. Genotypic relationships showed no genes encoding for D-sucrose degradation. As noted by Buron-Moles et al., these findings emphasize that phenotypes cannot be always inferred from the presence/absence of specific genes but require the explicit integration of phenotypic information [14].

Phenotypic features in carbohydrate metabolism appeared to reflect the functional diversity observed for each *L. lactis* isolate. Moreover, the percentage of carbohydrates utilized per isolate was slightly higher (33% to 45%) in comparison with *L. lactis* type strains (29%). Their metabolism potential exhibited differences and seemed to be richer and diversified. Interestingly, the presence of some genes did not always imply a degradation of the associated carbohydrate in the API 50 CH systems. It could be explained by a progressive pseudogenization process (gene inactivation) linked to their current niche specification. The *L. lactis* isolates came from a farm environment before being introduced into raw milk and then into the UMSC. A comparative study revealed that an *L. lactis* gene set related to degradation products such as arabinose and xylose originated from plants [54]. Genomes of subspecies *lactis* KF147 and KF282, which were of plant origin, possessed genes set for the uptake and conversion of typical cell wall degradation products such as arabinose and xylose. All of the analyzed genomes possessed a key enzyme for L-arabinose degradation even though phenotypic profiles showed an inability to grow on arabinose (Table 6). Moreover, B group isolates possessed the *araD* gene, encoding an L-ribulose-5-phosphate 4-epimerase (EC 5.1.3.4), but it was annotated as a pseudogene. A mutation may have divided the gene into two-parts, and thus the *araD* gene became inactive. These results suggested that the human-made environment triggered microbial strains to genomic function adaptations to their dairy environment. UMSC isolates progressively appeared to be domesticated since they no longer could grow in a variety of plant-derived carbohydrates due to pseudogenization or transposon insertions in genes [56]. These genomic changes could be worsened
due to the "backslopping" process in which artisans are continuously transferring microbial material from the previous fermentation product to start a new batch [15].

### 3.2.2. Antimicrobial Activities

The antimicrobial activity of bacterial isolates against foodborne, food-spoilage and pathogenic bacteria was evaluated by the agar well diffusion assay. All isolates that were part of the bacterial library (n = 1616) were tested. A total of 34/1616 isolates showed antimicrobial activity against up to five indicator strains (Table 7 and Supplementary Materials Table S1). Isolates 14 and 15 showed the largest antimicrobial spectrum inhibiting the growth of the five indicator strains (S. ser. Typhimurium CIP 104115, *S. aureus* DSMZ 13661, *S. aureus* CIP 76.25, *E. faecalis* CIP 103,015 and *L. innocua* CIP 80.11). *S. aureus* DSMZ 13,661 was sensitive and showed the largest inhibition zones. S. ser. Typhimurium CIP 104,115 was the only Gram-negative bacteria to be sensitive toward isolates 14 and 15. The antimicrobial activity of the *L. lactis*-isolate 14 was further studied in co-cultivation assays against the two indicator strains: *S. aureus* DSMZ 13,661 or S. ser. Typhimurium CIP 104115.

<table><caption>**Table 7.** Antimicrobial spectrum of CFS and NCFS produced by LAB isolates against related LAB, food-spoilage and foodborne pathogens by the well diffusion assay.</caption><thead><tr><th rowspan="2">Tested Isolates</th><th colspan="2">S. ser. Typhimurium CIP 104115</th><th colspan="2">S. aureus DSMZ 13661</th><th colspan="2">C. maltaromaticum CIP 103158</th><th colspan="2">S. ser. Enteritidis CIP 82.97</th><th colspan="2">S. aureus CIP 76.25</th><th colspan="2">E. faecalis CIP 103015</th><th colspan="2">L. innocua CIP 80.11</th></tr><tr><th>CFS</th><th>NCFS</th><th>CFS</th><th>NCFS</th><th>CFS</th><th>NCFS</th><th>CFS</th><th>NCFS</th><th>CFS</th><th>NCFS</th><th>CFS</th><th>NCFS</th><th>CFS</th><th>NCFS</th></tr></thead><tbody><tr><td>1</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>+</td><td>-</td></tr><tr><td>3</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>+</td><td>-</td></tr><tr><td>4</td><td>-</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>+</td><td>-</td><td>+</td><td>-</td></tr><tr><td>5</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>+</td><td>-</td><td>+</td></tr><tr><td>6</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>+</td><td>-</td><td>+</td><td>-</td></tr><tr><td>8</td><td>-</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>+</td><td>-</td><td>+</td><td>-</td><td>+</td><td>-</td></tr><tr><td>9</td><td>-</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>+</td><td>-</td></tr><tr><td>10</td><td>-</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>+</td><td>-</td></tr><tr><td>11</td><td>-</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>+</td><td>-</td><td>+</td><td>-</td></tr><tr><td>12</td><td>-</td><td>-</td><td>+</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>+</td><td>-</td></tr><tr><td>13</td><td>-</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>+</td><td>-</td><td>+</td><td>-</td></tr><tr><td>14</td><td>+</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>+</td><td>-</td><td>+</td><td>-</td><td>+</td><td>+</td></tr><tr><td>15</td><td>+</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>+</td><td>-</td><td>+</td><td>-</td><td>+</td><td>+</td></tr><tr><td>16</td><td>-</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>+</td><td>-</td><td>+</td><td>-</td></tr><tr><td>17</td><td>-</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>+</td><td>-</td></tr><tr><td>18</td><td>+</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>+</td><td>-</td><td>+</td><td>-</td></tr><tr><td>20</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>+</td><td>-</td><td>+</td><td>-</td></tr><tr><td>23</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>+</td><td>-</td></tr><tr><td>24</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>+</td><td>-</td></tr><tr><td>25</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>+</td><td>-</td></tr><tr><td>26</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>+</td><td>-</td></tr><tr><td>27</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>+</td><td>-</td></tr><tr><td>30</td><td>-</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>+</td><td>-</td></tr><tr><td>32</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>+</td><td>-</td></tr><tr><td>33</td><td>-</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>+</td><td>-</td></tr><tr><td>34</td><td>+</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>+</td><td>-</td><td>-</td><td>-</td><td>+</td><td>-</td></tr></tbody></table>

CFS: cell-free supernatant, NCFS: neutralized cell-free supernatant (adjusted to pH 7), +: presence of an inhibition halo (>5 mm), -: no inhibition halo.

During co-culture assays, the growth of each bacterium was monitored by bacterial counts on a specific medium (data not shown). Growth controls were obtained using bacterial counts performed on single cultures of each strain (indicators and isolate 14). Temporary growth slowdowns were observed during co-cultures and compared with their respective growth controls. Growth slowdowns in co-culture assays were highlighted in the exponential growth phase of S. ser. Typhimurium (between 8 and 20 h) and in the stationary growth phase of *S. aureus* DSMZ 13,661 (between 26 and 34 h). Their respective controls showed no growth disruption. The growth in the co-culture of isolate 14 remained similar to that of the growth control containing no indicator strain (data not shown).
To further examine isolate 14 antimicrobial activity, CFSs of each batch culture were collected during these growth slowdowns. Their antimicrobial activity was assessed by the agar well diffusion assay. Thus, CFSs were collected between 26 and 34 h of growth for *S. aureus* DSMZ 13,661 assays (mono and co-cultures). Their CFSs largest inhibition halo was obtained after 32 h of growth. Similarly, CFSs were collected between 8 and 20 h of growth for *S. ser.* Typhimurium CIP 104,115 assays (mono and co-cultures). Their CFSs largest inhibition halo was obtained after 16 h of growth (data not shown). Table 8 summarizes inhibition halo values (mm) of CFSs both obtained from mono and co-cultures after 32 h of growth for *S. aureus* assays and after 16 h of growth for *S. ser.* Typhimurium assays.

<table><caption>**Table 8.** Antimicrobial activities of CFSs collected from the culture of *L. lactis* subsp. *lactis*-isolate 14 cultivated with and without *S. ser.* Typhimurium CIP 104,115 or *S. aureus* DSMZ 13661. Plates were examined for clear inhibition halo (mm), evaluated by the CFS agar well diffusion assay. Largest inhibition halos are recorded below. CFSs were collected after 16 h of growth for *S. ser.* Typhimurium CIP 104,115 assays and after 32 h of growth for *S. aureus* DSMZ 13,661 assays.</caption><thead><tr><th>Tested CFS</th><th>S. ser.<br>Typhimurium<br>CIP 104115</th><th>Lactobacillus<br>Sakei CIP<br>104494</th><th>S. aureus<br>DSMZ 13661</th><th>L. innocua<br>CIP 80.11</th><th>L. lactis IL1403<sup>d</sup></th></tr></thead><tbody><tr><td><i>L. lactis</i> isolate 14</td><td>8.0 ± 0.7<sup>a</sup></td><td>12.0 ± 0.7</td><td>8.0 ± 0.0</td><td>24.0 ± 0.7</td><td>13.0 ± 1.4</td></tr><tr><td><i>S. ser.</i> Typhimurium<br>CIP 104115</td><td><sup>b</sup></td><td>/</td><td>/</td><td>/</td><td>/</td></tr><tr><td>Isolate 14 + <i>S. ser.</i><br>Typhimurium CIP<br>104115</td><td>11.0 ± 0.0</td><td>14.0 ± 1.4</td><td>ND<sup>c</sup></td><td>22.0 ± 1.4</td><td>15.0 ± 2.12</td></tr><tr><td><i>S. aureus</i> DSMZ 13661</td><td>/</td><td>10.0 ± 2.12</td><td>/</td><td>/</td><td>/</td></tr><tr><td>Isolate 14 + <i>S. aureus</i><br>DSMZ 13661</td><td>ND</td><td>11.0 ± 1.4</td><td>14.0 ± 0.0</td><td>18.0 ± 0.7</td><td>17.0 ± 0.7</td></tr></tbody></table>

CFSs of isolate 14, alone or in the presence of indicator strains, showed antimicrobial activity against *S. ser.* Typhimurium CIP 104115, *Lactobacillus sakei* CIP 104494, *S. aureus* DSMZ 13661, *L. innocua* CIP 80.11 and *L. lactis* subsp. *lactis* IL1403. CFSs of indicator strains showed no inhibition halo against themselves. These results highlighted an increase in inhibition halo sizes when isolate 14 was co-cultivated with *S. aureus* or *S. ser.* Typhimurium. The recorded pH values between 8 and 20 h and 26 and 34 h of incubation time ranging from 6.6 to 6.8 for each batch culture. Comparative studies described the growth of *S. aureus* under much more stringent conditions, at pH values of 5.25 and 4.48 [58] or 5.6 in milk [59]. The inhibition of *S. aureus* in buffered media by *Lactococcus* strains was also observed at pH 6.8 and cannot, therefore, be attributed to a pH drop since *S. aureus* can grow at lower pH values [59]. Mufandaedza et al. also reported the survival and the growth of *Salmonella* sp. in fermented product associated with low pH [60].

The web-based tool BAGEL4 was used to mine secondary metabolite gene clusters within the whole genome *L. lactis* L3 (isolate 14). A coding sequence sharing 51.35% identity with the *lcnB* gene was recovered. The gene is known to be involved in the production of the lactococcin B bacteriocin (LcnB) [61]. A sequence sharing 69.23% identity with the *lciB* gene was also recovered and is involved in the immunity toward LcnB. A putative export system encoded by analogs of *lcnC* and *lcnD* genes was also identified with 98.37% identities. This analysis highlighted the presence of four genes sharing identities with the lactococcin operon [62]. The LcnB affects membrane permeability, dissipating its potential and its pH gradient caused by leakage of ions. In *L. lactis* strains carrying the immunity protein, the LcnB is inactive [63]. The LcnB antimicrobial activity was detected against close-related Lactococci strains. None of the Lactococci strains LcnB producers could inhibit the growth of other Gram-positive or Gram-negative bacteria [64,65]. Isolate 14 inhibited the growth of *L. lactis* IL1403, a plasmid-free strain sensitive to bacteriocin. These results suggested that isolate 14 could produce an antimicrobial peptide close to the LcnB. Şahİngİl et al. also reported inhibition of *Salmonella* sp. by a lactococcin-like bacteriocin, produced by an *L. lactis* subsp. *lactis* strain [66]. That lactococcin activity

<footnote><sup>a</sup> Mean ± SD from duplicate determinations, <sup>b</sup>: no inhibition halo, <sup>c</sup> ND: not determined, <sup>d</sup> plasmid-free strain; indicator strain for lactococcin B-like bacteriocin [57].</footnote>
was observed during the early logarithmic phase (6–7 h of growth) and had maximal activity at 18 h of incubation time, as with the results shown here for isolate 14.

**4. Conclusions**

This study described the microbial community of an undefined starter culture with its richness in *Lactococcus* strains, naturally present in the raw milk microbiota. Five *L. lactis* variants were identified, and five carbohydrate metabolism profiles were highlighted in the *Lactococcus* genus. These *Lactococcus* strains differ in their phenotypic properties from the strains commonly used as industrial starters. Moreover, some *L. lactis* isolates presented an antimicrobial activity detected by cultural methods. These bacteria could act as a biocontrol agent when combined with other procedures in the so-called hurdle technology. Their GRAS status, their origin from an ecological niche and their role in flavor formation support the use of the UMSC in cheese manufacture procedures. Furthermore, considering the culture characteristics and their relationship with product quality, further studies could be helpful to determine their effects on technological aptitude and sensory properties of the cheese.

**Supplementary Materials:** The following are available online at http://www.mdpi.com/2304-8158/9/5/622/s1, Figure S1. Carbohydrate metabolism tested using the API 50 CH system (BioMérieux, Marcy l'Etoile, France). A blue dot highlights a positive response to carbohydrate degradation. Figure S2. Pairwise average nucleotide identity (ANI) values across six *Lactococcus* genomes and the reference genome *L. lactis* subsp. *lactis* IL1403. Colors in the heatmap represent pairwise ANI values, with a gradient from blue (low identity) to green (high identity). Table S1. Antimicrobial spectrum of LAB candidate isolates against related LAB, food-spoilage and foodborne pathogens evaluated by the agar well diffusion assay. Plates were recorded for inhibition halos (mm).

**Author Contributions:** Conceptualization, S.S. (Sabrina Saltaji); Data curation, S.S. (Sabrina Saltaji); Formal analysis, S.S. (Sabrina Saltaji) and O.R.; Funding acquisition, S.S. (Sophie Sablé), F.T., S.D. and R.C.; Investigation, S.S. (Sabrina Saltaji) and R.C.; Methodology, S.S. (Sabrina Saltaji); Project administration, S.S. (Sophie Sablé) and S.D.; Resources, V.S. and F.T.; Software, O.R.; Validation, S.S. (Sabrina Saltaji) and O.R.; Visualization, S.S. (Sabrina Saltaji); Writing—original draft, S.S. (Sabrina Saltaji); Writing—review and editing, R.C. All authors read and approved the final manuscript.

**Funding:** This research received no external funding.

**Acknowledgments:** We thank Stéphane Aymerich for the *L. lactis* IL1403 strain. The INRAE MIGALE bioinformatics platform is thanked for access to computational resources and data storage facilities. The LABGeM (CEA/Genoscope and CNRS UMR8030), the France Génomique and French Bioinformatics Institute national infrastructures (funded as part of Investissement d'Avenir program managed by Agence Nationale pour la Recherche, contracts ANR-10-INBS-09 and ANR-11-INBS-0013) are acknowledged for support within the MicroScope annotation platform. This work was financed by Flores de terroirs company and "Association Nationale Recherche Technologie" (convention CIFRE N° 2017/0745).

**Conflicts of Interest:** The authors declare no conflict of interest.

**References**

1. 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]

2. De Angelis, M.; de Candia, S.; Calasso, M.P.; Faccia, M.; Guinee, T.P.; Simonetti, M.C.; Gobbetti, M. Selection and use of autochthonous multiple strain cultures for the manufacture of high-moisture traditional Mozzarella cheese. *Int. J. Food Microbiol.* **2008**, *125*, 123–132. [CrossRef]

3. De Pasquale, I.; Di Cagno, R.; Buchin, S.; de Angelis, M.; Gobbetti, M. Use of autochthonous mesophilic lactic acid bacteria as starter cultures for making Pecorino Crotonese cheese: Effect on compositional, microbiological and biochemical attributes. *Food Res. Int.* **2019**, *116*, 1344–1356. [CrossRef]

4. Kupiec, B.; Revell, B. Speciality and artisanal cheeses today: The product and the consumer. *Br. Food J.* **1998**, *100*, 236–243. [CrossRef]

5. Antonsson, M.; Molin, G.; Ardö, Y. *Lactobacillus* strains isolated from Danbo cheese as adjunct cultures in a cheese model system. *Int. J. Food Microbiol.* **2003**, *85*, 159–169. [CrossRef]
6. Bissonnette, F.; Labrie, S.; Deveau, H.; Lamoureux, M.; Moineau, S. Characterization of Mesophilic Mixed Starter Cultures Used for the Manufacture of Aged Cheddar Cheese. *J. Dairy Sci.* **2000**, *83*, 620–627. [CrossRef]

7. Frantzen, C.A.; Kot, W.; Pedersen, T.B.; Ardö, Y.M.; Broadbent, J.R.; Neve, H.; Hansen, L.H.; Dal Bello, F.; Østlie, H.M.; Kleppen, H.P.; et al. Genomic Characterization of Dairy Associated *Leuconostoc* Species and Diversity of *Leuconostocs* in Undefined Mixed Mesophilic Starter Cultures. *Front. Microbiol.* **2017**, *8*, 132. [CrossRef] [PubMed]

8. Carafa, I.; Clementi, F.; Tuohy, K.; Franciosi, E. Microbial evolution of traditional mountain cheese and characterization of early fermentation cocci for selection of autochtonous dairy starter strains. *Food Microbiol.* **2016**, *53*, 94–103. [CrossRef] [PubMed]

9. 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*, 1–12. [CrossRef]

10. Koczura, M.; Martin, B.; Turille, G.; De Marchi, M.; Kreuzer, M.; Berard, J. Milk composition, but not cheese properties, are impaired the day after transhumance to alpine pastures. *Int. Dairy J.* **2019**, *99*, 1–9. [CrossRef]

11. Wouters, J.T.M.; Ayad, E.H.E.; Hugenholtz, J.; Smit, G. Microbes from raw milk for fermented dairy products. *Int. Dairy J.* **2002**, *12*, 91–109. [CrossRef]

12. Taïbi, A.; Dabour, N.; Lamoureux, M.; Roy, D.; LaPointe, G. Evaluation of the genetic polymorphism among *Lactococcus lactis* subsp. *cremoris* strains using comparative genomic hybridization and multilocus sequence analysis. *Int. J. Food Microbiol.* **2010**, *144*, 20–28.

13. Smid, E.J.; Kleerebezem, M. Production of aroma compounds in lactic fermentations. *Annu. Rev. Food Sci. Technol.* **2014**, *5*, 313–326. [CrossRef] [PubMed]

14. Buron-Moles, G.; Chailyan, A.; Dolejs, I.; Forster, J.; Mikš, M.H. Uncovering carbohydrate metabolism through a genotype-phenotype association study of 56 lactic acid bacteria genomes. *Appl. Microbiol. Biotechnol.* **2019**, *103*, 3135–3152. [CrossRef] [PubMed]

15. Steensels, J.; Gallone, B.; Voordeckers, K.; Verstrepen, K.J. Domestication of Industrial Microbes. *Curr. Biol.* **2019**, *29*, R381–R393. [CrossRef]

16. Torres Manno, M.; Zuljan, F.; Alarcón, S.; Esteban, L.; Blancato, V.; Espariz, M.; Magni, C. Genetic and phenotypic features defining industrial relevant *Lactococcus lactis*, *L. cremoris* and *L. lactis* biovar. *diacetylactis* strains. *J. Biotechnol.* **2018**, *282*, 25–31. [CrossRef]

17. Wu, C.; Huang, J.; Zhou, R. Genomics of lactic acid bacteria: Current status and potential applications. *Crit. Rev. Microbiol.* **2017**, *43*, 393–404. [CrossRef]

18. Bayjanov, J.R.; Starrenburg, M.J.; van der Sijde, M.R.; Siezen, R.J.; van Hijum, S.A. Genotype-phenotype matching analysis of 38 *Lactococcus lactis* strains using random forest methods. *BMC Microbiol.* **2013**, *13*, 1–11. [CrossRef]

19. Ceapa, C.; Lambert, J.; van Limpt, K.; Wels, M.; Smokvina, T.; Knol, J.; Kleerebezem, M. Correlation of *Lactobacillus rhamnosus* Genotypes and Carbohydrate Utilization Signatures Determined by Phenotype Profiling. *Appl. Environ. Microbiol.* **2015**, *81*, 5458–5470. [CrossRef]

20. Park, S.-H.; Itoh, K.; Kikuchi, E.; Niwa, H.; Fujisawa, T. Identification and Characteristics of Nisin Z-Producing *Lactococcus lactis* subsp. *lactis* Isolated from Kimchi. *Curr. Microbiol.* **2003**, *46*, 0385–0388. [CrossRef]

21. Yang, E.; Fan, L.; Jiang, Y.; Doucette, C.; Fillmore, S. Antimicrobial activity of bacteriocin-producing lactic acid bacteria isolated from cheeses and yogurts. *AMB Express* **2012**, *2*, 1–12. [CrossRef] [PubMed]

22. O'Connor, P.M.; Kuniyoshi, T.M.; Oliveira, R.P.; Hill, C.; Ross, R.P.; Cotter, P.D. Antimicrobials for food and feed; a bacteriocin perspective. *Curr. Opin. Biotechnol.* **2020**, *61*, 160–167. [CrossRef] [PubMed]

23. Wilson, P.D.G.; Brocklehurst, T.F.; Arino, S.; Thuault, D.; Jakobsen, M.; Lange, M.; Farkas, J.; Wimpenny, J.W.T.; van Impe, J.F. Modelling microbial growth in structured foods: Towards a unified approach. *Int. J. Food Microbiol.* **2002**, *73*, 275–289. [CrossRef]

24. Dumarché, C.; Leclercq, S. Method for Preparing a Leaven From Unpasteurized Milk, Leaven Obtained by Means of This Method and Use of This Leaven for Producing a Milk Product. U.S. Patent 20090772755, 12 May 2011.

25. Leser, T.D.; Amenuvor, J.Z.; Jensen, T.K.; Lindecrona, R.H.; Boye, M.; Møller, K. Culture-Independent Analysis of Gut Bacteria: The Pig Gastrointestinal Tract Microbiota Revisited. *Appl. Environ. Microbiol.* **2002**, *68*, 673–690. [CrossRef]
26. Weisburg, W.G.; Barns, S.M.; Pelletier, D.A.; Lane, D.J. 16S ribosomal DNA amplification for phylogenetic study. *J. Bacteriol.* **1991**, *173*, 697–703. [CrossRef]

27. Escudié, F.; Auer, L.; Bernard, M.; Mariadassou, M.; Cauquil, L.; Vidal, K.; Maman, S.; Hernandez-Raquet, G.; Combes, S.; Pascal, G.; et al. FROGS: Find, Rapidly, OTUs with Galaxy Solution. *Bioinformatics* **2018**, *34*, 1287–1294. [CrossRef]

28. Mahé, F.; Rognes, T.; Quince, C.; de Vargas, C.; Dunthorn, M. Swarm: Robust and fast clustering method for amplicon-based studies. *PeerJ* **2014**, *2*, 1–13. [CrossRef]

29. Rognes, T.; Flouri, T.; Nichols, B.; Quince, C.; Mahé, F. VSEARCH: A versatile open source tool for metagenomics. *PeerJ* **2016**, *4*, 1–22. [CrossRef]

30. Pruesse, E.; Quast, C.; Knittel, K.; Fuchs, B.M.; Ludwig, W.; Peplies, J.; Glöckner, F.O. SILVA: A comprehensive online resource for quality checked and aligned ribosomal RNA sequence data compatible with ARB. *Nucleic Acids Res.* **2007**, *35*, 7188–7196. [CrossRef]

31. 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]

32. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. *Mol. Biol. Evol.* **2016**, *33*, 1870–1874. [CrossRef] [PubMed]

33. Altschul, S.F.; Madden, T.L.; Schäffer, A.A.; Zhang, J.; Zhang, Z.; Miller, W.; Lipman, D.J. Gapped BLAST and PSI-BLAST: A new generation of protein database search programs. *Nucleic Acids Res.* **1997**, *25*, 3389–3402. [CrossRef] [PubMed]

34. Chin, C.-S.; Alexander, D.H.; Marks, P.; Klammer, A.A.; Drake, J.; Heiner, C.; Clum, A.; Copeland, A.; Huddleston, J.; Eichler, E.E.; et al. Nonhybrid, finished microbial genome assemblies from long-read SMRT sequencing data. *Nat. Methods* **2013**, *10*, 563–569. [CrossRef] [PubMed]

35. Bankevich, A.; Nurk, S.; Antipov, D.; Gurevich, A.A.; Dvorkin, M.; Kulikov, A.S.; Lesin, V.M.; Nikolenko, S.I.; Pham, S.; Prjibelski, A.D.; et al. SPAdes: A New Genome Assembly Algorithm and Its Applications to Single-Cell Sequencing. *J. Comput. Biol.* **2012**, *19*, 455–477. [CrossRef]

36. Vallenet, D.; Calteau, A.; Dubois, M.; Amours, P.; Bazin, A.; Beuvin, M.; Burlot, L.; Bussell, X.; Fouteau, S.; Gautreau, G.; et al. MicroScope: An integrated platform for the annotation and exploration of microbial gene functions through genomic, pangenomic and metabolic comparative analysis. *Nucleic Acids Res.* **2020**, *48*, D579–D589. [CrossRef]

37. Lee, I.; Ouk Kim, Y.; Park, S.-C.; Chun, J. OrthoANI: An improved algorithm and software for calculating average nucleotide identity. *Int. J. Syst. Evol. Microbiol.* **2016**, *66*, 1100–1103. [CrossRef]

38. Kanehisa, M.; Sato, Y.; Kawashima, M.; Furumichi, M.; Tanabe, M. KEGG as a reference resource for gene and protein annotation. *Nucleic Acids Res.* **2016**, *44*, D457–D462. [CrossRef]

39. Van Heel, A.; Jong, A.; Song, C.; Viel, J.; Kok, J.; Kuipers, O. BAGEL4: A user-friendly web server to thoroughly mine RiPPs and bacteriocins. *Nucleic Acids Res.* **2018**, *46*, W278–W281. [CrossRef]

40. Bolotin, A.; Wincker, P.; Mauger, S.; Jaillon, O.; Malarne, K.; Weissenbach, J.; Ehrlich, S.D.; Sorokin, A. The Complete Genome Sequence of the Lactic Acid Bacterium *Lactococcus lactis* ssp. *lactis* IL1403. *Genome Res.* **2001**, *11*, 731–753. [CrossRef]

41. Dubourg, G.; Elsawi, Z.; Raoult, D. Assessment of the in vitro antimicrobial activity of *Lactobacillus* species for identifying new potential antibiotics. *Int. J. Antimicrob. Agents* **2015**, *46*, 590–593. [CrossRef]

42. Hwanhlem, N.; Ivanova, T.; Haertlé, T.; Jaffrès, E.; Dousset, X. Inhibition of food-spoilage and foodborne pathogenic bacteria by a nisin Z-producing *Lactococcus lactis* subsp. *lactis* KT2W2L. *LWT—Food Sci. Technol.* **2017**, *82*, 170–175. [CrossRef]

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

44. Turbes, G.; Linscott, T.D.; Tomasino, E.; Waite-Cusic, J.; Lim, J.; Meunier-Goddik, L. Evidence of terroir in milk sourcing and its influence on Cheddar cheese. *J. Dairy Sci.* **2016**, *99*, 5093–5103. [CrossRef] [PubMed]

45. Lamei, S.; Hu, Y.O.O.; Olofsson, T.C.; Andersson, A.F.; Forsgren, E.; Vásquez, A. Improvement of identification methods for honeybee specific Lactic Acid Bacteria; future approaches. *PLoS ONE* **2017**, *12*, 1–12. [CrossRef]
46. Dereeper, A.; Guignon, V.; Blanc, G.; Audic, S.; Buffet, S.; Chevenet, F.; Dufayard, J.-F.; Guindon, S.; Lefort, V.; Lescot, M.; et al. Phylogeny.fr: Robust phylogenetic analysis for the non-specialist. *Nucleic Acids Res.* **2008**, *36*, W465–W469. [CrossRef]

47. Callahan, B.J.; McMurdie, P.J.; Holmes, S.P. Exact sequence variants should replace operational taxonomic units in marker-gene data analysis. *ISME J.* **2017**, *11*, 2639–2643. [CrossRef]

48. Pillidge, C.J.; Sheehy, L.M.; Shihata, A.; Pu, Z.-Y.; Dobos, M.; Powell, I.B. Intragenomic 16S rRNA gene heterogeneity in *Lactococcus lactis* subsp. *cremoris*. *Int. Dairy J.* **2009**, *19*, 222–227. [CrossRef]

49. Ndoye, B.; Rasolofo, E.A.; LaPointe, G.; Roy, D. A review of the molecular approaches to investigate the diversity and activity of cheese microbiota. *Dairy Sci. Technol.* **2011**, *91*, 495–524. [CrossRef]

50. Sun, Z.; Zhang, W.; Guo, C.; Yang, X.; Liu, W.; Wu, Y.; Song, Y.; Kwok, L.Y.; Cui, Y.; Menghe, B.; et al. Comparative Genomic Analysis of 45 Type Strains of the Genus *Bifidobacterium*: A Snapshot of Its Genetic Diversity and Evolution. *PLoS ONE* **2015**, *10*, e1–14. [CrossRef]

51. Richter, M.; Rosselló-Móra, R. Shifting the genomic gold standard for the prokaryotic species definition. *Proc. Natl. Acad. Sci. USA* **2009**, *106*, 19126–19131. [CrossRef]

52. Plummer, E.; Twin, J.; Bulach, D.M.; Garland, S.M.; Tabrizi, S.N. A Comparison of Three Bioinformatics Pipelines for the Analysis of Preterm Gut Microbiota using 16S rRNA Gene Sequencing Data. *J. Proteom. Bioinform.* **2015**, *8*, 283–291. [CrossRef]

53. Yu, J.; Song, Y.; Ren, Y.; Qing, Y.; Liu, W.; Sun, Z. Genome-level comparisons provide insight into the phylogeny and metabolic diversity of species within the genus *Lactococcus*. *BMC Microbiol.* **2017**, *17*, 1–10. [CrossRef] [PubMed]

54. Siezen, R.J.; Starrenburg, M.J.C.; Boekhorst, J.; Renckens, B.; Molenaar, D.; van Hylckama Vlieg, J.E. Genome-Scale Genotype-Phenotype Matching of Two *Lactococcus lactis* Isolates from Plants Identifies Mechanisms of Adaptation to the Plant Niche. *Appl. Environ. Microbiol.* **2008**, *74*, 424–436. [CrossRef] [PubMed]

55. Benthin, S.; Nielsen, J.; Villadsen, J. Galactose Expulsion during Lactose Metabolism in *Lactococcus lactis* subsp. *cremoris* FD1 Due to Dephosphorylation of Intracellular Galactose 6-Phosphate. *Appl. Environ. Microbiol.* **1994**, *60*, 1254–1259. [CrossRef]

56. Kelly, W.J.; Ward, L.J.H.; Leahy, S.C. Chromosomal Diversity in *Lactococcus lactis* and the Origin of Dairy Starter Cultures. *Genome Biol. Evol.* **2010**, *2*, 729–744. [CrossRef]

57. Venema, K.; Dost, M.H.; Beun, P.A.; Haandrikman, A.J.; Venema, G.; Kok, J. The genes for secretion and maturation of lactococcins are located on the chromosome of *Lactococcus lactis* IL1403. *Appl. Environ. Microbiol.* **1996**, *62*, 1689–1692. [CrossRef]

58. Sutherland, J.P.; Bayliss, A.J.; Roberts, T.A. Predictive modelling of growth of *Staphylococcus aureus*: The effects of temperature, pH and sodium chloride. *Int. J. Food Microbiol.* **1994**, *21*, 217–236. [CrossRef]

59. Alomar, J.; Loubiere, P.; Delbes, C.; Nouaille, S.; Montel, M.C. Effect of *Lactococcus garvieae*, *Lactococcus lactis* and *Enterococcus faecalis* on the behaviour of *Staphylococcus aureus* in microfiltered milk. *Food Microbiol.* **2008**, *25*, 502–508. [CrossRef]

60. Mufandaedza, J.; Viljoen, B.C.; Feresu, S.B.; Gadaga, T.H. Antimicrobial properties of lactic acid bacteria and yeast-LAB cultures isolated from traditional fermented milk against pathogenic *Escherichia coli* and *Salmonella enteritidis* strains. *Int. J. Food Microbiol.* **2006**, *108*, 147–152. [CrossRef]

61. Boeckmann, B.; Bairoch, A.; Apweiler, R.; Blatter, M.-C.; Estreicher, A.; Gasteiger, E.; Martin, M.J.; Michoud, K.; O'Donovan, C.; Phan, I.; et al. The SWISS-PROT protein knowledgebase and its supplement TrEMBL in 2003. *Nucleic Acids Res.* **2003**, *31*, 365–370. [CrossRef]

62. Fath, M.J.; Kolter, R. ABC transporters: Bacterial exporters. *Microbiol. Mol. Biol. Rev.* **1993**, *57*, 995–1017. [CrossRef]

63. Venema, K.; Abee, T.; Haandrikman, A.J.; Leenhouts, K.J.; Kok, J.; Konings, W.N.; Venema, G. Mode of Action of Lactococcin B, a Thiol-Activated Bacteriocin from *Lactococcus lactis*. *Appl. Environ. Microbiol.* **1993**, *59*, 1041–1048. [CrossRef] [PubMed]

64. Kojic, M.; Strahinic, I.; Fira, D.; Jovcic, B.; Topisirovic, L. Plasmid content and bacteriocin production by five strains of *Lactococcus lactis* isolated from semi-hard homemade cheese. *Can. J. Microbiol.* **2006**, *52*, 1110–1120. [CrossRef] [PubMed]
65. Morgan, S.; Ross, R.P.; Hill, C. Bacteriolytic activity caused by the presence of a novel lactococcal plasmid encoding lactococcins A, B, and M. *Appl. Environ. Microbiol.* **1995**, *61*, 2995–3001. [CrossRef]

66. Şahingil, D.; İşleroğlu, H.; Yildirim, Z.; Akçelik, M.; Yildirim, M. Characterization of lactococcin BZ produced by *Lactococcus lactis* subsp. *lactis* BZ isolated from boza. *Turk. J. Biol.* **2011**, *35*, 21–33.

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

© 2020 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 (http://creativecommons.org/licenses/by/4.0/).