Preprint
Article

Antibiotic Resistance/Susceptibility Profiles of Staphylococcus equorum Strains from Cheese, and Genome Analysis for Antibiotic Resistance Genes

Altmetrics

Downloads

193

Views

85

Comments

0

A peer-reviewed article of this preprint also exists.

Submitted:

06 June 2023

Posted:

06 June 2023

You are already at the latest version

Alerts
Abstract
In this work, the resistance/susceptibility (R/S) profile of Staphylococcus equorum strains (n=30) from cheese to 16 antibiotics was determined by broth microdilution. The minimum inhibitory concentration (MIC) for all antibiotics was low in most strains, although higher MICs compatible with acquired genes were also noted. Genome analysis of 13 strains showed the S. equorum resistome to be composed of intrinsic mechanisms, acquired mutations, and acquired genes. However, the genetic data did not always correlate with the phenotype. As such, a cat gene providing resistance to chloramphenicol was found on a plasmid in one strain; this was able to provide resistance to Staphylococcus aureus after electroporation. An msr(A) polymorphic gene was identified in five strains. The Mrs(A) variants were associated with variable resistance to erythromycin. All strains harboured a polymorphic fosB/fosD gene, although only one acquired copy was associated with strong resistance to fosfomycin. Similarly, a plasmid-associated blaR1-blaZI operon encoding a penicillinase was identified in five ampicillin- and penicillin G-susceptible strains. Identified genes not associated with resistance further included mph(C) in two strains and norA in all strains. The antibiotic R/S status and gene content of S. equorum strains intended to be employed in food systems should be carefully determined.
Keywords: 
Subject: Biology and Life Sciences  -   Food Science and Technology

1. Introduction

The discovery of antibiotics is among the 20th centuryʼs most important achievements in human and veterinary medicine [1]. Their promise has, however, been tarnished by the appearance of resistance and the subsequent transfer of the associated genetic determinants to pathogenic bacteria, hindering the treatment of infections [2]. Although the presence of antibiotic resistance in beneficial and commensal bacteria poses no direct risk to human or animal health, populations of such microorganisms are potential reservoirs of resistance genes that pathogens may acquire [3,4]. The food chain is pivotal in the transmission of such genes, a consequence of the high cell densities and stress conditions the associated microbiota may experience [5,6]. Transfer can occur either during food processing or intestinal transit [7,8].
Milk and fermented dairy products are key players in the transmission of antibiotic resistance in food [9]. Large gene loads for antimicrobial resistance have been repeatedly reported in milk and dairy products [10,11,12,13]. Interest in antibiotic resistance in the dairy setting has largely focused on lactic acid bacteria (LAB) [3,14,15], but this clashes with more recent postulates of the One Health concept, which understands all environments -and thus their microbiota- to be connected [16]. Given the huge microbial diversity seen in some cheeses [17], plus the fact that the population sizes of several taxa may exceed those of LAB [18,19,20], focusing antibiotic resistance only on the latter type of microorganisms is insufficient. However, distinguishing intrinsic from acquired resistance in food-associated bacteria can be challenging, as it is distinguishing acquired resistance due to mutations from resistance due to acquired genes [14]; as regards transference, only the latter resistance is of great concern [21]. Indeed, resistance/susceptibility (R/S) cut-offs for many species have yet to be established [22].
Staphylococcus equorum is a member of the coagulase-negative staphylococci (CNS), a bacterial group whose members are frequently detected in food-processing environments and fermented foods [23,24,25,26,27]; certainly, S. equorum is a majority commensal population in some cheeses [28,29,30]. CNS species help maintain the colour of fermented meats [31], and contribute to the flavour of all fermented foods via the production of low molecular weight aromatic compounds [30]. They also inhibit the growth of undesirable microorganisms [32]. Therefore, although they do not enjoy Generally Regarded as Safe (GRAS) status, and carry no European Qualified Presumption of Safety (QPS) label, selected strains of S. equorum have been proposed as starters or adjunct cultures for some fermented foods [33,34,35]. No incidents of food poisoning or infections with S. equorum (or any other CNS species) have ever been reported [36]. However, in order of not to spread antibiotic resistance genes through the food chain, starters, and adjunct cultures should contain no transferable ARGs. Work on the antibiotic R/S profiles of S. equorum strains has been particularly scant [24,37,38,39]. Further, these works have mostly involved disk diffusion assays, which do not allow MICs to be established, which are needed for the establishment of reliable R/S cut-offs.
The present work reports a broth microdilution survey of the MICs for 16 antibiotics in S. equorum strains (n=30) isolated from traditional blue-veined cheeses made from raw milk. To examine the links between strain phenotype and genetics, 13 strains were then subjected to genome sequencing. This prompted further phenotypic testing with respect to additional antibiotics, and allowed acquired resistance caused by mutations to be distinguished from resistance afforded by acquired genes.

2. Results

2.1. Antimicrobial testing

The genotyping results revealed wide genetic diversity among the studied S. equorum strains. Indeed, among the 53 isolates, 30 different strains were identified with <90% similarity (i.e., below the repeatability level of the techniques used) (Figure S1). Representative isolates of the 30 different typing profiles were then tested for antibiotic resistance by broth microdilution.
The antibiotic R/S survey of S. equorum strains returned no intrinsic resistance to the set of antibiotics tested, except for moderate resistance to chloramphenicol (range 4-16 µg mL-1; mode 8 µg mL-1) (Table 1).
The majority of the antibiotics showed very low MICs in most strains, suggesting the latter to be highly susceptible. For gentamicin, kanamycin, and neomycin, the MICs in all strains were equal to or lower than the lowest antibiotic concentrations assayed. Low MIC values just one or two dilutions above the weakest concentrations tested were also obtained for streptomycin, vancomycin, and ciprofloxacin. Since S. equorum-specific cut-offs are yet to be established, the breakpoints proposed by the EUCAST and CLSI for catalase-negative staphylococci and Staphylococcus spp. were used (Table 1). This determined that all the S. equorum strains were susceptible to the above six antibiotics. The MICs for the remaining antibiotics varied more strongly among the strains. Nonetheless, the MICs for tetracycline, linezolid, ampicillin, and trimethoprim were still low and, when established, below the cut-offs; consequently, the strains were again considered susceptible. For clindamycin, quinupristin-dalfopristin, and rifampicin, the MICs were slightly above the EUCAST clinical breakpoints, but still below those of the CLSI for some strains; these strains were also deemed susceptible. The MICs of penicillin G were above the CLSI cut-off in seven strains (ranging from 0.5 to 16 µg mL-1). The same resistance range was obtained by an independent penicillin G microdilution test. However, variability in the strains giving lower or higher MICs was noted. Penicillin testing by an MTS assay proved the seven strains to resist higher antibiotic concentrations than all others, although the MICs were much lower than those obtained with the Sensititre plates (Table 1). In addition, the MTS assay showed the occasional presence of colonies growing on the verge of the inhibition halo, suggesting the presence of low-resistant mutants. Finally, clear-cut MICs compatible with acquired resistance were scored in only four strains: for erythromycin in S. equorum 8A3C, 16A1C, and 50A2C (MICs 64, 32, and 24 μg mL-1, respectively), and for chloramphenicol in S. equorum 35A3C (MIC 64 μg mL-1) (Table S1).
In order to establish S. equorum species-specific R/S cut-offs, the distribution of the MICs obtained in this work, and those already recorded in the literature, were analysed using ECOFFinder and NRI worksheet programs. The distribution curves obtained were almost identical with the two programs for most antibiotics or differed by just one dilution. For the sake of simplicity, only those obtained in the NRI analysis are reported (Figure S2). Narrow-peaked curves were seen for most antibiotics, while the majority of the remainder had normal distribution curves, except for those of ampicillin and penicillin G. Tentative S. equorum-specific R/S cut-offs could therefore be proposed for all tested antibiotics (Table 1).

2.2. Genome sequencing

Thirteen strains phenotypically either susceptible to all the tested antibiotics or resistant to any of them were selected for genome sequencing. On average, the recorded genome sequences contained 2,937,190 bp and had a GC content of 32.86%. Genome sequences were distributed in a number of contigs, ranging from 24 (in S. equorum 23A3C) to 107 (in S. equorum 50A2C). Pairwise digital DNA-DNA hybridization (dDDH) of the genome of these strains with those of selected type strains for Staphylococcus species returned an identity value of >95% to S. equorum subsp. equorum NCTC 12414T for seven strains (T17, 2A3C, 8A3C, 16A1C, 23A3C, 35A3C, and 50A2C). Three other strains (CL10P, 48A3I, and 1BCExtra) showed the strongest identity to S. equorum subsp. linens DSM 15097T. These results confirmed their previous identification based on the partial amplification and analysis of their 16S rRNA genes, the manganese-dependent superoxide dismutase gene (sodA), and the glyceraldehyde-3-phosphate dehydrogenase gene (gap) (Vázquez et al., unpublished). The same assignments were inferred from the orthoANI values (Table S2) and the phylogenomic analysis (Figure S3) results for the isolated strains and related staphylococcal-type strains. Although closely related -and treated here as S. equorum- three related isolates (5A3I, 11A1I, and 30A2I) showed dDDH and orthoANI values of 58.8-66.0% and 94.74-95.96%, respectively, with respect to the two S. equorum subspecies type strains. These values are lower than, or on the verge of, the most accepted thresholds (70 and 95% for dDDH and orthoANI, respectively), suggesting they may indicate a different taxonomic unit. The genetic diversity of the 13 sequenced strains was also examined by comparing their genome sequences with those in public databases. The results largely agreed with those obtained by PCR-based genotyping. The scattered distribution of the database strains over the phylogenomic tree (Figure S4) strongly suggests that most of the sequenced strains are genetically unrelated.

