Front. Microbiol.Frontiers in MicrobiologyFront. Microbiol.1664-302XFrontiers Media S.A.10.3389/fmicb.2023.1216850MicrobiologyOriginal ResearchWhole genome sequence-based analysis of Staphylococcus aureus isolated from bovine mastitis in Thuringia, GermanyMoawadAmira A.12*El-AdawyHosny13LindeJörg1JostInes4TanjaGärtner4KatjaHruschka4KarstenDonat4NeubauerHeinrich1MoneckeStefan567TomasoHerbert11Institute of Bacterial Infections and Zoonoses, Friedrich-Loeffler-Institut, Jena, Germany2Animal Health Research Institute, Agriculture Research Center (ARC), Giza, Egypt3Faculty of Veterinary Medicine, Kafrelsheikh University, Kafr El-Sheikh, Egypt4Tiergesundheitsdienst der Thueringer Tierseuchenkasse, Jena, Germany5Leibniz Institute of Photonic Technology (IPHT), Jena, Germany6InfectoGnostics Research Campus Jena e.V., Jena, Germany7Institute for Medical Microbiology and Virology, Dresden University Hospital, Dresden, Germany
Edited by: Rustam Aminov, University of Aberdeen, United Kingdom
Reviewed by: Mohammad Gharaibeh, Jordan University of Science and Technology, Jordan; Catherine M. Logue, University of Georgia, United States
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
Background
Bovine mastitis is a common disease of dairy cattle causing major economic losses due to reduced yield and poor quality of milk worldwide. The current investigation aimed to gain insight into the genetic diversity, antimicrobial resistance profiles and virulence associated factors of Staphylococcus (S.) aureus isolated from clinical bovine mastitis in dairy farms in Thuringia, Germany.
Methods
Forty Staphylococcus aureus isolates collected from clinical bovine mastitis cases from 17 Thuringian dairy farms were phenotyped and genetically characterized using whole genome sequencing.
Results
Out of 40 S. aureus, 30 (75%) were confirmed as methicillin resistant isolates. The isolates showed elevated antimicrobial resistance against penicillin, tetracycline and oxacillin, i.e., 77.5, 77.5, and 75%, respectively. Lower resistance rates were found against moxifloxacin, ciprofloxacin, gentamicin and trimethoprim/sulfamethoxazole, i.e., 35, 35, 30, and 22.5%, respectively. While resistance against clindamycin and erythromycin was rarely found (5 and 2.5%, respectively). All isolates were susceptible to linezolid, teicoplanin, vancomycin, tigecycline, fosfomycin, fusidic acid and rifampicin. These isolates were further allocated into five different sequence types: ST398 (n = 31), ST1074 (n = 4), ST504 (n = 3), ST582 (CC15) (n = 1) and ST479 (n = 1). These isolates were also assigned to seven clusters with up to 100 SNP which has facilitated geographical mapping and epidemiological distribution in Thuringia. Strains belonging to ST398 were classified into clusters 1, 2, 3, 4 and 7. The isolates of ST504 were of cluster 5, those of ST1074 were belonging to cluster 6. Resistance genes blaZ, blaI and blaR associated with penicillin resistance were found in 32 (80%) strains, all except one were belonging to ST398. Methicillin resistance associated mecA was identified in 30 (96.8%) isolates of ST398. All tetracycline and erythromycin resistant isolates were of ST398, and all harbored both tetM and ermA. About 90.3% of tetracycline resistant isolates assigned to ST398 were also carrying tetK gene. The point mutations parC_S80F, gyrA_S84L and parC_S80Y in gyrA and parC associated with quinolone resistance were found in all phenotypically resistant isolates to ciprofloxacin and moxifloxacin (n = 14). Sixty-eight virulence genes were identified among isolates. Both lukD/E and lukM/F-PV-P83 were identified in 22.5% of isolates, all were non-ST398.
Conclusion
In this study, ST398 had the highest potential to cause disease and had a massive prevalence in bovine mastitis cases. Five different sequence types and seven clusters were identified in the federal state of Thuringia. The circulation of some clusters in the same region over several years shows the persistence of cluster-associated infection despite the intensive veterinary care. On the other hand, some regions had different clusters at the same year or in different consecutive years. Different sequence types and associated different clusters of S. aureus were geographically widely distributed among dairy farms in Thuringia. The findings of this study show that various clusters have the potential to spread over a large geographical scale. The detection of LA-MRSA on dairy farms, which is known for cabapility to widely spread among different groups of animals, humans and their environment urges for the implementation of national wide strategic programs. The identification of CA-MRSA among the isolates such as ST398 poses a significant risk for the transmission of such strains between animals and humans on dairy farms.
bovine mastitisStaphylococcus aureusMRSAWGSST398antimicrobial resistanceThuringiasection-at-acceptanceAntimicrobials, Resistance and ChemotherapyBackground
Staphylococcus (S.) aureus is among the most common pathogens associated with bovine clinical and subclinical mastitis in Germany and can cause difficult-to-treat cases (Tenhagen et al., 2018; Kadlec et al., 2019). Staphylococcus aureus infections primarily occur during lactation and rarely in the dry period causing increased somatic cell counts (SCC) in bulk tank milk (BTM) and severe losses of milk due to poor cure rates (Schnitt and Tenhagen, 2020).
Bovine mastitis is usually treated with β-lactams, fluoroquinolones, cephalosporins, macrolides, and aminoglycosides (Krömker and Leimbach, 2017; TÄHAV, 2018; Ziesch et al., 2018; Bolte et al., 2020; Preine et al., 2022). The first line therapy for udder intramammary infections is the administration of benzylpenicillin/cloxacillin. An alternative treatment is the intramammary (IMM) or intramuscular (IM) combination of cefalexin and kanamycin (Bolte et al., 2020). Additional systemic antimicrobial treatment depends on the veterinarian’s decision and is restricted to severe cases including methicillin resistant S. aureus (MRSA) infections.
Recent allelic variations of MRSA strains show the independent development and spread of these strains to the German federal states. Transmission between dairy farms within one federal state is more likely attributed to local races, trading networks and similar raising methods (Lienen et al., 2020). Additionally, MRSA can be transmitted between animals and farm personnel causing sporadic severe infections in humans such as otitis, wound infections, sepsis, dermatitis, pneumonia, and endocarditis (Goerge et al., 2017).
Methicillin-resistance is primarily encoded by mecA and/or its homolog mecC. These genes are part of the mobile genetic element called “staphylococcal cassette chromosome mec” (SCCmec) (Wendlandt et al., 2013; Schnitt and Tenhagen, 2020).
The cassette often harbors additional resistance genes against several classes of antibiotics including aminoglycosides, macrolides, fluoroquinolones and the last resort drugs such as linezolid. Therefore, MRSA are mostly considered ‘multi-drug resistant’ superbugs (Cuny et al., 2017). Bovine mastitis caused by S. aureus is also fortified by numerous virulence factors such as toxins, leukocidins, superantigens (SAs) such as enterotoxins seA–seQ, and toxic shock syndrome toxin 1 (tsst-1). The accessory genome apart from the core genome, comprises approximately 15%–25% of the total S. aureus genome, leading to great variations in genomic content. The accessory genome includes mobile genetic elements (MGEs), in which the virulence and antimicrobial resistance (AMR) associated genes are included (Alibayov et al., 2014; Feßler et al., 2018; Hoekstra et al., 2020).
Population genetic studies have shown that the genotypes of S. aureus are often correlated with host species-specificity. The bovine associated sequence types (STs) 97, 133, 151, 479, and 771 have been confirmed as the predominant sequence types incriminated in bovine mastitis in some previous studies (Schlotter et al., 2012; Hoekstra et al., 2020).
The livestock associated MRSA (LA-MRSA) are less likely to cause severe human infections due to a lack of Immune invasion cluster (IEC) genes, Panton-Valentine leucocidin (PVL) or genes encoding the toxic shock syndrome toxin (TSST) (Cuny et al., 2015). Nonetheless, LA-MRSA was identified in bovine mastitis in Germany with the likelihood of transmission to human (Kadlec et al., 2019; Lienen et al., 2020).
The LA-MRSA of clonal complex CC398 has a broad spectrum of hosts including humans and is increasingly reported as the predominant sequence type found in German dairy farms (Tenhagen et al., 2018; Lienen et al., 2020; Schnitt et al., 2020).
The aim of this study was to identify and compare the genotypes, antimicrobial resistance profiles, virulence factors and geographic mapping of clinical S. aureus/MRSA strains isolated from milk samples on dairy farms in Thuringia, Germany. Whole-genome sequencing (WGS) of strains and sequence data analysis were applied to assess and investigate the epidemiologic situation of clinical mastitis and to evaluate the arising potential of public health risk.
MethodsIdentification and phenotypic characterization of <italic>Staphylococcus aureus</italic> isolates
Forty S. aureus isolates from milk were collected from cows with clinical mastitis from May 2020 to July 2022 from different 17 dairy herds in different locations in Thuringia (Supplementary Table S1). The samples were cultured and tested by the Animal Health Service of Thuringia (TGD) according to the guidelines of the German Veterinary Association (G.V.A, 2018) and international standards (NMC 2017) (Schnitt et al., 2020). All isolates were further investigated at the Institute of Bacterial Infections and Zoonoses, Friedrich-Loeffler-Institut, Germany. In order to identify and differentiate MRSA, all isolates were cultivated on Chromogenic agar medium (Mast Diagnostica GmbH, Reinfeld, Germany) at 37°C for 24 h under aerobic conditions. All isolates were confirmed as S. aureus using MALDI-TOF MS (Manukumar and Umesha, 2017; Yu et al., 2022). The interpretation of MALDI-TOF MS results was performed according to the manufacturer’s recommendation. The reference database was provided by Bruker (MBT-BDAL-8468). Multiplex real-time PCR assays were also performed for genus and species identification (Velasco et al., 2014; Wang et al., 2014).
Phenotypic antimicrobial susceptibility testing
The minimum inhibitory concentration (MIC) testing for S. aureus isolates was performed using a VITEK-2 system (VITEK-2, bioMérieux Deutschland GmbH, Nürtingen, Germany) with VITEK 2 AST-P592 Group B Staphylococci, Enterococci, and Streptococci Susceptibility Testing Cards (bioMérieux Deutschland GmbH, Nürtingen, Germany) according to the manufacturer’s instructions. The tested antimicrobial agents included in the VITEK 2 AST cards were penicillin, ampicillin, oxacillin, imipenem, gentamicin, ciprofloxacin, moxifloxacin, erythromycin, clindamycin, linezolid, teicoplanin, vancomycin, tetracycline, tigecycline, fosfomycin, fusidic acid, rifampicin, and trimethoprim/sulfamethoxazole. The interpretation of results was carried out according to the EUCAST standard (Version 13.0, 2023/NAK Deutschland).
Phenotypic methicillin/oxacillin resistance was also confirmed using E-TEST oxacillin (0.016–256 μg/mL) gradient strips (bioMérieux, Baden-Württemberg, Germany) (CLSI, 2020). For quality control of susceptibility testing, the S. aureus (DSM 2569) and E. coli (ATTC® 352,189) reference strains were used (CLSI, 2018).
DNA extraction and purification
A loopful from each fresh culture was suspended in 1 mL sterilized phosphate buffer saline (PBS) and heat inactivated at 96°C for 30 min. Genomic DNA was extracted and purified from bacterial cultures using QIAGEN Genomic-tip 20/G Kit (Qiagen GmbH, Hilden, Germany) according to the instructions of the manufacturer with a prior in-house modification step, i.e., adding of a lysis mixture (10–20 μL lysostaphin, 20 μL lysozyme, 2 μL ribonuclease A (2 μL of 10 mg/mL) and 45 μL proteinase K) followed by an incubation at 37°C for 2 h with slight shaking.
The concentration and quality of eluted DNA was determined photometrically using a Colibri spectrophotometer (Titertek, Berthold Technologies GmbH & Co. KG, Germany) and additionally measured using a Qubit 3 fluorometer (Fisher Scientific GmbH, Dreieich, Germany). The prepared DNA was preserved at −20°C for further investigations.
Determination of the presence of <italic>mec</italic>A
The PCR was carried out in a 25 μL reaction mix including 10 μL of LC480 Probes Master PCR Reaction Mix (Roche Diagnostics, Mannheim, Germany), 3.25 μL of PCR water (Roche Diagnostics, Mannheim, Germany), 1 μL of each primer, 0.25 μL of each probe, 5 μL of template and H2O for the non-template control. Amplification reactions were performed in a CFX96 Touch Real-Time PCR Detection Thermocycler System (Bio-Rad Laboratories GmbH, Feldkirchen, Germany).
PCR conditions were as follows; 1 cycle (50°C, 2 min.), 1 cycle (95°C, 10 min.) and 50 cycles (95°C, 20 s; 60°C, 40 s) (Velasco et al., 2014; Wang et al., 2014). S. aureus (DSM 2569) was used as positive control in each reaction.
Whole genome sequencing
The extracted DNA was sequenced using an Illumina MiSeq2000 platform. Sequencing libraries were created using the Nextera XT DNA Library Preparation Kit (Illumina Inc., San Diego, CA). Paired end sequencing producing 300 bp long reads was performed on an Illumina MiSeq instrument according to the manufacturer’s instructions (Illumina Inc., San Diego, CA). Raw sequencing data were deposited at the European Nucleotide Archive (ENA) as BioProject PRJEB61659. The bioinformatic analysis started with quality control of the raw paired end reads. The Linux based bioinformatics pipeline WGSBAC v. 2.2.01 was used for data analysis. Unless other settings are mentioned, all tools were used in their default standard settings. For quality control of raw reads, FastQC v. 0.11.7 (Andrews, 2018) was used and the coverage was calculated. Based on raw reads, WGSBAC performed assembly using Shovill v. 1.0.4.2 The quality of the assembled genomes was checked viaQUAST v. 5.0.2 (Bankevich et al., 2012). In order to check for potential contamination on both reads and assemblies, the pipeline uses Kraken 2 v. 1.1 (Wood et al., 2019) with the database Kraken2DB. For the detection of genetic markers for antimicrobial resistance and virulence, WGSBAC uses the software ABRicate (v. 0.8.10)3 and the databases ResFinder (Florensa et al., 2022), NCBI (Sayers et al., 2020), and Virulence Factor Database (VFDB) (Chen et al., 2005). Moreover, WGSBAC used AMRFinderPlus (v. 3.6.10) (Feldgarden et al., 2019) which detects genes and point mutations leading to AMR. AMRFinderPlus is used in organism-specific settings (i.e., Staphylococcus aureus). Platon is used for detection of plasmid borne contigs (Schwengers et al., 2020).
For genotyping, WGSBAC uses classical multilocus sequence typing (MLST) on assembled genomes using the software mlst v. 2.16.1 that incorporates the PubMLST database for the seven housekeeping genes of S. aureus.4 For high resolution genotyping, WGSBAC performed mapping-based genotyping using core genome single nucleotide polymorphisms (cgSNPs) identified by Snippy v. 4.3.6.5 The genome of S. aureus NCTC 8325 (accession NC_007795.1) was used as reference. RAxML (Randomized Axelerated Maximum Likelihood) v. 8 (Stamatakis, 2014) was used for phylogenetic tree construction based on cgSNP analysis. The tree was rooted to the reference genome and visualized using the interactive Tree of Life (iTOL) v. 4 web tool.6 Hierarchical clustering within the statistical language R was performed based on the pairwise SNP distances.
Results
All 40 isolates were identified as S. aureus using MALDI-TOF MS and confirmed by real-time PCR (Supplementary Table S1). Out of 40 S. aureus, 30 (77.5%) were confirmed as methicillin resistant S. aureus (MRSA) isolates by growing on chromogenic agar media and amplification of the mecA gene by real-time PCR (Supplementary Table S1).
1. Phenotypic antimicrobial resistance profiles
Thirty (75%), 31 (77.5%) and 31 (77.5%) S. aureus were phenotypically resistant to oxacillin, penicillin and tetracycline, respectively. Antimicrobial resistance against moxifloxacin, ciprofloxacin, gentamicin and trimethoprim/sulfamethoxazole was detected in 14 (35%), 14 (35%), 12 (30%) and 9 (22.5%) strains, respectively. Few isolates were resistant against clindamycin (2/40, 5%) and erythromycin (1/40, 2.5%). All isolates were sensitive to linezolid, teicoplanin, vancomycin, tigecycline, fosfomycin, fusidic acid and rifampicin. All oxacillin-resistant isolates were also resistant to penicillin and tetracycline (Table 1). Twenty-four strains (60%) were resistant to at least three classes of antibiotics and confirmed as multi-drug resistant (MDR) (Supplementary Table S1).
Phenotypic resistance profiles of 40 S. aureus isolated from clinical bovine mastitis in Thuringia, Germany.
Antimicrobial class
Antimicrobial agent
Resistance
n.
%
Beta-lactam antibiotics
Penicillin
31
77.5
Oxacillin
30
75
Aminoglycosides
Gentamicin
12
30
Fluoroquinolones
Ciprofloxacin
14
35
Moxifloxacin
14
35
Macrolides
Erythromycin
1
2.5
Lincomycin
Clindamycin
2
5
Oxazolidinones
Linezolid
0
0
Glycopeptides
Teicoplanin
0
0
Vancomycin
0
0
Tetracyclines
Tetracycline
31
77.5
Tigecycline
0
0
Phosphonic antibiotics
Fosfomycin
0
0
Fusidane class
Fusidic acid
0
0
Rifampicin
0
0
Sulfonamides plus trimethoprim
Trimethoprim/Sulfamethoxazole
9
22.5
2. Genetic characterization of Staphylococcus aureus isolates
Five different sequence types (ST) were identified. The majority of isolates (31/40) were assigned to ST398. Of note, chimerism with CC9, as frequently observed in livestock ST398 (Fetsch et al., 2017), was ruled out, based on an absence of Q2G1R6-cstB as well as on a presence of C5Q1F1-ycjY and G7ZTC1 (Monecke et al., 2018).
Four isolates were identified as ST1074 and three isolates as ST504. Both of these two sequence types belong to CC705. One isolate each was assigned to ST479 and ST582 (CC15) (Supplementary Table S1 and Figures 1, 2).
Phylogenetic analysis tree of 40 Staphylococcus aureus isolated from clinical mastitis cases from dairy farms in Thuringia using WGS analysis. The tree was rooted to the reference genome and visualized using the interactive Tree of Life (iTOL) v. 4 web tool (https://itol.embl.de/login.cgi). Hierarchical clustering within the statistical language R was performed based on the pairwise SNP distances. Data included are; Survey year, ST types, city, hclust maxSNPs_20, 50, and 100. Phenotypic resistance to penicillin, methicillin and tetracycline. Resistance genes; mecA/C, blaZ, tetM and tetK. Virulence genes lukD, E and lukF-PV.
Geographical distribution of different sequence types (STs) of S. aureus among clinical mastitis cases from dairy farms in Thuringia.
3. Genetic characterization of antimicrobial resistance
Thirty isolates were mecA positive, fully in accordance to phenotypic test results. All MRSA isolates identified in this study belonged to ST398. The most abundant SCCmec element was the SCC [mec VT + czrC] composite element as previously sequenced in SO385; GenBank AM990992.1. This was identified in twenty-six isolates. Three isolates harbored other variants of SCCmec VT. In one of them, an additional ccrA-1 recombinase gene was detected. Finally, a single isolate carried a SCCmec IVa element.
The penicillin resistance causative gene blaZ was was identified in 32 (80%) of isolates. All except one belonged to ST398 (one isolate was on ST582/CC15).
A simultaneous presence of both tetracycline resistance genes, tetM and tetK, was identified in 28 isolates that were all belonging to ST398, while tetM alone was identified in all ST398 (Chen et al., 2005) isolates. Data for all other resistance associated genes are shown in Supplementary Table S1.
4. Relatedness between phenotypic and genotypic resistance profiles
The isolates showed high similarity in phenotypic and genotypic resistance profiles for most of the tested drugs (Figure 1 and Supplementary Table S1).
The resistance gene blaZ gene was identified in all isolates that were phenotypically resistant to penicillin (100%). All phenotypically methicillin resistant isolates harbored mecA (100%). All tetracycline resistant isolates carried either tetM alone or both, tetK and tetM. All trimethoprim/sulfamethoxazole and erythromycin resistant isolates carried either dfrG or dfrK and ermA, respectively.
Additionally, a concordance of 100% was observed in case of ciprofloxacin and moxifloxacin resistance in relation to the carriage of point mutations related to the drug resistance.
5. Geographical distribution of genetic variants
In order to study the phylogeographic variation of S. aureus in Thuringia, clustering of cgSNP general genomic clusters (100 cgSNPs) and more closely related strains (20 cgSNP) were applied. Thirty out of 40 strains (75%) were clustered into seven clusters. The most prominent clusters were; 1 (22.5%), 2 (12.5%) and 3 (12.5%). The distribution of all clusters is shown in Figure 1 and Supplementary Table S1.
Clustering based on 100 cgSNPs has higher discriminatory power than MLST. Therefore, strains belonging to ST398 were separated into clusters 1, 2, 3, 4 and 7. The isolates of cluster 5 were identified as ST504. Isolates of cluster 6 were belonged to ST1074. The predominant sequence type (ST398) was confirmed in all MRSA of 5 different clusters and circulated in 14 locations (Figures 1, 2).
To further refine clustering of the closely related strains, 20cgSNPs were used as a cutoff. This clustered 15 out of 40 strains (37.5%) into four clusters. Cluster 1 contained eight strains which were all isolated in the district Unstrut-Hainich-Kreis between 2020 and 2021. Three strains all isolated in the district Nesse-Apfelstädt between 2021 and 2022 were grouped within cluster 2. Cluster 3 consisted of two strains collected in Neuhaus-Schierschnitz and Sonneberg in 2020 and 2021. Within cluster 4 two strains each from districts Altenburger Land and Unstrut-Hainich-Kreis in 2021 were grouped.
6. Geographical distribution of antimicrobial resistance
The geographical distribution of phenotypic antimicrobial resistance of S. aureus isolates based on 100 cgSNPs is illustrated in Figure 1 and Supplementary Table S1. Interestingly, all S. aureus of cluster 1 isolated in 2020 and 2021 in Unstrut-Hainich-Kreis had the same phenotypic antimicrobial resistance profiles. The same effect was noticed for isolates of clusters 3, 4, 6, and 7, that were isolated in 2020 and 2021 from different locations (Supplementary Table S1).
On the contrary, isolates of clusters 2 and 5 did not share the same phenotypic resistance profile among the same cluster.
7. Virulence associated genes
Sixty-eight virulence genes were identified in 40 S. aureus isolates using WGS (Supplementary Table S1). They were assigned to five functional categories (adherence related virulence factors (n = 14), toxins; including toxin system genes (n = 17), enterotoxins and enterotoxin-like genes (n = 9), iron acquisition and metabolism-related genes (n = 7), host immune evasion genes (n = 8) and exoenzymes (n = 13)).
Regarding leukotoxins, the hlg-locus comprising of lukS, lukF and hlgA were present in all isolates, as they were the chromosomal genes encoding a bi-component leukocidin lukA/B (=lukG/H). The genomic island associated leukotoxins lukD/E were present in nine out of forty isolates, i.e., all CC15, CC479 and CC705 isolates, but absent from all CC398. The phage-borne leukocidin genes lukM/lukF-P83 were identified in the CC479 and CC705 sequences, i.e., in 9 out of 40 isolates. The genes encoding Panton-Valentine leukocidin, lukF/S-PV were absent from all sequences.
The enterotoxin gene cluster egc, comprising seg, sei, selm, seln, selo and selu was detected in the CC479 and CC705 sequences, i.e., in 8 out of 40 isolates. The CC705 isolates additionally harbored an enterotoxin gene homologue, ORF CM14. Only a single isolate (the CC398 MRSA with the composite SCCmec VT + ccrA1 element) isolate carried one of the “classical” enterotoxin genes, sea. This sequence also contained genes sak, chp and scn (encoding staphylokinase, chemotaxis-inhibiting protein CHIPS and a staphylococcal complement inhibitor) that are frequently co-localized with sea on hemolysin-beta-converting phages in S. aureus strains from humans.
Genes encoding epidermal cell differentiation inhibitors or exfoliative toxins were not identified in any genome sequence examined.
Discussion
Staphylococcus (S.) aureus is one of the most prevalent bacterial species incriminated in clinical and subclinical bovine mastitis infections as reported by different studies that have focused on S. aureus in dairy herds (Tenhagen et al., 2006; Minst et al., 2012; Tenhagen et al., 2014; Bolte et al., 2020). The pathogen has been isolated very often from bovine mastitis outbreaks in Germany and has developed multiple resistances (Tenhagen et al., 2006; Minst et al., 2012; Bundestierärztekammer, 2015; Zimmermann, 2017; Bolte et al., 2020).
In this study, we examined the presence of both methicillin resistant and methicillin sensitive S. aureus (MRSA/MSSA) associated with clinical mastitis in dairy cattle in addition to associated antimicrobial resistance and virulence factors in 40 isolates from Thuringian dairy farms.
The most conspicuous result was the high prevalence of MRSA. Previous studies from the same region indicated an MRSA prevalence of a mere 2% among bovine mastitis isolates collected 10 years ago in the same geographical region (Schlotter et al., 2012). Although different sampling procedures and low isolate number hinders a substantiated epidemiological assessment, it is assumed that CC398 MRSA among cattle mastitis cases is not that rare anymore.
The MLST (generated using WGS data) in this study has clustered S. aureus isolates into sequence types (STs). The majority of isolates were assigned to ST398 (77.5%), ST1074 (10%), ST504 (7.5%), ST582 (2.5%) and ST479 (2.5%). ST504 and ST1074 are related to each other and both belong to CC705, which is a well-known bovine lineage. ST582 belongs to CC15 which is common among humans, suggesting anthropo-zoonotic transmission. The detection of these CC479, CC705 and CC398 strains from bovine milk is not surprising as it is consistent with most studies nationally and internationally (Hoekstra et al., 2018, 2020; Lienen et al., 2020; Schnitt et al., 2020; Campos et al., 2022). Except in Asia, where ST9 is more prevalent (Annamanedi et al., 2021), ST398 is the most common LA-MRSA lineage worldwide, which is responsible for a wide range of different infections, from minor skin lesions to invasive infections and even fatalities in humans. However, the presence of distinct strains within CC398 in this study should be noted, out of which some are adapted to humans and carry IEC genes, sometimes even PVL while others are associated only with livestock (Kashif et al., 2019).
Previous molecular examination of many LA-MRSA ST398 strains confirmed the low pathogenicity and low number of associated virulence factors in these strains (Kraushaar et al., 2017; Ji et al., 2021). These findings are in accordance with the low number of virulence factors identified in ST398 in the current study.
According to the national reference center for staphylococci in Germany, about 3.7% of human clinical MRSA infections are associated with LA-MRSA ST398 (Schnitt and Tenhagen, 2020). The transmission of LA-MRSA between humans and animals, including dairy cows has only been reported on a few occasions (Cuny et al., 2015; Schnitt and Tenhagen, 2020).
Many mastitis outbreaks caused by LA-MRSA ST398 have been reported in Germany (Tenhagen et al., 2018; Hansen et al., 2019). These findings agree with our finding that ST398 is the predominant identified ST in 77.5% of our isolates.
In the current study, we also found only one (2.5%) isolate assigned to ST479 which is known for its dominant lineage associated with bovine mastitis that carries variable resistance and virulence genes in European countries, including Germany (Hoekstra et al., 2020).
Moreover, health care acquired methicillin resistant S. aureus (HA-MRSA) was found in dairy cows in Germany (Tenhagen et al., 2018) and indeed HA-MRSA ST582 was identified once (2.5%) among isolates in this study.
Antimicrobial resistance
Antimicrobial resistance (AMR) usually spreads among S. aureus isolates due to the selective pressure as a result of extensive use of antimicrobials to treat bovine mastitis. The European Network for Optimization of Veterinary Antimicrobial Treatment (ENOVAT), COST Action CA 18217, has listed bovine mastitis as a priority infectious disease that contributes significantly to antimicrobial consumption in animals in European countries (Preine et al., 2022).
In the current study, resistance to the antibiotic classes, i.e., aminoglycosides, β-lactams, trimethoprim-, tetracyclines, macrolides, streptogramin, lincosamides, quinolones and phosphonic antibiotics was confirmed showing notable differences between the dairy farms of the current study.
The phenotypic resistance was mostly in agreement with bioinformatic predictions made from sequence data mined for AMR genes. Resistance to oxacillin was mediated in all MRSA strains in this study by the mecA gene while the variant mecC gene was not detected in this study. This gene is usually restricted to a few clonal complexes, including CC130, CC599 and CC425. The former ones occur in small wild mammals such as hedgehogs, also locally in Thuringia (Monecke et al., 2016). The latter has also been found in cattle (García-Álvarez et al., 2011), but CC425-MRSA-XI have, to the best of our knowledge, not yet been identified in Thuringia.
All oxacillin resistant isolates in this study were also resistant to penicillin and tetracycline. This finding might be attributed to the fact that hydrophilic molecules such as β-lactams, tetracyclines and some fluoroquinolones are particularly affected by the same changes in permeability of the outer membrane as a mechanism of bacterial resistance (Munita and Arias, 2016). It is a fact that all CC398 MRSA carry at least one tet gene, and usually even several.
The blanket dry cow therapy (DCT) with long acting β-lactam antimicrobials, including penicillin and cloxacillin, is still often used to treat intramammary S. aureus infections. Hence, it could be considered as a factor that might pose a selective pressure favoring the dissemination of MRSA strains (Schnitt et al., 2020).
Antimicrobial resistance against penicillin is primarily attributed to enzymatic inactivation of β-lactams and target site modification. Not unexpectedly, phenotypic resistance to penicillin was identified in 77.5% of all S. aureus isolates and in 100% of MRSA isolates in this study. This resistance was confirmed by identification of β-lactamase encoding genes blaI, blaZ, blaR and blaPC. Tenhagen et al., have accordingly earlier reported a 100% penicillin resistance in MRSA isolates in German dairy herds (Tenhagen et al., 2018). Another study recorded that blaZ was found in 87.9% of MRSA isolates from dairy farms (Lienen et al., 2020).
Tetracyclines are extensively used in veterinary medicine, thus promoting the occurrence of tetracycline resistant strains. Tetracycline resistance is often mediated by tetL, tetK and tetM, which code for membrane associated efflux proteins of the Major Facilitator Superfamily. In this study, phenotypic resistance to tetracycline was confirmed in 77.5% of all isolates, which was mediated by tetK and/or tetM. These findings were in agreement with results of other studies where all MRSA strains from dairy farms were resistant to tetracycline (Fessler et al., 2012; Kadlec et al., 2019; Lienen et al., 2020). Besides, tetracycline resistance was detected in 95.1% of the MRSA in BTM from German dairy farms (Tenhagen et al., 2018).
Macrolide, lincosamide, and streptogramin resistance in staphylococci is based on the presence of one or more erm genes that modify the target site(s) in the 23S rRNA and thereby inhibit the binding of these drugs. Resistance to erythromycin was rare among strains examined (2.5%) in this study. This finding fits well with reports that macrolide resistance in S. aureus retrieved from bovine mastitis is low (Lienen et al., 2020). The ermA gene, known to confer resistance to macrolides, lincosamides and streptogramin B (Leclercq, 2002), was identified in the sole erythromycin resistant isolate in this study. The isolate was also resistant to clindamycin.
In contrast to most mastitis therapeutics, which are administered by direct injection into the infected mammary gland, fluoroquinolones are usually administrated systemically in several Scandinavian countries. These countries have used a common strategy for treating clinical mastitis since several years. According to the “Nordic Guidelines for Mastitis Therapy” a restrictive use of penicillin is followed by a reduced application of cephalosporins and quinolones. Both cephalosporins and quinolones are considered as critically important antimicrobials (CIAs) belonging to human reserved list (Landin et al., 2011; Preine et al., 2022). Fluoroquinolones are used in dairy cattle and calves in Germany. However, published data on use and resistance against these drugs are limited. In this study, phenotypic resistance to quinolones, i.e., ciprofloxacin and moxifloxacin was found in 35% of all isolates. These isolates were all identified as MRSA. The resistance was reflected by the presence of point mutations parC_S80F, gyrA_S84L and parC_S80Y in all quinolone resistant isolates. Ciprofloxacin resistance in S. aureus is mediated by point mutations that alter ciprofloxacin targets, including topoisomerase IV (grlA, grlB) and DNA gyrase (gyrA, gyrB). The most common single point mutations, such as grlA E84K and grlA S80F/Y, are often found in ciprofloxacin resistant clinical isolates (Papkou et al., 2020). Amongst MRSA isolates from German dairy herds, quinolone resistance was either not found (Spohr et al., 2011; Kadlec et al., 2019), or at a low prevalence (8.3%) (Kreausukon, 2011).
Resistance to trimethoprim was detected in 22.5% of all isolates and 30% of MRSA isolates in the current study. That was quite lower than results obtained in previous studies in 2012 and 2021, where 43.2% or more than half of MRSA isolates from dairy farms were resistant to trimethoprim (Fessler et al., 2012; Lienen et al., 2020). The resistance to trimethoprim in this study was mediated by the dfrG and/or dfrK genes that were identified in all phenotypically resistant isolates. Both genes were only identified in CC398 isolates, which was in agreement with other studies that identified dfrK only in CC398-MRSA isolates in German dairy herds (Fessler et al., 2010; Lienen et al., 2020).
Staphylococci are non-target bacteria to pleuromutilin. However, the use of these drugs (instead of their use) in pig farming selects for multi resistant MRSA (Lienen et al., 2020). In this study, pleuromutilin resistance was confirmed by the detection of vgaA or vgaE genes. Thus, a hypothesis of pig-transmission should be taken in consideration. The vgaA gene was reported in MRSA strains from several farms located in different German federal states (Schwendener and Perreten, 2011; Feßler et al., 2018). Additionally, Hauschild et al, detected the vgaE gene in CC398 in dairy cattle (Hauschild et al., 2012). Although the vgaA gene was reported as the most widespread among the vga genes (Feßler et al., 2018), in the current study the vgaE (12.5%) gene was more common than vgaA (2.5%) across the MRSA strains.
Aminoglycosides are widely used in veterinary medicine (EMA, 2017). Enzymatic inactivation (acetylation and phosphorylation) are the main mechanisms for the resistance against them. Phenotypic resistance against gentamicin was identified in 36.7% of MRSA isolates in this study. All gentamicin resistant isolates belonged to ST398 MRSA strains and carried the resistance associated genes ant(6)-Ia, aac(6′)-Ie-aph(2″), aadD1, ant(9)-Ia, ant(4′)-Ia and aph(2″)-Ih. In another study, the resistance against gentamicin among MRSA CC398 isolated from German dairy farms was very similar to our results (38.4%) (Fessler et al., 2012).
An earlier study by Schlotter et al. in German dairy herds recorded phenotypic resistance to both linezolid and teicoplanin (38.9 and 5.6%) in MRSA-CC130 isolates, while resistance to vancomycin, tigecycline, fosfomycin, or rifampicin was not observed amongst the isolates (Schlotter et al., 2014). In this study, none of the isolates showed resistance to linezolid, tigecycline, vancomycin, fosfomycin or teicoplanin.
The presence of MDR efflux pumps norA, LmrS, mepA and mepR genes was observed in all S. aureus isolates in this study. However, not all of them were associated with expressed phenotypic drug resistance. Of note, other studies reported the high prevalence of the same elements in Staphylococcus spp. (Antiabong et al., 2017).
The data of earlier studies could be confirmed in the current study. Lower values of resistance against various agents might be caused by the limited number of strains investigated. There are obvious indications in this study for transmission of S. aureus strains from carrier host species like human or pig to cattle.
Virulence factors
Virulence factors (VFs) such as enterotoxins, leukocidins and haemolysins produced by staphylococci isolated from milk, could generate potential public health implications for consumers. S. aureus has the ability to survive inside mammary epithelial cells and immune cells after invading the teat canal (Mullarky et al., 2001; Barkema et al., 2006). Therefore, infections may persist over a long period and antimicrobial therapy is often hindered.
The distribution of a large number of VFs in isolates of this study was determined. These VFs were grouped into five functional categories.
Haemolysins α, β, γ and δ toxins of S. aureus play a major role in its pathogenicity. All 4 hemolysin genes were identified in all S. aureus isolates, which is consistent with other reports (Schmidt et al., 2017).
Leukocidins include Panton Valentine leukocidin (lukF/S-PV, lukM/lukFP83), lukA/B and lukD/E. They cause structural changes that facilitate membrane insertion and formation of octameric pores. In the current study, lukD/E and lukM/F-PV-P83 were found in non-ST398 isolates.
Enterotoxins are heat stable and usually persist after heat treatment in dairy products and remain active in the gastrointestinal system after ingestion. Enterotoxins are associated with staphylococcal food poisoning caused by cow milk and other dairy products (Benkerroum, 2018). Sea, seb, sec, sed and see represent the “classical types” of enterotoxins.
The distribution of enterotoxin genes here differed greatly among isolates and different STs. While among the most prevalent ST (398), only a single isolate contained the sea gene (n = 1), other enterotoxins sem, sei, seu, and sen were identified in other STs.
Conclusion
In the current study, S. aurus isolated from bovine clinical mastitis cases from farms of a single German federal state were investigated. Considering the limited number of isolates investigated in this study, a significant number of different sequence types and associated clusters was found and geographically widely distributed among dairy farms. Of note, the prevalence of MRSA was much higher than in previous studies on bovine mastitis in the same region. Persistance of clusters in the same region for different years shows that veterinary control measures do not have the impact to control disease. Detection of LA-MRSA in dairy farms, which is known for cabapility to spread easily among different groups of animals, humans, and in the environment suggests the need for urgent implementation of national wide strategic programs as already demanded in the animal health law (AHL) 2016 and recently confirmed by various EFSA reports. The identification of CA-MRSA such as ST398 among isolates poses a high risk for transmission and interplay of such strains between animals and humans on the dairy farms. The hygienic surveillance and awareness among the workers and veterinarians should be intensified in the future in order to expand the AMR and safeguard public health integrity.
The findings of this study show that various clusters have the potential to spread over a large geographical area (e.g., cluster 4), where some other clusters were more geographically restricted (e.g., cluster 1 and 3).
The mandatory antimicrobial susceptibility testing proves to be crucial in combination with prudent use of antimicrobials to avoid the emergence of resistance. The proposed action favours three of the current agricultural regulations, i.e., implementation of EU legislation, increase of income through avoiding losses due to reduced milk yield and decrease the veterinary costs. Most importantly is increasing the animal welfare through reducing suffering from clinical and subclinical mastitis infections in dairy farms in Germany.
As bovine mastitis caused by S. aureus especially MRSA is a world-wide problem, we aimed in this study to spot light on the situation in Thuringia as one of German federal states in Germany having relatively high number of cattle farms for milk production and suffering from recurrent mastitis problems of different causes.
This study aimed also to update the situation of the bovine mastitis caused by S. aureus/MRSA and using WGS for full genetic characterization in addition to antibiotic resistance prediction. These are among the main objectives of a joint project to innovate a fast robust tool for rapid detection of AMR in field cases of mastitis.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary material.
Author contributions
AAM, HE-A, HT, and HN: conceptualization. AAM, JL, HT, HE-A, IJ, GT, HK, and DK: methodology. JL and SM: software and validation. AAM, JL, HE-A, and SM: formal analysis. AAM, IJ, GT, HK, and DK: investigation. AAM, IJ, GT, HK, DK, and HE-A: resources. AAM, JL, HE-A, HT, SM, and HN: data curation. AAM: writing—original draft preparation. AAM, HE-A, JL, SM, HT, and HN: writing—review and editing, visualization, and supervision. All authors contributed to the article and approved the submitted version.
Funding
This research was funded by the German Federal Ministry of Education and Research (BMBF) and the German Federal Ministry for Economic Affairs and Energy, within the framework of joint project “Adaptable decentralized diagnostics for veterinary and human medicine (ADA)”; Subproject: Adaptable diagnostics for prudent antibiosis of mastitis in livestock (ADA-M) (Funding No.: 13GW0456D: BMBF).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
The authors thank Anja Hackbart and Peggy Methner at the Institute of Bacterial Infections and Zoonoses, Friedrich-Loeffler-Institut, Jena, Germany and the laboratory operation at the animal health service (TGD Thuringia) for their excellent technical assistance. The authors thank Christin Weber at Center for Applied Research, InfectoGnostics Research Campus Jena e.V.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2023.1216850/full#supplementary-material
ReferencesAlibayovB.Baba-MoussaL.SinaH.ZdeňkováK.DemnerováK. (2014). Staphylococcus aureus mobile genetic elements. Mol. Biol. Rep.41, 5005–5018. doi: 10.1007/s11033-014-3367-3AndrewsS., A quality control tool for high throughput sequence data UK (2018) UK, (2018). Available online at: https://www.bioinformatics.babraham.ac.uk/projects/fastqc/ (Accessed August 01, 2023).AnnamanediM.SheelaP.SundareshanS.IsloorS.GuptaP.JasmeenP.. (2021). Molecular fingerprinting of bovine mastitis-associated Staphylococcus aureus isolates from India. Sci. Rep.11:15228. doi: 10.1038/s41598-021-94760-x, PMID: 34315981AntiabongJ. F.KockM. M.BelleaN. M.EhlersM. M. (2017). Diversity of multidrug efflux genes and phenotypic evaluation of the in vitro resistance dynamics of clinical Staphylococcus aureus isolates using methicillin; a model β-lactam. Open Microbiol J11, 132–141. doi: 10.2174/1874285801711010132, PMID: 28839493BankevichA.NurkS.AntipovD.GurevichA. A.DvorkinM.KulikovA. S.. (2012). SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol.19, 455–477. doi: 10.1089/cmb.2012.0021, PMID: 22506599BarkemaH. W.SchukkenY. H.ZadoksR. N. (2006). Invited review: the role of cow, pathogen, and treatment regimen in the therapeutic success of bovine Staphylococcus aureus mastitis. J. Dairy Sci.89, 1877–1895. doi: 10.3168/jds.S0022-0302(06)72256-1, PMID: 16702252BenkerroumN. (2018). Staphylococcal enterotoxins and enterotoxin-like toxins with special reference to dairy products: an overview. Crit. Rev. Food Sci. Nutr.58, 1943–1970. doi: 10.1080/10408398.2017.1289149BolteJ.ZhangY.WenteN.MahmmodY. S.SvennesenL.KrömkerV. (2020). Comparison of phenotypic and genotypic antimicrobial resistance patterns associated with Staphylococcus aureus mastitis in German and Danish dairy cows. J. Dairy Sci.103, 3554–3564. doi: 10.3168/jds.2019-17765, PMID: 32089308BundestierärztekammerLeitlinien für den sorgfältigen Umgang mit antibakteriell wirksamen Tierarzneimitteln mit Erläuterungen –Überarbeitete Fassung Beilage zum Deutschen Tierärzteblatt 3/2015 Schluetersche Verlag, Hannover (2015). 1–24.CamposB.PickeringA. C.RochaL. S.AguilarA. P.Fabres-KleinM. H.de Oliveira MendesT. A.. (2022). Diversity and pathogenesis of Staphylococcus aureus from bovine mastitis: current understanding and future perspectives. BMC Vet. Res.18:115. doi: 10.1186/s12917-022-03197-5, PMID: 35331225ChenL.YangJ.YuJ.YaoZ.SunL.ShenY.. (2005). VFDB: a reference database for bacterial virulence factors. Nucleic Acids Res.33, D325–D328. doi: 10.1093/nar/gki008, PMID: 15608208CLSI. Performance standards for antimicrobial susceptibility testing for Bacteria isolated from animals, 4th edition, CLSI supplement VET08. Wayne, PA, USA (2018). 2018.CLSI. Performance standards for antimicrobial susceptibility testing; CLSI supplement(ISBN 978-1-68440-066-9 [print]; ISBN electronic). Clinical and Laboratory Standards Institute, Pennsylvania, USA, (2020). 58–67CunyC.AbdelbaryM.LayerF.WernerG.WitteW. (2015). Prevalence of the immune evasion gene cluster in Staphylococcus aureus CC398. Vet. Microbiol.177, 219–223. doi: 10.1016/j.vetmic.2015.02.031, PMID: 25778546CunyC.ArnoldP.HermesJ.EckmannsT.MehrajJ.SchoenfelderS.. (2017). Occurrence of CFR-mediated multiresistance in staphylococci from veal calves and pigs, from humans at the corresponding farms, and from veterinarians and their family members. Vet. Microbiol.200, 88–94. doi: 10.1016/j.vetmic.2016.04.002, PMID: 27102205EMA. Reflection paper on use of aminoglycosides in animals in the European Union: Development of resistance and 5 impact on human and animal health. European medicines agencyLondon, UK, (2017). Available at: https://www.ema.europa.eu/en/documents/scientific-guideline/reflection-paper-use-aminoglycosides-animals-european-union-development-resistance-impact-human_en.pdf (Accessed August 01, 2023).FeldgardenM.BroverV.HaftD. H.PrasadA. B.SlottaD. J.TolstoyI.. (2019). Validating the AMRFinder tool and resistance gene database by using antimicrobial resistance genotype-phenotype correlations in a collection of isolates. Antimicrob. Agents Chemother.63:e00483. doi: 10.1128/AAC.00483-19, PMID: 31427293FeßlerA.KadlecK.WangY.ZhangW. J.WuC.ShenJ.. (2018). Small antimicrobial resistance plasmids in livestock-associated methicillin-resistant Staphylococcus aureus CC398. Front. Microbiol.9:2063. doi: 10.3389/fmicb.2018.02063, PMID: 30283407FesslerA. T.Olde RiekerinkR. G.RothkampA.KadlecK.SampimonO. C.LamT. J.. (2012). Characterization of methicillin-resistant Staphylococcus aureus CC398 obtained from humans and animals on dairy farms. Vet. Microbiol.160, 77–84. doi: 10.1016/j.vetmic.2012.05.005, PMID: 22655975FesslerA.ScottC.KadlecK.EhrichtR.MoneckeS.SchwarzS. (2010). Characterization of methicillin-resistant Staphylococcus aureus ST398 from cases of bovine mastitis. J. Antimicrob. Chemother.65, 619–625. doi: 10.1093/jac/dkq021FetschA.KraushaarB.KäsbohrerA.HammerlJ. A. (2017). Turkey meat as source of CC9/CC398 methicillin-resistant Staphylococcus aureus in humans?Clin. Infect. Dis.64, 102–103. doi: 10.1093/cid/ciw687 PMID: 27794020FlorensaA. F.KaasR. S.ClausenP.Aytan-AktugD.AarestrupF. M. (2022). ResFinder – an open online resource for identification of antimicrobial resistance genes in next-generation sequencing data and prediction of phenotypes from genotypes. Microb Genom8:748. doi: 10.1099/mgen.0.000748G.V.A. guidelines for antiseptic Milk sampling and guidelines to isolate and identify mastitis pathogens, Guidelines for antiseptic Milk sampling and guidelines to isolate and identify mastitis pathogens, Verlag der Deutschen Veterinärmedizinischen Gesellschaft e.V., Gießen, Germany, (2018)García-ÁlvarezL.HoldenM. T.LindsayH.WebbC. R.BrownD. F.CurranM. D.. (2011). Meticillin-resistant Staphylococcus aureus with a novel mecA homologue in human and bovine populations in the UK and Denmark: a descriptive study. Lancet Infect. Dis.11, 595–603. doi: 10.1016/S1473-3099(11)70126-8, PMID: 21641281GoergeT.LorenzM. B.van AlenS.HübnerN.-O.BeckerK.KöckR. (2017). MRSA colonization and infection among persons with occupational livestock exposure in Europe: prevalence, preventive options and evidence. Vet. Microbiol.200, 6–12. doi: 10.1016/j.vetmic.2015.10.027, PMID: 26658156HansenJ. E.RoncoT.SteggerM.SieberR. N.FertnerM. E.MartinH. L.. (2019). LA-MRSA CC398 in dairy cattle and veal calf farms indicates spillover from pig. Production10:02733. doi: 10.3389/fmicb.2019.02733HauschildT.FesslerA. T.KadlecK.BillerbeckC.SchwarzS. (2012). Detection of the novel VGA(E) gene in methicillin-resistant Staphylococcus aureus CC398 isolates from cattle and poultry. J. Antimicrob. Chemother.67, 503–504. doi: 10.1093/jac/dkr446, PMID: 22013012HoekstraJ.RuttenV.SommelingL.van WervenT.SpaninksM.DuimB.. (2018). High production of LukMF in Staphylococcus aureus field strains is associated with clinical bovine mastitis. Toxins10:200. doi: 10.3390/toxins10050200, PMID: 29762488HoekstraJ.ZomerA. L.RuttenV.BenedictusL.StegemanA.SpaninksM. P.. (2020). Genomic analysis of European bovine Staphylococcus aureus from clinical versus subclinical mastitis. Sci. Rep.10:18172. doi: 10.1038/s41598-020-75179-2, PMID: 33097797JiX.KrügerH.TaoJ.WangY.FeßlerA. T.BaiR.. (2021). Comparative analysis of genomic characteristics, fitness and virulence of MRSA ST398 and ST9 isolated from China and Germany. Emerg Microbes Infect10, 1481–1494. doi: 10.1080/22221751.2021.1951125, PMID: 34210245KadlecK.EntorfM.PetersT. (2019). Occurrence and characteristics of livestock-associated methicillin-resistant Staphylococcus aureus in quarter milk samples from dairy cows in Germany. Front. Microbiol.10:1295. doi: 10.3389/fmicb.2019.01295, PMID: 31244807KashifA.McClureJ. A.LakhundiS.PhamM.ChenS.ConlyJ. M.. (2019). Staphylococcus aureus ST398 virulence is associated with factors carried on prophage ϕSa3. Front. Microbiol.10:2219. doi: 10.3389/fmicb.2019.02219, PMID: 31608039KraushaarB.HammerlJ. A.KienölM.HeinigM. L.SperlingN.Dinh ThanhM.. (2017). Acquisition of virulence factors in livestock-associated MRSA: lysogenic conversion of CC398 strains by virulence gene-containing phages, vol. 7, 1–13.KreausukonK., Usage of antimicrobials on 60 dairy farms in northern Germany and characterization of methicillin-resistant Staphylococcus aureus (MRSA) and extended spectrum beta-lactamases producing Escherichia coli (ESBLs-producing E. coli) isolated from bulk tank milk samples, Free University of Berlin, (2011) Mensch und Buch VerlagKrömkerV.LeimbachS. (2017). Mastitis treatment-reduction in antibiotic usage in dairy cows. Reproduction in domestic animals. Zuchthygiene52, 21–29. doi: 10.1111/rda.13032, PMID: 28815847LandinH.CarlssonJ.KatholmJ.KulkasL.Plym-ForshellK.RattenborgE.. (2011). Nordic guidelines for mastitis therapyThe Swedish Dairy Association, Lund, Sweden.LeclercqR. (2002). Mechanisms of resistance to macrolides and lincosamides: nature of the resistance elements and their clinical implications. Clin. Infect. Dis.34, 482–492. doi: 10.1086/324626, PMID: 11797175LienenT.SchnittA.HammerlJ. A.MaurischatS.TenhagenB. A. (2020). Genomic distinctions of LA-MRSA ST398 on dairy farms from different German Federal States with a low risk of severe human infections. Front. Microbiol.11:575321. doi: 10.3389/fmicb.2020.575321ManukumarH. M.UmeshaS. (2017). MALDI-TOF-MS based identification and molecular characterization of food associated methicillin-resistant Staphylococcus aureus. Sci. Rep.7:11414. doi: 10.1038/s41598-017-11597-zMinstK.MärtlbauerE.MillerT.MeyerC. (2012). Short communication: Streptococcus species isolated from mastitis milk samples in Germany and their resistance to antimicrobial agents. J. Dairy Sci.95, 6957–6962. doi: 10.3168/jds.2012-5852, PMID: 22999286MoneckeS.Gavier-WidénD.HotzelH.PetersM.GuentherS.LazarisA.. (2016). Diversity of Staphylococcus aureus isolates in European wildlife. PLoS One11:e0168433. doi: 10.1371/journal.pone.0168433, PMID: 27992523MoneckeS.SlickersP.GawlikD.MüllerE.ReissigA.Ruppelt-LorzA.. (2018). Variability of SCCmec elements in livestock-associated CC398 MRSA. Vet. Microbiol.217, 36–46. doi: 10.1016/j.vetmic.2018.02.024, PMID: 29615254MullarkyI. K.SuC.FriezeN.ParkY. H.SordilloL. M. (2001). Staphylococcus aureus AGR genotypes with enterotoxin production capabilities can resist neutrophil bactericidal activity. Infect. Immun.69, 45–51. doi: 10.1128/IAI.69.1.45-51.2001, PMID: 11119487MunitaJ. M.AriasC. A. (2016). Mechanisms of antibiotic resistance. Microbiol Spectr4, 1–24. doi: 10.1128/microbiolspec.VMBF-0016-2015NMC (2017). Current Concepts of Bovine Mastitis. Fifth Edition, New Prague, USA: National Mastitis Council.PapkouA.HedgeJ.KapelN.YoungB.MacLeanR. C. (2020). Efflux pump activity potentiates the evolution of antibiotic resistance across S. aureus isolates. Nat. Commun.11:3970. doi: 10.1038/s41467-020-17735-y, PMID: 32769975PreineF.HerreraD.ScherpenzeelC.KalmusP.McCoyF.SmulskiS.. (2022). Different European perspectives on the treatment of clinical mastitis in lactation. Antibiotics11:1107. doi: 10.3390/antibiotics11081107SayersE. W.BeckJ.BristerJ. R.BoltonE. E.CaneseK.ComeauD. C.. (2020). Database resources of the National Center for biotechnology information. Nucleic Acids Res.48, D9–D16. doi: 10.1093/nar/gkz899, PMID: 31602479SchlotterK.EhrichtR.HotzelH.MoneckeS.PfefferM.DonatK. (2012). Leukocidin genes lukF-P83 and lukM are associated with Staphylococcus aureus clonal complexes 151, 479 and 133 isolated from bovine udder infections in Thuringia Germany. Vet. Res.43:42. doi: 10.1186/1297-9716-43-42, PMID: 22587484SchlotterK.Huber-SchlenstedtR.GanglA.HotzelH.MoneckeS.MüllerE.. (2014). Multiple cases of methicillin-resistant CC130 Staphylococcus aureus harboring MECC in milk and swab samples from a Bavarian dairy herd. J. Dairy Sci.97, 2782–2788. doi: 10.3168/jds.2013-7378, PMID: 24630668SchmidtT.KockM. M.EhlersM. M. (2017). Molecular characterization of Staphylococcus aureus isolated from bovine mastitis and close human contacts in south African dairy herds: genetic diversity and inter-species host transmission. Front. Microbiol.8:511. doi: 10.3389/fmicb.2017.00511SchnittA.LienenT.Wichmann-SchauerH.CunyC.TenhagenB. A. (2020). The occurrence and distribution of livestock-associated methicillin-resistant Staphylococcus aureus ST398 on German dairy farms. J. Dairy Sci.103, 11806–11819. doi: 10.3168/jds.2020-18958, PMID: 33041041SchnittA.TenhagenB. A. (2020). Risk factors for the occurrence of methicillin-resistant Staphylococcus aureus in dairy herds: an update. Foodborne Pathog. Dis.17, 585–596. doi: 10.1089/fpd.2019.2638, PMID: 31433237SchwendenerS.PerretenV. (2011). New transposon Tn6133 in methicillin-resistant Staphylococcus aureus ST398 contains vga(E), a novel streptogramin a, pleuromutilin, and lincosamide resistance gene. Antimicrob. Agents Chemother.55, 4900–4904. doi: 10.1128/AAC.00528-11, PMID: 21768510SchwengersO.BarthP.FalgenhauerL.HainT.ChakrabortyT.GoesmannA. (2020). Platon: identification and characterization of bacterial plasmid contigs in short-read draft assemblies exploiting protein sequence-based replicon distribution scores. Microb Genom6:398. doi: 10.1099/mgen.0.000398SpohrM.RauJ.FriedrichA.KlittichG.FetschA.GuerraB.. (2011). Methicillin-resistant Staphylococcus aureus (MRSA) in three dairy herds in Southwest Germany. Zoonoses Public Health58, 252–261. doi: 10.1111/j.1863-2378.2010.01344.x, PMID: 20630047StamatakisA. (2014). RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics30, 1312–1313. doi: 10.1093/bioinformatics/btu033, PMID: 24451623TÄHAV, Verordnung über tierärztliche Hausapotheken. Available at: https://www.gesetze-im-internet.de/t_hav/BJNR021150975.html (Accessed November 25, 2022). (2018).TenhagenB. A.AltK.PfefferkornB.WiehleL.KäsbohrerA.FetschA. (2018). Short communication: methicillin-resistant Staphylococcus aureus in conventional and organic dairy herds in Germany. J. Dairy Sci.101, 3380–3386. doi: 10.3168/jds.2017-12939, PMID: 29395146TenhagenB. A.KösterG.WallmannJ.HeuwieserW. (2006). Prevalence of mastitis pathogens and their resistance against antimicrobial agents in dairy cows in Brandenburg. J. Dairy Sci.89, 2542–2551. doi: 10.3168/jds.S0022-0302(06)72330-X, PMID: 16772573TenhagenB. A.VossenkuhlB.KäsbohrerA.AltK.KraushaarB.GuerraB.. (2014). Methicillin-resistant Staphylococcus aureus in cattle food chains - prevalence, diversity, and antimicrobial resistance in Germany. J. Anim. Sci.92, 2741–2751. doi: 10.2527/jas.2014-7665, PMID: 24778337VelascoV.SherwoodJ. S.Rojas-GarcíaP. P.LogueC. M. (2014). Multiplex real-time PCR for detection of Staphylococcus aureus, mecA and Panton-Valentine Leukocidin (PVL) genes from selective enrichments from animals and retail meat. PLoS One9:e97617. doi: 10.1371/journal.pone.0097617, PMID: 24849624WangH. Y.KimS.KimJ.ParkS. D.UhY.LeeH. (2014). Multiplex real-time PCR assay for rapid detection of methicillin-resistant staphylococci directly from positive blood cultures. J. Clin. Microbiol.52, 1911–1920. doi: 10.1128/JCM.00389-14, PMID: 24648566WendlandtS.FeßlerA. T.MoneckeS.EhrichtR.SchwarzS.KadlecK. (2013). The diversity of antimicrobial resistance genes among staphylococci of animal origin. Int. J. Med. Microbiol.303, 338–349. doi: 10.1016/j.ijmm.2013.02.006, PMID: 23499306WoodD. E.LuJ.LangmeadB. (2019). Improved metagenomic analysis with kraken 2. Genome Biol.20:257. doi: 10.1186/s13059-019-1891-0, PMID: 31779668YuJ.TienN.LiuY. C.ChoD. Y.ChenJ. W.TsaiY. T.. (2022). Rapid identification of methicillin-resistant Staphylococcus aureus using MALDI-TOF MS and machine learning from over 20,000 clinical isolates. Microbiol. Spectr.10:e0048322. doi: 10.1128/spectrum.00483-22, PMID: 35293803ZieschM.WenteN.ZhangY.ZarembaW.EnglS.KrömkerV. (2018). Noninferiority trial investigating the efficacy of a nonantibiotic intramammary therapy in the treatment of mild-to-moderate clinical mastitis in dairy cows with longer lasting udder diseases. J. Vet. Pharmacol. Ther.41, 11–21. doi: 10.1111/jvp.12415, PMID: 28449183Zimmermann, Beziehungen zwischen der Innerherdenprävalenz subklinischer Streptokokkenmastitiden und der Haltungs- und Melkhygiene in Thüringer Milcherzeugerbetrieben, (2017). Available at: https://nbn-resolving.org/urn:nbn:de:bsz:15-qucosa-224871 (Accessed August 01, 2023).