Next Article in Journal
Streptococcus pyogenes Lineage ST62/emm87: The International Spread of This Potentially Invasive Lineage
Previous Article in Journal
Occurrence, Typing, and Resistance Genes of ESBL/AmpC-Producing Enterobacterales in Fresh Vegetables Purchased in Central Israel
Previous Article in Special Issue
Emergence and Genomic Features of a mcr-1 Escherichia coli from Duck in Hungary
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genomic Insights into Methicillin-Resistant Staphylococci and Mammaliicocci from Bulk Tank Milk of Dairy Farms in Serbia

1
Directorate for National Reference Laboratories, Ministry of Agriculture, Forestry and Water Management, Batajnički drum 7, 11186 Belgrade, Serbia
2
Institute of Microbiology and Parasitology, Veterinary Faculty, University of Ljubljana, Gerbičeva 60, 1000 Ljubljana, Slovenia
3
Department of Food Hygiene and Technology, Faculty of Veterinary Medicine, University of Belgrade, Bulevar Oslobođenja 18, 11000 Belgrade, Serbia
*
Author to whom correspondence should be addressed.
Antibiotics 2023, 12(10), 1529; https://doi.org/10.3390/antibiotics12101529
Submission received: 15 September 2023 / Revised: 29 September 2023 / Accepted: 7 October 2023 / Published: 11 October 2023

Abstract

:
The potential risk to human and animal health provides a rationale for research on methicillin-resistant staphylococci (MRS) and mammaliicocci (MRM) in dairy herds. Here, we aimed to estimate their occurrence in the bulk tank milk (BTM) samples collected in 2019–2021 from 283 bovine dairy farms in the Belgrade district. We used whole-genome sequencing to characterize the obtained isolates and assess their genetic relatedness. A total of 70 MRS/MRM were recovered, most frequently Staphylococcus haemolyticus and Mammaliicoccus sciuri. Five clusters of 2–4 genetically related isolates were identified and epidemiological data indicated transmission through, e.g., farm visits by personnel or milk collection trucks. Most MRSA isolates belonged to the typical livestock-associated lineage ST398-t034. One MRSA isolate (ST152-t355) harbored the PVL-encoding genes. Since MRS/MRM isolates obtained in this study frequently harbored genes conferring multidrug resistance (MDR), this argues for their role as reservoirs for the spread of antimicrobial resistance genes. The pipeline milking system and total bacterial count >100,000 CFU/mL were significantly associated with higher occurrences of MRS/MRM. Our study confirms that BTM can be a zoonotic source of MRS, including MDR strains. This highlights the urgent need for good agricultural practices and the continuous monitoring of MRS/MRM in dairy farms.

1. Introduction

Staphylococcus aureus is a well-known pathogen of many animal species and humans, and is associated with a variety of manifestations ranging from benign skin infections to life-threatening conditions [1]. In dairy herds, mastitis caused by S. aureus remains a major disease burden despite numerous herd management programs [2]. Although many different S. aureus genotypes circulate in dairy herds worldwide, a limited number of clones are responsible for most S. aureus mastitis cases, with clonal complexes (CCs) 1, 5, 8 and 97 being the most prevalent [2]. From a food safety perspective, S. aureus is considered one of the most important causative agents of food intoxication worldwide [3]. Conversely, coagulase-negative staphylococci, now more often referred to as non-S. aureus staphylococci (NAS) [4], have traditionally been considered non-pathogenic due to their established commensal relationship with humans and animals, although the pathogenicity of different NAS species and their involvement in bovine mastitis remain to be fully evaluated [4].
Methicillin-resistant staphylococci (MRS) have been repeatedly isolated from dairy herds, albeit at low prevalence rates [5]. A meta-analysis of published data showed that the estimated global prevalence of methicillin-resistant S. aureus (MRSA) from various sources on dairy cattle farms was 3.2% [6]. A slightly higher global prevalence of MRSA (4.3%) was reported for clinical and subclinical bovine mastitis cases [7]. Methicillin-resistant NAS, which are often isolated from dairy farms, serve as a potential reservoir of antimicrobial resistance (AMR) and virulence genes for mastitis-associated pathogens [8,9,10]. Since AMR genes are often located on mobile genetic elements, they can be horizontally transferred between different species of the family Staphylococcaceae, and the resistant mammaliicocci can also play a major role as AMR reservoirs for pathogenic staphylococci [4,11,12]. Namely, a novel genus Mammaliicoccus with five species (M. sciuri, M. fleurettii, M. lentus, M. vitulinus and M. stepanovicii) was recently reclassified from the genus Staphylococcus to comprise members of the former Staphylococcus sciuri group [13], and was reported to harbor AMR genes encoding resistance to several antimicrobial classes [12,14] and unusual SCCmec elements [12,15,16]. The presence of mammaliicocci on dairy farms has also been confirmed [10,12].
Agriculture has a long tradition in Serbia and accounts for 6.3% of the gross domestic product [17]. According to the official statistics for 2021, there are a total of 408,000 dairy cows in Serbia [17]. Milk production reaches 1473 × 106 L/year, with an average milk yield of 3626 L/cow, and is heavily dependent on small, family-run farms. Despite the great importance of milk production in Serbia and the existence of a national mastitis control program under the national Regulation on Animal Health Program [18], its implementation in practice has not been successful due to lack of resources and defined protocols. In addition, reliable data on the farm use of antimicrobials are lacking [19]. Moreover, data on the prevalence of major mastitis-associated pathogens, average somatic cell counts (SCC) and total bacterial counts (TBC) in the bulk tank milk (BTM), are generally not available either at national or regional level.
Studies on MRSA and other staphylococci from Serbia are scarce and mainly focused on the prevalence and characterization of MRSA from humans [20,21,22,23], and to a lesser extent on MRSA from farmed [24,25,26] and companion [27] animals. Since the first reported case of MRSA-associated subclinical bovine mastitis [26], efforts have been made to characterize dairy-related MRSA [28,29], but specific data on the prevalence and molecular epidemiology of MRS in the Serbian dairy supply chain are still missing.
Several studies used BTM samples to determine MRSA prevalence in dairy herds [30,31,32,33,34,35]. The objectives of the present study were to (i) estimate the occurrence of MRS and methicillin-resistant mammaliicocci (MRM) in BTM samples, (ii) perform an in-depth genetic characterization of the obtained MRS and MRM isolates using whole-genome sequencing (WGS) and (iii) identify risk factors associated with the occurrence of MRS/MRM in BTM.

2. Results

2.1. The Occurrence of MRS/MRM in BTM Samples

A total of 283 BTM samples were collected from a different dairy farm each to evaluate the occurrence of MRS/MRM. The results showed that 68/283 (24.0%) farms were positive for MRS and/or MRM, with a total of 70 MRS/MRM isolates identified (Tables S1 and S2). Of the 68 farms positive for MRS and/or MRM, 66 were positive for either MRS or MRM. In addition, the farm with code 255 was positive for two different MRS species (S. aureus isolate MRSA4 and S. haemolyticus isolate MRS7), and the farm with code 250 was positive for MRS and MRM species (M. fleurettii isolate MRS36a_Sf and S. haemolyticus isolate MRS36_Sh) (Tables S1 and S2). In total, 41/283 (14.5%) farms were positive for MRS and 28/283 (9.9%) farms for MRM; the percentage of MRSA-positive farms was 3.5% (10/283) (Table S2).

2.2. Species Identification of MRS and MRM

The performance of the VITEK 2 and matrix-assisted laser desorption/ionization–time of flight (MALDI–TOF) VITEK mass spectrometry (MS) systems for MRS/MRM identification was assessed using WGS as the reference method. Species identification with three different methods was discordant for 10/70 isolates: 9 M. fleuretii isolates (misidentified by VITEK 2 and MALDI-TOF VITEK MS systems) and 1 S. aureus isolate (misidentified by VITEK 2 system) (Table S2). The studied 70 MRS/MRM isolates belonged to the following species according to WGS: four Staphylococcus species (S. aureus (n = 10), S. epidermidis (n = 10), S. haemolyticus (n = 17) and S. saprophyticus (n = 5)) and three Mammaliicoccus species (M. fleurettii (n = 9), M. lentus (n = 2) and M. sciuri (n = 17)). All isolates were mecA-positive by PCR.

2.3. WGS Characterization of MRS and MRM

Regarding multilocus sequence typing (MLST), spa typing and SCCmec typing, the 10 studied S. aureus isolates belonged to four different sequence types (STs): ST398 (six isolates; all spa type t034 and SCCmec type V), ST7882 (two isolates; t127 and t693, both SCCmec type IVa), ST152 (one isolate; t355 and SCCmec type V) and ST7 (one isolate; t091 and SCCmec type V). The basic characteristics of the constructed ad hoc (species-specific) cgMLST schemes are shown in Table S3 and the typing results for other MRS and MRM species in Table S4.
In total, 17 virulence genes were detected in MRSA isolates; each MRSA isolate harbored at least 4 of them. One isolate belonging to ST152 harbored the Panton–Valentine leucocidin (PVL)-encoding genes (Figure 1, Table S5). Other MRS and MRM species harbored no virulence genes, except for the arginine catabolic mobile element (ACME) present in 5/10 S. epidermidis isolates (Figure 1, Table S5).
In total, 31 AMR genes were detected in MRS/MRM isolates, conferring resistance to eight antimicrobial groups (aminoglycosides, beta-lactams, macrolide/lincosamide/streptogramin (MLS) group, tetracyclines, trimethoprim, phenicols, fosfomycin and fusidic acid) (Tables S6 and S7). Additionally, in S. aureus isolates (6/10), two chromosomal point mutations conferring fluoroquinolone resistance were identified: S80Y mutation in the grlA/parC gene and S84L mutation in the gyrA gene. AMR genes identified in each MRS/MRM isolate are shown in Figure 1 and Figure 2. The following AMR genes were identified in two M. lentus isolates: (i) mecA, mph(C), lnu(A), erm(B), tet(K), fosD, str, aac(6′)-aph(2″), fexA and dfrG in MRS44, and (ii) mecA, mph(C), aac(6′)-aph(2″), cat(pC221) in MRS50. In total, 61.4% isolates harbored genes conferring multidrug resistance (MDR; resistance to ≥3 antimicrobial classes [36]). The MDR rate was significantly higher in MRS isolates than in MRM isolates (36/42 vs. 7/28, respectively; p < 0.0001, Fisher’s exact test).
Clustering based on the core genome multilocus sequence typing (cgMLST) was performed for all species with more than two representative isolates, and the constructed phylogenetic trees showed high genetic diversity within the species (Figure 1 and Figure 2). The following cgMLST clusters (≤24 cgMLST allele differences) were identified: two S. aureus clusters, two M. fleurettii clusters and one M. sciuri cluster (Table S4, Figure 1 and Figure 2). Clusters 1 and 2 encompassed isolates from the farms in close proximity to each other (same village; ~2 km radius), whereas other clusters encompassed isolates from geographically more distant farms (20–50 km). A single veterinarian supervised both farms associated with cluster 1. Farms associated with clusters 2 and 5 belonged to the milk-shed area of the same dairy (same raw milk supply line).

2.4. Risk Factors for the Occurrence of MRS/MRM in Dairy Farms

The milking system, housing system, herd size, breed, TBC and SCC were included into the logistic regression model to assess their association with the occurrence of MRS/MRM in dairy farms. The milking system and TBC were significantly associated with the occurrence of MRS/MRM in dairy farms, as assessed by the univariate risk factor analysis (Table 1), and were included in the final multivariate logistic regression model (Table 2). The Hosmer–Lemeshow test showed that the regression model was a good fit (p = 0.056). The regression model showed that farms with a pipeline milking system had 2.5-fold higher odds of harboring MRS/MRM than farms with a bucket milking system (p = 0.003). In addition, farms with TBC > 100,000 CFU (colony forming units)/mL had 2.8-fold higher odds of harboring MRS/MRM compared with farms with TBC < 100,000 CFU/mL (p = 0.001).