2.3. Genome analysis for ARGs

The genome sequences of the 13 strains were examined for genes involved in antimicrobial resistance by comparison with the sequences in the dedicated antimicrobial resistance databases CARD, NCBI-RGC, and ResFinder. Nucleotide sequences and deduced proteins of interest were further sought in PATRIC annotations and compared with sequences in the NCBI database by BLAST analysis. As a result, genes associated with resistance to six classes of antibiotics -penam (β-lactams), macrolides, lincosamides, phenicols, fluoroquinolones, and phosphonic acids (fosfomycin)- were detected (Table 2). Only in a minority of strains was the presence of well-known acquired genes associated with a concomitant phenotypic resistance.
Despite all strains being susceptible to ampicillin and penicillin G, genes involved in the production of β-lactamase enzymes were possessed by all strains (Table 2). The blaR1-blaZI operon described in Staphylococcus aureus [42], which encodes a class A betalactamase (BlaZ, a penicillinase), was identified in five strains (5A3I, 11A1I, 30A2I, 48A3I and 50A2C). The operon was located in small contigs of 6.5-8.9 kbp within the same genetic context and surrounded by plasmid-associated genes. Electrotransformation of competent S. equorum cells with plasmid DNA from all five strains produced no viable ampicillin- or penicillin-resistant transformants. In addition to blaZ, CARD, NCBI-RGC, and ResFinder searches identified another gene in S. equorum 35A3C coding for another putative class A betalactamase: bla. Inspection of PATRIC annotations revealed bla to be present in all the studied genomes. The gene shared the identical genetic context in all strains and was flanked by chromosomally encoded genes. Nonetheless, bla gene was polymorphic and the protein variants (Figure S5A) may have differential activity.
An mph(C) gene was detected for two strains (8A3C and 16A1C) in database comparisons. Both strains also carried an msr(A) gene also found in five other strains (T17, 2A3C, 23A3C, 35A3C, and 50A2C). mph(C) and msr(A) code for, respectively, a macrolide 2ʼ-phosphotransferase (inactivating the antibiotic) and an ABC-F type ribosomal protection protein (protecting the target). Nucleotide identity of the mph(C) gene lead to a single Mph(C) protein, while Mrs(A) showed several variants, of which two (those of 8A3C and 16A1C, and of 50A2C) were associated with high erythromycin resistance (MIC 24-64 μg mL-1; Figure S5B). All these strains, however, proved to be clindamycin-susceptible, even though two strains (1BCExtra and 2A3C) also carried an additional plasmid-encoded gene homologous to the lnu(A) gene from Staphylococcus species. lnu(A) codes for lincosamide nucleotidyltransferase and has been shown involved in resistance to lincosamides [42,43]. The presence of lnu(A) leads to phenotypic examination for lincomycin resistance in all 13 strains. The MICs for this antibiotic varied widely (from 1 to 64 µg mL-1), with no correlation seen with respect to the presence or absence of lnu(A). Finally, although the MICs for ciprofloxacin were low, database searches identified a polymorphic norA gene in all strains (Figure S5C). This chromosomal gene codes for a housekeeping efflux pump of the major facilitator superfamily (MFS) known to be involved in fluoroquinolone resistance in the Staphylococcus genus [44].
Although fosfomycin was not initially tested, a fosB/fosD polymorphic gene coding for a fosfomycin-inactivating enzyme was identified in the chromosome of all strains; the open reading frame (ORF) was disrupted in some of them. Phenotypic testing of the strains led by this finding showed a surprisingly wide range of resistance to fosfomycin (MICs from 2 to >1024 µg mL-1). Further, the strain showing the greatest resistance (1BCExtra) contained two copies of the gene, although one was disrupted. In fact, the complete copy of fosB/fosD in 1BCExtra (1BCExtra-1) was identified by database searches to be an acquired gene. The copy was deemed to be chromosomally encoded, but it was flanked on one end by Tn552-associated ORFs. Further, the protein variant of strain 1BCExtra proved to be the most dissimilar to all others (Figure S5D).
In contrast to the complex situation of most antibiotics, a cat gene matching phenotypic resistance to chloramphenicol was detected in strain 35A3C. Genes supporting moderate phenotypic resistance to chloramphenicol (MIC 8-16 µg mL-1) were not identified. The cat gene encodes a type A-7 chloramphenicol o-acetyltransferase. It was located in a contig of 4.6 kbp harbouring genes (repB, mobCAB) and signals (dso, sso, oriT) involved in plasmid replication and transfer (Figure S6A). The Cat protein was identical to that encoded on pC221 and very similar to that of pC223, both plasmids from S. aureus (Figure S6B). Further, the contig and the whole pC223 (GenBank accession no. NC_005243.1) were very similar at the nucleotide level (92% identity and 98% length coverage). By transforming the plasmid complement of 35A3C into S. aureus cells, chloramphenicol-resistant colonies were obtained. All transformants contained a plasmid with a restriction enzyme digestion profile compatible with that determined in silico for the cat contig. These results indicate that the contig contained the complete sequence of the cat plasmid of 35A3C (pCAT). The MIC for chloramphenicol increased from 4 µg mL-1 to 32 µg mL-1 in S. aureus pCAT-containing transformants. In contrast, transformants were never recovered when the plasmid DNA of this strain was electroporated into Lactococcus lactis, Enterococcus faecalis, or Escherichia coli cells, suggesting pCAT does not replicate in these hosts.
Finally, the use of antibiotic resistance databases identified several genes (qacJ, norC, sepA, sdrM, and qacJ) coding for multidrug efflux transporters similar to proteins involved in resistance to disinfectants and antiseptics. Of these, qacJ from CL10P returned identity and length coverage percentages to genes in databases of 98.13-100% and 70.14-95.54%, respectively. qacJ was located in a small contig of (2323 bp) containing an ORF coding for a replication protein with a nucleotide identity of 93% to the S. aureus plasmid pNVH01 (NC_004562.1).

3. Discussion

