Next Article in Journal
β-Lactamase Producing Escherichia coli Encoding blaCTX-M and blaCMY Genes in Chicken Carcasses from Egypt
Next Article in Special Issue
Involvement of Versatile Bacteria Belonging to the Genus Arthrobacter in Milk and Dairy Products
Previous Article in Journal
Effect of Palm-Based Shortenings of Various Melting Ranges as Animal Fat Replacers on the Physicochemical Properties and Emulsion Stability of Chicken Meat Emulsion
Previous Article in Special Issue
Changes in Key Aroma Compounds and Esterase Activity of Monascus-Fermented Cheese across a 30-Day Ripening Period
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Foods
Volume 12
Issue 3
10.3390/foods12030599
foods-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

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

Comparative Genomic Analysis Reveals the Functional Traits and Safety Status of Lactic Acid Bacteria Retrieved from Artisanal Cheeses and Raw Sheep Milk

by
Ilias Apostolakos
1,
Spiros Paramithiotis
2 and
Marios Mataragas
1,*
1
Department of Dairy Research, Institution of Technology of Agricultural Products, Hellenic Agricultural Organization “DIMITRA”, 3 Ethnikis Antistaseos St., 45221 Ioannina, Greece
2
Department of Food Science and Human Nutrition, Agricultural University of Athens, 75 Iera Odos St., 11855 Athens, Greece
*
Author to whom correspondence should be addressed.
Foods 2023, 12(3), 599; https://doi.org/10.3390/foods12030599
Submission received: 20 December 2022 / Revised: 21 January 2023 / Accepted: 28 January 2023 / Published: 1 February 2023
(This article belongs to the Special Issue Cheese: Chemistry, Physics and Microbiology)
Browse Figures
Versions Notes

Abstract

:
Lactic acid bacteria (LAB) are valuable for the production of fermented dairy products. We investigated the functional traits of LAB isolated from artisanal cheeses and raw sheep milk, assessed their safety status, and explored the genetic processes underlying the fermentation of carbohydrates. Lactiplantibacillus plantarum had the largest and more functional genome compared to all other LAB, while most of its protein-encoding genes had unknown functions. A key finding of our analysis was the overall absence of acquired resistance genes (RGs), virulence genes (VGs), and prophages, denoting that all LAB isolates fulfill safety criteria and can be used as starter or adjunct cultures. In this regard, the identified mobile genetic elements found in LAB, rather than enabling the integration of RGs or VGs, they likely facilitate the uptake of genes involved in beneficial functions and in the adaptation of LAB in dairy matrices. Another important finding of our study was that bacteriocins and CAZymes were abundant in LAB though each species was associated with specific genes, which in turn had different activity spectrums and identified applications. Additionally, all isolates were able to metabolize glucose, lactose, maltose, and sucrose, but Lactiplantibacillus plantarum was strongly associated with the fermentation of rhamnose, mannose, cellobiose, and trehalose whereas Levilactobacillus brevis with the utilization of arabinose and xylose. Altogether these results suggest that to fully exploit the beneficial properties of LAB, a combination of strains as food additives may be necessary. Interestingly, biological processes involved in the metabolism of carbohydrates that are not of direct interest for the dairy industry may yield valuable metabolites or activate pathways associated with beneficial health effects. Our results provide useful information for the development of new probiotic artisanal cheeses and probiotic starter cultures.

1. Introduction

Lactic acid bacteria (LAB) are a diverse group of Gram-positive, non-sporulating, cocci or rods, which are microaerophilic or anaerobic. They are catalase-negative and produce lactic acid as the major end-product of carbohydrate fermentation [1]. LAB form two groups based on the end fermentation products, namely the homo- and hetero-fermentative groups. Hetero-fermentative LAB are further subdivided into facultatively and obligate fermentative species [2]. LAB play an important role in fermentation, producing a variety of metabolites and substances that are responsible for the sensorial properties and preservation of the final products [3]. They are of great economic importance for the production of fermented foods and for their use as starter cultures.
The use of lactic acid bacteria in the dairy industry is a well-established approach to enhance the organoleptic characteristics of fermented milk products. Cheese is one of the most popular fermented dairy products and is the most important product of the dairy industry. It is produced from milk by the action of rennet and LAB, which are added as starter cultures. The rennet breaks down the casein into peptides and amino acids, which are further metabolized by the LAB to produce energy and the main end products, lactic acid and carbon dioxide [4]. The lactic acid causes a decrease in the pH of the cheese, which is the main factor responsible for its preservation. LAB are also responsible for the production of a variety of volatile compounds and other substances that contribute to the flavor of the cheese [4,5,6,7]. Artisanal cheeses are traditional dairy products that are produced by small-scale farmers. Artisanal cheeses are usually produced by using a limited number of starter cultures and therefore are characterized by a high level of heterogeneity in their microbial composition and diversity. Artisanal cheeses constitute a complex system of microbial interactions between milk bacteria, starter cultures, and indigenous bacteria of the dairy environment [5].
With regard to their microbiological safety, LAB have the Qualified Presumption of Safety (QPS) status by the European Food Safety Authority (EFSA) and are generally recognized as safe (GRAS) in the USA [8]. In addition to technological properties, the use of LAB as food additives has increased in recent years due to their beneficial effects on human health. Some LAB strains have been used as probiotics, which are defined as live microorganisms that, when administered in adequate amounts, confer a health benefit to the host [3,9].
Whole genome sequencing (WGS) is a powerful tool for the investigation of bacterial species, as it provides a high-resolution view of the genetic diversity of bacteria [10]. The WGS technology has been applied to the study of the diversity of LAB species, including Lactobacillus acidophilus, Lactiplantibacillus plantarum, and Lacticaseibacillus casei [11,12]. The genomes of LAB species have been sequenced and analyzed in several studies, and the functions of the genes in these genomes have been characterized. However, the underlying genetic basis of LAB species that drives their functional traits and fermentation patterns is still not fully explored. Additionally, few studies have investigated the LAB of Greek artisanal cheeses and raw sheep milk [13,14], which limits the possibility to exploit the beneficial properties of these bacteria. In this study, we used WGS technology to conduct a comparative genomic analysis of 51 LAB strains, belonging to six different LAB species and retrieved from traditional Greek artisanal cheeses and raw sheep milk. The scope of the study was to investigate the functional traits of the LAB collection, assess their safety status, and explore the genetic processes involved in and related to the fermentation of several carbohydrates.

2. Materials and Methods

2.1. Bacterial Strains Used in the Study

The LAB used throughout this study (n = 51) were isolated from sheep milk, artisanal Feta, and Kefalograviera cheeses [13,14], and deposited in the culture collection of the Dairy Research Department (DRD) of Hellenic Agricultural Organization “DIMITRA” (ELGO-DIMITRA). Storage and culture conditions were performed as described by Tsigkrimani et al. [14]. The collection included the following LAB species: Lactiplantibacillus plantarum (Lb. plantarum, n = 16), Levilactobacillus brevis (Lb. brevis, n = 11), Lactococcus lactis subsp. lactis (Lc. lactis lactis, n = 9), Pediococcus pentosaceus (Pd. pentosaceus, n = 5), Leuconostoc mesenteroides (Ln. mesenteroides, n = 5), Latilactobacillus curvatus (Lb. curvatus, n = 3), and Lactococcus lactis subsp. cremoris (Lc. lactis cremoris, n = 2).

2.2. Whole Genome Sequencing, Assembly, and Quality Control