3. Discussion

Dairy production is the largest sector of Serbian agriculture, accounting for approximately 8% of the agriculture production value [37]. It depends on family-run farms and lacks national mastitis control programs. To this aim, we have estimated the occurrence of MRS and MRM in BTM samples from 283 bovine dairy farms located in the Belgrade district. The selected farms were representative of the Serbian dairy sector, since the vast majority (99.95%) of national agricultural holdings (AHs) specialized in milk production are family-run [38]. In addition, the analyzed farms were characterized by the tie-stall housing system and domestic Simmental dairy breed, which are also prevalent at the national level [39,40].
One quarter of the analyzed farms were positive for MRS and/or MRM, with 3.5% of farms positive for MRSA. To date, several studies have investigated MRSA prevalence in dairy farms [10,30,32,33,34,35,41,42,43,44,45,46,47,48]. In general, the prevalence of MRSA in dairy herds has been reported to be low [5,6], although an increase is expected [48]. However, differences in study designs (sample types, inoculum sizes, isolation procedures, detection methods and farm management systems) make comparisons difficult, in addition to the impact of other factors (e.g., geographic location, sanitation practices, on-farm biosecurity, antimicrobial use and mastitis control strategy) that may also influence MRSA prevalence. The low MRSA occurrence rate observed in this study is consistent with data from similar studies that used BTM cultures and a selective enrichment protocol [30,32,33,44]. Higher MRSA occurrence rates in BTM reported in other studies could be explained by pre-selection of dairy farms based on a history of MRSA-associated clinical/subclinical mastitis [42] or previous reports of MRSA from BTM [49,50].
Most previous studies on staphylococci in dairy herds have focused on MRSA, although methicillin resistance has been reported in several other Staphylococcus and Mammaliicoccus species commonly found on dairy farms [10,51,52,53,54]. Methicillin resistance rates in NAS are generally higher than in S. aureus [5]. In this study, non-S. aureus MRS/MRM were detected in 20.8% of BTM samples. This rate is lower than those previously reported in Switzerland (62.0%; [51]) and Germany (42.1%; [10]). Of note, in the German study, farms were selected based on previous MRSA detection. In contrast, a low incidence of non-S. aureus MRS/MRM in BTM has been reported in the United States (2.4%; [53]) and the United Kingdom (~5%; [54]). The low methicillin resistance rate in NAS and mammaliicocci noted in the above-mentioned studies could be due to the controlled and restrictive use of antimicrobials in dairy farming, high levels of hygiene on farms (especially during milking) and implementation of biosecurity measures [53,54]. M. sciuri and S. haemolyticus were the most frequently isolated species in this study. M. sciuri is widespread in nature and can be found in various animal hosts [14]. It also represents the most likely evolutionary origin of mecA [55]. M. sciuri is most commonly isolated from animal teats, whereas S. haemolyticus is mainly associated with milk samples [56,57] and it has a higher prevalence of the mecA gene compared to other staphylococci species [58]. Previous studies on the occurrence of MRS/MRM in cattle also reported a high rate of methicillin resistance in M. sciuri [8,50]. M. sciuri has also been reported as the most abundant MRS/MRM species from BTM in England and Wales [54].
In this study, TBC (>100,000 CFU/mL) and the milking system (pipeline milking) were significantly associated with the occurrence of MRS/MRM in dairy farms. Elevated TBC is an indicator of unsanitary practices on the farm, including an absence of precautionary measures for milking operation, a low milking cow cleanliness score, poor equipment hygiene and water quality, as well as high California Mastitis Test (CMT) scores and inadequate raw milk cooling [59]. Effective milking preparation and housing hygiene are paramount in preventing contagious and environmental mastitis [60]. S. aureus spreads from cow to cow and udder to udder, primarily during the milking process. The milking hygiene score is negatively correlated with MRSA prevalence [49]. Inadequate hygiene during milking has been observed previously in all MRSA-affected farms [50]. Biofilms in the milking equipment are also a potential source of contamination of BTM with S. aureus and other pathogens [61,62]. Although hygiene during milking and the presence of biofilms in the milking systems were not evaluated in this study, we assume that the lack of pre- and post-dipping protocols, improper handling and hygiene of milking unit clusters, lack of mastitis control measures such as rapid segregation and culling of infected cows are positively correlated with the MRS/MRM-positive status of the dairy farms investigated. On the contrary, the housing system, herd size, breed and SCC were not found to be significant risk factors for the MRS/MRM-positive status of dairy herds. This contrasts with previous studies that identified herd size as a risk factor for MRS/MRM occurrence, with larger herds more likely to be positive than smaller herds [32,34]. Of note, most (76.3%) of the analyzed farms were typically in the range of 1–19 cows and often kept multiple animal species including pigs, which are considered a risk factor for MRSA occurrence in dairy herds [5]. A previous study on udder pathogens identified milking and housing systems as important risk factors for the occurrence of intramammary infections and showed that these risk factors are pathogen specific [63].
MALDI-TOF MS has been successfully used to identify staphylococci from bovine milk [64,65,66]. In this study, a comparative evaluation of the VITEK 2 and MALDI-TOF VITEK MS systems in parallel with WGS as a reference method showed that nine out of ten misidentified isolates belonged to M. fleurettii. This is likely due to limitations of the reference databases used (VITEK 2 GP and SARAMIS), which did not include M. fleurettii. As previously reported [67,68], curated and updated databases are essential for optimal performance of VITEK MS automated systems. An advantage of WGS as the reference method for species identification compared with other methods is that it provides the most accurate phylogeny and the analyzed strains can be compared with any strain/species in the database, including recently described species. An additional advantage of WGS used here was the determination of AMR and virulence genes, which eliminated the need for conventional phenotypic and PCR testing.
In this study, five cgMLST clusters of MRS/MRM isolates were observed. Clusters 1 and 2 included isolates from farms located in the same village, suggesting MRSA transmission through personnel visits (e.g., farm workers, veterinarians or professional visitors). Other possible transmission vehicles include animal vectors (e.g., flies or rodents) or dust [50]. Because the farms associated with clusters 2 and 5 were part of the milk-shed area of the same dairy, milk trucks represent a likely transmission route. These trucks routinely visit multiple dairy herds on the same day, and no disinfection baths were present at the entrance to each farm. One MRSA-positive farm associated with cluster 2 (farm code 224) kept pigs in close proximity (<50 m) to dairy cows. Thus, pigs could also have been involved in MRSA transmission to dairy cows. In all cases, MRS/MRM transmission could also have occurred through local livestock trade, but such data were not available. Additional epidemiological data on the farms studied would be needed to identify the role of the above possible transmission routes in MRS/MRM transmission between herds.
Several MRSA clones have been identified in milk and dairy products worldwide, with CC398 and CC1 being the most frequently reported MLST CCs in Europe [69]. Most (6/10) MRSA isolates analyzed in this study belonged to the livestock-associated (LA) lineage ST398-SCCmec V-t034. Previous studies from Europe reported ST398 as the predominant clone in dairy herds [32,34,70,71,72,73]. Two MRSA isolates from our study were of CC1, one belonging to spa type t127. MRSA CC1 has been recognized as a successful hospital- (HA) and community-acquired (CA) MRSA lineage in humans and has also been responsible for invasive infections [74,75]. MRSA CC1-t127 has also been reported in cattle from Germany, Italy and Switzerland [34,74,76,77]. In Hungary, MRSA CC1-t127 was found in cattle with subclinical mastitis, and transmission to a farm worker in close contact with cattle was also demonstrated [78]. In this study, one of the identified MRSA isolates was of genotype ST152-t355 and harbored the PVL-encoding genes. This genotype has been linked to community cases throughout Europe and is prevalent in the Balkan region [79]. This genotype was the most common (67.5%) MRSA genotype isolated from milk and dairy products in southern Italy, and most isolates were PVL-positive [80]. It is worth noting that the ST152-t355 genotype has been reported in Serbia since 2008 [20] and it is still the dominant PVL-positive human genotype in this region [23]. Due to the zoonotic potential of MRSA [73,81], the introduction of the human-adapted MRSA genotype ST152-t355 into the analyzed farm by colonized humans cannot be excluded.
Virulence genes were distributed with different frequencies among the MRSA isolates studied, which may indicate evolution and adaptation to different niches or environments. In agreement with previous results [82], the most frequently detected virulence genes were hlgA and hlgB, followed by aur and hlgC. This is not surprising since the aur gene is highly conserved in S. aureus [83]. Most (8/10) MRSA isolates examined had the same pattern of virulence genes encoding hemolysins (hlgA, hlgB and hlgC) and aureolysin (aur) as previously described in two MRSA ST398 from BTM originating from healthy cows in China [84]. Additionally, the human immune evasion cluster (IEC) genes scn and sak were detected in three MRSA isolates of genotype ST398-t034. The re-acquisition of IEC by LA-MRSA ST398 suggests the re-adaptation of LA-MRSA ST398 to the human host and favors the (re)emergence of MRSA ST398 in humans [85,86]. Consistent with previous findings [74,87,88,89,90], both MRSA CC1 isolates examined harbored seh, sak and scn. Because staphylococcal enterotoxin H (SEH)-producing strains have been involved in staphylococcal food poisoning [91], the risk of such strains to cause foodborne intoxication needs to be evaluated. The presence of IEC in the MRSA CC1 isolates studied suggests that this clone may be human adapted. In addition, the CC1 isolates examined carried other protease genes (splA, splB) and lukE/D genes encoding leukotoxins. In addition, the two sea-positive MRSA ST398-t034 isolates from BTM samples could also pose a serious threat to human health. Among the classical staphylococcal enterotoxins (SE), staphylococcal enterotoxin A (SEA) is most commonly associated with staphylococcal food poisoning worldwide [92,93]. In contrast, staphylococcal enterotoxin C (SEC) is the most common staphylococcal enterotoxin produced by S. aureus from dairy animals [94,95], and is particularly common in strains from bovine mastitis [96], but was not detected in this study. Although the involvement of MRSA ST398 in food poisoning has not been previously reported [97], the evolving nature of MRSA ST398 and its ability to re-acquire IEC variants containing the sea gene [98], as observed in this study, suggest that special attention needs to be paid to the early detection of such strains.
A comprehensive WGS-based analysis of 441 bovine NAS isolates revealed a diverse repertoire of virulence genes, but a clear link between them and mastitis was not established [99]. In a recent study, the virulence profile of methicillin-resistant NAS and mammaliicocci (MRNASM) was investigated using WGS in the staphylococcal population of the nasal microbiota of dromedary camels from Algeria [100]. In agreement with the present results, among the MRNASM isolates, only S. epidermidis harbored virulence factors. In this study, half of the S. epidermidis isolates harbored ACME, which is involved in an increased ability to colonize skin and mucous membranes [101]. Consistent with this finding, nearly half of the S. epidermidis isolates examined from bovine mastitis milk in China carried ACME [102].
The most recent surveillance report [103] shows that 71,440 kg of antimicrobials were sold for veterinary use in Serbia in 2020, of which 218 kg were used for intramammary treatment. However, comprehensive national programs for AMR surveillance in veterinary medicine are still lacking. According to the survey conducted among Serbian livestock veterinarians, enrofloxacin was the most frequently prescribed antibiotic in mastitis therapy, followed by amoxicillin (alone or in combination with clavulanic acid), penicillin, ceftriaxone and tetracycline [19]. Although a linear relationship between antibiotic use and resistance patterns or the presence of AMR genes has not been demonstrated [104], it can be hypothesized that the frequent prescription of enrofloxacin in mastitis therapy favors the emergence of fluoroquinolone resistance, which was observed in the MRSA ST398 isolates studied. In addition to resistance to fluoroquinolones and beta-lactams, all MRSA ST398-t034 isolates examined harbored genes conferring resistance to additional antimicrobial classes (lincosamides, aminoglycosides, tetracyclines and trimethoprim). MLS resistance encoded by erm(C) was observed in one isolate. This is of concern since MRSA spa type t034 has been shown to be involved in clinical mastitis [105]. CA-MRSA clones generally have lower MDR rates and are susceptible to various non-beta-lactam antimicrobials [106,107], which is consistent with the present results where the CA-MRSA ST152-t355 isolate harbored only the aac(6′)-aph(2″) gene conferring resistance to aminoglycosides (besides the genes conferring resistance to beta-lactams).
Nobrega and colleagues [58] demonstrated that the prevalence of AMR genes in NAS from dairy herds (clinical and subclinical mastitis) was species-specific. Few studies have examined AMR genes in NAS from the milk of healthy cows [54,108]. Fosfomycin has gained attention as a combination therapy for MRSA infections. In this study, fosB was observed only in methicillin-resistant NAS isolates, with the highest occurrence rate observed in S. epidermidis (10/10). A similar occurrence of fosB has been reported for NAS isolated from bovine milk samples and from retail ground meat in Japan [109,110]. The MLS resistance gene erm(C) was identified in at least one representative of the MRS species analyzed, which contrasts with the results of a previous study in which erm(C) was detected exclusively in NAS species [109] or the erm(C) gene was not observed [108]. In this study, the lincosamide resistance gene Inu(A) was rarely identified in NAS isolates, namely in two S. haemolyticus and three S. saprophyticus isolates. This contrasts with the results of a previous study where Inu(A) was frequently identified in NAS [109]. In S. haemolyticus, high rates of MLS resistance genes msr(A) and mph(C) were observed, followed by erm(C), whereas Inu(A), vga(A) and vga(A)LC were less abundant. It should be noted that mph(C) confers high-level resistance to macrolides in S. aureus only when msr(A) is also present [111]. Contrary to the present findings, the msr(A) gene was not identified in S. haemolyticus from milk of dairy cows in China, although mph(C) was the most prevalent AMR gene [108]. Based on WGS-predicted susceptibility, most (36/42) MRS isolates analyzed were classified as MDR, supporting the hypothesis that MRS from cow milk are a potential reservoir for AMR genes [54,108].
The number of different AMR genes varied among Mammaliicoccus species: M. fleurettii isolates mostly harbored only mecA, whereas M. sciuri and M. lentus harbored AMR genes conferring resistance to six and seven antimicrobial groups, respectively. This result is consistent with previous studies [12,112]. Compared with a large number of blaZ-positive MRS isolates, blaZ was present in only one M. sciuri isolate (MRS60). In addition, blaZ was also absent in all S. saprophyticus isolates examined, which is in agreement with previous studies [54,109]. Although the carriage of tetracycline resistance genes is common in M. sciuri [12,112], tet genes were absent in most (12/17) M. sciuri isolates examined. In M. sciuri, three different tet genes were identified. M. sciuri isolates with multiple AMR genes conferring MDR clustered together in the cgMLST tree and were phylogenetically distant from the remaining M. sciuri isolates harboring only mecA, mecA1 and sal(A). Although previously reported in M. lentus from dairy farms [12], the cfr gene, which confers resistance to a last-resort antibiotic linezolid, was not detected in the MRM isolates examined.
Two limitations of the WGS-based prediction of AMR (and other) genes used in this study should be mentioned. First, the AMR gene’s presence does not always equate to phenotypic resistance, and phenotypic susceptibility was not tested in this study. Second, the screening of AMR genes in draft genomes does not necessarily recover all AMR genes present in the genome due to the fragmented nature of the assemblies.