Commensal and beneficial food-borne bacteria may play a role in disseminating antibiotic resistance to pathogens by acting as reservoirs of genetic determinants [45]. A comprehensive understanding of the whole pool of genes available and the transfer process is essential if we are to combat this resistance [35]. S. equorum populations can become large in a number of fermented foods [21,23,24,25,26,28], yet the antibiotic resistance properties of this speciesʼ strains have been little studied [38,46,47]. In this work, strains of S. equorum isolated from a Spanish traditional blue-veined cheese were subjected to phenotypic, genetic, and genomic analyses in order to determine their antibiotic R/S profiles, to identify the genetic basis of the phenotypic resistances noted, and to search their genomes for antibiotic resistance-associated genes. As in other bacterial populations from cheese [48,49,50], the combined fingerprinting profiles revealed wide genetic diversity, which was confirmed by phylogenomic analysis of the sequenced strains. The examination of unrelated, non-clonal strains is important if the results of R/S assays results are to be reliable and extendible to new S. equorum strains.
In agreement with previous studies [37,39,51], low MICs were returned for most antibiotics in the present work; indeed, all were considered susceptible. Some MICs, however, fell above the cut-offs for Staphylococcus spp. established by the EUCAST [40], although they were below those established by the CLSI [41]; these strains were, therefore, also considered susceptible to the antibiotics concerned. Further genetic analyses of the strains showing the highest MICs might be advisable to exclude the presence of low-resistance conferring systems [52,53]. The MICs above the CLSI cut-off to penicillin G correlated well with an enhanced MIC for ampicillin, although a cut-off for this latter antibiotic has yet to be established. EUCAST has not developed cut-offs for ampicillin and penicillin G due to widespread resistance to these antibiotics among staphylococcal isolates, which makes their clinical use impractical [40]. The MIC values for ampicillin and penicillin G obtained with the microdilution assays were not reproduced by the MTS system. Despite of this, the largest MICs were always exhibited by the same strains, suggesting a differential resistance as compared to more susceptible ones. This resistance level, however, was much lower than that recorded in clinically-resistant S. aureus strains (16-256 µg mL-1) [54,55]. Low resistance might be provided by dissimilar cell structures (cell-wall permeability, membrane charges, penicillin binding protein variants, etc.) or the differential activity of unspecific transporters (multidrug systems, efflux pumps) [56]. Anyway, the bimodal topology of the distribution curves for these antibiotics (Figure S2), strongly suggests that a part of the population has got acquired resistance. The presence of colonies within the inhibition halos points towards undefined mutations accounting for such resistance. In contrast to β-lactams, some strains showed different (and quite high) MIC values for erythromycin and chloramphenicol. These values were considered compatible with the presence of dedicated resistance systems and the involvement of acquired (and possibly transferable) genes.
In the absence of species-specific R/S cut-offs, authorized agencies recommend the use of those of closely related species [40,57]. However, knowing whether a given strain is susceptible or resistant to an antibiotic is still challenging, particularly if the cut-offs cited by different sources (such as those of EUCAST and CLSI) are dissimilar [40,41]. In addition to clinical R/S cut-offs, microbiological (MCOFFs) or ecological cut-offs (ECOFFs) may also serve to distinguish susceptible from resistant strains [36,58]. MCOFFs and ECOFFs describe the MIC above which bacterial isolates have phenotypically detectable acquired resistance mechanisms [59]. This should better allow the identification of strains carrying transferable resistances, and thus help prevent their use in food (and feed) systems [60]. Examining the data in the literature has the potential to help establish S. equorum such microbiological or ecological cut-offs, but differences between surveys in terms of the antibiotics tested, the concentration ranges contemplated, the methodologies used (disk diffusion, plate, and broth microdilution, Etest, etc.), and the media and culture conditions imposed, etc., hamper any direct comparisons of results [24,37,38,39,51]. Analysis of the distribution of MICs from large sets of antibiotic resistance assays using software programs such as ECOFFinder [61] or NRI [62] were utilized in this study to propose S. equorum-specific cut-offs. The incorporation of new MIC data by surveying more strains from different sources can help establish robust and reliable cut-offs for this species.
Genome sequencing and comparative analysis are currently the gold standards for characterizing the genetic potential of microorganisms, including the detection of ARGs in bacteria [63], the prediction of phenotypic resistance [64,65], and the assessment of its transferability [66]. Comparing the genomes of antibiotic-resistant and susceptible bacteria would further allow distinctions to be made between genes known to be spread through bacterial species and housekeeping genes involved in antibiotic resistance [67]. The distribution of the resistome in chromosomal and plasmid contigs has recently been reported following genome analysis of different Staphylococcus species [68,69,70], including S. equorum strains [71]. Chromosomal genes are found on large contigs harbouring well-recognized housekeeping genes, while genes on plasmids are usually found in small contigs and in the vicinity of genes coding for replication proteins or proteins involved in plasmid maintenance and mobilization. The genome analysis undertaken in the present work led to further phenotypic testing for antibiotics that were not included in the Sensititre panel. It also allowed for acquired but silent genes to be detected. In the absence of antibiotic pressure, the acquisition of tightly regulated genes (silent genes) that do not lead to the phenotypic expression of resistance have been abundantly described [47], of these some might easily activate under certain conditions [72]. The transfer of silent genes to clinically important pathogens could lead to therapy failure [73].
Genes of the bla and mec families have been repeatedly reported to mediate resistance in Staphylococcus spp. to β-lactam antibiotics [74]. In the S. equorum strains of this work, only genes of the bla family were identified: an identical blaR1-blaZI operon on plasmid contigs in five strains, and a bla gene in all other strains on the bacterial chromosome. Whether the contigs harbour the complete plasmid molecules is not currently known. Neither blaZ nor bla, however, afforded moderate or strong resistance to ampicillin or penicillin G. Indeed, the strains carrying the genes for both betalactamases were among the most susceptible to both antibiotics. These results lead us to foresee that the phenotypic resistance to penicillin G is not due to acquired genes. Even though large plasmids have been recently introduced in Gram-positive bacteria by electroporation [75], no colonies were obtained by transforming plasmid DNA from the five strains into S. aureus, suggesting the S. equorum plasmidic BlaZ system is not functional. As reported for S. aureus, sequence variation in the blaR1-blaZI region [47,76], or mutations in the promoter region [77], may also account for the discrepancies observed between genotype and phenotype. As pointed out above, despite being silent, acquired genes such as blaZ and lnu(A) may still represent a hazard.
The msr(A) and fosB/fosD genes proved to be polymorphic. The protein variants encoded by the different ORFs were thought to be responsible for the different MICs recorded for erythromycin and fosfomycin. As regards erythromycin resistance, the simultaneous presence of msr(A) and mph(C) in a single strain, as might be the case in some of the present strains, has previously been associated with enhanced MICs [43]. norA encodes an efflux pump of the major facilitator superfamily that has been associated with resistance to fluoroquinolones in Staphylococcus [44]. However, this gene in the S. equorum strains was not associated with ciprofloxacin resistance. Lüthje and Schwarz [78] propose norA to be a core gene in staphylococci and suggest fusaric acid and siderophores be its natural substrates.
Plasmids are the major route of dissemination of resistance determinants; therefore, resistance genes harboured on plasmids are the most likely to be transferred [79]. Small plasmids conferring chloramphenicol resistance via cat genes coding for chloramphenicol o-acetyl transferases have been abundantly characterized in Staphylococcus species of different origins [80,81]. Plasmids such as pC221 and pC223 belong to the rolling-circle replicating (RCR) plasmids of the pT181 family and have been shown to replicate in several Gram-positive species, including staphylococci and Bacillus subtilis [82]. In addition to cat, blaR1-blaZI, and lnu(A) determinants were also found to be plasmid-encoded in this work. The in vitro conjugal transfer of lnuA-containing plasmids, and their lincomycin-associated resistance, from S. equorum to Staphylococcus spp., E. faecalis and Tetragenococcus halophilus, has already been reported [46]. lnuA-containing plasmids from S. equorum cheese strains are larger (contigs of 32-34 kbp) than those in S. equorum strains from fermented seafood (2.6-2.8 kbp; these last plasmids also replicate by an RCR mechanism) [46]. Differences in the size and mode of replication can influence the transfer and host range of plasmids, which may limit their spread capability. In the present work, pCAT was successfully transferred to S. aureus by electroporation. The presence of an oriT sequence followed by ORFs encoding Mob proteins further suggests pCAT has a capacity for mobilization by other means (e.g., conjugation). Certainly, the 1BCExtra-1 fosB/fosD gene, which is thought to lie on the chromosome, was flanked by ORFs encoding Tn552-associated proteins, strongly suggesting horizontal acquisition by transposon-mediated integration. Indeed, Tn552 has long been recognized for its capacity to integrate plasmids on the chromosome [83]. Although posing a lesser risk, this fosfomycin resistance gene is thought to be easily transmitted via horizontal transfer events. Together, the above results suggest that S. equorum is already a reservoir of easily transferable ARGs. Therefore, full phenotypic and genetic characterizations of any S. equorum strains being contemplated for use as starters or adjunct cultures in food fermentations are advisable.

4. Materials and Methods

4.1. Bacterial strains and culture conditions

A total of 53 S. equorum isolates from traditional blue-veined Cabrales cheese were initially investigated. These isolates were routinely cultured in Brain Heart Infusion (BHI) broth (VWR International, Darmstadt, Germany) supplemented with 2.5% (w/v) NaCl (3% in total), at 32ºC for 48 h. Staphylococcus aureus RN4220 was grown without agitation at 37ºC in Tryptone Sosy Broth (TSB; Scharlab, Barcelona, Spain). Enterococcus faecalis 52c and Lactococcus lactis NZ9000 were grown at 32ºC in M17 broth (Formedium, Swaffham, UK) supplemented with 1% (w/v) glucose (GM17). Escherichia coli DH10B was grown in 2xTY broth at 37ºC with shaking. When required, agar (2% w/v) was added to the broth media for the preparation of solid plates.

4.2. Typing of the strains

Total genomic DNA was isolated using the GenElute Bacterial Genomic DNA Kit (Sigma-Aldrich, St. Louis, CA, USA) following the manufacturerʼs recommendations for Gram-positive bacteria. Genotyping of the strains was performed by combining the fingerprinting profiles obtained with primers BoxA2R (5ʼ-ACGTGGTTTGAAGAGATTTTCG-3ʼ), OPA18 (5ʼ-AGGTGACCGT-3ʼ) and M13 (5ʼ-GAGGGTGGCGGTTCT-3ʼ), as previously reported [49]. Clustering of the profiles was performed using the Unweighted Pair Group Method with Arithmetic Means algorithm (UPGMA) and Jaccard similarity coefficients, using MVSP v.3.21 software (Kovach Computing Services, Pentraeth, UK).

4.3. Antibiotic testing