Extraction of DNA took place with the aid of the GenElute Bacterial Genomic DNA Kit (Sigma-Aldrich, St. Louis, MO, USA) according to the instructions of the manufacturer. Quality of the DNA was assessed by agarose gel electrophoresis [15] and quantification by Qubit 2.0 (ThermoFisher Scientific, Waltham, MA, USA). Library preparation was performed as follows: DNA was randomly fragmented by sonication; the ends were polished, A-tailed and ligated with Illumina’s sequencing adapters before PCR amplification using P5 and P7 primers. The AMPure XP system (Beckman Coulter, Brea, CA, USA) was employed for purification and size selection of the PCR products. The Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA) was used for assessment of library size; quantification was performed by qPCR. Sequencing of the qualified libraries was performed by Novogene Genomics Service (Novogene Co., Cambridge, UK) on the Illumina Novaseq 6000 platform (Illumina, San Diego, CA, USA) (2 × 150 bp). The FastQC v.0.11 software (Babraham Bioinformatics, Cambridge, UK) that is available in the KBase platform [16,17] was used for quality assessment of the adapter-free raw reads. Polishing of the raw reads and de novo assembling into contigs were performed with the Unicycler assembler and Pilon, available on the PATRIC v3.6.8 web platform (accessed on 1 January 2022) [18,19,20]. Organization of the contigs into scaffolds was performed by the Multi-Draft based Scaffolder (MeDuSa) v1.6 [21]. Then, scaffolds were ordered and oriented on the basis of the complete reference genomes present in the NCBI database (https://www.ncbi.nlm.nih.gov/ accessed on 10 January 2022); namely Lb. plantarum SK151, Lb. brevis NPS-QW-145, Lc. lactis LAC460, Pd. pentosaceus ATCC 25745, Ln. mesenteroides SRCM102733, Lb. curvatus JCM 1096 = DSM 20019, and Lc. cremoris subsp. cremoris KW2. Quality assessment of the contigs and scaffolds, in terms of completeness (≥95%) and contamination (≤5%), was performed with the CheckM tool v1.21 [22], available on the PATRIC v3.6.8 web platform. The Skew Index Test (SkweIT) v1.0 [23] was used to assess possible mis-assemblies after scaffolding.

2.3. In Silico Typing and Comparative Genomic Analysis

QUAST [24] was used to assess the assembled genomes’ quality, whereas Kraken2 [25] and TYGS [26] were used to identify LAB species. The genomes were annotated using PROKKA [27]. Functional annotation and subsystem analysis of ORFs was performed with the COG database [28] via the eggnog-mapper tool [29]. The presence of clustered Regularly Interspaced Short Palindromic Repeats (CRISPRs) was assessed with the CRISPRCasFinder tool [30], while prophages integrated in LAB genomes were detected with PHASTER [31]. Abricate [32] was used to assess the existence of resistance genes (RGs), virulence genes (VGs), mobile genetic elements (MGEs), and plasmids using the Resfinder [33], VFDB [34], MobileElementFinder [35], and PlasmidFinder [36] databases, respectively. Bacteriocins were identified with BAGEL4 [34]. Presence of Carbohydrate-active enzymes (CAZymes) was assessed with the Run_dbcan V3 tool [37], whereas Traitar [38] was used to predict the phenotypic traits of LAB, in addition to those validated in a previous study performing laboratory experiments [14]. To reveal important differences in LAB pertaining to the gene-content of CAZymes, bacteriocins, and MGEs, as well as to their phenotypic traits, cluster heatmaps were computed using as input the presence/absence patterns of these features. Validation of the heatmap clusters was done with statistical analysis for overrepresentation of these features in specific LAB species. Moreover, we assessed overrepresented Gene Ontology (GO)-terms as follows. The LAB pangenome was annotated with the respective GO terms with eggnog [29]. Next, we explored enriched GO-terms with Scoary [39] and visualized the output with GoFigure! [40]. The pangenome and core-genome alignment were computed with Roary [41]. Proteins were clustered together if their amino acid sequence identity was ≥75%. The threshold for the core-genome was set at 90%. We used FastTree [42] to generate a phylogenetic tree, which was visualized with iTOL [43]. Galaxy platform [44] was used for running some of the tools.

2.4. Statistical Analysis

To establish which COG categories, MGEs, CAZymes, bacteriocins, phenotypic traits, and GO terms were significantly overrepresented in each LAB species, we used the respective presence/absence data matrices as input to Scoary [39]. The significance level (alpha) was set at 0.05. The p-values were adjusted with Benjamini–Hochberg’s method for multiple comparisons correction. Statistical comparison in the R programming language [45] was also conducted for the basic genomic features (genome length, GC content, and number of CDS) of LAB strains with One-way ANOVA (ANalysis Of VAriance), followed by post-hoc Tukey HSD (Honestly Significant Difference). Throughout the study, LAB strains were analyzed at the species level, thus the two Lc. lactis cremoris strains were analyzed together with the nine Lc. lactis lactis strains as one species (Lc. lactis).

3. Results and Discussion

3.1. Assembly Statistics and Subsystem Analysis

The strain identification presented by Tsigkrimani et al. [13] was verified by the taxonomic classification with Kraken2 and TYGS. Details of species identification, isolation source, and assembly statistics for each genome are presented in Table 1. The basic characteristics of LAB genomes were evaluated statistically. The average genome length of our collection was 2.65 Mb, with Lb. plantarum having the largest genome (2.98 Mb, on average). Genome lengths were overall not significantly different in pairwise comparisons except for Lb. plantarum, which had a significantly larger genome (0.85 Mb, on average) than Lb. brevis (2.20 Mb on average). GC content was 44.13% on average and we found that Lb. plantarum and Lc. lactis had richer GC content than Pd. pentosaceus but no difference among them. Lastly, Lc. lactis and Lb. curvatus had the highest number of CDS on average, respectively, 3217 and 3087, and both had significantly more CDS than the other LAB species, especially compared to Pd. pentosaceus which had the smallest CDS count (1935, on average).
A set of CDS that collectively implement a specific biological process or structural complex is defined as a subsystem [46]. Subsystem analysis revealed the presence of 21 enriched COG categories (Figure 1). The COG category of unknown function (S) was the most enriched one, with an average of 451 CDS across LAB genomes. In this category as well as in 16/21 of COG categories overall, Lb. plantarum was significantly more enriched compared to all other LAB genomes (Figure 1). These findings indicate that Lb. plantarum has a highly functional genome and at the same time, a large part of its CDS have unknown functions that warrant further characterization [47]. Other enriched COG categories were those of transcription (K), carbohydrate metabolism and transport (G), replication and recombination (L), translation (J), and amino acid metabolism (E) having on average 209, 174, 168, 164, and 161 CDS, respectively (Figure 1). The least enriched categories were those of extracellular structures (W), present only in Ln. mesenteroides and Pd. pentosaceus (each with two CDS on average), followed by RNA processing (A), which was present only in Lc. lactis genomes (one CDS on average).

3.2. Phylogenetic Analysis and Assessment of RGs and VGs

We found 15,585 COG clusters in the pangenome, whereas 38 COGs comprised the core-genome. The phylogenetic analysis and resulting tree indicated the genetic relatedness of LAB (Figure 2). Clusters were formed according to the species of LAB with no observable overlaps. Moreover, formation of subclusters was observed within the species clusters, which were formed according to their originating food (Feta cheese, Kefalograviera cheese, and raw sheep milk), especially for Ln. mesenteroides and Lc. lactis.
Analysis for presence of resistance and virulence genes (RGs, VGs) showed absence of VGs in all LAB, except for Lc. lactis genomes (n = 11), which contained the lmrD gene, an inherent, chromosomally encoded efflux pump that encodes for resistance to lincosamides in Streptomyces lincolnensis and Lc. lactis. In addition, it has been found that this gene confers resistance to bile acids in Lc. lactis [7,48].

3.3. Comparative Genomics

A clustered heatmap was generated from the presence/absence patterns of the bacteriocins, CAZymes, and MGEs. Three major clusters were formed: Lb. plantarum and Lc. lactis formed the largest cluster, whereas Lb. curvatus and Lb. brevis clustered together with 3/5 isolates of Ln. mesenteroides (Figure 3). The smallest cluster contained the remaining (2/5) Ln. mesenteroides isolates and Pd. pentosaceus. This finding suggests that although phylogenetically distinct, LAB isolates share similar functional traits, most probably due to their adaptation to the specific environment of dairy products [49].
Bacteriocins are small, ribosomally synthesized peptides that are produced by bacteria as a means of inhibiting the growth of other bacteria. They are a form of antimicrobial peptide and are often specific to a particular group of bacteria or even to a single species. Bacteriocins are typically active against closely related strains but can also be active against more distantly related microorganisms [50]. They can be classified based on their mode of action, targeting the bacterial cell wall, the cytoplasmic membrane, or even the DNA/RNA of the microbes. Bacteriocins have been tested in foods as natural preservatives to inhibit pathogenic and spoilage bacteria [51].
Bacteriocins were present in all LAB species except for Lb. brevis. Pediococcus pentosaceus was the most bacteriocin-enriched species; its isolates had at least one bacteriocin-encoding gene, an average of 35 genes, and 37 unique bacteriocins (Figure 3). Overrepresentation analysis with Scoary showed strong association [Odds ratio (OR) >> 1, p << 0.05)] of the species with Bac43, Bacteriocin_31 (bacA), Bavaricin_MN, and Divercin_RV41. These bacteriocins are active against E. faecalis, E. faecium, E. hirae, E. durans, and L. monocytogenes [52,53]. The second most enriched species was Ln. mesenteroides, followed by Lb. plantarum with an average of 18 and 10 bacteriocins, respectively. Interestingly, we found that Ln. mesenteroides was not significantly associated with particular bacteriocins as most of its bacteriocin-encoding genes (n = 27 unique) were shared and more frequently found in Pd. pentosaceus (Figure 3). In contrast, Lb. plantarum was strongly associated (OR >> 1, p << 0.05) with Plantaricins (plnA, plnE, plnF, and plnJ) and Acidocin_B (acdB). Plantaricins have been extensively studied as natural preservatives in dairy and meat products, for their role in the inhibition of persistent bacterial infections as well their protective effect in irritable bowel syndrome disease (IBS) and urinary tract infections (UTIs) [54]. Acidocin_B (acdB) was first purified from Lb. acidophilus strain M46 and shows strong activity against L. monocytogenes, Clostridium sporogenes, and Brochothrix thermosphacta but is generally inactive against LAB species, hence it does not compromise the microflora stability during food fermentation [55]. Furthermore, Lc. lactis and Lb. curvatus had the least bacteriocins on average (three and one, respectively), with Lb. curvatus being significantly associated with carnobacteriocin B2 (cbnB2), a peptide active against Listeria and Enterococcus spp. [56,57].
CAZymes, short for carbohydrate-active enzymes, are a group of enzymes that catalyze the degradation, modification, or synthesis of complex carbohydrates. They are classified based on the type of carbohydrate they act on as well as the type of reaction they catalyze, such as hydrolysis, transfer, or oxidation. CAZymes play a vital role in the breakdown and utilization of carbohydrates by microorganisms. They also have potential applications in industrial processes, such as biofuel production and the production of enzymes used in food and beverage production [3]. Moreover, combining pre- and probiotics results in beneficial effects, such as the inhibition of inflammatory processes and the reduction of cholesterol levels [58]. CAZymes were ubiquitous in LAB. Overall, LAB isolates contained 17 CAZymes on average, with Lc. lactis being the most enriched species (n = 19), followed by Lb. plantarum and Lb. brevis (n = 18, each). In total, 28 unique CAZymes belonging to five major groups were identified. The majority of the CAZymes belonged to the Glycosyltransferase (GT) group with GT2, GT24, and GT28 being present in all LAB isolates (n = 51). Moreover, glycoside hydrolase (GH) 51 was present in 50/51 isolates whereas Carbohydrate-binding module (CBM) 50 in 46/51. Scoary analysis showed that Lb. plantarum and Lc. lactis were significantly associated with GT5, GT26, GT35, as well as CBM48. Of note, CBMs catalyze long CAZymes, such as glycoside hydrolases (GHs), and this process plays a key role the catabolism of complex carbohydrates such as lactose and starch [3]. In contrast, Lb. brevis was exclusively associated with several CAZymes (n = 8) of the GH group. Collectively, the aforementioned CAZymes cooperatively contribute to the dietary carbohydrate deconstruction [59].
A prophage is a bacteriophage genome that integrates in the bacterial chromosome or exists as an extrachromosomal element within the bacterial cell. Prophages that integrate in bacterial genomes often harbor resistance or virulence genes that can be transferred to the host bacterium [60]. With regard to the prophage content, none of the isolates had intact prophage regions in their genomes. The existence of CRISPR/Cas systems can help to protect bacterial genomes from prophage integration [61]. In this regard, none of the LAB isolates had robust evidence (evidence level = 4) of CRISPR sequences and cas genes in their genomes.
Presence of plasmids was evaluated with plasmidFinder, the detailed results are presented in Table 2. Lactiplantibacillus plantarum and Lb. curvatus had the highest plasmid content with three plasmids on average, followed by Lc. lactis, Lb. brevis, and Ln. mesenteroides (two plasmids on average, each). The most abundant plasmid type in these strains was a rep28 plasmid, initially reported in Lb. plantarum strain P-8 (NCBI accession No: CP005948) [62]. A similar distribution pattern was observed for MGEs (Figure 3). ISLsa1 (IS30 family), the most abundant insertion sequence (IS) element, was significantly associated with Lb. brevis and Lb. plantarum. Interestingly, this IS was first reported in the meat-borne Latilactobacillus sakei strain 23K and is presumably involved in the translocation of genes responsible for resilience during the harsh conditions of food processing, such as fluctuating redox and oxygen levels [63]. In contrast, Lc. lactis was associated with ISLL6 (IS3 family), an inherent IS element that together with IS982 make up 1 in 10 genes encoded in on lactococcal plasmids. These IS elements have been linked with essential functional roles in dairy products, such as the production of bacteriocins, carbohydrate metabolism (CAZymes), and resistance to bacteriophage infection [64]. The predisposition of bacteria to harbor MGEs is essential in the evolution of bacterial pathogens. The presence of MGEs in the genome signifies the ability of the species to persistently harbor pathogenicity and microbial resistance factors [36,65]. Using robust alignment criteria, we showed overall absence of acquired VGs and RGs in our collection. This absence may indicate that MGEs in LAB from traditional Greek artisanal cheeses and raw sheep milk play a role in the adaptation of the bacteria in different dairy matrices [66] and are involved in beneficial functional roles such as the production of bioactive peptides and enzymes [64]. However, both aspects warrant further investigations to fully understand the function of MGEs.