4. Materials and Methods

4.1. Study Design and Sample Collection

During November 2019 and April 2021, BTM samples were collected from 283 bovine dairy farms located in the Belgrade district, covering 3234 km2 surface area, in the scope of quality control scheduled by the national reference laboratory for milk quality and safety. The selection criteria for participation in the study were as follows: farms that had a refrigerated bulk tank for storing and cooling of milk, cows were not milked by hand and the farmers were willing to participate in the study. The farms included in the study represented 1.4% of the total AHs specialized in milk production in Serbia. Most farms (279/283; 98.6%) were classified as family AHs. Four farms (4/283; 1.4%) were AHs with a legal entity status (AHLE). General farm data were collected through a survey that included data on farm location, herd size, breed, housing and milking system, herd average daily milk production, and SCC and TBC in BTM. Farm location data were anonymized in accordance with the general data protection regulation (farm codes 1–283). Metadata of the dairy farms under study are shown in Table S1.
The average number of cows per farm was 21 (range: 2–505). The Simmental breed was kept in 55.5% of the farms studied, whereas the Holstein-Friesian breed was in 29.3% of the farms; both breeds were kept in 15.2% of the farms. Dairy cows were kept exclusively indoors, without grazing. Tie-stall housing was the predominant housing type (92.9%), whereas on the remaining farms (7.1%) cows were kept in a loose (free-stall) system. Four types of milking systems were used: bucket (44.2%), pipeline (48.4%), milking parlor (6.4%) and automatic milking systems (1.1%) (Table S1).
The reported average daily milk yield per farm was 430 L (range: 30–13,064) (Table S1). The average SCC, expressed as a log10 value of somatic cells/mL (mean ± SD), was 5.40 ± 0.28 (range: 4.71–6.26). The overall average TBC, expressed as a log10 value of CFU/mL (mean ± SD), was 4.87 ± 0.46 (range: 3.90–6.45) (Table S1). The analyzed farms delivered raw milk to the two regional dairy plants situated in the Belgrade district. One of them collected milk from more than 14,000 farms, the other from more than 500 farms.
At each farm, a single BTM sample was collected in a sterile 50 mL container directly from the BTM cooler using a clean and sanitized dipper. Samples were immediately refrigerated and transported to the laboratory, where they were processed on the day of sampling.

4.2. Isolation and Confirmation of Presumptive MRS and MRM

The isolation procedure followed a two-step enrichment protocol. For primary enrichment, milk samples (25 mL) were mixed with 225 mL of Mueller–Hinton broth (Oxoid, Wesel, Germany) supplemented with 6.5% NaCl. After incubation at 37 °C for 20 h, 1 mL of the suspension was transferred to 9 mL of tryptic soy broth (Merck, Darmstadt, Germany) supplemented with 3.5 mg/L of cefoxitin (Sigma-Aldrich Chemie, Munich, Germany) and 75 mg/L of aztreonam (Sigma-Aldrich Chemie, Munich, Germany), and incubated at 37 °C for another 20 h.
After the second enrichment, one loopful (10 µL) of suspension was inoculated in parallel onto Baird-Parker (BP) agar plates (Oxoid, Basingstoke, UK) and ChromoBio MRSA chromogenic agar plates selective for MRSA (Biolab, Budapest, Hungary). Suspect MRS/MRM were collected from BP and ChromoBio MRSA agar plates, respectively. A single suspect MRS/MRM colony from each positive BTM sample was selected for further examination, except in two cases where two presumptive MRS/MRM isolates (after examining cultures for purity with subculturing on sheep blood agar plates), were recovered (farm codes 225 and 250; Table S1) and included into the analysis. All suspect MRS/MRM colonies underwent Gram staining and catalase testing.
All Gram-positive, catalase-positive colonies selected from both isolation media (n = 211) were screened for methicillin susceptibility by the disk diffusion test with 30 μg cefoxitin disks on Mueller–Hinton agar plates (both from Oxoid, Wesel, Germany) according to the European Committee on Antimicrobial Susceptibility Testing (EUCAST) guidelines [113]. S. aureus ATCC 29213 was used as a quality control strain.
For molecular confirmation of methicillin resistance, DNA was extracted from bacterial colonies by the boiling method as described previously [114]. The presence of mecA and mecC genes was tested for by a multiplex PCR according to the instructions provided by the EU reference laboratory for antimicrobial resistance [115].

4.3. Species Identification of MRS and MRM Using VITEK 2 and MALDI-TOF VITEK MS Systems

Species identification by VITEK 2 (bioMérieux, Marcy l’Etoile, France) was performed using GP ID cards according to the manufacturer’s instructions. The obtained results were interpreted according to the ID-GP database. S. aureus ATCC 29213 and S. epidermidis ATCC 14990 were used as quality control strains.
Species identification was also performed by MALDI-TOF MS analysis with a VITEK MS system (bioMérieux, Marcy l’Etoile, France). Before the analysis, isolates were streaked onto tryptic soy agar plates (Merck, Darmstadt, Germany) and incubated at 37 °C overnight. One colony was smeared onto target slides covered with 1 µL of VITEK MS matrix and air-dried. The acquired spectra were analyzed using the VITEK MS Plus SARAMIS Knowledge Base v4.12 database.

4.4. WGS and Bioinformatics Analyses