The MICs for 16 antibiotics in S. equorum strains was determined by broth microdilution using Sensititre EULACBI1 and EULACBI2 plates (Trek Diagnostic Systems, East Grinstead, UK). Briefly, colonies grown for 48 h on Mueller-Hinton (M-H) (Oxoid, Basingstoke, UK) agar plates were used to prepare cell suspensions in 0.9% NaCl solution (density corresponding to McFarland standard 1; i.e., ~3 x 108 cfu mL-1). The suspensions were then further diluted 1:750 in M-H to achieve a final concentration of about 4 x 105 cfu mL-1. One hundred microlitres were then inoculated into the wells of the Sensititre plates and incubated at 32ºC for 48 h. MICs for ampicillin, penicillin G, and fosfomycin (Sigma-Aldrich), and lincomycin (ThermoFisher) were also determined by microdilution in M-H broth using two-fold antibiotic dilutions. MICs were defined as the lowest antibiotic concentration at which no growth was observed.
When the concentration range of an antibiotic on the Sensititre plates was unable to pinpoint the MIC, and when MIC values had to be confirmed, MIC Test Strips (MTS; Liofilchem, Roseto degli Abruzzi, Italy) were used. This involved dipping a sterile cotton swab into a cell suspension corresponding to McFarland standard 0.5 (~1.5 x 108 cfu mL-1) and using it to inoculate the surface of an M-H agar plate. After drying for 15 min, the antibiotic-containing strip was placed on the plate and incubated at 32ºC for 48 h. MICs for ampicillin and penicillin G were further confirmed with MTS.
Strains were considered susceptible or resistant based on the breakpoints recommended for Staphylococcus spp. by the European Committee on Antimicrobial Susceptibility Testing [40], and by the American Clinical and Laboratory Standards Institute, as recovered from the M100-ED32:2023 supplement [41]. A strain was classified as resistant when the MIC exceed the breakpoint in more than one dilution.

4.4. MIC analysis and tentative R/S cut-offs