3.4. Analysis of Phenotypic Traits

The phenotypic traits evaluated experimentally in a previous work [14] were coupled with Traitar analysis for additional phenotypic characteristics and showed that, irrespective of their species, all isolates can utilize sugars such as glucose, lactose, maltose, and sucrose (Figure 4). The clustered heatmap of phenotypic traits indicated an overlap of species clusters, suggestive of shared phenotypic profiles between the studied LAB species. As with the clustered heatmap of genomic features (Figure 3), three major phenotype clusters were formed, albeit with different composition (Figure 4). The largest cluster comprised isolates of Ln. mesenteroides and Lb. plantarum. Scoary analysis indicated that Lb. plantarum isolates were strongly associated (OR >> 1, p << 0.05) with the fermentation of multiple carbohydrates (rhamnose, mannose, cellobiose, and trehalose), as well as with the sorbitol and mannitol sugar alcohols (Figure 4). In contrast, Ln. mesenteroides was not significantly linked to any of the assessed phenotypic traits. The second largest cluster was formed by Lc. lactis, Pd. pentosaceus, and Lb. curvatus. All three species were not significantly associated with specific carbohydrate fermentation. Of note, Lc. lactis and Pd. pentosaceus showed significant association with utilization of starch and malonate, suggesting that these species could be used in both dairy and vegetable fermentation [67]. The last cluster was exclusively formed by Lb. brevis isolates, which were significantly associated with utilization of arabinose and xylose.
In order to explore the genomic processes involved in and related to the fermentation of particular carbohydrates, we analyzed LAB genomes to identify significantly enriched Gene Ontology (GO) terms, associated with the metabolism of carbohydrates (Supplementary File S1). For cellobiose, the most overrepresented GO term clusters were (a) the positive regulation of ribosome biogenesis (GO:0090070) and (b) the establishment of competence for transformation (GO:0030420). These findings support the notion that the metabolism of cellobiose is predominantly controlled at the transcriptional level [68] while uptake of DNA from the bacterial environment (transformation) might also play a role, e.g., by the uptake of MGEs that facilitate the exchange of genes involved in carbohydrate fermentation [69]. Pathways involved in the fermentation of cellobiose might not be of direct interest for dairy products, however, it was recently reported that the same pathways drive the production of epilactose, a novel and promising prebiotic sugar which occurs from the metabolism of lactose [70].
Moreover, we found that the lipid A biosynthesis (GO:0009245) was overrepresented in LAB and able to catabolize raffinose (Supplementary File S1). Of note, it was recently reported that this pathway is upregulated in the gut microbiome of obese adults that receive whole grain wheat (WGW) products rich in oligosaccharides. In the presence of oligosaccharides such as raffinose, the lipid A biosynthesis of bacteria is upregulated. This effect, together with other modulations in the gut microbiome, lead to the amelioration of nonalcoholic fatty liver disease in obese subjects [71]. The most abundant gene in raffinose-fermenting strains involved in this pathway was acpP (Acyl carrier protein), encoding for a protein-carrier of the growing fatty acid chain during fatty acid biosynthesis [72].
For ribose, we found that the most overrepresented GO cluster was the tRNA modification (GO:0006400) comprising 20 unique GO terms (Supplementary File S1). This finding corroborates the results of McLeod et al. [73] that uptake and metabolism of ribose in Lb. sakei and other LAB is predominantly regulated at the transcription level. Indeed, we found that that the enzyme-encoding genes dtd (D-tyrosyl-tRNA(tyr) deacylase) and truA (tRNA pseudouridylate synthase A (pseudouridylate synthase I) were abundant in ribose-fermenting LAB strains.
Furthermore, the pathway for breakdown of arabinose to xylulose 5-phosphate and other compounds (GO:0019569) was significantly overrepresented in strains able to ferment arabinose (Supplementary File S1). The typical glycolytic routes, e.g., for glucose, in bacteria, are the Embden-Meyerhof-Parnas, the Entner-Doudoroff (ED), and the oxidative pentose phosphate pathways. However, the metabolism pathway of pentoses, such as xylose and arabinose, involve isomerases and kinases to produce xylulose 5-phosphate [74]. In this context, the gene involved in this pathway was the adolase araD (L-ribulose-5-phosphate 4-epimerase), a key enzyme that links the arabinose metabolic pathway to the pentose phosphate pathway and allows the bacteria to use arabinose as an energy source [75]. This gene was abundant in Lb. brevis strains, which were strongly associated with the fermentation of arabinose and xylose.
Lastly, the extracellular polysaccharide biosynthetic process (GO:0045226) was enriched in LAB strains, mainly belonging to Lb. plantarum, that were able to metabolize mannitol (Supplementary File S1). Interestingly, LAB can divert carbon resources away from glycolysis and into production of extracellular polysaccharides (EPS) [76]. EPSs have a significant role in food industrial applications; they have been employed as thickeners, stabilizers, and gelling agents in food products, like yogurt, in which they enhance the sensorial properties and “mouthfeel”. Moreover, EPSs increase the gastrointestinal transit time, thus promoting the rate of colonization of probiotic bacteria [77].

4. Conclusions

Our analysis demonstrated that Lb. plantarum had the largest genome and at the same time had most of the COG categories significantly more enriched compared to all other LAB genomes. Of note, most COGs of Lb. plantarum, and to a lesser extent of the other LAB, had unknown functions. These findings indicate that Lb. plantarum has a highly functional genome and at the same time a large part of its CDS is unexplored and warrant further characterization [78]. Furthermore, we showed overall absence of acquired RGs, VGs, and prophages in all LAB. According to EFSA, all of the analyzed LAB species have the QPS status as long as “the strains do not harbor any acquired antimicrobial resistance genes to clinically relevant antimicrobials” [79]. In this context, our results suggest that LAB isolated from artisanal Feta, Kefalograviera, and raw sheep milk are expected to be safe for use as starter, adjunct cultures, or as food additives in general. In this context, the identified MGEs, rather than enabling the integration of pathogenic and virulence determinants, they are likely to facilitate the uptake of genes involved in beneficial functional roles, such as the production of bioactive peptides and enzymes as well as for the adaptation of LAB in various dairy matrices [64]. Bacteriocins and CAZymes were largely ubiquitous in LAB. However, each LAB species was significantly associated with specific genes, which in turn had different activity spectrums and identified applications. Altogether these findings suggest that to fully exploit the beneficial functional properties of LAB, a combination of strains as food additives may be necessary [80].
The ability to ferment multiple carbohydrates is a key trait of LAB that is essential for their industrial application. All isolates were able to metabolize important dairy carbohydrates. Of note, Lb. plantarum and Lb. brevis were strongly associated with the fermentation processes of specific carbohydrates that were less abundant in the other LAB. GO analysis shed light in biological pathways involved in and related to the fermentation these carbohydrates. Interestingly, biological processes involved in the metabolism of carbohydrates that are not of direct interest for the dairy industry (e.g., cellobiose), may yield valuable metabolites or activate pathways associated with beneficial health effects.
Our results provide useful information for the further development of new probiotic artisanal cheeses and for the design of probiotic starter cultures.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/foods12030599/s1, Supplementary File S1. Clusters of Gene Ontology (GO) biological processes that were overrepresented in isolates able to catabolize particular carbohydrates (cellobiose, raffinose, ribose, arabinose, and mannitol). The representative GO terms from each cluster are plotted in a two-dimensional semantic space where similar terms are placed closer to each other on the scatterplot. The size of each circle represents the number of unique GO terms found in the respective cluster.

Author Contributions

Conceptualization, M.M.; methodology, I.A. and M.M.; software, I.A.; formal analysis, I.A. and M.M.; investigation, I.A. and S.P.; resources, M.M.; data curation, I.A. and M.M.; writing—original draft, I.A.; writing—review and editing, I.A., S.P. and M.M.; supervision, M.M.; project administration, M.M.; funding acquisition, M.M. All authors have read and agreed to the published version of the manuscript.

Funding

European Union and Greek national funds; RESEARCH–CREATE–INNOVATE (T1EDK-02087).

Data Availability Statement

This Whole Genome Shotgun project (BioProject number PRJNA846557, PRJNA847016, PRJNA847025, PRJNA847015, PRJNA847019, PRJNA847022, PRJNA916921) has been deposited at DDBJ/ENA/GenBank under the accession numbers JAMQJD000000000 to JAMQJS000000000 (Lb. plantarum), JAMRWH000000000 to JAMRWR000000000 (Lb. brevis), JAMRVU000000000 to JAMRVY000000000 (Lc. lactis subsp. lactis strains DRD-85, DRD-89, DRD-122, DRD-132, DRD-134), JAMRWS000000000 to JAMRWW000000000 (Pd. pentosaceus), JAMRWC000000000 to JAMRWG000000000 (Ln. mesenteroides), JAMRVZ000000000 to JAMRWB000000000 (Lb. curvatus), JAQGDA000000000 to JAQGDD000000000 (Lc. lactis subsp. lactis strains DRD-203, DRD-205, DRD-206, DRD-207) and JAQGDE000000000 to JAQGDF000000000 (Lc. lactis subsp. cremoris strains DRD-208, DRD-210). The version described in this paper is version JAMQJD010000000 to JAMQJS010000000 (Lb. plantarum), JAMRWH010000000 to JAMRWR010000000 (Lb. brevis), JAMRVU010000000 to JAMRVY010000000 (Lc. lactis subsp. lactis; strains DRD-85, DRD-89, DRD-122, DRD-132, DRD-134), JAMRWS010000000 to JAMRWW010000000 (Pd. pentosaceus), JAMRWC010000000 to JAMRWG010000000 (Ln. mesenteroides) and JAMRVZ010000000 to JAMRWB010000000 (Lb. curvatus), JAQGDA010000000 to JAQGDD010000000 (Lc. lactis subsp. lactis strains DRD-203, DRD-205, DRD-206, DRD-207) and JAQGDE010000000 to JAQGDF010000000 (Lc. lactis subsp. cremoris strains DRD-208, DRD-210).

Acknowledgments

This research has been co-financed by the European Union and Greek national funds through the Operational Program Competitiveness, Entrepreneurship, and Innovation, under the call RESEARCH–CREATE–INNOVATE (T1EDK-02087).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cano-Garrido, O.; Seras-Franzoso, J.; Garcia-Fruitós, E. Lactic Acid Bacteria: Reviewing the Potential of a Promising Delivery Live Vector for Biomedical Purposes. Microb. Cell Fact. 2015, 14, 137. [Google Scholar] [CrossRef] [PubMed]
  2. Wang, Y.; Wu, J.; Lv, M.; Shao, Z.; Hungwe, M.; Wang, J.; Bai, X.; Xie, J.; Wang, Y.; Geng, W. Metabolism Characteristics of Lactic Acid Bacteria and the Expanding Applications in Food Industry. Front. Bioeng. Biotechnol. 2021, 9, 612285. [Google Scholar] [CrossRef] [PubMed]
  3. Sun, Z.; Harris, H.M.B.; McCann, A.; Guo, C.; Argimón, S.; Zhang, W.; Yang, X.; Jeffery, I.B.; Cooney, J.C.; Kagawa, T.F.; et al. Expanding the Biotechnology Potential of Lactobacilli through Comparative Genomics of 213 Strains and Associated Genera. Nat. Commun. 2015, 6, 8322. [Google Scholar] [CrossRef] [PubMed]
  4. Mayo, B.; Rodríguez, J.; Vázquez, L.; Flórez, A.B. Microbial Interactions within the Cheese Ecosystem and Their Application to Improve Quality and Safety. Foods 2021, 10, 602. [Google Scholar] [CrossRef] [PubMed]
  5. Kelleher, P.; Murphy, J.; Mahony, J.; van Sinderen, D. Next-Generation Sequencing as an Approach to Dairy Starter Selection. Dairy Sci. Technol. 2015, 95, 545–568. [Google Scholar] [CrossRef]
  6. Tofalo, R.; Perpetuini, G.; Battistelli, N.; Pepe, A.; Ianni, A.; Martino, G.; Suzzi, G. Accumulation γ-Aminobutyric Acid and Biogenic Amines in a Traditional Raw Milk Ewe’s Cheese. Foods 2019, 8, 401. [Google Scholar] [CrossRef]
  7. Tofalo, R.; Schirone, M.; Fasoli, G.; Perpetuini, G.; Patrignani, F.; Manetta, A.C.; Lanciotti, R.; Corsetti, A.; Martino, G.; Suzzi, G. Influence of Pig Rennet on Proteolysis, Organic Acids Content and Microbiota of Pecorino Di Farindola, a Traditional Italian Ewe’s Raw Milk Cheese. Food Chem. 2015, 175, 121–127. [Google Scholar] [CrossRef]
  8. 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. [Google Scholar] [CrossRef]
  9. Sałański, P.; Kowalczyk, M.; Bardowski, J.K.; Szczepankowska, A.K. Health-Promoting Nature of Lactococcus lactis IBB109 and Lactococcus lactis IBB417 Strains Exhibiting Proliferation Inhibition and Stimulation of Interleukin-18 Expression in Colorectal Cancer Cells. Front. Microbiol. 2022, 13, 822912. [Google Scholar] [CrossRef]
  10. Pightling, A.W.; Pettengill, J.B.; Luo, Y.; Baugher, J.D.; Rand, H.; Strain, E. Interpreting Whole-Genome Sequence Analyses of Foodborne Bacteria for Regulatory Applications and Outbreak Investigations. Front. Microbiol. 2018, 9, 01482. [Google Scholar] [CrossRef] [Green Version]
  11. Bull, M.J.; Jolley, K.A.; Bray, J.E.; Aerts, M.; Vandamme, P.; Maiden, M.C.J.; Marchesi, J.R.; Mahenthiralingam, E. The Domestication of the Probiotic Bacterium Lactobacillus acidophilus. Sci. Rep. 2015, 4, 7202. [Google Scholar] [CrossRef]
  12. Yu, J.; Sun, Z.; Liu, W.; Bao, Q.; Zhang, J.; Zhang, H. Phylogenetic Study of Lactobacillus acidophilus Group, L. casei Group and L. plantarum Group Based on Partial hsp60, pheS and tuf Gene Sequences. Eur. Food Res. Technol. 2012, 234, 927–934. [Google Scholar] [CrossRef]
  13. Tsigkrimani, M.; Bakogianni, M.; Paramithiotis, S.; Bosnea, L.; Pappa, E.; Drosinos, E.H.; Skandamis, P.N.; Mataragas, M. Microbial Ecology of Artisanal Feta and Kefalograviera Cheeses, Part I: Bacterial Community and Its Functional Characteristics with Focus on Lactic Acid Bacteria as Determined by Culture-Dependent Methods and Phenotype Microarrays. Microorganisms 2022, 10, 161. [Google Scholar] [CrossRef]
  14. Tsigkrimani, M.; Panagiotarea, K.; Paramithiotis, S.; Bosnea, L.; Pappa, E.; Drosinos, E.H.; Skandamis, P.N.; Mataragas, M. Microbial Ecology of Sheep Milk, Artisanal Feta, and Kefalograviera Cheeses. Part II: Technological, Safety, and Probiotic Attributes of Lactic Acid Bacteria Isolates. Foods 2022, 11, 459. [Google Scholar] [CrossRef]
  15. Syrokou, M.K.; Themeli, C.; Paramithiotis, S.; Mataragas, M.; Bosnea, L.; Argyri, A.A.; Chorianopoulos, N.G.; Skandamis, P.N.; Drosinos, E.H. Microbial Ecology of Greek Wheat Sourdoughs, Identified by a Culture-Dependent and a Culture-Independent Approach. Foods 2020, 9, 1603. [Google Scholar] [CrossRef]
  16. Andrews, S. FastQC: A Quality Control Tool for High Throughput Sequence Data; Babraham Institute: Cambridge, UK, 2019. [Google Scholar]
  17. Arkin, A.P.; Cottingham, R.W.; Henry, C.S.; Harris, N.L.; Stevens, R.L.; Maslov, S.; Dehal, P.; Ware, D.; Perez, F.; Canon, S.; et al. KBase: The United States Department of Energy Systems Biology Knowledgebase. Nat. Biotechnol. 2018, 36, 566–569. [Google Scholar] [CrossRef]
  18. Wick, R.R.; Judd, L.M.; Gorrie, C.L.; Holt, K.E. Completing Bacterial Genome Assemblies with Multiplex MinION Sequencing. Microb. Genomics 2017, 3. [Google Scholar] [CrossRef]
  19. Walker, B.J.; Abeel, T.; Shea, T.; Priest, M.; Abouelliel, A.; Sakthikumar, S.; Cuomo, C.A.; Zeng, Q.; Wortman, J.; Young, S.K.; et al. Pilon: An Integrated Tool for Comprehensive Microbial Variant Detection and Genome Assembly Improvement. PLoS ONE 2014, 9, e112963. [Google Scholar] [CrossRef]
  20. Davis, J.J.; Wattam, A.R.; Aziz, R.K.; Brettin, T.; Butler, R.; Butler, R.M.; Chlenski, P.; Conrad, N.; Dickerman, A.; Dietrich, E.M.; et al. The PATRIC Bioinformatics Resource Center: Expanding Data and Analysis Capabilities. Nucleic Acids Res. 2020, 48, D606–D612. [Google Scholar] [CrossRef]
  21. Bosi, E.; Donati, B.; Galardini, M.; Brunetti, S.; Sagot, M.F.; Lió, P.; Crescenzi, P.; Fani, R.; Fondi, M. MeDuSa: A Multi-Draft Based Scaffolder. Bioinformatics 2015, 31, 2443–2451. [Google Scholar] [CrossRef] [Green Version]
  22. Parks, D.H.; Imelfort, M.; Skennerton, C.T.; Hugenholtz, P.; Tyson, G.W. CheckM: Assessing the Quality of Microbial Genomes Recovered from Isolates, Single Cells, and Metagenomes. Genome Res. 2015, 25, 1043. [Google Scholar] [CrossRef] [PubMed]
  23. Lu, J.; Salzberg, S.L. SkewIT: The Skew Index Test for Large-Scale GC Skew Analysis of Bacterial Genomes. PLoS Comput. Biol. 2020, 16, e1008439. [Google Scholar] [CrossRef] [PubMed]
  24. Gurevich, A.; Saveliev, V.; Vyahhi, N.; Tesler, G. QUAST: Quality Assessment Tool for Genome Assemblies. Bioinformatics 2013, 29, 1072–1075. [Google Scholar] [CrossRef] [PubMed]
  25. Wood, D.E.; Lu, J.; Langmead, B. Improved Metagenomic Analysis with Kraken 2. Genome Biol. 2019, 20, 257. [Google Scholar] [CrossRef] [PubMed]
  26. Meier-Kolthoff, J.P.; Göker, M. TYGS Is an Automated High-Throughput Platform for State-of-the-Art Genome-Based Taxonomy. Nat. Commun. 2019, 10, 2182. [Google Scholar] [CrossRef]
  27. Seemann, T. Prokka: Rapid Prokaryotic Genome Annotation. Bioinformatics 2014, 30, 2068–2069. [Google Scholar] [CrossRef]
  28. Tatusov, R.L.; Fedorova, N.D.; Jackson, J.D.; Jacobs, A.R.; Kiryutin, B.; Koonin, E.V.; Krylov, D.M.; Mazumder, R.; Smirnov, S.; Nikolskaya, A.N.; et al. The COG Database: An Updated Vesion Includes Eukaryotes. BMC Bioinform. 2003, 4, 1–14. [Google Scholar] [CrossRef]
  29. Cantalapiedra, C.P.; Hernández-Plaza, A.; Letunic, I.; Bork, P.; Huerta-Cepas, J. EggNOG-Mapper v2: Functional Annotation, Orthology Assignments, and Domain Prediction at the Metagenomic Scale. Mol. Biol. Evol. 2021, 38, 5825–5829. [Google Scholar] [CrossRef]
  30. Couvin, D.; Bernheim, A.; Toffano-Nioche, C.; Touchon, M.; Michalik, J.; Néron, B.; Rocha, E.P.C.; Vergnaud, G.; Gautheret, D.; Pourcel, C. CRISPRCasFinder, an Update of CRISRFinder, Includes a Portable Version, Enhanced Performance and Integrates Search for Cas Proteins. Nucleic Acids Res. 2018, 46, W246–W251. [Google Scholar] [CrossRef]
  31. Arndt, D.; Marcu, A.; Liang, Y.; Wishart, D.S. PHAST, PHASTER and PHASTEST: Tools for Finding Prophage in Bacterial Genomes. Brief. Bioinform. 2019, 20, 1560–1567. [Google Scholar] [CrossRef]
  32. Seemann, T. Abricate; Github: San Francisco, CA, USA, 2020. [Google Scholar]
  33. Zankari, E.; Hasman, H.; Cosentino, S.; Vestergaard, M.; Rasmussen, S.; Lund, O.; Aarestrup, F.M.; Larsen, M.V. Identification of Acquired Antimicrobial Resistance Genes. J. Antimicrob. Chemother. 2012, 67, 2640–2644. [Google Scholar] [CrossRef]
  34. Chen, L.; Zheng, D.; Liu, B.; Yang, J.; Jin, Q. VFDB 2016: Hierarchical and Refined Dataset for Big Data Analysis—10 Years On. Nucleic Acids Res. 2016, 44, D694–D697. [Google Scholar] [CrossRef]
  35. Johansson, M.H.K.; Bortolaia, V.; Tansirichaiya, S.; Aarestrup, F.M.; Roberts, A.P.; Petersen, T.N. Detection of Mobile Genetic Elements Associated with Antibiotic Resistance in Salmonella enterica Using a Newly Developed Web Tool: MobileElementFinder. J. Antimicrob. Chemother. 2021, 76, 101–109. [Google Scholar] [CrossRef]
  36. Carattoli, A.; Zankari, E.; Garciá-Fernández, A.; Larsen, M.V.; Lund, O.; Villa, L.; Aarestrup, F.M.; Hasman, H. In Silico Detection and Typing of Plasmids Using Plasmidfinder and Plasmid Multilocus Sequence Typing. Antimicrob. Agents Chemother. 2014, 58, 3895–3903. [Google Scholar] [CrossRef]
  37. Zhang, H.; Yohe, T.; Huang, L.; Entwistle, S.; Wu, P.; Yang, Z.; Busk, P.K.; Xu, Y.; Yin, Y. DbCAN2: A Meta Server for Automated Carbohydrate-Active Enzyme Annotation. Nucleic Acids Res. 2018, 46, W95–W101. [Google Scholar] [CrossRef]
  38. Weimann, A.; Mooren, K.; Frank, J.; Pope, P.B.; Bremges, A.; McHardy, A.C. From Genomes to Phenotypes: Traitar, the Microbial Trait Analyzer. mSystems 2016, 1. [Google Scholar] [CrossRef]
  39. Brynildsrud, O.; Bohlin, J.; Scheffer, L.; Eldholm, V. Rapid Scoring of Genes in Microbial Pan-Genome-Wide Association Studies with Scoary. Genome Biol. 2016, 17, 238. [Google Scholar] [CrossRef]
  40. Reijnders, M.J.M.F.; Waterhouse, R.M. Summary Visualizations of Gene Ontology Terms with GO-Figure! Front. Bioinform. 2021, 1, 638255. [Google Scholar] [CrossRef]
  41. Page, A.J.; Cummins, C.A.; Hunt, M.; Wong, V.K.; Reuter, S.; Holden, M.T.G.; Fookes, M.; Falush, D.; Keane, J.A.; Parkhill, J. Roary: Rapid Large-Scale Prokaryote Pan Genome Analysis. Bioinformatics 2015, 31, 3691–3693. [Google Scholar] [CrossRef]
  42. Price, M.N.; Dehal, P.S.; Arkin, A.P. FastTree: Computing Large Minimum Evolution Trees with Profiles Instead of a Distance Matrix. Mol. Biol. Evol. 2009, 26, 1641–1650. [Google Scholar] [CrossRef]
  43. Letunic, I.; Bork, P. Interactive Tree Of Life (ITOL) v5: An Online Tool for Phylogenetic Tree Display and Annotation. Nucleic Acids Res. 2021, 49, W293–W296. [Google Scholar] [CrossRef] [PubMed]
  44. Afgan, E.; Baker, D.; Batut, B.; van den Beek, M.; Bouvier, D.; Cech, M.; Chilton, J.; Clements, D.; Coraor, N.; Grüning, B.A.; et al. The Galaxy Platform for Accessible, Reproducible and Collaborative Biomedical Analyses: 2018 Update. Nucleic Acids Res. 2018, 46, W537–W544. [Google Scholar] [CrossRef] [PubMed]
  45. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2020. [Google Scholar]
  46. Overbeek, R.; Begley, T.; Butler, R.M.; Choudhuri, J.V.; Chuang, H.Y.; Cohoon, M.; de Crécy-Lagard, V.; Diaz, N.; Disz, T.; Edwards, R.; et al. The Subsystems Approach to Genome Annotation and Its Use in the Project to Annotate 1000 Genomes. Nucleic Acids Res. 2005, 33, 5691–5702. [Google Scholar] [CrossRef] [PubMed]
  47. Tenea, G.N.; Ortega, C. Genome Characterization of Lactiplantibacillus Plantarum Strain UTNGt2 Originated from Theobroma grandiflorum (White Cacao) of Ecuadorian Amazon: Antimicrobial Peptides from Safety to Potential Applications. Antibiotics 2021, 10, 383. [Google Scholar] [CrossRef] [PubMed]
  48. Flórez, A.B.; de Los Reyes-Gavilán, C.G.; Wind, A.; Mayo, B.; Margolles, A. Ubiquity and Diversity of Multidrug Resistance Genes in Lactococcus lactis Strains Isolated between 1936 and 1995. FEMS Microbiol. Lett. 2006, 263, 21–25. [Google Scholar] [CrossRef]
  49. Almeida, M.; Hébert, A.; Abraham, A.L.; Rasmussen, S.; Monnet, C.; Pons, N.; Delbès, C.; Loux, V.; Batto, J.M.; Leonard, P.; et al. Construction of a Dairy Microbial Genome Catalog Opens New Perspectives for the Metagenomic Analysis of Dairy Fermented Products. BMC Genomics 2014, 15. [Google Scholar] [CrossRef]
  50. Silva, C.C.G.; Silva, S.P.M.; Ribeiro, S.C. Application of Bacteriocins and Protective Cultures in Dairy Food Preservation. Front. Microbiol. 2018, 9, 594. [Google Scholar] [CrossRef]
  51. Henning, C.; Gautam, D.; Muriana, P. Identification of Multiple Bacteriocins in Enterococcus Spp. Using an Enterococcus-Specific Bacteriocin PCR Array. Microorganisms 2015, 3, 1–16. [Google Scholar] [CrossRef]
  52. Franz, C.M.A.P.; Van Belkum, M.J.; Holzapfel, W.H.; Abriouel, H.; Gálvez, A. Diversity of Enterococcal Bacteriocins and Their Grouping in a New Classification Scheme. FEMS Microbiol. Rev. 2007, 31, 293–310. [Google Scholar] [CrossRef]
  53. Todokoro, D.; Tomita, H.; Inoue, T.; Ike, Y. Genetic Analysis of Bacteriocin 43 of Vancomycin-Resistant Enterococcus faecium. Appl. Environ. Microbiol. 2006, 72, 6955–6964. [Google Scholar] [CrossRef] [Green Version]
  54. Abdulhussain Kareem, R.; Razavi, S.H. Plantaricin Bacteriocins: As Safe Alternative Antimicrobial Peptides in Food Preservation—A Review. J. Food Saf. 2020, 40. [Google Scholar] [CrossRef]
  55. Leer, R.J.; van der Vossen, J.M.B.M.; van Giezen, M.; van Noort Johannes, M.; Pouwels, P.H. Genetic Analysis of Acidocin B, a Novel Bacteriocin Produced by Lactobacillus acidophilus. Microbiology 1995, 141, 1629–1635. [Google Scholar] [CrossRef]
  56. Hammami, R.; Zouhir, A.; Le Lay, C.; Ben Hamida, J.; Fliss, I. BACTIBASE Second Release: A Database and Tool Platform for Bacteriocin Characterization. BMC Microbiol. 2010, 10, 22. [Google Scholar] [CrossRef]
  57. Markusková, B.; Lichvariková, A.; Szemes, T.; Koreňová, J.; Kuchta, T.; Drahovská, H. Genome Analysis of Lactic Acid Bacterial Strains Selected as Potential Starters for Traditional Slovakian Bryndza Cheese. FEMS Microbiol. Lett. 2018. [Google Scholar] [CrossRef]
  58. Tarrah, A.; Pakroo, S.; Lemos Junior, W.J.F.; Guerra, A.F.; Corich, V.; Giacomini, A. Complete Genome Sequence and Carbohydrates-Active EnZymes (CAZymes) Analysis of Lactobacillus paracasei DTA72, a Potential Probiotic Strain with Strong Capability to Use Inulin. Curr. Microbiol. 2020, 77, 2867–2875. [Google Scholar] [CrossRef]
  59. Shen, J.; Zheng, L.; Chen, X.; Han, X.; Cao, Y.; Yao, J. Metagenomic Analyses of Microbial and Carbohydrate-Active Enzymes in the Rumen of Dairy Goats Fed Different Rumen Degradable Starch. Front. Microbiol. 2020, 11, 1003. [Google Scholar] [CrossRef]
  60. He, Q.; Hou, Q.; Wang, Y.; Li, J.; Li, W.; Kwok, L.-Y.; Sun, Z.; Zhang, H.; Zhong, Z. Comparative Genomic Analysis of Enterococcus faecalis: Insights into Their Environmental Adaptations. BMC Genomics 2018, 19, 527. [Google Scholar] [CrossRef]
  61. Ghattargi, V.C.; Gaikwad, M.A.; Meti, B.S.; Nimonkar, Y.S.; Dixit, K.; Prakash, O.; Shouche, Y.S.; Pawar, S.P.; Dhotre, D.P. Comparative Genome Analysis Reveals Key Genetic Factors Associated with Probiotic Property in Enterococcus faecium Strains. BMC Genomics 2018, 19, 652. [Google Scholar] [CrossRef]
  62. Zhang, W.; Sun, Z.; Menghe, B.; Zhang, H. Short Communication: Single Molecule, Real-Time Sequencing Technology Revealed Species- and Strain-Specific Methylation Patterns of 2 Lactobacillus Strains. J. Dairy Sci. 2015, 98, 3020–3024. [Google Scholar] [CrossRef]
  63. Chaillou, S.; Champomier-Vergès, M.-C.; Cornet, M.; Crutz-Le Coq, A.-M.; Dudez, A.-M.; Martin, V.; Beaufils, S.; Darbon-Rongère, E.; Bossy, R.; Loux, V.; et al. The Complete Genome Sequence of the Meat-Borne Lactic Acid Bacterium Lactobacillus sakei 23K. Nat. Biotechnol. 2005, 23, 1527–1533. [Google Scholar] [CrossRef] [Green Version]
  64. Ainsworth, S.; Stockdale, S.; Bottacini, F.; Mahony, J.; van Sinderen, D. The Lactococcus lactis Plasmidome: Much Learnt, yet Still Lots to Discover. FEMS Microbiol. Rev. 2014, 38, 1066–1088. [Google Scholar] [CrossRef] [PubMed]
  65. Apostolakos, I.; Paramithiotis, S.; Mataragas, M. Functional and Safety Characterization of Weissella paramesenteroides Strains Isolated from Dairy Products through Whole-Genome Sequencing and Comparative Genomics. Dairy 2022, 3, 799–813. [Google Scholar] [CrossRef]
  66. Durrant, M.G.; Li, M.M.; Siranosian, B.A.; Montgomery, S.B.; Bhatt, A.S. A Bioinformatic Analysis of Integrative Mobile Genetic Elements Highlights Their Role in Bacterial Adaptation. Cell Host Microbe 2020, 27, 140–153.e9. [Google Scholar] [CrossRef] [PubMed]
  67. Kiousi, D.E.; Efstathiou, C.; Tegopoulos, K.; Mantzourani, I.; Alexopoulos, A.; Plessas, S.; Kolovos, P.; Koffa, M.; Galanis, A. Genomic Insight Into Lacticaseibacillus paracasei SP5, Reveals Genes and Gene Clusters of Probiotic Interest and Biotechnological Potential. Front. Microbiol. 2022, 13, 2038. [Google Scholar] [CrossRef]
  68. Liu, D.; Liu, Y.; Zhang, D.; Chen, X.; Liu, Q.; Xiong, B.; Zhang, L.; Wei, L.; Wang, Y.; Fang, H.; et al. Quantitative Proteome Profiling Reveals Cellobiose-Dependent Protein Processing and Export Pathways for the Lignocellulolytic Response in Neurospora crassa. Appl. Environ. Microbiol. 2020, 86. [Google Scholar] [CrossRef]
  69. Allen, J.R.; Torres-Acosta, M.A.; Mohan, N.; Lye, G.J.; Ward, J.M. Segregationally Stabilised Plasmids Improve Production of Commodity Chemicals in Glucose-Limited Continuous Fermentation. Microb. Cell Fact. 2022, 21, 229. [Google Scholar] [CrossRef]
  70. Krewinkel, M.; Gosch, M.; Rentschler, E.; Fischer, L. Epilactose Production by 2 Cellobiose 2-Epimerases in Natural Milk. J. Dairy Sci. 2014, 97, 155–161. [Google Scholar] [CrossRef]
  71. Van Trijp, M.P.H.; Schutte, S.; Esser, D.; Wopereis, S.; Hoevenaars, F.P.M.; Hooiveld, G.J.E.J.; Afman, L.A. Minor Changes in the Composition and Function of the Gut Microbiota During a 12-Week Whole Grain Wheat or Refined Wheat Intervention Correlate with Liver Fat in Overweight and Obese Adults. J. Nutr. 2021, 151, 491–502. [Google Scholar] [CrossRef]
  72. Lai, C.-Y.; Cronan, J.E. β-Ketoacyl-Acyl Carrier Protein Synthase III (FabH) Is Essential for Bacterial Fatty Acid Synthesis. J. Biol. Chem. 2003, 278, 51494–51503. [Google Scholar] [CrossRef]
  73. McLeod, A.; Snipen, L.; Naterstad, K.; Axelsson, L. Global Transcriptome Response in Lactobacillus sakei during Growth on Ribose. BMC Microbiol. 2011, 11, 145. [Google Scholar] [CrossRef] [Green Version]
  74. Watanabe, S.; Fukumori, F.; Nishiwaki, H.; Sakurai, Y.; Tajima, K.; Watanabe, Y. Novel Non-Phosphorylative Pathway of Pentose Metabolism from Bacteria. Sci. Rep. 2019, 9, 155. [Google Scholar] [CrossRef]
  75. Okano, K.; Yoshida, S.; Tanaka, T.; Ogino, C.; Fukuda, H.; Kondo, A. Homo-D-Lactic Acid Fermentation from Arabinose by Redirection of the Phosphoketolase Pathway to the Pentose Phosphate Pathway in L-Lactate Dehydrogenase Gene-Deficient Lactobacillus plantarum. Appl. Environ. Microbiol. 2009, 75, 5175–5178. [Google Scholar] [CrossRef]
  76. Welman, A.D.; Maddox, I.S. Exopolysaccharides from Lactic Acid Bacteria: Perspectives and Challenges. Trends Biotechnol. 2003, 21, 269–274. [Google Scholar] [CrossRef]
  77. Angelin, J.; Kavitha, M. Exopolysaccharides from Probiotic Bacteria and Their Health Potential. Int. J. Biol. Macromol. 2020, 162, 853–865. [Google Scholar] [CrossRef]
  78. Yilmaz, B.; Bangar, S.P.; Echegaray, N.; Suri, S.; Tomasevic, I.; Manuel Lorenzo, J.; Melekoglu, E.; Rocha, J.M.; Ozogul, F. The Impacts of Lactiplantibacillus plantarum on the Functional Properties of Fermented Foods: A Review of Current Knowledge. Microorganisms 2022, 10, 826. [Google Scholar] [CrossRef]
  79. Koutsoumanis, K.; Allende, A.; Alvarez-Ordóñez, A.; Bolton, D.; Bover-Cid, S.; Chemaly, M.; Davies, R.; De Cesare, A.; Hilbert, F.; Lindqvist, R.; et al. Scientific Opinion on the Update of the List of QPS-recommended Biological Agents Intentionally Added to Food or Feed as Notified to EFSA (2017–2019). EFSA J. 2020, 18. [Google Scholar] [CrossRef]
  80. Teixeira, C.G.; da Silva, R.R.; Fusieger, A.; Martins, E.; de Freitas, R.; de Carvalho, A.F. Weissella Na Indústria de Alimentos: Uma Revisão. Res. Soc. Dev. 2021, 10, e8310514557. [Google Scholar] [CrossRef]
Foods 12 00599 g001 550
Figure 1. Bar chart plot showing the average number of genes per COG category for each LAB species. Abbreviation symbols for COG categories are according to the legend on the right. Species’ bars annotated with an asterisk indicate that the COG category was significantly more enriched for that species compared to all other LAB species. No genes belonging to the COG category “Mobilome: Prophages & Transposons” were identified in the isolates of our collection.
Figure 1. Bar chart plot showing the average number of genes per COG category for each LAB species. Abbreviation symbols for COG categories are according to the legend on the right. Species’ bars annotated with an asterisk indicate that the COG category was significantly more enriched for that species compared to all other LAB species. No genes belonging to the COG category “Mobilome: Prophages & Transposons” were identified in the isolates of our collection.
Foods 12 00599 g001
Foods 12 00599 g002 550
Figure 2. Phylogenetic tree of the 51 LAB isolates. The colored ring indicates the isolates’ species according to the legend. Symbols at the tree’s leaves denote the isolation source of each isolate according to the legend.
Figure 2. Phylogenetic tree of the 51 LAB isolates. The colored ring indicates the isolates’ species according to the legend. Symbols at the tree’s leaves denote the isolation source of each isolate according to the legend.
Foods 12 00599 g002
Foods 12 00599 g003 550
Figure 3. Clustered heatmap generated using a presence-absence data matrix of MGEs, CAZymes, and bacteriocins of all LAB isolates (n = 51). Each column represents one isolate and each row one genomic element. Isolates and genomic elements are color-coded according to the legend. Presence (dark orange) or absence (light orange) of the element in the respective isolates is shown.
Figure 3. Clustered heatmap generated using a presence-absence data matrix of MGEs, CAZymes, and bacteriocins of all LAB isolates (n = 51). Each column represents one isolate and each row one genomic element. Isolates and genomic elements are color-coded according to the legend. Presence (dark orange) or absence (light orange) of the element in the respective isolates is shown.
Foods 12 00599 g003
Foods 12 00599 g004 550
Figure 4. Clustered heatmap generated using a presence-absence data matrix of LAB phenotypic traits. Traits with a double asterisk (**) have been experimentally confirmed in the work of Tsigkrimani et al. (2022) [14]. Each column represents one isolate and each row one phenotypic trait. Isolates and traits are color-coded according to the legend. Presence (dark orange) or absence (light orange) of the trait in the respective isolates is shown.
Figure 4. Clustered heatmap generated using a presence-absence data matrix of LAB phenotypic traits. Traits with a double asterisk (**) have been experimentally confirmed in the work of Tsigkrimani et al. (2022) [14]. Each column represents one isolate and each row one phenotypic trait. Isolates and traits are color-coded according to the legend. Presence (dark orange) or absence (light orange) of the trait in the respective isolates is shown.
Foods 12 00599 g004
Table
Table 1. Species identification and assembly statistics for all LAB genomes.
Table 1. Species identification and assembly statistics for all LAB genomes.
Strain ID Microorganism SourceGenome Size (Mb) GC Content (%) No. of Scaffolds N50 (Mb) No. of CDSs
DRD-10Lb. plantarumKefalograviera1.738.04.01.72445
DRD-15Lb. plantarumFeta3.044.922.01.32435
DRD-16Lb. plantarumFeta2.545.823.02.32449
DRD-31Lb. plantarumFeta3.544.2126.03.22170
DRD-32Lb. plantarumFeta3.544.2125.03.12450
DRD-34Lb. plantarumFeta2.037.624.01.02449
DRD-36Lb. plantarumFeta3.044.920.02.32453
DRD-38Lb. plantarumFeta3.444.295.02.92449
DRD-41Lb. plantarumFeta3.544.2105.02.12450
DRD-44Lb. plantarumFeta2.545.819.02.32454
DRD-46Lb. plantarumKefalograviera3.544.2103.02.92450
DRD-63Lb. plantarumFeta3.544.2106.03.22143
DRD-65Lb. plantarumFeta3.544.2109.03.22147
DRD-67Lb. plantarumFeta1.737.610.01.32147
DRD-76Lb. plantarumSheep milk3.544.2115.03.22741
DRD-124Lb. plantarumSheep milk3.444.287.03.23213
DRD-208Lc. lactis cremorisKefalograviera1.937.114.01.43253
DRD-210Lc. lactis cremorisKefalograviera2.545.821.02.33254
DRD-85Lc. lactis lactisSheep milk2.545.822.02.32753
DRD-89Lc. lactis lactisSheep milk2.545.821.02.33232
DRD-122Lc. lactis lactisSheep milk2.545.822.02.33253
DRD-132Lc. lactis lactisSheep milk1.937.19.01.93264
DRD-134Lc. lactis lactisSheep milk3.644.1227.03.13260
DRD-203Lc. lactis lactisKefalograviera3.044.918.01.03262
DRD-205Lc. lactis lactisKefalograviera3.544.2115.03.03251
DRD-206Lc. lactis lactisKefalograviera3.544.296.03.23243
DRD-207Lc. lactis lactisKefalograviera2.434.919.01.73364
DRD-164Lb. curvatusFeta2.534.926.02.32754
DRD-170Lb. curvatusFeta2.534.930.00.73255
DRD-171Lb. curvatusFeta3.444.292.03.23254
DRD-2Ln. mesenteroidesKefalograviera2.635.027.02.22562
DRD-30Ln. mesenteroidesFeta2.534.931.01.92560
DRD-138Ln. mesenteroidesSheep milk2.545.820.02.32507
DRD-140Ln. mesenteroidesSheep milk2.237.515.02.12630
DRD-141Ln. mesenteroidesSheep milk2.545.821.02.32458
DRD-12Lb. brevisFeta2.037.714.01.12486
DRD-35Lb. brevisFeta2.137.442.01.22555
DRD-51Lb. brevisKefalograviera2.037.217.01.82551
DRD-52Lb. brevisKefalograviera2.141.741.01.92506
DRD-59Lb. brevisKefalograviera2.141.741.01.92398
DRD-60Lb. brevisKefalograviera2.141.741.01.92504
DRD-136Lb. brevisSheep milk2.037.27.01.82253
DRD-139Lb. brevisSheep milk2.346.412.02.22020
DRD-195Lb. brevisKefalograviera2.545.820.02.32113
DRD-198Lb. brevisKefalograviera2.545.821.02.31710
DRD-201Lb. brevisKefalograviera2.634.926.00.52016
DRD-42Pd. pentosaceusFeta2.534.925.01.42008
DRD-48Pd. pentosaceusKefalograviera2.534.926.01.31998
DRD-61Pd. pentosaceusKefalograviera2.534.915.02.31662
DRD-144Pd. pentosaceusFeta2.535.441.01.82019
DRD-185Pd. pentosaceusFeta2.535.443.01.81986
Table
Table 2. Plasmids identified in LAB strains by the plasmidFinder tool. The number of isolates per species having a plasmid (based on the NCBI Accession column) is shown on the first column.
Table 2. Plasmids identified in LAB strains by the plasmidFinder tool. The number of isolates per species having a plasmid (based on the NCBI Accession column) is shown on the first column.
Species
(Number of Isolates with Plasmid/Total Number of Isolates)		PlasmidsIdentity (%)Length (bp)NoteNCBI
Accession
Lc. lactis (5/11)repUS4901108repA(pCI2000)AF178424
Lc. lactis (8/11)
Ln. mesenteroides (1/5)		rep32971151pli0023(pLI100)AL592102
Lc. lactis (7/11)repUS331001352repA(pGdh442)AY849557
Ln. mesenteroides (5/5)rep31871132LKI10596(LkipL4719)CP001755
Lb. plantarum (11/16)repUS73941100rep(pLBUC02)CP002654
Lb. plantarum (12/16)rep38811031rep(pLBUC03)CP002655
Lb. brevis (1/11)
Pd. pentosaceus (1/5)		rep2892932repA(pCIS4)CP003162
Ln. mesenteroides (1/5)repUS72981036C27008541(pKLC4)CP003855
Lb. brevis (9/11)
Lb. curvatus (3/3)
Lb. plantarum (12/16)		rep3898885repA(LBPp1)CP005943
Lb. brevis (10/11)
Lb. curvatus (3/3)
Lb. plantarum (13/16)		rep2899915LBPp6g007(LBPp6)CP005948
Ln. mesenteroides (5/5)rep31881133LCKp400005(pLCK4)DQ489739
Lb. curvatus (3/3)repUS5191662rep(pCPS49)FN806792
Lc. lactis (4/11)rep33831131rep(pSMA198)HE613570
Lc. lactis (1/11)repUS421001157repB(pVF18)JN172910
Lc. lactis (1/11)rep32811168repB(pVF22)JN172912
Lb. plantarum (1/16)
Pd. pentosaceus (1/5)		repUS6493956repA(pR18)JN601038
Lc. lactis (1/11)rep33861144rep(pK214)X92946
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Apostolakos, I.; Paramithiotis, S.; Mataragas, M. Comparative Genomic Analysis Reveals the Functional Traits and Safety Status of Lactic Acid Bacteria Retrieved from Artisanal Cheeses and Raw Sheep Milk. Foods 2023, 12, 599. https://doi.org/10.3390/foods12030599

AMA Style

Apostolakos I, Paramithiotis S, Mataragas M. Comparative Genomic Analysis Reveals the Functional Traits and Safety Status of Lactic Acid Bacteria Retrieved from Artisanal Cheeses and Raw Sheep Milk. Foods. 2023; 12(3):599. https://doi.org/10.3390/foods12030599

Chicago/Turabian Style

Apostolakos, Ilias, Spiros Paramithiotis, and Marios Mataragas. 2023. "Comparative Genomic Analysis Reveals the Functional Traits and Safety Status of Lactic Acid Bacteria Retrieved from Artisanal Cheeses and Raw Sheep Milk" Foods 12, no. 3: 599. https://doi.org/10.3390/foods12030599

APA Style

Apostolakos, I., Paramithiotis, S., & Mataragas, M. (2023). Comparative Genomic Analysis Reveals the Functional Traits and Safety Status of Lactic Acid Bacteria Retrieved from Artisanal Cheeses and Raw Sheep Milk. Foods, 12(3), 599. https://doi.org/10.3390/foods12030599

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

Article Metrics

Back to TopTop