A total of 70 MRS/MRM isolates underwent WGS (Table S2). Genomic DNA was extracted using a DNA Blood & Tissue Kit (Qiagen, Hilden, Germany) according to the manufacturer’s instructions for Gram-positive bacteria. DNA libraries were prepared with the Illumina TruSeq DNA Nano Library Prep Kit (Illumina, San Diego, CA, USA). Sequencing was performed on the NextSeq 500 System using the 2 × 150 bp chemistry (Illumina, San Diego, CA, USA) to a minimum coverage of 150×. Sequencing data were submitted to the NCBI Sequence Read Archive (SRA) database under the BioProject accession number PRJNA1016076.
Raw reads were assembled with Shovill v1.18 (https://github.com/tseemann/shovill, accessed on 5 July 2023) using the --trim option and SPAdes v3.13.1 [116] as the underlying assembler with default parameters. Quality of the assemblies was assessed using Quast v5.0.2 [117]; only the genomes with N50 > 20,000 bp, number of contigs (longer than 1000 bp) less than 500 and total assembly length of ~2.8 Mbp were included into further analyses.
WGS-based species identification based on read mapping was performed using Kraken2 v2.1.2 [118], which was run with default parameters and the MiniKraken database v2. In cases of unreliable species identification using Kraken2, additional species identification based on average nucleotide identity (ANI) values was performed. To this aim, the ANI values based on the MUMmer algorithm were calculated using JSpeciesWS [119], and the 95–96% ANI threshold [120] was used for species delineation.
All MRS/MRM isolates were also screened for the presence of virulence and AMR genes. AMR genes were identified using ResFinder v4.1 [121] with default parameters. Virulence genes were identified using VirulenceFinder v2.0 (https://cge.food.dtu.dk/services/VirulenceFinder/, accessed on 11 July 2023) with default parameters and the S. aureus database. In addition, non-S. aureus draft genomes were screened for virulence genes against the Virulence factor database (VFDB; [122]) using ABRicate v1.0.1 (https://github.com/tseemann/abricate, accessed on 11 July 2023) with 60% coverage and 90% identity cut-offs. spa typing of S. aureus was performed using spatyper v1.0 [123]. MLST of the species with available PubMLST databases (S. aureus, S. epidermidis and S. haemolyticus) was performed using the Sequence query tool implemented in PubMLST [124]. SCCmec typing of all MRS/MRM isolates was performed using SCCmecFinder v1.2 [125] with default parameters. cgMLST was performed using chewBBACA v2.8.5 [126]; to maximize the discriminatory power, an ad hoc cgMLST scheme was constructed for each species with ≥2 representatives and core loci were defined as loci present in 100% genomes under study (Table S3). A neighbor-joining tree based on cgMLST allele profiles was constructed using Grapetree v1.5 [127] and was further annotated using iTol v6.7.3 [128]. A threshold of 24 cgMLST allele differences was used for cluster delineation [129].

4.5. Statistical Analysis

To determine the risk factors significantly associated with the occurrence of MRS/MRM in dairy farms, statistical analysis was performed using SPSS software v21 (SPSS Inc., Chicago, IL, USA). TBC and SCC thresholds were selected according to criteria set up in Regulation No 853/2004 [130]. An arbitrary value was used for herd size classification, taking into account that 87% of national AHs specialized in milk production have fewer than 10 cows [131]. For univariate logistic analysis, milking system, housing system, herd size, breed, TBC and SCC were included as independent categorical variables to identify the risk factors significantly associated with the occurrence of MRS/MRM in dairy farms (binary dependent variable). Subsequently, the variables with a significance level of p ≤ 0.15 in the univariate model were included in the final multivariate logistic regression model to determine the combined (adjusted) effects of the risk factors. The goodness-of-fit of the logistic regression model was assessed using the Hosmer–Lemeshow test. Results were reported as odds ratios (ORs) and their 95% confidence intervals (Cls).

5. Conclusions

This study provides insights into the occurrence and genomic characteristics of MRS/MRM from dairy farms in Serbia and highlights the importance of WGS to study such zoonotic opportunistic pathogens. In total, 24.0% farms were positive for MRS and/or MRM, with 70 MRS/MRM isolates identified. A high within-species genetic diversity was observed. Non-S. aureus MRS/MRM isolates harboring genes conferring multidrug resistance were frequently observed, highlighting their role as potential reservoirs of AMR genes. Virulence genes were rarely detected in MRS/MRM. Most MRSA isolates belonged to the typical LA-MRSA lineage ST398-t034, and one MRSA isolate of genotype ST152-t355 harbored PVL-encoding genes. Despite the low prevalence of MRSA, the presence of IEC-, SE-encoding and AMR genes in MRSA isolates from BTM poses a risk to animal and human health. Considering the zoonotic potential of MRSA, the continuous surveillance of MRSA in agriculture settings is essential. The pipeline milking system and TBC > 100,000 CFU/mL were identified as risk factors for the occurrence of MRS/MRM in BTM. These findings support the implementation of new and targeted strategies for MRS/MRM control and prevention.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/antibiotics12101529/s1, Table S1: Farm metadata; Table S2: Isolate metadata; Table S3: cgMLST schemes; Table S4: Strain characteristics; Table S5: Presence/absence profiles of virulence genes; Table S6: Presence/absence profiles of antimicrobial resistance genes and resistance-associated mutations (staphylococci); Table S7: Presence/absence profiles of antimicrobial resistance genes (mammaliicocci).

Author Contributions

Conceptualization, A.K. and S.B.; investigation, A.K., T.L., J.A., M.G., D.K., J.Đ. and B.P.; resources, A.K., M.G. and S.B.; data curation, A.K., J.A., M.G., D.K., B.P. and S.B.; writing—original draft preparation, J.A., M.G., D.K., B.P. and S.B.; writing—review and editing, J.A., M.G., D.K., B.P. and S.B.; visualization, B.P.; supervision, S.B.; funding acquisition, A.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was financially supported by the Ministry of Science, Technological Development and Innovation of the Republic of Serbia (Contract number 451-03-47/2023-01/200143), by the Directorate for National Reference Laboratories, Belgrade, Serbia and by the Slovenian Research and Innovation Agency (research core funding no. P4-0092 ‘Animal health, environment and food safety’ and Slovenia-Serbia bilateral project grant no. BI-RS/23-25-007 ‘Raw milk as a potential source of multidrug resistant bacteria and resistance genes’).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Raw sequencing data are available under the NCBI BioProject accession number PRJNA1016076.

Acknowledgments

Aleksandra Popović and Sreten Urošević are warmly acknowledged for supporting the identification of study farms. Dara Jovanović is warmly acknowledged for performing VITEK MS analysis.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Howden, B.P.; Giulieri, S.G.; Wong Fok Lung, T.; Baines, S.L.; Sharkey, L.K.; Lee, J.Y.H.; Hachani, A.; Monk, I.R.; Stinear, T.P. Staphylococcus aureus host interactions and adaptation. Nat. Rev. Microbiol. 2023, 21, 380–395. [Google Scholar] [CrossRef]
  2. Campos, B.; Pickering, A.C.; Souza Rocha, L.; Pereira Aguilar, A.; Fabres-Klein, M.H.; de Oliveira Mendes, T.A.; Ross Fitzgerald, J.; de Oliveira Barros Ribon, A. Diversity and pathogenesis of Staphylococcus aureus from bovine mastitis: Current understanding and future perspectives. BMC Vet. Res. 2022, 18, 115. [Google Scholar] [CrossRef]
  3. Fetsch, A.; Johler, S. Staphylococcus aureus as a foodborne pathogen. Curr. Clin. Microbiol. Rep. 2018, 5, 88–96. [Google Scholar] [CrossRef]
  4. De Buck, J.; Ha, V.; Naushad, S.; Nobrega, D.B.; Luby, C.; Middleton, J.R.; De Vliegher, S.; Barkema, H.W. Non-aureus staphylococci and bovine udder health: Current understanding and knowledge gaps. Front. Vet. Sci. 2021, 8, 658031. [Google Scholar] [CrossRef]
  5. Schnitt, A.; Tenhagen, B.A. Risk factors for the occurrence of methicillin-resistant Staphylococcus aureus in dairy herds: An update. Foodborne Pathog. Dis. 2020, 17, 585–596. [Google Scholar] [CrossRef]
  6. Khanal, S.; Boonyayatra, S.; Awaiwanont, N. Prevalence of methicillin-resistant Staphylococcus aureus in dairy farms: A systematic review and meta-analysis. Front. Vet. Sci. 2022, 9, 947154. [Google Scholar] [CrossRef]
  7. Zaatout, N.; Hezil, D. A meta-analysis of the global prevalence of methicillin-resistant Staphylococcus aureus (MRSA) isolated from clinical and subclinical bovine mastitis. J. Appl. Microbiol. 2022, 132, 140–154. [Google Scholar] [CrossRef]
  8. Vanderhaeghen, W.; Vandendriessche, S.; Crombé, F.; Nemeghaire, S.; Dispas, M.; Denis, O.; Hermans, K.; Haesebrouck, F.; Butaye, P. Characterization of methicillin-resistant non-Staphylococcus aureus staphylococci carriage isolates from different bovine populations. J. Antimicrob. Chemother. 2013, 68, 300–307. [Google Scholar] [CrossRef]
  9. Argudin, M.A.; Vanderhaeghen, W.; Butaye, P. Diversity of antimicrobial resistance and virulence genes in methicillin-resistant non-Staphylococcus aureus staphylococci from veal calves. Res. Vet. Sci. 2015, 99, 10–16. [Google Scholar] [CrossRef]
  10. Schnitt, A.; Lienen, T.; Wichmann-Schauer, H.; Tenhagen, B.A. The occurrence of methicillin-resistant non-aureus staphylococci in samples from cows, young stock, and the environment on German dairy farms. J. Dairy Sci. 2021, 104, 4604–4614. [Google Scholar] [CrossRef]
  11. Saraiva, M.M.S.; de Leon, C.; Silva, N.; Raso, T.F.; Serafini, P.P.; Givisiez, P.E.N.; Gebreyes, W.A.; Oliveira, C.J.B. Staphylococcus sciuri as a reservoir of mecA to Staphylococcus aureus in non-migratory seabirds from a remote oceanic island. Microb. Drug Resist. 2021, 27, 553–561. [Google Scholar] [CrossRef] [PubMed]
  12. Lienen, T.; Schnitt, A.; Hammerl, J.A.; Maurischat, S.; Tenhagen, B.-A. Mammaliicoccus spp. from German dairy farms exhibit a wide range of antimicrobial resistance genes and non-wildtype phenotypes to several antibiotic classes. Biology 2022, 11, 152. [Google Scholar] [CrossRef] [PubMed]
  13. Madhaiyan, M.; Wirth, J.S.; Saravanan, V.S. Phylogenomic analyses of the Staphylococcaceae family suggest the reclassification of five species within the genus Staphylococcus as heterotypic synonyms, the promotion of five subspecies to novel species, the taxonomic reassignment of five Staphylococcus species to Mammaliicoccus gen. nov., and the formal assignment of Nosocomiicoccus to the family Staphylococcaceae. Int. J. Syst. Evol. Microbiol. 2020, 70, 5926–5936. [Google Scholar] [CrossRef] [PubMed]
  14. Nemeghaire, S.; Argudín, M.A.; Feßler, A.T.; Hauschild, T.; Schwarz, S.; Butaye, P. The ecological importance of the Staphylococcus sciuri species group as a reservoir for resistance and virulence genes. Vet. Microbiol. 2014, 171, 342–356. [Google Scholar] [CrossRef]
  15. Harrison, E.M.; Paterson, G.K.; Holden, M.T.G.; Ba, X.; Rolo, J.; Morgan, F.J.E.; Pichon, B.; Kearns, A.; Zadoks, R.N.; Peacock, S.J.; et al. A novel hybrid SCCmec-mecC region in Staphylococcus sciuri. J. Antimicrob. Chemother. 2013, 69, 911–918. [Google Scholar] [CrossRef]
  16. Paterson, G.K. Genomic epidemiology of methicillin-resistant Staphylococcus sciuri carrying a SCCmec-mecC hybrid element. Infect. Genet. Evol. 2020, 79, 104148. [Google Scholar] [CrossRef]
  17. SORS (Statistical Office of the Republic of Serbia). Statistical Yearbook of the Republic of Serbia: Agriculture; Publication of SORS: Belgrade, Serbia, 2022; pp. 211–250. Available online: https://publikacije.stat.gov.rs/G2022/PdfE/G20222055.pdf (accessed on 27 February 2023).
  18. OGSR (Official Gazette of the Republic of Serbia). OGRS No. 23/2023: Regulation on Animal Health Program for the Year 2023. in Serbian. Available online: http://demo.paragraf.rs/demo/combined/Old/t/t2023_03/SG_023_2023_008.htm (accessed on 29 September 2023).
  19. Vidović, J.; Stojanović, D.; Cagnardi, P.; Kladar, N.; Horvat, O.; Ćirković, I.; Bijelić, K.; Stojanac, N.; Kovačević, Z. Farm animal veterinarians’ knowledge and attitudes toward antimicrobial resistance and antimicrobial use in the Republic of Serbia. Antibiotics 2022, 11, 64. [Google Scholar] [CrossRef]
  20. Ćirković, I.; Sørum, M.; Radenković, D.; Švabić Vlahović, M.; Larsen, A.R. National surveillance reveals findings of Panton-Valentine leukocidin positive methicillin-resistant Staphylococcus aureus in Serbia. J. Med. Microbiol. 2013, 62, 342–344. [Google Scholar] [CrossRef]
  21. Ćirković, I.; Đukić, S.; Carević, B.; Mazić, N.; Mioljević, V.; Stepanović, S. Methicillin-resistant Staphylococcus aureus nasal carriage among hospitalized patients and health care workers in medical center of Serbia. Arch. Biol. Sci. Belgrade 2014, 66, 87–92. [Google Scholar] [CrossRef]
  22. Cirkovic, I.; Stepanovic, S.; Skov, R.; Trajkovic, J.; Grgurevic, A.; Larsen, A.R. Carriage and genetic diversity of methicillin-resistant Staphylococcus aureus among patients and healthcare workers in a Serbian university hospital. PLoS ONE 2015, 10, e0127347. [Google Scholar] [CrossRef]
  23. Rakonjac, B.; Lepšanović, Z.; Šuljagić, V.; Jovčić, B.; Kojić, M.; Larsen, A.R.; Urić, M.; Ćirković, I. Predominance of t355/ST152/SCCmec V clonal type among PVL-positive MRSA isolates in a tertiary care hospital in Belgrade, Serbia. PLoS ONE 2022, 17, e0273474. [Google Scholar] [CrossRef] [PubMed]
  24. Velebit, B.; Fetsch, A.; Mirilovic, M.; Teodorovic, V.; Jovanovic, M. MRSA in pigs in Serbia. Vet. Rec. 2010, 167, 183–184. [Google Scholar] [CrossRef]
  25. Zutic, M.; Cirkovic, I.; Pavlovic, L.; Asanin, J.; Jovanovic, S.; Zutic, J.; Asanin, R. First isolation of methicillin-resistant Staphylococcus aureus from pigs’ clinical samples in Serbia. Acta Vet. Brno 2012, 81, 225–227. [Google Scholar] [CrossRef]
  26. Zutic, M.; Cirkovic, I.; Pavlovic, L.; Zutic, J.; Asanin, J.; Radanovic, O.; Pavlovic, N. Occurrence of methicillin-resistant Staphylococcus aureus in milk samples from Serbian cows with subclinical mastitis. Afr. J. Microbiol. Res. 2012, 6, 5887–5889. [Google Scholar] [CrossRef]
  27. Asanin, J.; Misic, D.; Aksentijevic, K.; Tambur, Z.; Rakonjac, B.; Kovacevic, I.; Spergser, J.; Loncaric, I. Genetic profiling and comparison of human and animal methicillin-resistant Staphylococcus aureus (MRSA) isolates from Serbia. Antibiotics 2019, 8, 26. [Google Scholar] [CrossRef] [PubMed]
  28. Pajić, M.J.; Rašić, Z.B.; Velebit, B.M.; Boboš, S.F.; Mihajlović-Ukropina, M.M.; Radinović, M.Ž.; Galfi, A.L.; Petković, J.M.; Trojačanec, S.I. The Prevalence of Methicillin Resistance and Panton-Valentine leukocidin Synthesis Genes in Staphylococcus aureus Isolates of Bovine and Human Origin. Vet. Arhiv 2014, 84, 205–214. Available online: https://hrcak.srce.hr/121942 (accessed on 23 March 2023).
  29. Bulajic, S.; Colovic, S.; Misic, D.; Djordjevic, J.; Savic-Radovanovic, R.; Asanin, J.; Ledina, T. Enterotoxin production and antimicrobial susceptibility in Staphylococci isolated from traditional raw milk cheeses in Serbia. J. Environ. Sci. Health B 2017, 52, 864–870. [Google Scholar] [CrossRef]
  30. Kreausukon, K.; Fetsch, A.; Kraushaar, B.; Alt, K.; Müller, K.; Krömker, V.; Zessin, K.-H.; Käsbohrer, A.; Tenhagen, B.-A. Prevalence, antimicrobial resistance, and molecular characterization of methicillin-resistant Staphylococcus aureus from bulk tank milk of dairy herds. J. Dairy Sci. 2012, 95, 4382–4388. [Google Scholar] [CrossRef]
  31. Paterson, G.K.; Morgan, F.J.E.; Harrison, E.M.; Peacock, S.J.; Parkhill, J.; Zadoks, R.N.; Holmes, M.A. Prevalence and properties of mecC methicillin-resistant Staphylococcus aureus (MRSA) in bovine bulk tank milk in Great Britain. J. Antimicrob. Chemother. 2014, 69, 598–602. [Google Scholar] [CrossRef]
  32. Cortimiglia, C.; Luini, M.; Bianchini, V.; Marzagalli, L.; Vezzoli, F.; Avisani, D.; Bertoletti, M.; Ianzano, A.; Franco, A.; Battisti, A. Prevalence of Staphylococcus aureus and of methicillin-resistant S. aureus clonal complexes in bulk tank milk from dairy cattle herds in Lombardy Region (Northern Italy). Epidemiol. Infect. 2016, 144, 3046–3051. [Google Scholar] [CrossRef]
  33. Locatelli, C.; Cremonesi, P.; Bertocchi, L.; Zanoni, M.G.; Barberio, A.; Drigo, I.; Varisco, G.; Castiglioni, B.; Bronzo, V.; Moroni, P. Short communication: Methicillin-resistant Staphylococcus aureus in bulk tank milk of dairy cows and effect of swine population density. J. Dairy Sci. 2016, 99, 2151–2156. [Google Scholar] [CrossRef]
  34. Tenhagen, B.-A.; Alt, K.; Pfefferkorn, B.; Wiehle, L.; Käsbohrer, A.; Fetsch, A. Short communication: Methicillin-resistant Staphylococcus aureus in conventional and organic dairy herds in Germany. J. Dairy Sci. 2018, 101, 3380–3386. [Google Scholar] [CrossRef] [PubMed]
  35. Cui, C.; Ba, X.; Holmes, M.A. Prevalence and characterization of mecC MRSA in bovine bulk tank milk in Great Britain, 2017–2018. JAC Antimicrob. Resist. 2021, 3, dlaa125. [Google Scholar] [CrossRef] [PubMed]
  36. Magiorakos, A.-P.; Srinivasan, A.; Carey, R.B.; Carmeli, Y.; Falagas, M.E.; Giske, C.G.; Harbarth, S.; Hindler, J.F.; Kahlmeter, G.; Olsson-Liljequist, B.; et al. Multidrug-resistant, extensively drug-resistant and pandrug-resistant bacteria: An international expert proposal for interim standard definitions for acquired resistance. Clin. Microbiol. Infect. 2012, 18, 268–281. [Google Scholar] [CrossRef] [PubMed]
  37. MAEP-RS (Ministry of Agriculture and Environmental Protection Republic of Serbia). Republic of Serbia IPARD Programme for 2014–2020. 2019. Available online: https://uap.gov.rs/wp-content/uploads/2017/11/IPARDII-final-III-modification-ENG27062019.pdf (accessed on 12 September 2023).
  38. SORS (Statistical Office of the Republic of Serbia). Area, Livestock (LSU), Labour Force and Standard Output (SO) by Legal Status of Holding and Agricultural Size of Farm (UAA). STAT Database (Agriculture, Forestry and Fishery, Census of Agriculture, Farm Structure 2018, Key Farm Variables); SORS: Belgrade, Serbia, 2018. Available online: https://data.stat.gov.rs/Home/Result/1300020101?languageCode=en-US&displayMode=table (accessed on 27 February 2023).
  39. Ostojić Andrić, D.; Hristov, S.; Petrović, M.M.; Pantelić, V.; Bojkovski, J.; Novaković, Ž.; Lazarević, M.; Nikšić, D. Housing Conditions and Welfare of Dairy Cows in Serbia. In Proceedings of the 4th International Congress ‘New Perspectives and Challenges of Sustainable Livestock Production’, Belgrade, Serbia, 7–9 October 2015; pp. 62–73. Available online: https://r.istocar.bg.ac.rs/bitstream/handle/123456789/591/Proceedings-2015-Dusica-Ostojic-Andric.pdf?sequence=1&isAllowed=y (accessed on 12 September 2023).
  40. Relic, R.; Lakic, N.; Jankovic, L.; Davidovic, V.; Staric, J.; Jezek, J. Factors affecting rearing practices and health of calves on family farms. Span. J. Agric. Res. 2021, 19, e0501. [Google Scholar] [CrossRef]
  41. Virgin, J.E.; Van Slyke, T.M.; Lombard, J.E.; Zadoks, R.N. Short Communication: Methicillin-resistant Staphylococcus aureus detection in US bulk tank milk. J. Dairy Sci. 2009, 92, 4988–4991. [Google Scholar] [CrossRef]
  42. Spohr, M.; Rau, J.; Friedrich, A.; Klittich, G.; Fetsch, A.; Guerra, B.; Hammerl, J.A.; Tenhagen, B.-A. Methicillin-resistant Staphylococcus aureus (MRSA) in three dairy herds in southwest Germany. Zoonoses Public Health 2011, 58, 252–261. [Google Scholar] [CrossRef]
  43. Haran, K.P.; Godden, S.M.; Boxrud, D.; Jawahir, S.; Bender, J.B.; Sreevatsan, S. Prevalence and characterization of Staphylococcus aureus, including methicillin-resistant Staphylococcus aureus, isolated from bulk tank milk from Minnesota dairy farms. J. Clin. Microbiol. 2012, 50, 688–695. [Google Scholar] [CrossRef]
  44. Parisi, A.; Caruso, M.; Normanno, G.; Latorre, L.; Sottili, R.; Miccolupo, A.; Fraccalvieri, R.; Santagada, G. Prevalence, antimicrobial susceptibility and molecular typing of methicillin-resistant Staphylococcus aureus (MRSA) in bulk tank milk from southern Italy. Food Microbiol. 2016, 58, 36–42. [Google Scholar] [CrossRef]
  45. Locatelli, C.; Cremonesi, P.; Caprioli, A.; Carfora, V.; Ianzano, A.; Barberio, A.; Morandi, S.; Casula, A.; Castiglioni, B.; Bronzo, V.; et al. Occurrence of methicillin-resistant Staphylococcus aureus in dairy cattle herds, related swine farms, and humans in contact with herds. J. Dairy Sci. 2017, 100, 608–619. [Google Scholar] [CrossRef]
  46. Papadopoulos, P.; Papadopoulos, T.; Angelidis, A.S.; Boukouvala, E.; Zdragas, A.; Papa, A.; Hadjichristodoulou, C.; Sergelidis, D. Prevalence of Staphylococcus aureus and of methicillin-resistant S. aureus (MRSA) along the production chain of dairy products in north-western Greece. Food Microbiol. 2018, 69, 43–50. [Google Scholar] [CrossRef] [PubMed]
  47. Selim, A.; Kelis, K.; AlKahtani, M.D.F.; Albohairy, F.M.; Attia, K.A. Prevalence, antimicrobial susceptibilities and risk factors of Methicillin resistant Staphylococcus aureus (MRSA) in dairy bovines. BMC Vet. Res. 2022, 18, 293. [Google Scholar] [CrossRef] [PubMed]
  48. Tenhagen, B.A.; Alt, K.; Grobbel, M.; Maurischat, S. MRSA in bulk tank milk of dairy herds in Germany—Changes over time. Tierarztl. Prax. Ausg. G Grosstiere Nutztiere. 2023, 51, 63–69. [Google Scholar] [CrossRef] [PubMed]
  49. Antoci, E.; Pinzone, M.R.; Nunnari, G.; Stefani, S.; Cacopardo, B. Prevalence and Molecular Characteristics of Methicillin-Resistant Staphylococcus aureus (MRSA) among Subjects Working on Bovine Dairy Farms. Infez. Med. 2013, 21, 125–129. Available online: https://www.infezmed.it/media/journal/Vol_21_2_2013_5.pdf (accessed on 12 September 2023). [PubMed]
  50. Schnitt, A.; Lienen, T.; Wichmann-Schauer, H.; Cuny, C.; Tenhagen, B.A. The occurrence and distribution of livestock-associated methicillin-resistant Staphylococcus aureus ST398 on German dairy farms. J. Dairy Sci. 2020, 103, 11806–11819. [Google Scholar] [CrossRef] [PubMed]
  51. Huber, H.; Ziegler, D.; Pflüger, V.; Vogel, G.; Zweifel, C.R.; Stephan, R. Prevalence and characteristics of methicillin-resistant coagulase-negative staphylococci from livestock, chicken carcasses, bulk tank milk, minced meat, and contact persons. BMC Vet. Res. 2011, 7, 6. [Google Scholar] [CrossRef]
  52. Frey, Y.; Rodriguez, J.P.; Thomann, A.; Schwendener, S.; Perreten, V. Genetic characterization of antimicrobial resistance in coagulase-negative staphylococci from bovine mastitis milk. J. Dairy Sci. 2013, 96, 2247–2257. [Google Scholar] [CrossRef]
  53. Cicconi-Hogan, K.M.; Belomestnykh, N.; Gamroth, M.; Ruegg, P.L.; Tikofsky, L.; Schukken, Y.H. Short communication: Prevalence of methicillin resistance in coagulase-negative staphylococci and Staphylococcus aureus isolated from bulk milk on organic and conventional dairy farms in the United States. J. Dairy Sci. 2014, 97, 2959–2964. [Google Scholar] [CrossRef]
  54. Fisher, E.A.; Paterson, G.K. Prevalence and characterisation of methicillin-resistant staphylococci from bovine bulk tank milk in England and Wales. J. Glob. Antimicrob. Resist. 2020, 22, 139–144. [Google Scholar] [CrossRef]
  55. Rolo, J.; Worning, P.; Boye Nielsen, J.; Sobral, R.; Bowden, R.; Bouchami, O.; Damborg, P.; Guardabassi, L.; Perreten, V.; Westh, H.; et al. Evidence for the evolutionary steps leading to mecA-mediated β-lactam resistance in staphylococci. PLoS Genet. 2017, 13, e1006674. [Google Scholar] [CrossRef]
  56. Traversari, J.; van den Borne, B.H.P.; Dolder, C.; Thomann, A.; Perreten, V.; Bodmer, M. Non-aureus staphylococci species in the teat canal and milk in four commercial Swiss dairy herds. Front. Vet. Sci. 2019, 6, 186. [Google Scholar] [CrossRef] [PubMed]
  57. Rosa, N.M.; Penati, M.; Fusar-Poli, S.; Addis, M.F.; Tola, S. Species identification by MALDI-TOF MS and gap PCR-RFLP of non-aureus Staphylococcus, Mammaliicoccus, and Streptococcus spp. associated with sheep and goat mastitis. Vet. Res. 2022, 53, 84. [Google Scholar] [CrossRef]
  58. Nobrega, D.B.; Naushad, S.; Naqvi, S.A.; Condas, L.A.Z.; Saini, V.; Kastelic, J.P.; Luby, C.; De Buck, J.; Barkema, H.W. Prevalence and genetic basis of antimicrobial resistance in non-aureus staphylococci isolated from Canadian dairy herds. Front. Microbiol. 2018, 9, 256. [Google Scholar] [CrossRef] [PubMed]
  59. Martin, N.H.; Evanowski, R.L.; Wiedmann, M. Invited review: Redefining raw milk quality—Evaluation of raw milk microbiological parameters to ensure high-quality processed dairy products. J. Dairy Sci. 2023, 106, 1502–1517. [Google Scholar] [CrossRef]
  60. Zigo, F.; Vasil’, M.; Ondrašovicová, S.; Výrostková, J.; Bujok, J.; Pecka-Kielb, E. Maintaining optimal mammary gland health and prevention of mastitis. Front. Vet. Sci. 2021, 8, 607311. [Google Scholar] [CrossRef] [PubMed]
  61. Latorre, A.A.; Pachá, P.A.; González-Rocha, G.; San Martín, I.; Quezada-Aguiluz, M.; Aguayo-Reyes, A.; Bello-Toledo, H.; Oliva, R.; Estay, A.; Pugin, J.; et al. On-farm surfaces in contact with milk: The role of Staphylococcus aureus-containing biofilms for udder health and milk quality. Foodborne Pathog. Dis. 2020, 17, 44–51. [Google Scholar] [CrossRef]
  62. Pacha, P.A.; Munoz, M.A.; González-Rocha, G.; San Martín, I.; Quezada-Aguiluz, M.; Aguayo-Reyes, A.; Bello-Toledo, H.; Latorre, A.A. Molecular diversity of Staphylococcus aureus and the role of milking equipment adherences or biofilm as a source for bulk tank milk contamination. J. Dairy Sci. 2021, 104, 3522–3531. [Google Scholar] [CrossRef]
  63. Smistad, M.; Bakka, H.C.; Sølverød, L.; Jørgensen, H.J.; Wolff, C. Prevalence of udder pathogens in milk samples from Norwegian dairy cows recorded in a national database in 2019 and 2020. Acta Vet. Scand. 2023, 65, 19. [Google Scholar] [CrossRef]
  64. Nonnemann, B.; Lyhs, U.; Svennesen, L.; Kristensen, K.A.; Klaas, I.C.; Pedersen, K. Bovine mastitis bacteria resolved by MALDI-TOF mass spectrometry. J. Dairy Sci. 2019, 102, 2515–2524. [Google Scholar] [CrossRef]
  65. Conesa, A.; Dieser, S.; Barberis, C.; Bonetto, C.; Lasagno, M.; Vay, C.; Odierno, L.; Porporatto, C.; Raspanti, C. Differentiation of non-aureus staphylococci species isolated from bovine mastitis by PCR-RFLP of groEL and gap genes in comparison to MALDI-TOF mass spectrometry. Microb. Pathog. 2020, 149, 104489. [Google Scholar] [CrossRef]
  66. Hamel, J.; Zhang, Y.; Went, N.; Krömker, V. Non-S. aureus staphylococci (NAS) in milk samples: Infection or contamination? Vet. Microbiol. 2020, 242, 108594. [Google Scholar] [CrossRef] [PubMed]
  67. Bergeron, M.; Dauwalder, O.; Gouy, M.; Freydiere, A.M.; Bes, M.; Meugnier, H.; Benito, Y.; Etienne, J.; Lina, G.; Vandenesch, F.; et al. Species identification of staphylococci by amplification and sequencing of the tuf gene compared to the gap gene and by matrix-assisted laser desorption ionization time-of-flight mass spectrometry. Eur. J. Clin. Microbiol. Infect. Dis. 2011, 30, 343–354. [Google Scholar] [CrossRef] [PubMed]
  68. Mahmmod, Y.S.; Nonnemann, B.; Svennesen, L.; Pedersen, K.; Klaas, I.C. Typeability of MALDI-TOF assay for identification of non-aureus staphylococci associated with bovine intramammary infections and teat apex colonization. J. Dairy Sci. 2018, 101, 9430–9438. [Google Scholar] [CrossRef] [PubMed]
  69. Titouche, Y.; Akkou, M.; Houali, K.; Auvray, F.; Hennekinne, J.-A. Role of milk and milk products in the spread of methicillin-resistant Staphylococcus aureus in the dairy production chain. J. Food Sci. 2022, 87, 3699–3723. [Google Scholar] [CrossRef]
  70. Tegegne, H.A.; Florianová, M.; Gelbíčová, T.; Karpíšková, R.; Koláčková, I. Detection and molecular characterization of methicillin-resistant Staphylococcus aureus isolated from bulk tank milk of cows, sheep, and goats. Foodborne Pathog. Dis. 2018, 16, 68–73. [Google Scholar] [CrossRef]
  71. Hansen, J.E.; Ronco, T.; Stegger, M.; Sieber, R.N.; Fertner, M.E.; Martin, H.L.; Farre, M.; Toft, N.; Larsen, A.R.; Pedersen, K. LA-MRSA CC398 in dairy cattle and veal calf farms indicates spillover from pig production. Front. Microbiol. 2019, 10, 2733. [Google Scholar] [CrossRef]
  72. Kadlec, K.; Entorf, M.; Peters, T. Occurrence and characteristics of livestock-associated methicillin-resistant Staphylococcus aureus in quarter milk samples from dairy cows in Germany. Front. Microbiol. 2019, 10, 1295. [Google Scholar] [CrossRef]
  73. Lienen, T.; Schnitt, A.; Cuny, C.; Maurischat, S.; Tenhagen, B.-A. Phylogenetic tracking of LA-MRSA ST398 intra-farm transmission among animals, humans and the environment on German dairy farms. Microorganisms 2021, 9, 1119. [Google Scholar] [CrossRef]
  74. Alba, P.; Feltrin, F.; Cordaro, G.; Porrero, M.C.; Kraushaar, B.; Argudín, M.A.; Nykäsenoja, S.; Monaco, M.; Stegger, M.; Aarestrup, F.M.; et al. Livestock-associated methicillin resistant and methicillin susceptible Staphylococcus aureus sequence type (CC)1 in European farmed animals: High genetic relatedness of isolates from Italian cattle herds and humans. PLoS ONE 2015, 10, e0137143. [Google Scholar] [CrossRef]
  75. Monaco, M.; Pedroni, P.; Sanchini, A.; Bonomini, A.; Indelicato, A.; Pantosti, A. Livestock-associated methicillin-resistant Staphylococcus aureus responsible for human colonization and infection in an area of Italy with high density of pig farming. BMC Infect. Dis. 2013, 13, 258. [Google Scholar] [CrossRef]
  76. Huber, H.; Koller, S.; Giezendanner, N.; Stephan, R.; Zweifel, C. Prevalence and characteristics of meticillin-resistant Staphylococcus aureus in humans in contact with farm animals, in livestock, and in food of animal origin, Switzerland, 2009. Euro Surveill. 2010, 15, 19542. [Google Scholar] [CrossRef] [PubMed]
  77. Pilla, R.; Castiglioni, V.; Gelain, M.E.; Scanziani, E.; Lorenzi, V.; Anjum, M.; Piccinini, R. Long-term study of MRSA ST1, t127 mastitis in a dairy cow. Vet. Rec. 2012, 170, 312. [Google Scholar] [CrossRef]
  78. Juhász-Kaszanyitzky, E.; Jánosi, S.; Somogyi, P.; Dán, A.; van der Graaf-van Bloois, L.; van Duijkeren, E.; Wagenaar, J.A. MRSA transmission between cows and humans. Emerg. Infect. Dis. 2007, 13, 630–632. [Google Scholar] [CrossRef] [PubMed]
  79. Baig, S.; Rhod Larsen, A.; Martins Simões, P.; Laurent, F.; Johannesen, T.B.; Lilje, B.; Tristan, A.; Schaumburg, F.; Egyir, B.; Cirkovic, I.; et al. Evolution and population dynamics of clonal complex 152 community-associated methicillin-resistant Staphylococcus aureus. mSphere 2020, 5, e00226-20. [Google Scholar] [CrossRef] [PubMed]
  80. Basanisi, M.G.; La Bella, G.; Nobili, G.; Franconieri, I.; La Salandra, G. Genotyping of methicillin-resistant Staphylococcus aureus (MRSA) isolated from milk and dairy products in South Italy. Food Microbiol. 2017, 62, 141–146. [Google Scholar] [CrossRef]
  81. Crespo-Piazuelo, D.; Lawlor, P.G. Livestock associated methicillin resistant Staphylococcus aureus (LA MRSA) prevalence in humans in close contact with animals and measures to reduce on farm colonization. Ir. Vet. J. 2021, 74, 21. [Google Scholar] [CrossRef]
  82. Patel, K.; Godden, S.M.; Royster, E.E.; Crooker, B.A.; Johnson, T.J.; Smith, E.A.; Sreevatsan, S. Prevalence, antibiotic resistance, virulence and genetic diversity of Staphylococcus aureus isolated from bulk tank milk samples of U.S. dairy herds. BMC Genom. 2021, 22, 367. [Google Scholar] [CrossRef]
  83. Sabat, A.; Kosowska, K.; Poulsen, K.; Kasprowicz, A.; Sekowska, A.; van Den Burg, B.; Travis, J.; Potempa, J. Two allelic forms of the aureolysin gene (aur) within Staphylococcus aureus. Infect. Immun. 2000, 68, 973–976. [Google Scholar] [CrossRef]
  84. Cui, M.; Li, J.; Ali, T.; Kalim, K.; Wang, H.; Song, L.; Li, Z.; Ren, X.; Ma, F.; Zou, M.; et al. Emergence of livestock-associated MRSA ST398 from bulk tank milk, China. J. Antimicrob. Chemother. 2020, 75, 3471–3474. [Google Scholar] [CrossRef]
  85. Sieber, R.N.; Urth, T.R.; Petersen, A.; Moller, C.H.; Price, L.B.; Skov, R.L.; Larsen, A.R.; Stegger, M.; Larsen, J. Phage-mediated immune evasion and transmission of livestock-associated methicillin-resistant Staphylococcus aureus in humans. Emerg. Infect. Dis. 2020, 26, 2578–2585. [Google Scholar] [CrossRef]
  86. Chaguza, C.; Smith, J.T.; Bruce, S.A.; Gibson, R.; Martin, I.W.; Andam, C.P. Prophage-encoded immune evasion factors are critical for Staphylococcus aureus host infection, switching, and adaptation. Cell Genom. 2022, 2, 100194. [Google Scholar] [CrossRef] [PubMed]
  87. Monecke, S.; Müller, E.; Dorneanu, O.S.; Vremerǎ, T.; Ehricht, R. Molecular typing of MRSA and of clinical Staphylococcus aureus isolates from Iaşi, Romania. PLoS ONE 2014, 9, e97833. [Google Scholar] [CrossRef] [PubMed]
  88. Benito, D.; Lozano, C.; Rezusta, A.; Ferrer, I.; Vasquez, M.A.; Ceballos, S.; Zarazaga, M.; Revillo, M.J.; Torres, C. Characterization of tetracycline and methicillin resistant Staphylococcus aureus strains in a Spanish hospital: Is livestock-contact a risk factor in infections caused by MRSA CC398? Int. J. Med. Microbiol. 2014, 304, 1226–1232. [Google Scholar] [CrossRef] [PubMed]
  89. Hummerjohann, J.; Naskova, J.; Baumgartner, A.; Graber, H.U. Enterotoxin-producing Staphylococcus aureus genotype B as a major contaminant in Swiss raw milk cheese. J. Dairy Sci. 2014, 97, 1305–1312. [Google Scholar] [CrossRef]
  90. Earls, M.R.; Shore, A.C.; Brennan, G.I.; Simbeck, A.; Schneider-Brachert, W.; Vremerǎ, T.; Dorneanu, O.S.; Slickers, P.; Ehricht, R.; Monecke, S.; et al. A novel multidrug-resistant PVL-negative CC1-MRSA-IV clone emerging in Ireland and Germany likely originated in South-Eastern Europe. Infect. Genet. Evol. 2019, 69, 117–126. [Google Scholar] [CrossRef]
  91. Ikeda, T.; Tamate, N.; Yamaguchi, K.; Makino, S. Mass outbreak of food poisoning disease caused by small amounts of staphylococcal enterotoxins A and H. Appl. Environ. Microbiol. 2005, 71, 2793–2795. [Google Scholar] [CrossRef] [PubMed]
  92. Kadariya, J.; Smith, T.C.; Thapaliya, D. Staphylococcus aureus and staphylococcal food-borne disease: An ongoing challenge in public health. Biomed. Res. Int. 2014, 2014, 827965. [Google Scholar] [CrossRef]
  93. Bastos, C.P.; Bassani, M.T.; Mata, M.M.; Lopes, G.V.; da Silva, W.P. Prevalence and expression of staphylococcal enterotoxin genes in Staphylococcus aureus isolated from food poisoning outbreaks. Can. J. Microbiol. 2017, 63, 834–840. [Google Scholar] [CrossRef]
  94. Linage, B.; Rodríguez-Calleja, J.M.; Otero, A.; García-López, M.L.; Santos, J.A. Characterization of coagulase-positive staphylococci isolated from tank and silo ewe milk. J. Dairy Sci. 2012, 95, 1639–1644. [Google Scholar] [CrossRef]
  95. Vitale, M.; Gaglio, S.; Galluzzo, P.; Cascone, G.; Piraino, C.; Di Marco Lo Presti, V.; Alduina, R. Antibiotic resistance profiling, analysis of virulence aspects and molecular genotyping of Staphylococcus aureus isolated in Sicily, Italy. Foodborne Pathog. Dis. 2018, 15, 177–185. [Google Scholar] [CrossRef]
  96. Fang, R.; Cui, J.; Cui, T.; Guo, H.; Ono, H.K.; Park, C.-H.; Okamura, M.; Nakane, A.; Hu, D.-L. Staphylococcal enterotoxin C is an important virulence factor for mastitis. Toxins 2019, 11, 141. [Google Scholar] [CrossRef] [PubMed]
  97. Cuny, C.; Wieler, L.H.; Witte, W. Livestock-associated MRSA: The impact on humans. Antibiotics 2015, 4, 521–543. [Google Scholar] [CrossRef] [PubMed]
  98. van Wamel, W.J.B.; Rooijakkers, S.H.M.; Ruyken, M.; van Kessel, K.P.M.; van Strijp, J.A.G. The innate immune modulators staphylococcal complement inhibitor and chemotaxis inhibitory protein of Staphylococcus aureus are located on beta-hemolysin-converting bacteriophages. J. Bacteriol. 2006, 188, 1310–1315. [Google Scholar] [CrossRef] [PubMed]
  99. Naushad, S.; Naqvi, S.A.; Nobrega, D.; Luby, C.; Kastelic, J.P.; Barkema, H.W.; De Buck, J. Comprehensive virulence gene profiling of bovine non-aureus staphylococci based on whole-genome sequencing data. mSystems 2019, 4, e00098-18. [Google Scholar] [CrossRef]
  100. Belhout, C.; Boyen, F.; Vereecke, N.; Theuns, S.; Taibi, N.; Stegger, M.; de la Fé-Rodríguez, P.Y.; Bouayad, L.; Elgroud, R.; Butaye, P. Prevalence and molecular characterization of methicillin-resistant staphylococci (MRS) and mammaliicocci (MRM) in dromedary camels from Algeria: First detection of SCCmec-mecC hybrid in methicillin-resistant Mammaliicoccus lentus. Antibiotics 2023, 12, 674. [Google Scholar] [CrossRef]
  101. Diep, B.A.; Stone, G.G.; Basuino, L.; Graber, C.J.; Miller, A.; des Etages, S.A.; Jones, A.; Palazzolo-Ballance, A.M.; Perdreau-Remington, F.; Sensabaugh, G.F.; et al. The arginine catabolic mobile element and staphylococcal chromosomal cassette mec linkage: Convergence of virulence and resistance in the USA300 clone of methicillin-resistant Staphylococcus aureus. J. Infect. Dis. 2008, 197, 1523–1530. [Google Scholar] [CrossRef]
  102. Tong, C.; Wu, Z.; Zhao, X.; Xue, H. Arginine catabolic mobile elements in livestock-associated methicillin-resistant staphylococcal isolates from bovine mastitic milk in China. Front. Microbiol. 2018, 9, 1031. [Google Scholar] [CrossRef]
  103. ALIMS (Medicines and Medical Devices Agency of Serbia). Promet Veterinarskih Lekova 2019–2020; ALIMS (Agencija za Lekova I Medicinska Sredstva Srbije): Belgrade, Serbia, 2021. Available online: https://www.alims.gov.rs/wp-content/uploads/2022/03/PPL-VET-2020.pdf (accessed on 4 May 2023). (In Serbian)
  104. Schwarz, S.; Feßler, A.T.; Loncaric, I.; Wu, C.; Kadlec, K.; Wang, Y.; Shen, J. Antimicrobial resistance among staphylococci of animal origin. Microbiol. Spectr. 2018, 6, ARBA-0010-2017. [Google Scholar] [CrossRef]
  105. Monecke, S.; Kuhnert, P.; Hotzel, H.; Slickers, P.; Ehricht, R. Microarray based study on virulence-associated genes and resistance determinants of Staphylococcus aureus isolates from cattle. Vet. Microbiol. 2007, 125, 128–140. [Google Scholar] [CrossRef]
  106. Chua, K.; Laurent, F.; Coombs, G.; Grayson, M.L.; Howden, B.P. Not community-associated methicillin-resistant Staphylococcus aureus (CA-MRSA)! A clinician’s guide to community MRSA—Its evolving antimicrobial resistance and implications for therapy. Clin. Infect. Dis. 2011, 52, 99–114. [Google Scholar] [CrossRef]
  107. Said-Salim, B.; Dunman, P.M.; McAleese, F.M.; Macapagal, D.; Murphy, E.; McNamara, P.J.; Arvidson, S.; Foster, T.J.; Projan, S.J.; Kreiswirth, B.N. Global regulation of Staphylococcus aureus genes by Rot. J. Bacteriol. 2003, 185, 610–619. [Google Scholar] [CrossRef] [PubMed]
  108. Shoaib, M.; Xu, J.; Meng, X.; Wu, Z.; Hou, X.; He, Z.; Shang, R.; Zhang, H.; Pu, W. Molecular epidemiology and characterization of antimicrobial-resistant Staphylococcus haemolyticus strains isolated from dairy cattle milk in Northwest, China. Front. Cell. Infect. Microbiol. 2023, 13, 1183390. [Google Scholar] [CrossRef] [PubMed]
  109. Fergestad, M.E.; Touzain, F.; De Vliegher, S.; De Visscher, A.; Thiry, D.; Tchamba, C.N.; Mainil, J.G.; L’Abee-Lund, T.; Blanchard, Y.; Wasteson, Y. Whole genome sequencing of staphylococci isolated from bovine milk samples. Front. Microbiol. 2021, 12, 715851. [Google Scholar] [CrossRef]
  110. Osada, M.; Aung, M.S.; Urushibara, N.; Kawaguchiya, M.; Ohashi, N.; Hirose, M.; Kobayashi, N. Prevalence and antimicrobial resistance of Staphylococcus aureus and coagulase-negative Staphylococcus/Mammaliicoccus from retail ground meat: Identification of broad genetic diversity in fosfomycin resistance gene fosB. Pathogens 2022, 11, 469. [Google Scholar] [CrossRef] [PubMed]
  111. Lüthje, P.; Schwarz, S. Antimicrobial resistance of coagulase-negative staphylococci from bovine subclinical mastitis with particular reference to macrolide-lincosamide resistance phenotypes and genotypes. J. Antimicrob. Chemother. 2006, 57, 966–969. [Google Scholar] [CrossRef]
  112. Schauer, B.; Szostak, M.P.; Ehricht, R.; Monecke, S.; Feßler, A.T.; Schwarz, S.; Spergser, J.; Krametter-Frötscher, R.; Loncaric, I. Diversity of methicillin-resistant coagulase-negative Staphylococcus spp. and methicillin-resistant Mammaliicoccus spp. isolated from ruminants and New World camelids. Vet. Microbiol. 2021, 254, 109005. [Google Scholar] [CrossRef]
  113. EUCAST (European Committee on Antimicrobial Susceptibility Testing). Breakpoint Tables for Interpretation of MICs and Zone Diameters; Version 11.0; EUCAST: Växjö, Sweden, 2021; Available online: https://www.eucast.org/fileadmin/src/media/PDFs/EUCAST_files/Breakpoint_tables/v_11.0_Breakpoint_Tables.pdf (accessed on 4 May 2023).
  114. Schouls, L.M.; Spalburg, E.C.; van Luit, M.; Huijsdens, X.W.; Pluister, G.N.; van Santen-Verheuvel, M.G.; van der Heide, H.G.J.; Grundmann, H.; Heck, M.E.O.C.; de Neeling, A.J. Multiple-locus variable number tandem repeat analysis of Staphylococcus aureus: Comparison with pulsed-field gel electrophoresis and spa-typing. PLoS ONE 2009, 4, e5082. [Google Scholar] [CrossRef]
  115. Stegger, M.; Andersen, P.S.; Kearns, A.; Pichon, B.; Holmes, M.A.; Edwards, G.; Laurent, F.; Teale, C.; Skov, R.; Larsen, A.R. Rapid detection, differentiation and typing of methicillin-resistant Staphylococcus aureus harbouring either mecA or the new mecA homologue mecALGA251. Clin. Microbiol. Infect. 2012, 18, 395–400. [Google Scholar] [CrossRef]
  116. Bankevich, A.; Nurk, S.; Antipov, D.; Gurevich, A.A.; Dvorkin, M.; Kulikov, A.S.; Lesin, V.M.; Nikolenko, S.I.; Pham, S.; Prjibelski, A.D.; et al. SPAdes: A new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. 2012, 19, 455–477. [Google Scholar] [CrossRef]
  117. Gurevich, A.; Saveliev, V.; Vyahhi, N.; Tesler, G. QUAST: Quality assessment tool for genome assemblies. Bioinformatics 2013, 29, 1072–1075. [Google Scholar] [CrossRef]
  118. Wood, D.E.; Lu, J.; Langmead, B. Improved metagenomics analysis with Kraken 2. Genome Biol. 2019, 20, e257. [Google Scholar] [CrossRef] [PubMed]
  119. Richter, M.; Rosselló-Móra, R.; Glöckner, F.O.; Peplies, J. JSpeciesWS: A web server for prokaryotic species circumscription based on pairwise genome comparison. Bioinformatics 2016, 32, 929–931. [Google Scholar] [CrossRef] [PubMed]
  120. Richter, M.; Rosselló-Móra, R. Shifting the genomic gold standard for the prokaryotic species definition. Proc. Natl. Acad. Sci. USA 2009, 106, e19126-31. [Google Scholar] [CrossRef] [PubMed]
  121. 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.; et al. ResFinder 4.0 for predictions of phenotypes from genotypes. J. Antimicrob. Chemother. 2020, 75, 3491–3500. [Google Scholar] [CrossRef] [PubMed]
  122. Liu, B.; Zheng, D.D.; Zhou, S.Y.; Chen, L.H.; Yang, J. VFDB 2022: A general classification scheme for bacterial virulence factors. Nucleic Acids Res. 2022, 50, D912–D917. [Google Scholar] [CrossRef]
  123. Bartels, M.D.; Petersen, A.; Worning, P.; Nielsen, J.B.; Larner-Svensson, H.; Johansen, H.K.; Andersen, L.P.; Jarløv, J.O.; Boye, K.; Larsen, A.R.; et al. Comparing whole-genome sequencing with Sanger sequencing for spa typing of methicillin-resistant Staphylococcus aureus. J. Clin. Microbiol. 2014, 52, 4305–4308. [Google Scholar] [CrossRef]
  124. Jolley, K.A.; Maiden, M.C. BIGSdb: Scalable analysis of bacterial genome variation at the population level. BMC Bioinformatics 2010, 11, e595. [Google Scholar] [CrossRef]
  125. Kaya, H.; Hasman, H.; Larsen, J.; Stegger, M.; Johannesen, T.B.; Allesøe, R.L.; Lemvigh, C.K.; Aarestrup, F.M.; Lund, O.; Larsen, A.R. SCCmecFinder, a web-based tool for typing of Staphylococcal Cassette Chromosome mec in Staphylococcus aureus using whole-genome sequence data. mSphere 2018, 3, e00612-17. [Google Scholar] [CrossRef]
  126. Silva, M.; Machado, M.P.; Silva, D.N.; Rossi, M.; Moran-Gilad, J.; Santos, S.; Ramirez, M.; Carriço, J.A. chewBBACA: A complete suite for gene-by-gene schema creation and strain identification. Microb. Genom. 2018, 4, e000166. [Google Scholar] [CrossRef]
  127. Zhou, Z.; Alikhan, N.F.; Sergeant, M.J.; Luhmann, N.; Vaz, C.; Francisco, A.P.; Carriço, J.A.; Achtman, M. GrapeTree: Visualization of core genomic relationships among 100,000 bacterial pathogens. Genome Res. 2018, 28, 1395–1404. [Google Scholar] [CrossRef]
  128. 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]
  129. Schürch, A.C.; Arredondo-Alonso, S.; Willems, R.J.L.; Goering, R.V. Whole genome sequencing options for bacterial strain typing and epidemiologic analysis based on single nucleotide polymorphism versus gene-by-gene-based approaches. Clin. Microbiol. Infect. 2018, 24, 350–354. [Google Scholar] [CrossRef] [PubMed]
  130. EU (European Union). Corrigendum to Regulation (EC) No 853/2004 of the European Parliament and of the Council of 29 April 2004 Laying down Specific Hygiene Rules for Food of Animal Origin. Off. J. Eur. Union 2004, L 226, 22–82. Available online: https://eur-lex.europa.eu/eli/reg/2004/853/corrigendum/2004-06-25/oj (accessed on 11 September 2023).
  131. SORS (Statistical Office of the Republic of Serbia). Farms and Heads of Animals by Livestock Units (LSU). STAT Database (Agriculture, Census of Agriculture, Census 2012, Livestock and Bees); SORS: Belgrade, Serbia, 2012. Available online: https://data.stat.gov.rs/Home/Result/1300010204?languageCode=en-US (accessed on 6 June 2023).
Figure 1. cgMLST phylogenetic tree of methicillin-resistant Staphylococcus spp. isolates. Violet and orange squares next to the tree denote the presence of antimicrobial resistance and virulence genes, respectively. Light colored squares represent genes with <100% identity and <100% coverage to the reference gene. Clusters of isolates differing in ≤24 cgMLST alleles are highlighted in red rectangles.
Figure 1. cgMLST phylogenetic tree of methicillin-resistant Staphylococcus spp. isolates. Violet and orange squares next to the tree denote the presence of antimicrobial resistance and virulence genes, respectively. Light colored squares represent genes with <100% identity and <100% coverage to the reference gene. Clusters of isolates differing in ≤24 cgMLST alleles are highlighted in red rectangles.
Antibiotics 12 01529 g001
Figure 2. cgMLST phylogenetic tree of methicillin-resistant Mammaliicoccus spp. isolates. Violet squares next to the tree represent the presence of antimicrobial resistance genes. Light violet squares represent genes with <100% identity and <100% coverage to the reference gene. Clusters of isolates differing in ≤24 cgMLST alleles are highlighted in red rectangles.
Figure 2. cgMLST phylogenetic tree of methicillin-resistant Mammaliicoccus spp. isolates. Violet squares next to the tree represent the presence of antimicrobial resistance genes. Light violet squares represent genes with <100% identity and <100% coverage to the reference gene. Clusters of isolates differing in ≤24 cgMLST alleles are highlighted in red rectangles.
Antibiotics 12 01529 g002
Table 1. Univariate risk factor analysis for the occurrence of methicillin-resistant staphylococci (MRS) and mammaliicocci (MRM) in dairy farms. Statistically significant p values are highlighted in bold.
Table 1. Univariate risk factor analysis for the occurrence of methicillin-resistant staphylococci (MRS) and mammaliicocci (MRM) in dairy farms. Statistically significant p values are highlighted in bold.
VariablenMRS/MRM-
Positive Farms (%)
OR95% CIp
Milking system
Bucket 12516.81.00 0.008
Pipeline13732.82.421.34–4.37
Robot333.32.480.26–28.57
Milking parlor185.60.290.04–2.31
Housing system
Tie-stall26325.11.00 0.146
Free-stall2010.00.330.07–1.47
Herd size (no. of dairy cows)
>1019424.21.00 0.908
<108923.60.970.54–1.74
Breed
Simmental15728.01.00 0.216
Holstein Friesian8319.30.610.32–1.17
Mix4318.60.590.25–1.36
TBC (CFU/mL)
<100,00019317.61.00 0.001
>100,0009037.82.841.61–4.99
SCC per mL
<400,00022722.51.00 0.218
>400,0005630.41.500.79–2.88
Total28324.0
CFU, colony forming units; CI, confidence interval; OR, odds ratio; SCC, somatic cell count; TBC, total bacterial count.
Table 2. Adjusted odds ratios (OR) and 95% confidence intervals (CI) for the risk factors significantly associated with the occurrence of methicillin-resistant staphylococci (MRS) and mammaliicocci (MRM) in dairy farms as assessed by multivariate logistic regression. Statistically significant p values are highlighted in bold.
Table 2. Adjusted odds ratios (OR) and 95% confidence intervals (CI) for the risk factors significantly associated with the occurrence of methicillin-resistant staphylococci (MRS) and mammaliicocci (MRM) in dairy farms as assessed by multivariate logistic regression. Statistically significant p values are highlighted in bold.
VariableAdjusted OR95% CIp
Milking system
Bucket1.00
Pipeline2.511.37–4.590.003
Robot2.570.21–31.790.461
Milking parlor0.370.05–3.000.353
TBC (CFU/mL)
<100,0001.00
>100,0002.781.55–4.970.001
CFU, colony forming units; CI, confidence interval; OR, odds ratio; TBC, total bacterial count.
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

Kos, A.; Papić, B.; Golob, M.; Avberšek, J.; Kušar, D.; Ledina, T.; Đorđević, J.; Bulajić, S. Genomic Insights into Methicillin-Resistant Staphylococci and Mammaliicocci from Bulk Tank Milk of Dairy Farms in Serbia. Antibiotics 2023, 12, 1529. https://doi.org/10.3390/antibiotics12101529

AMA Style

Kos A, Papić B, Golob M, Avberšek J, Kušar D, Ledina T, Đorđević J, Bulajić S. Genomic Insights into Methicillin-Resistant Staphylococci and Mammaliicocci from Bulk Tank Milk of Dairy Farms in Serbia. Antibiotics. 2023; 12(10):1529. https://doi.org/10.3390/antibiotics12101529

Chicago/Turabian Style

Kos, Andrea, Bojan Papić, Majda Golob, Jana Avberšek, Darja Kušar, Tijana Ledina, Jasna Đorđević, and Snežana Bulajić. 2023. "Genomic Insights into Methicillin-Resistant Staphylococci and Mammaliicocci from Bulk Tank Milk of Dairy Farms in Serbia" Antibiotics 12, no. 10: 1529. https://doi.org/10.3390/antibiotics12101529

APA Style

Kos, A., Papić, B., Golob, M., Avberšek, J., Kušar, D., Ledina, T., Đorđević, J., & Bulajić, S. (2023). Genomic Insights into Methicillin-Resistant Staphylococci and Mammaliicocci from Bulk Tank Milk of Dairy Farms in Serbia. Antibiotics, 12(10), 1529. https://doi.org/10.3390/antibiotics12101529

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