To suggest some tentative cut-offs, the distribution of the presently determined and literature-reported MICs was analysed using ECOFFinder software (https://clsi.org/meetings/microbiology/ecoffinder/) and the Normalized Resistance Interpretation (NRI) spreadsheet (http://www.bioscand.se/nri/). Both methods are based on statistics involving non-linear regression curves, as reported by Turnidge et al. [61], and Kronvall [62], respectively.

4.5. Whole genome sequencing and analysis

Genomic DNA from overnight cultures of S. equorum was extracted using the QIAmp DNA Mini Kit (Qiagen, Hilden, Germany) following the manufacturerʼs protocol. Next-generation sequencing libraries were prepared using standard methods and paired-end (2*150 bp) sequenced in a NovaSeq 6000 sequencer at GATC (Eurofins, Ebersberg, Germany). Genomes were annotated using PATRIC services (https://www.patricbrc.org/). For this, reads were first checked for quality with FastQC software, and assembled in contigs using the Unicycler program after comparing different variables using the Quality Assessment Tool for Genome Assemblies (QUAST). Errors were polished using Pilon and Racon software. Annotation with PATRIC was performed using the RAST tool kit (RASTtk). Antibiotic resistance was investigated by genome analysis and comparison against sequences in the CARD database (https://card.mcmaster.ca/) (80% identity, 70% length coverage), the NCBI AMR Reference Gene Catalogue (NCBI-RGC) (using AMRFinderPlus: https://www.ncbi.nlm.nih.gov/pathogens/antimicrobial-resistance/AMRFinder/), and the ResFinder database (https://cge. cbs.dtu.dk/services/ResFinder/) (80% protein identity, 60% length coverage). PATRIC annotation data were also examined to seek out ARGs. Using BLAST tools (https://blast.ncbi.nlm.nih.gov/Blast.cgi), DNA and deduced protein sequences of interest were individually compared (at the nucleotide and amino acid levels, respectively) against information in the NCBI database. Whole-genome sequence data were used to ascertain the phylogenetic relationships between the sequenced strains and the type strains of related staphylococcal species by means of digital DNA-DNA hybridization (dDDH) and orthologous average nucleotide identity (orthoANI), as reported by Meier-Kolthoff and Göker [84] and Yoon et al. [85], respectively.

4.6. Phylogenetic and phylogenomic analyses

Phylogenomic analysis of the strains with type strains of Streptococcus species were carried out at the DSMZ Type Strain Genome Server (http://ggdc.dsmz.de/). Phylogenomic comparison of the sequenced strains with six selected S. equorum strains from NCBI database was carried out at the BV-BRC Server (https://www.bv.brc.org) using 500 genes. The phylogenetic tree for Cat proteins was constructed using sequences retrieved from the NCBI database (https://www.ncbi.nlm.nih.gov/protein/) aligned with MUSCLE and using the maximum likelihood estimation test. Alignment and phylogenetic reconstructions of other proteins were performed using the function “build” of ETE3 3.1.2 as implemented on GenomeNet (https://www.genome.jp/tools/ete/).

4.7. Isolation and transformation of plasmid DNA

Plasmid DNA from S. equorum strains was extracted and purified as described by OʼSullivan and Klaenhammer [86] with minor modifications. Briefly, 4 µL of mutanolysin (5 U µL-1), 5 µL of lysostaphin (2 mg mL-1), 20 µL of proteinase K (30 mg mL-1), and 4 µL of RNase (20 mg mL-1) were added to the lysis buffer. Electrocompetent cells of E. coli, S. aureus, L. lactis, and E. faecalis were prepared as reported by Sambrook and Russell [87], Schneewind and Missiakas [88], Holo and Nes [89], and Pérez et al. [90], respectively. Electrotransformation (electroporation) was performed using a Gene Pulser device (Bio-Rad, Richmond, CA, USA) following, as required, standard protocols for Gram-positive or Gram-negative bacteria. Transformants were selected on appropriate agar plates supplemented with antibiotics (chloramphenicol 10 μg mL-1; ampicillin and penicillin G 2, 4, and 8 μg mL-1), and, if recovered, analysed to determine their plasmid content.

5. Conclusions

S. equorum strains from cheese phenotypically showed little clear-cut resistance to the tested antibiotics, and the concomitant presence of phenotypic resistance and supporting genes was only observed in a minority of strains. Acquired resistance by mutations was thought to provide low resistance to ampicillin and penicillin G and high resistance to erythromycin. Acquired genes -some of which were silent- known to be spread across diverse bacterial groups, involved in resistance to β-lactams, chloramphenicol, and fosfomycin were found in the genome of some strains. Multi-resistant strains or strains carrying more than one acquired gene were not detected. Strains showing genes on large contigs containing well-recognized chromosomally-encoded genes were identified. Others were found in smaller contigs and in the vicinity of genes coding for plasmid-replication proteins or proteins involved in plasmid maintenance or mobilization. These were deemed to be plasmid-encoded and thought to pose the greatest risk of horizontal transfer. Indeed, both the origin of replication of pCAT and its associated cat gene were found to be functional in S. aureus. Expressed or silent, starter candidates of S. equorum should be free of acquired ARGs, demanding their thorough phenotypic and genomic characterization before use.

Supplementary Materials

The following supporting information can be downloaded at the website of this paper posted on Preprints.org, Figure S1: Dendrogram of similarity for the S. equorum isolates from cheese; Figure S2: ECOFFinder analysis of the MICs distribution and proposed R/S cut-offs for S. equorum; Figure S3: Phylogenomic analysis of the S. equorum strains of this study with the type strains of several staphylococci species.; Figure S4: Phylogenomic analysis of the S. equorum strains from cheese and others from databases; Figure S5: Phylogenetic trees of proteins encoded by the antibiotic resistance genes; Figure S6: Genetic map of pCAT and phylogenetic tree of Cat proteins. Table S1: MIC values of 16 antibiotics to 30 S equorum strains from cheese; Table S2: Genomic Distance Analysis of Staphylococcus species type strains with the sequenced S. equorum strains of this study.

Author Contributions

Conceptualization, B.M. and A.B.F.; Methodology, L.V., A.B.F., and B.M. Investigation, L.V., J.R.; Data analysis, L.V., M.E.S., and B.M.; Original draft preparation, L.V. and B.M. Writing, review and editing, L.V., M.E.S., J.R., A.B.F., and B.M.; Funding acquisition, B.M.

Funding

This research was funded by projects from the Spanish Ministry of Science and Innovation (PID2019-110549RB-I00/AEI/10.13039/501100011033) and Asturias Principality (AYUD/2021/50916). L.V. was awarded a postdoctoral research contract from the Health Research Institute of Asturias Principality (ISPA-FINBA). A Ph.D. grant of the Severo Ochoa Program from Asturias Principality was awarded to J.R. (BP19-098).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The complete genome sequences of the 13 S. equorum strains examined were deposited in the GenBank database under the BioProject and BioSample accession numbers PRJNA940711, and SAMN33577425 through SAMN33577437.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hutchings, M.I.; Truman, A.W.; Wilkinson, B. Antibiotics: past, present and future. Curr. Opin. Microbiol. 2019, 51, 72–80. [Google Scholar] [CrossRef] [PubMed]
  2. Hossain, A.; Habibullah-Al-Mamun, M.; Nagano, I.; Masunaga, S.; Kitazawa, D.; Matsuda, H. Antibiotics, antibiotic-resistant bacteria, and resistance genes in aquaculture: risks, current concern, and future thinking. Environ. Sci. Pollut. Res. Int. 2022, 29, 11054–11075. [Google Scholar] [CrossRef] [PubMed]
  3. Nunziata, L.; Brasca, M.; Morandi, S.; Silvetti, T. Antibiotic resistance in wild and commercial non-enterococcal lactic acid bacteria and bifidobacteria strains of dairy origin: An update. Food Microbiol. 2022, 104, 103999. [Google Scholar] [CrossRef] [PubMed]
  4. Penders, J.; Stobberingh, E.E.; Savelkoul, P.H.; Wolffs, P.F. The human microbiome as a reservoir of antimicrobial resistance. Front. Microbiol. 2013, 4, 87. [Google Scholar] [CrossRef] [PubMed]
  5. Flórez-Cuadrado, D.; Moreno, M.A.; Ugarte-Ruiz, M.; Domínguez, L. Antimicrobial resistance in the food chain in the European Union. Adv. Food Nutr. Res. 2018, 86, 115–136. [Google Scholar] [CrossRef] [PubMed]
  6. Wong, A.; Matijasic, B.B.; Ibana, J.A.; Lim, R.L.H. Editorial: Antimicrobial resistance along the food chain: Are We What We Eat? Front. Microbiol. 2022, 13, 881882. [Google Scholar] [CrossRef] [PubMed]
  7. Crits-Christoph, A.; Hallowell, H.A.; Koutouvalis, K.; Suez, J. Good microbes, bad genes? The dissemination of antimicrobial resistance in the human microbiome. Gut Microbes 2022, 14, 2055944. [Google Scholar] [CrossRef] [PubMed]
  8. Kumar, S.B.; Arnipalli, S.R.; Ziouzenkova, O. Antibiotics in food chain: The consequences for antibiotic resistance. Antibiotics (Basel) 2020, 9, 688. [Google Scholar] [CrossRef]
  9. Lambraki, I.A.; Cousins, M.; Graells, T.; Léger, A.; Henriksson, P.; Harbarth, S.; Troell, M.; Wernli, D.; Søgaard Jørgensen, P.; Desbois, A.P.; Carson, C.A.; Parmley, E.J.; Majowicz, S.E. Factors influencing antimicrobial resistance in the European food system and potential leverage points for intervention: A participatory, One Health study. PLoS One 2022, 17, e0263914. [Google Scholar] [CrossRef]
  10. Alexa Oniciuc, E.A.; Walsh, C.J.; Coughlan, L.M.; Awad, A.; Simon, C.A.; Ruiz, L.; Crispie, F.; Cotter, P.D.; Alvarez-Ordóñez, A. Dairy products and dairy-processing environments as a reservoir of antibiotic resistance and quorum-quenching determinants as revealed through functional metagenomics. mSystems 2020, 5, e00723–19. [Google Scholar] [CrossRef]
  11. Blanco-Picazo, P.; Gómez-Gómez, C.; Morales-Cortes, S.; Muniesa, M.; Rodríguez-Rubio, L. Antibiotic resistance in the viral fraction of dairy products and a nut-based milk. Int. J. Food Microbiol. 2022, 367, 109590. [Google Scholar] [CrossRef] [PubMed]
  12. Flórez, A.B.; Alegría, A.; Rossi, F.; Delgado, S.; Felis, G.E.; Torriani, S.; Mayo, B. Molecular identification and quantification of tetracycline and erythromycin resistance genes in Spanish and Italian retail cheeses. Biomed. Res. Int. 2014, 2014, 746859. [Google Scholar] [CrossRef] [PubMed]
  13. Wang, W.; Wei, X.; Wu, L.; Shang, X.; Cheng, F.; Li, B.; Zhou, X.; Zhang, J. The occurrence of antibiotic resistance genes in the microbiota of yak, beef and dairy cattle characterized by a metagenomic approach. J. Antibiot. (Tokyo) 2021, 74, 508–518. [Google Scholar] [CrossRef] [PubMed]
  14. Ammor, M.S.; Flórez, A.B.; Mayo, B. Antibiotic resistance in non-enterococcal lactic acid bacteria and bifidobacteria. Food Microbiol. 2007, 24, 559–570. [Google Scholar] [CrossRef] [PubMed]
  15. Álvarez-Cisneros, Y.M.; Ponce-Alquicira, E. Antibiotic resistance in lactic acid bacteria. In Antimicrobial Resistance. A global threat; Kumar, Y., Ed.; IntechOpen: London, UK, 2018; Chapter 4; pp. 53–73. [Google Scholar] [CrossRef]
  16. Sagar, P.; Aseem, A.; Banjara, S.K.; Veleri, S. The role of food chain in antimicrobial resistance spread and One Health approach to reduce risks. Int. J. Food Microbiol. 2023, 391-393, 110148. [Google Scholar] [CrossRef] [PubMed]
  17. 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]
  18. Quigley, L.; OʼSullivan, O.; Beresford, T.P.; Ross, R.P.; Fitzgerald, G.F.; Cotter, P.D. Molecular approaches to analysing the microbial composition of raw milk and raw milk cheese. Int. J. Food Microbiol. 2011, 150, 81–94. [Google Scholar] [CrossRef] [PubMed]
  19. Dugat-Bony, E.; Garnier, L.; Denonfoux, J.; Ferreira, S.; Sarthou, A.S.; Bonnarme, P.; Irlinger, F. Highlighting the microbial diversity of 12 French cheese varieties. Int. J. Food Microbiol. 2016, 238, 265–273. [Google Scholar] [CrossRef]
  20. Yeluri Jonnala, B.R.; McSweeney, P.L.H.; Sheehan, J.J.; Cotter, P.D. Sequencing of the cheese microbiome and its relevance to industry. Front. Microbiol. 2018, 9, 1020. [Google Scholar] [CrossRef]
  21. Jian, Z.; Zeng, L.; Xu, T.; Sun, S.; Yan, S.; Yang, L.; Huang, Y.; Jia, J.; Dou, T. Antibiotic resistance genes in bacteria: Occurrence, spread, and control. J. Basic Microbiol. 2021, 61, 1049–1070. [Google Scholar] [CrossRef]
  22. EUCAST. MIC distributions and epidemiological cut-off value (ECOFF) setting. EUCAST SOP 10.2, 2021 (http://www.eucast.org).
  23. Leroy, S.; Lebert, I.; Chacornac, J.P.; Chavant, P.; Bernardi, T.; Talon, R. Genetic diversity and biofilm formation of Staphylococcus equorum isolated from naturally fermented sausages and their manufacturing environment. Int. J. Food Microbiol. 2009, 134, 46–51. [Google Scholar] [CrossRef] [PubMed]
  24. Cachaldora, A.; Fonseca, S.; Franco, I.; Carballo. J. Technological and safety characteristics of Staphylococcaceae isolated from Spanish traditional dry-cured sausages. Food Microbiol. 2013, 33, 61–8. [Google Scholar] [CrossRef] [PubMed]
  25. Coton, E.; Desmonts, M.H.; Leroy, S.; Coton, M.; Jamet, E.; Christieans, S.; Donnio, P.Y.; Lebert, I.; Talon, R. Biodiversity of coagulase-negative staphylococci in French cheeses, dry fermented sausages, processing environments and clinical samples. Int. J. Food Microbiol. 2010, 137, 221–229. [Google Scholar] [CrossRef] [PubMed]
  26. Jeong, D.W.; Kim, H.R.; Lee, J.H. Genetic diversity of Staphylococcus equorum isolates from Saeu-jeotgal evaluated by multilocus sequence typing. Antonie van Leeuwenhoek 2014, 106, 795–808. [Google Scholar] [CrossRef] [PubMed]
  27. Meugnier, H.; Bes, M.; Vernozy-Rozand, C.; Mazuy, C.; Brun, Y.; Freney, J.; Fleurette, J. Identification and ribotyping of Staphylococcus xylosus and Staphylococcus equorum strains isolated from goat milk and cheese. Int. J. Food Microbiol. 1996, 31, 325–31. [Google Scholar] [CrossRef] [PubMed]
  28. Deetae, P.; Bonnarme, P.; Spinnler, H.E.; Helinck, S. Production of volatile aroma compounds by bacterial strains isolated from different surface-ripened French cheeses. Appl. Microbiol. Biotechnol. 2007, 76, 1161–1171. [Google Scholar] [CrossRef] [PubMed]
  29. Schornsteiner, E.; Mann, E.; Bereuter, O.; Wagner, M.; Schmitz-Esser, S. Cultivation-independent analysis of microbial communities on Austrian raw milk hard cheese rinds. Int. J. Food Microbiol. 2014, 180, 88–97. [Google Scholar] [CrossRef]
  30. Unno, R.; Matsutani, M.; Suzuki, T.; Kodama, K.; Matsushita, H.; Yamasato, K.; Koizumi, Y.; Ishikawa, M. Lactic acid bacterial diversity in Brie cheese focusing on salt concentration and pH of isolation medium and characterisation of halophilic and alkaliphilic lactic acid bacterial isolates. Int. Dairy J. 2020, 109, 104757. [Google Scholar] [CrossRef]
  31. Landete, G.; Curiel, J.A.; Carrascosa, A.V.; Muñoz, R.; de las Rivas, B. Characterization of coagulase-negative staphylococci isolated from Spanish dry cured meat products. Meat Sci. 2013, 93, 387–96. [Google Scholar] [CrossRef]
  32. Stavropoulou, D.A.; De Vuyst, L.; Leroy, F. Nonconventional starter cultures of coagulase-negative staphylococci to produce animal-derived fermented foods, a SWOT analysis. J. Appl. Microbiol. 2018, 25, 1570–1586. [Google Scholar] [CrossRef]
  33. Lee, J.-H.; Heo, S.; Jeong, D.-W. Genomic insights into Staphylococcus equorum KS1039 as a potential starter culture for the fermentation of high-salt foods. BMC Genomics 2018, 19, 136. [Google Scholar] [CrossRef] [PubMed]
  34. Place, R.B.; Hiestand, D.; Gallmann, H.R.; Teuber, M. Staphylococcus equorum subsp. linens, subsp. nov., a starter culture component for surface ripened semi-hard cheeses. Syst. Appl. Microbiol. 2003, 26, 30–7. [Google Scholar] [CrossRef]
  35. Rodríguez-González, M.; Fonseca, S.; Centeno, J.A.; Carballo, J. Biochemical changes during the manufacture of Galician chorizo sausage as affected by the addition of autochthonous starter cultures. Foods 2020, 9, 1813. [Google Scholar] [CrossRef] [PubMed]
  36. Irlinger, F. Safety assessment of dairy microorganisms: coagulase-negative staphylococci. Int. J. Food Microbiol. 2008, 126, 302–10. [Google Scholar] [CrossRef] [PubMed]
  37. Marty, E.; Bodenmann, C.; Buchs, J.; Hadorn, R.; Eugster-Meier, E.; Lacroix, C.; Meile, L. Prevalence of antibiotic resistance in coagulase-negative staphylococci from spontaneously fermented meat products and safety assessment for new starters. Int. J. Food Microbiol. 2012, 159, 74–83. [Google Scholar] [CrossRef] [PubMed]
  38. Mikulášová, M.; Valáriková, J.; Dušinský, R.; Chovanová, R.; Belicová, A. Multiresistance of Staphylococcus xylosus and Staphylococcus equorum from Slovak Bryndza cheese. Folia Microbiol. (Praha) 2014, 59, 223–227. [Google Scholar] [CrossRef] [PubMed]
  39. Resch, M.; Nagel, V.; Hertel, C. Antibiotic resistance of coagulase-negative staphylococci associated with food and used in starter cultures. Int. J. Food Microbiol. 2008, 127, 99–104. [Google Scholar] [CrossRef]
  40. EUCAST. Breakpoint tables for interpretation of MICs and zone diameters. Version 13.0, valid from 2023-01-01. https://www.eucast.org/clinical_breakpoints/, 2023.
  41. CLSI. Performance standards for antimicrobial susceptibility testing, 33rd ed. CLSI supplement M100. Clinical and Laboratory Standards Institute; 2023.
  42. Olsen, J.E.; Christensen, H.; Aarestrup, F.M. Diversity and evolution of blaZ from Staphylococcus aureus and coagulase-negative staphylococci. J. Antimicrob. Chemother. 2006, 57, 450–60. [Google Scholar] [CrossRef]
  43. Lüthje, P.; von Kockritz-Blickwede, M.; Schwarz, S. Identification and characterization of small staphylococcal plasmids carrying the lincosamide nucleotidyltransferase gene lnu(A). J. Antimicrob. Chemother. 2007, 59, 600–606. [Google Scholar] [CrossRef]
  44. Ferreira, C.; Abrantes, P.; Costa, S.S.; Viveiros, M.; Couto, I. Occurrence and variability of the efflux pump gene norA across the Staphylococcus genus. Int. J. Mol. Sci. 2022, 23, 15306. [Google Scholar] [CrossRef]
  45. Wang, H.H.; Manuzon, M.; Lehman, M.; Wan, K.; Luo, H.; Wittum, T.E.; Yousef, A.; Bakaletz, L.O. Food commensal microbes as a potentially important avenue in transmitting antibiotic resistance genes. FEMS Microbiol. Lett. 2006, 254, 226–31. [Google Scholar] [CrossRef] [PubMed]
  46. Lee, J.-H.; Jeong, D.-W. Characterization of mobile Staphylococcus equorum plasmids isolated from fermented seafood that confer lincomycin resistance. PLoS One 2015, 10, e140190. [Google Scholar] [CrossRef]
  47. Seitter, M.; Nerz, C.; Rosenstein, R.; Götz, F.; Hertel, C. DNA microarray-based detection of genes involved in safety and technologically relevant properties of food associated coagulase-negative staphylococci. Int. J. Food Microbiol. 2011, 145, 449–58. [Google Scholar] [CrossRef] [PubMed]
  48. Frantzen, C.A. , Kleppen, H.P., Holo, H. Lactococcus lactis diversity in undefined mixed dairy starter cultures as revealed by comparative genome analyses and targeted amplicon sequencing of epsD. Appl. Environ. Microbiol. 2018, 84, e02199-17. [Google Scholar] [CrossRef] [PubMed]
  49. Rodríguez, J.; González-Guerra, A.; Vázquez, L.; Fernández-López, R.; Flórez, A.B.; de la Cruz, F.; Mayo, B. Isolation and phenotypic and genomic characterization of Tetragenococcus spp. from two Spanish traditional blue-veined cheeses made of raw milk. Int. J. Food Microbiol. 2022, 16, 371:109670. [Google Scholar] [CrossRef] [PubMed]
  50. Siezen, R.J.; Tzeneva, V.A.; Castioni, A.; Wels, M.; Phan, H.T.; Rademaker, J.L.; Starrenburg, M.J.; Kleerebezem, M.; Molenaar, D.; van Hylckama Vlieg, J.E. Phenotypic and genomic diversity of Lactobacillus plantarum strains isolated from various environmental niches. Environ. Microbiol. 2010, 12, 758–773. [Google Scholar] [CrossRef]
  51. Even, S.; Leroy, S.; Charlier, C.; Zakour, N.B.; Chacornac, J.P.; Lebert, I.; Jamet, E.; Desmonts, M.H.; Coton, E.; Pochet, S.; Donnio, P.Y.; Gautier, M.; Talon, R.; Le Loir, Y. Low occurrence of safety hazards in coagulase negative staphylococci isolated from fermented foodstuffs. Int. J. Food Microbiol. 2010, 139, 87–95. [Google Scholar] [CrossRef]
  52. Baquero, F. Low-level antibacterial resistance: a gateway to clinical resistance. Drug Resist. Updat. 2001, 4, 93–105. [Google Scholar] [CrossRef]
  53. Frimodt-Møller, J.; Rossi, E.; Haagensen, J.A.J.; Falcone, M.; Molin, S.; Johansen, H.K. Mutations causing low level antibiotic resistance ensure bacterial survival in antibiotic-treated hosts. Sci. Rep. 2018, 8, 12512. [Google Scholar] [CrossRef]
  54. Horvath, A.; Rozgonyi, F.; Pesti, N.; Kocsis, E.; Malmos, G.; Kristof, K.; Nagy, K.; Lagler, H.; Presterl, E.; Stich, K.; Gattringer, R.; Kotolacsi, G.; Cekovska, Z.; Graninger, W. Quantitative differences in antibiotic resistance between methicillin-resistant and methicillin-susceptible Staphylococcus aureus strains isolated in Hungary, Austria and Macedonia. J. Chemother. 2010, 22, 246–53. [Google Scholar] [CrossRef]
  55. Tamendjari, S.; Bouzebda, F.A.; Chaib, L.; Aggad, H.; Ramdani, M.; Bouzebda, Z. Antibiotic resistance of Staphylococcus aureus isolated from raw cow and goat milk produced in the Tiaret and Souk Ahras areas of Algeria. Vet. World. 2021, 14, 1929–1934. [Google Scholar] [CrossRef] [PubMed]
  56. Kapoor, G.; Saigal, S.; Elongavan, A. Action and resistance mechanisms of antibiotics: A guide for clinicians. J. Anaesthesiol. Clin. Pharmacol. 2017, 33, 300–305. [Google Scholar] [CrossRef] [PubMed]
  57. European Commission. Guidance on the approval and low-risk criteria linked to “antimicrobial resistance” applicable to microorganisms used for plant protection in accordance with regulation EC no. 1107/2009, 2020.
  58. Olsson-Liljequist, B.; Larsson, P.; Walder, M.; Miorner, H. Antimicrobial susceptibility testing in Sweden. III. Methodology for susceptibility testing. Scand. J. Infect. Dis. Suppl. 1997, 105, 13–23. [Google Scholar] [PubMed]
  59. Kahlmeter, G.; Turnidge, J. How to: ECOFFs-the why, the how, and the donʼts of EUCAST epidemiological cutoff values. Clin. Microbiol. Inf. 2022, 28, 952–954. [Google Scholar] [CrossRef] [PubMed]
  60. Founou, L.L.; Founou, R.C.; Essack, S.Y. Antibiotic resistance in the food chain: A developing country-perspective. Front. Microbiol. 2016, 7, 1881. [Google Scholar] [CrossRef] [PubMed]
  61. Turnidge, J.; Kahlmeter, G.; Kronvall, G. Statistical characterisation of bacterial wild-type MIC value distributions and the determination of epidemiological cut-off values. Clin. Microbiol. Infect. 2006, 12, 418–25. [Google Scholar] [CrossRef]
  62. Kronvall, G. Normalized resistance interpretation as a tool for establishing epidemiological MIC susceptibility breakpoints. J. Clin. Microbiol. 2010, 48, 4445–52. [Google Scholar] [CrossRef]
  63. Campedelli, I.; Mathur, H.; Salvetti, E.; Clarke, S.; Rea, M.C.; Torriani, S.; Ross, R.P.; Hill, C.; OʼToole, P.W. Genus-wide assessment of antibiotic resistance in Lactobacillus spp. Appl. Environ. Microbiol. 2018, 85, e01738–18. [Google Scholar] [CrossRef]
  64. Bortolaia, V.; Kaas, R.S.; Ruppe, E.; Roberts, M.C.; Schwarz, S.; Cattoir, V.; Philippon, A.; Allesoe, R.L.; Rebelo, A.R.; Florensa, A.F.; Fagelhauer, L.; Chakraborty, T.; Neumann, B.; Werner, G.; Bender, J.K.; Stingl, K.; Nguyen, M.; Coppens, J.; Xavier, B.B.; Malhotra-Kumar, S.; Westh, H.; Pinholt, M.; Anjum, M.F.; Duggett, N.A.; Kempf, I.; Nykäsenoja, S.; Olkkola, S.; Wieczorek, K.; Amaro, A.; Clemente, L.; Mossong, J.; Losch, S.; Ragimbeau, C.; Lund, O.; Aarestrup, F.M. ResFinder 4.0 for predictions of phenotypes from genotypes. J. Antimicrob. Chemother. 2020, 75, 3491–3500. [Google Scholar] [CrossRef]
  65. Su, M.; Satola, S.W.; Read, T.D. Genome-based prediction of bacterial antibiotic resistance. J. Clin. Microbiol. 2019, 57, e01405–18. [Google Scholar] [CrossRef]
  66. Schürch, A.C.; van Schaik, W. Challenges and opportunities for whole-genome sequencing-based surveillance of antibiotic resistance. Ann. N.Y. Acad. Sci. 2017, 1388, 108–120. [Google Scholar] [CrossRef] [PubMed]
  67. WHO. GLASS whole genome sequencing for surveillance of antimicrobial resistance. Geneva: World Health Organization. Licence: CC BY-NC-SA- 3.0 IGO, 2020.
  68. Altayb, H.N.; Elbadawi, H.S.; Baothman, O.; Kazmi, I.; Alzahrani, F.A.; Nadeem, M.S.; Hosawi, S.; Chaieb, K. Whole-genome sequence of multidrug-resistant methicillin-resistant Staphylococcus epidermidis carrying biofilm-associated genes and a unique composite of SCCmec. Antibiotics (Basel) 2022, 11, 861. [Google Scholar] [CrossRef] [PubMed]
  69. Pennone, V.; Prieto, M.; Álvarez-Ordóñez, A.; Cobo-Diaz, J.F. Antimicrobial resistance genes analysis of publicly available Staphylococcus aureus genomes. Antibiotics (Basel) 2022, 11, 1632. [Google Scholar] [CrossRef] [PubMed]
  70. Witteveen, S.; Hendrickx, A.P.A.; de Haan, A.; Notermans, D.W.; Landman, F.; van Santen-Verheuvel, M.G.; de Greeff, S.C.; Kuijper, E.J.; van Maarseveen, N.M.; Vainio, S.; Schouls, L.M. Genetic characteristics of methicillin-resistant Staphylococcus argenteus isolates collected in the Dutch national MRSA surveillance from 2008 to 2021. Microbiol, Spectr. 2022, 10, e0103522. [Google Scholar] [CrossRef]
  71. Jeong, D.-W.; Heo, S.; Ryu, S.; Blom, J.; Lee, J.-H. Genomic insights into the virulence and salt tolerance of Staphylococcus equorum. Sci. Rep. 2018, 7, 5383. [Google Scholar] [CrossRef] [PubMed]
  72. Wagner, T.M.; Howden, B.P.; Sundsfjord, A.; Hegstad, K. Transiently silent acquired antimicrobial resistance: an emerging challenge in susceptibility testing. J. Antimicrob. Chemother. 2023, 78, 586–598. [Google Scholar] [CrossRef]
  73. Kime, L.; Randall, C.P.; Banda, F.I.; Coll, F.; Wright, J.; Richardson, J.; Empel, J.; Parkhill, J.; OʼNeill, A.J. Transient silencing of antibiotic resistance by mutation represents a significant potential source of unanticipated therapeutic failure. mBio 2019, 10, e01755–19. [Google Scholar] [CrossRef]
  74. Schwendener, S.; Perreten, V. The bla and mec families of β-lactam resistance genes in the genera Macrococcus, Mammaliicoccus and Staphylococcus: an in-depth analysis with emphasis on Macrococcus. J. Antimicrob. Chemother. 2022, 77, 1796–1827. [Google Scholar] [CrossRef]
  75. Flórez, A.B.; Vázquez, L.; Rodríguez, J.; Mayo, B. Directed recovery and molecular characterization of antibiotic resistance plasmids from cheese bacteria. Int. J. Mol. Sci. 2021, 22, 7801. [Google Scholar] [CrossRef]
  76. Rocha, G.D.; Nogueira, J.F.; Gomes Dos Santos, M.V.; Boaventura, J.A.; Nunes Soares, R.A.; José de Simoni Gouveia, J.; Matiuzzi da Costa, M.; Gouveia, G.V. Impact of polymorphisms in blaZ, blaR1 and blaI genes and their relationship with β-lactam resistance in S. aureus strains isolated from bovine mastitis. Microb. Pathog. 2022, 165, 105453. [Google Scholar] [CrossRef]
  77. Ivanovic, I.; Boss, R.; Romanò, A.; Guédon, E.; Le-Loir, Y.; Luini, M.; Graber, H.U. Penicillin resistance in bovine Staphylococcus aureus: Genomic evaluation of the discrepancy between phenotypic and molecular test methods. J. Dairy Sci. 2023, 106, 462–475. [Google Scholar] [CrossRef]
  78. Lüthje, P.; Schwarz, S. Molecular basis of resistance to macrolides and lincosamides among staphylococci and streptococci from various animal sources collected in the resistance monitoring program BfT-GermVet. Int. J. Antimicrob. Agents. 2007, 29, 528–35. [Google Scholar] [CrossRef] [PubMed]
  79. Bukowski, M.; Piwowarczyk, R.; Madry, A.; Zagorski-Przybylo, R.; Hydzik, M.; Wladyka, B. Prevalence of antibiotic and heavy metal resistance determinants and virulence-related genetic elements in plasmids of Staphylococcus aureus. Front. Microbiol. 2019, 10, 805. [Google Scholar] [CrossRef] [PubMed]
  80. Cardoso, M.; Schwarz, S. Chloramphenicol resistance plasmids in Staphylococcus aureus isolated from bovine subclinical mastitis. Vet. Microbiol. 1992, 30, 223–32. [Google Scholar] [CrossRef] [PubMed]
  81. Tennent, J.M.; May, J.W.; Skurray, R.A. Characterisation of chloramphenicol resistance plasmids of Staphylococcus aureus and S. epidermidis by restriction enzyme mapping techniques. J. Med. Microbiol. 1986, 22, 79–84. [Google Scholar] [CrossRef] [PubMed]
  82. te Riele, H.; Michel, B.; Ehrlich, S.D. Single-stranded plasmid DNA in Bacillus subtilis and Staphylococcus aureus. Proc. Natl. Acad. Sci. USA 1986, 83, 2541–2545. [Google Scholar] [CrossRef]
  83. Sohail, M.; Dyke, K.G. Suppression of the thermosensitive replication phenotype of the derivative plasmid of pI9789::Tn552 in Staphylococcus aureus may involve integration of the plasmid into the host chromosome. FEMS Microbiol. Lett. 1996, 136, 129–136. [Google Scholar] [CrossRef]
  84. 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]
  85. Yoon, S.H.; Ha, S.M.; Lim, J.M.; Kwon, S.J.; Chun, J. A large-scale evaluation of algorithms to calculate average nucleotide identity. Antonie van Leeuwenhoek 2017, 110, 1281–1286. [Google Scholar] [CrossRef]
  86. Oʼsullivan, D.J.; Klaenhammer, T.R. Rapid mini-prep isolation of high-quality plasmid DNA from Lactococcus and Lactobacillus spp. Appl. Environ. Microbiol. 1993, 59, 2730–2733. [Google Scholar] [CrossRef]
  87. Sambrook, J.; Russell, D.W. Molecular cloning: a laboratory manual, 3rd ed.; Cold Spring Harbor Laboratory Press: New York, USA, 2001. [Google Scholar]
  88. Schneewind, O.; Missiakas, D. Genetic manipulation of Staphylococcus aureus. Curr. Protoc. Microbiol. 2014, 32, Unit–9C.3. [Google Scholar] [CrossRef] [PubMed]
  89. Holo, H.; Nes, I.F. High-frequency transformation by electroporation of Lactococcus lactis subsp. cremoris grown with glycine in osmotically stabilized media. Appl. Environ. Microbiol. 1989, 55, 3119–3123. [Google Scholar] [CrossRef] [PubMed]
  90. Pérez, M.; Calles-Enríquez, M.; Nes, I.; Martín, M.C.; Férnandez, M.; Ladero, V.; Álvarez, M.A. Tyramine biosynthesis is transcriptionally induced at low pH and improves the fitness of Enterococcus faecalis in acidic environments. Appl. Microbiol. Biotechnol. 2015, 99, 3547–58. [Google Scholar] [CrossRef] [PubMed]
Table 1. - Distribution of the Minimum Inhibitory Concentration (MIC) values of 16 antibiotics to 30 Staphylococcus equorum strains isolated from cheese obtained with the Sensititre EULACBI1 and EULACBI2 plates (Trek Diagnostic Systems).
Table 1. - Distribution of the Minimum Inhibitory Concentration (MIC) values of 16 antibiotics to 30 Staphylococcus equorum strains isolated from cheese obtained with the Sensititre EULACBI1 and EULACBI2 plates (Trek Diagnostic Systems).
Antibiotics Number of isolates with a MIC value (µg mL-1) Staphylococcus spp. cut-offsa S. equorum cut-offs
0.03 0.06 0.12 0.25 0.5 1 2 4 8 16 32 64 128 256 EUCAST CLSI This workb
S (≤) R (>) S (≤) I (=) R (≥) R (≥)
Gentamicin 30c 2 2 4 8 16 1
Kanamycin 30c 8 8 (-) 4
Streptomycin 12c 16 2 (-) (-) 4
Neomycin 30c (-) (-) 0.25
Tetracycline 14 14 2 1 1 4 8 16 2
Erythromycin 3 14 6 3 1 2d 1d 1 1 0.5 1-4 8 2
Clindamycin 2 8 7 7 5 1 0.25 0.25 0.5 1-2 4 4
Chloramphenicol 3 23 3 1 (-) 8 16 32 32
Ampicillin 4c 8 8 3 2d 3d 2d (-) (-) 0.5
Penicillin G 7c 4 12 2d 4d 1d (-) 0.12 (-) 0.25 0.5
Vancomycin 21 9 4 4 4 8-16 32 2
Quinupristin-dalfopristin 2 11 16 1 1 1 1 2 4 4
Linezolid 5 19 6 4 4 4 (-) 8 8
Trimethoprim 2 10 8 10 4 4 8 (-) 16 8
Ciprofloxacin 21c 9 0.001 1 (-) 2 4 1
Rifampicin 27c 3 0.06 0.06 1 2 4 0.5
aCut-offs (in µg mL-1) established by the European Committee on Antimicrobial Susceptibility Testing (EUCAST) [40] and the Clinical and Laboratory Standards Institute (CLSI) [41]. Output: S, susceptible; I, intermediate; R, resistant; (-), breakpoint not established. In grey, strains with an MIC value considered resistant by EUCAST, CLSI or both. bThe cut-offs proposed included the results of this work and those in Marty et al. [37]. cStrains that do not grow on the lowest concentration of the antibiotic assayed (MIC ≤ to this value). dCurrent MIC established with the MIC Test System (MTS; Liofilchem).
Table 2. - Genes related to antibiotic resistance identified in the genome of the Staphylococcus equorum strains isolated from cheese of this study.
Table 2. - Genes related to antibiotic resistance identified in the genome of the Staphylococcus equorum strains isolated from cheese of this study.
Antibiotic class/gene Activity/Resistance mechanism Strain(s) Identified by database and/or pipeline % identity/% length coveragea Amino acid (aa) identity/total aa Locationb (size kbp) Maximum homology to protein
Penams
blaR1-blaZI Class A beta-lactamase/antibiotic inactivation (AI) 5A3I, 11A1I, 30A2I, 48A3I, 50A2C CARD, NCBI-RGC, PATRIC, ResFinder 100/100 281/281 Plasmid (6.5-8.9) WP_069819195.1
bla Class A beta-lactamase CL10P, 1BCExtra, 5A3I, 8A3C,16A1C, 50A2C Manual revision 99-100/100 279-282/282 C WP_002508531.1
T17 99/100 279/282 WP_064783177.1
2A3C, 11A1I, 30A2I 99/100 282/282 WP_069813561.1
23A3C 100/100 282/282 WP_119627547.1
35A3C CARD, NCBI-RGC, ResFinder 100/100 282/282 WP_046465027.1
48A3I Manual revision 99/100 280/282 WP_197911012.1
Macrolides
mph(C) Macrolide 2ʼ-phosphotransferase/AI 8A3C, 16A1C CARD, NCBI-RGC, ResFinder 100/100 299/299 C WP_119544566.1
msr(A) ABC-F type ribosomal protection protein/target protection T17, 2A3C, 23A3C, 35A3C PATRIC, ResFinder 99/100 488/488 C WP_069813611.1
8A3C, 16A1C 99/100 487/488 WP_046465994.1
50A2C 100/100 488/488 WP_069854570.1
Lincosamides
lnu(A) Lincosamide nucleotidyltransferase/AI 1BCExtra, 2A3C Manual revision 100/100 161/161 Plasmid (32.0-34.6) WP_069813868.1
Phenicols
cat Type A chloramphenicol o-acetyl transferase/AI 35A3C CARD, NCBI-RGC, PATRIC, ResFinder 100/100 215/215 Plasmid (4.6) WP_053038759.1
Fluoroquinolones
norA CLP10, 1BCExtra, 48A3I PATRIC 99/100 385/386 C WP_002508336.1
Major facilitator superfamily of efflux pumps/antibiotic secretion T17, 23A3C,35A3C 100/100 386/386 WP_064783100.1
2A3C, 8A3C, 16A1C, 50A2C PATRIC, ResFinder 100/100 386/386 WP_021339414.1
5A3I*, 11A1I, 30A2I PATRIC 99*-100/100 385*-386/386 WP_069832674.1
Phosphonic acids
fosB/fosD Fosfomycin bacillithiol transferase/AI 1BCExtra-1 CARD, NCBI-RGC, PATRIC, ResFinder 100/100 139/139 C WP_000616116.1
T17, 23A3C, 35A3C PATRIC 100/100 139/139 WP_031266123.1
1BCExtra-2* CARD, PATRIC 84/51 70/139 WP_056935383.1
2A3C*,8A3C*, 16A1C*, 50A2C* CARD 84/32 45/139 WP_031266123.1
5A3I, 11A1I, 30A2I CARD, NCBI-RGC, PATRIC, ResFinder 100/100 139/139 WP_069833353.1
48A31* CARD, PATRIC 83/51 70/139 WP_031266123.1
aIdentity and coverage of DNA or protein sequence, depending on the database. Only identity and length coverage percentages of proteins are stated. bC, chromosome; P, plasmid. Plasmid location was considered when genes encoding plasmid-replicating, mobilization and/or maintenance proteins were found in the same contig. The size referrers to the size of the contig; that of the cat gene contained the whole plasmid. *Disrupted genes containing premature stop codons.
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.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

© 2024 MDPI (Basel, Switzerland) unless otherwise stated