Next Article in Journal
Severe West Nile Virus Neuroinvasive Disease: Clinical Characteristics, Short- and Long-Term Outcomes
Next Article in Special Issue
Canine Staphylococcus argenteus: Case Report from The Netherlands
Previous Article in Journal
Metabolic Capacity Differentiates Plenodomus lingam from P. biglobosus Subclade ‘brassicae’, the Causal Agents of Phoma Leaf Spotting and Stem Canker of Oilseed Rape (Brassica napus) in Agricultural Ecosystems
Previous Article in Special Issue
Bacteria Patterns on Tonsillar Surface and Tonsillar Core Tissue among Patients Scheduled for Tonsillectomy at Bugando Medical Centre, Mwanza, Tanzania
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Staphylococcus ratti sp. nov. Isolated from a Lab Rat

1
Department of Experimental Biology, Division of Genetics and Molecular Biology, Faculty of Science, Masaryk University, Kotlářská 2, 611 37 Brno, Czech Republic
2
Department of Experimental Biology, Czech Collection of Microorganisms, Faculty of Science, Masaryk University, Kamenice 5, 625 00 Brno, Czech Republic
3
Reference Laboratory for Staphylococci, National Institute of Public Health, Šrobárova 48, 100 42 Praha, Czech Republic
4
Leibniz Institute DSMZ-German Collection of Microorganisms and Cell Cultures, Inhoffenstraße 7 B, 38124 Braunschweig, Germany
5
Department of Public Health, Faculty of Medicine, Masaryk University, Kamenice 5, 625 00 Brno, Czech Republic
*
Author to whom correspondence should be addressed.
Pathogens 2022, 11(1), 51; https://doi.org/10.3390/pathogens11010051
Submission received: 30 November 2021 / Revised: 27 December 2021 / Accepted: 30 December 2021 / Published: 1 January 2022
(This article belongs to the Special Issue Staphylococcus Infections in Humans and Animals)

Abstract

:
Staphylococci from the Staphylococcus intermedius-Staphylococcus hyicus species group include numerous animal pathogens and are an important reservoir of virulence and antimicrobial resistance determinants. Due to their pathogenic potential, they are possible causative agents of zoonoses in humans; therefore, it is important to address the properties of these strains. Here we used a polyphasic taxonomic approach to characterize the coagulase-negative staphylococcal strain NRL/St 03/464T, isolated from the nostrils of a healthy laboratory rat during a microbiological screening of laboratory animals. The 16S rRNA sequence, MALDI-TOF mass spectrometry and positive urea hydrolysis and beta-glucuronidase tests clearly distinguished it from closely related Staphylococcus spp. All analyses have consistently shown that the closest relative is Staphylococcus chromogenes; however, values of digital DNA-DNA hybridization <35.3% and an average nucleotide identity <81.4% confirmed that the analyzed strain is a distinct Staphylococcus species. Whole-genome sequencing and expert annotation of the genome revealed the presence of novel variable genetic elements, including two plasmids named pSR9025A and pSR9025B, prophages, genomic islands and a composite transposon that may confer selective advantages to other bacteria and enhance their survival. Based on phenotypic, phylogenetic and genomic data obtained in this study, the strain NRL/St 03/464T (= CCM 9025T = LMG 31873T = DSM 111348T) represents a novel species with the suggested name Staphylococcus ratti sp. nov.

1. Introduction

Staphylococci are opportunistic pathogens widespread in nature; they are mainly isolated from the skin, skin glands, and mucous membranes of various animals [1], and less often from the blood of diseased animals [2]. The nasal cavity of many mammalian species is inhabited by distinctive staphylococcal species. Studies of Staphylococcus spp. distribution in wild rodents have shown the presence of predominantly coagulase-negative staphylococci, including Staphylococcus xylosus, Staphylococcus equorum, Staphylococcus succinus, Staphylococcus saprophyticus and Mammaliicoccus spp. [3,4]. The species composition of staphylococcal populations in wild and laboratory rats is largely unknown. The predominant species in laboratory rats is S. xylosus, followed by Staphylococcus aureus and Staphylococcus cohnii [5]. Recent studies have shown that rats are reservoir of various S. aureus clonal lineages, including methicillin-resistant strains [6]. On the other hand, murine S. aureus isolates exhibit host adaptation [7]. Both mice and rats are an animal infection model for studying S. aureus pathogenicity [8,9]. In wound infections, rats eliminated bacteria faster and more rapidly organized the inflammatory response than other rodents [10].
Similarly to the main pathogen S. aureus, some species in the Staphylococcus intermedius-Staphylococcus hyicus phylogenetic complex can also catalyze the polymerization of fibrinogen to fibrin through the enzyme coagulase, thereby contributing to the appearance of purulent foci—abscesses [11]. Besides coagulase, they also have other virulence factors, such as staphylococcal β-hemolysin (sphingomyelinase), a typical product of the species Staphylococcus pseudintermedius [12]. Staphylococcal species belonging to this complex are opportunistic pathogens of various animal species. S. pseudintermedius, Staphylococcus coagulans, and Staphylococcus canis occur as carriers in dogs and other carnivores [13,14,15,16], and when skin integrity is disrupted, they can cause skin and soft tissue infections or external ear otitis. Staphylococcus ursi was isolated from healthy black bears [17]. Staphylococcus felis is a common cat commensal and a potential urinary pathogen [18]. However, they can also be zoonotic pathogens, causing skin problems in humans; a typical example is a septic wound after a dog bite or cat scratch [19,20].
Staphylococcus hyicus forms a separate phylogenetic clade together with Staphylococcus agnetis and Staphylococcus chromogenes. S. hyicus causes exudative epidermitis in pigs, known as greasy pig disease, due to the expression of exfoliative toxins that selectively digest porcine desmoglein 1 [21]. Zoonotic infections of the bloodstream have been reported to be caused by S. hyicus, as has spondylodiscitis in animal keepers [22,23]. S. agnetis was primarily isolated from cows with mastitis [24], but recent findings demonstrate that S. agnetis may be associated with diseases and mortality in broiler chickens [25,26]. S. chromogenes is globally the most common cause of intramammary infections of dairy cows [27] and has also been reported in both healthy and diseased pigs [28], goats [29] and poultry [5]. There are also reports of nasal carriers of S. chromogenes among farm workers [30].
The prevalence of human infections caused by species of the S. intermedius/S. hyicus complex is low; however, due to their phenotypic similarity, the capture of some species may be underestimated in studies that did not use molecular techniques. In this work we provide a polyphasic characterization of an isolate of Staphylococcus sp. NRL/St 03/464T that was cultured during a microbiological examination of nasal swabs from healthy rats before a biological experiment at the Unit for Biomedicine and Welfare of Laboratory Animals, National Institute of Public Health (Prague, Czech Republic).

2. Results and Discussion

2.1. Phylogenetic Analyses

The isolate CCM 9025T (= NRL/St 03/464T) was assigned by a partial sequencing of the 16S rRNA gene to the Staphylococcus genus and previously defined Hyicus-Intermedius species group [31]; however, its biochemical profile, MALDI-TOF MS pattern and partial RNA polymerase subunit beta (rpoB) gene sequence did not allow for its classification into any known staphylococcal species. Therefore, a polyphasic taxonomic study was conducted which was focused on a detailed characterization of a new Staphylococcus species.
The obtained complete 16S rRNA gene sequence of strain CCM 9025T was compared to those of other taxa from the Hyicus-Intermedius species group of the Staphylococcus genus. The closest relatives were S. chromogenes (99.3% similarity), S. agnetis (99.2%) and S. hyicus (99.1%); other species were below 98% similarity. The topology of the maximum likelihood phylogenetic tree constructed with 16S rRNA gene sequences was similar to that of the neighbor-joining tree (Figure 1a).
Because the 16S rRNA analysis has limited discriminatory power for the identification of some staphylococcal species, the phylogenetic position of the new isolates was also assessed using the concatenated multilocus sequence data of six routinely used housekeeping genes (Figure 1b). The maximum likelihood phylogenetic trees for the housekeeping genes including the closest relatives had a very similar topology which corresponded to that of the 16S rRNA gene tree, clearly separated the novel isolate CCM 9025T from the established species and which confirmed S. chromogenes as the closest relative. The significant phylogenetic distance from the related staphylococcal species at the whole-genome level, with an average nucleotide identity (ANI) of <81.4% and digital DNA-DNA hybridization (dDDH) of <35.3%, were below the species delineation thresholds, which are 95–96% and 70%, respectively [32]. This confirmed that strain CCM 9025T represents a distinct Staphylococcus species named Staphylococcus ratti sp. nov. It forms a separate branch in the S. hyicus phylogenetic clade, as also shown by protein coding core genome analysis using the up-to-date bacterial core gene (UBCG) (Figure 1c,d).

2.2. Phenotypic, Genotypic and Chemotaxonomic Characteristics

Cells of strain CCM 9025T are irregular spherical cocci with diameter 705 ± 55 nm (n = 100) (Figure 2). The strain grew well on common media for staphylococci and was mesophilic and moderately halophilic, with the ability to hydrolyze biomacromolecules (gelatin, DNA). Carbohydrates seldom served as the source of carbon. The detailed phenotype data are subsequently mentioned in the species description in the text; here we only specify several notable results. Interestingly, the hyaluronidase test result was positive for CCM 9025T, which is a rare feature for coagulase-negative staphylococcal species. During biotyping, a few contradictory results depending on the tested conditions were found. Firstly, the Voges–Proskauer test (acetoin) was negative in a standard tube test, but positive with a commercial VPtest strip containing pyruvic acid instead of glucose. The second inconsistency was related to an enzymatic β-glucuronidase test included in the commercial kits, when β-glucuronidase was positive for CCM 9025T in STAPHYtest 24, but negative in API ZYM, probably due to different substrates being used for enzyme detection.
The tests distinguishing novel species from closely related staphylococci of the S. hyicus group are shown in Table 1, and these tests enable the correct identification of CCM 9025T at the species level.
Profiling analysis by matrix-assisted laser desorption/ionization—time-of-flight mass spectrometry (MALDI-TOF MS) is routinely used for bacterial identification. A distinct MALDI-TOF MS pattern is therefore a useful phenotypic feature in describing new species. Strain CCM 9025T generated a consistent MALDI-TOF MS profile containing signals in the mass range of 2–10 kDa, which was not related to any of those Staphylococcus species already represented in the commercial database at the time of testing. After the manual inclusion of profiles from 24 analyses of the strain CCM 9025T to the in-house database, the strain was re-analyzed and matched the new pattern with a score of 2.6, whereas the next closest species was S. schleiferi with a score of 1.2, which is far below the species identification threshold.
Staphylococci are easily differentiated from other Gram-positive cocci at the genus level by cellular fatty acid analysis, which is also one of the recommended tests at the species level. The cellular fatty acid profile of strain CCM 9025T revealed 4 major fatty acids (FAs) (≥10%), iso-C15:0 (35.0%), anteiso-C15:0 (24.1%), iso-C17:0 (10.8%) and anteiso-C17:0 (10.2%), similar to other validly named Staphylococcus spp. [33]. Comparison to the closest related Staphylococcus spp. showed qualitatively similar profiles of FAs with quantitative differences between the compared type strains (Table 2). Qualitative differences in iso-C15:0, iso-C17:0 and C20:0 clearly distinguish strain CCM 9025T from S. chromogenes CCM 3387T, and lower amounts of anteiso-C17:0 are specific for S. agnetis CCM 8869T. S. hyicus CCM 2368T has the most similar FA profile to CCM 9025T, but can be distinguished by lower amounts of branched C17:0 FAs and a higher amount of iso-C15:0.
The major respiratory quinone in strain CCM 9025T was MK-7 (95.7%). Menaquinones MK-6 and MK-8 were also detected as minor components of the electron transport system. The identification of MK-7 as the major component of the quinone system is in accordance with the genus characteristics, as members of the genus Staphylococcus reveal the presence of unsaturated menaquinones, typically with six, seven or eight isoprene units [33,34]. Analysis of the peptidoglycan structure revealed a cross-linkage type A structure of A3α L-Lys-Gly3-4 similar to the type A11.2 structure [35], with the molar amino acid ratio 2.0 Ala:3.4 Gly:1.0 Glu:0.8 Lys. The type A11.2 peptidoglycan structure was also identified in the closest related species S. chromogenes, S. agnetis and S. hyicus, as described by Schumann [35]. Unlike CCM 9025T and S. agnetis, S. chromogenes and S. hyicus were found to also contain minor amounts of serine, likely substituting for some glycine in the interpeptide bridge [24,33,36].
A DNA fingerprinting technique using repetitive sequence-based PCR (rep-PCR) fingerprinting with the (GTG)5 primer, previously shown to be suitable for the simultaneous detection and differentiation of Staphylococcus spp., was used to demonstrate the difference of strain CCM 9025T from related taxa. The rep-PCR clearly distinguished the analyzed strain from the type strains representing the phylogenetically close Staphylococcus spp. (Figure 3).

2.3. Whole Genome Characterization of Staphylococcus ratti sp. nov.

The genome of S. ratti sp. nov. type strain CCM 9025T was sequenced using Illumina and Oxford Nanopore platforms. The size of the complete chromosome assembly is 2.3 Mb with a mean coverage of 500-fold. Based on the NCBI automated annotation pipeline, a total of 2198 CDSs were identified in the genome, of which 2150 were protein-encoding genes. A total of 82 genes for RNAs were identified in the genome, including 59 tRNAs, 19 rRNAs including 7 (5S), 6 (16S) and 6 (23S), and 4 ncRNAs. Two plasmid sequences named pSR9025A (3311 bp) and pSR9025B (2455 bp) were assembled and annotated as separate extrachromosomal replicons. Plasmid pSR9025A encodes the gene for the Rep protein, which shares 100% amino acid identity with the Rep protein gene in p908 of S. agnetis [37] and a gsiB (glucose starvation-inducible protein B) gene homologue which is involved in a stress response. Further short similar regions were identified in the plasmids of many other coagulase-negative staphylococci. The plasmid pSR9025B is a cryptic plasmid similar (75.5% identity and 27% coverage) to the pLNU9 plasmid of S. chromogenes [38]. This is an indication of interspecies transfer of these types of variable genetic elements (VGEs).
Comparative genomic analysis of S. ratti with the two type strains S. chromogenes NCTC 10530T and S. hyicus ATCC 11249T revealed the presence of additional VGEs, including insertion sequence elements, a composite transposon, one prophage, and a genomic island (Figure 4).
Type I restriction-modification (RM) system genes were found downstream of orfX, but no evidence of staphylococcal cassette chromosome (SCC) integration was found. Type II RM system genes are localized near a cap operon. The clustered, regularly interspaced, short, palindromic repeats (CRISPR)/CRISPR-associated gene (Cas) system was identified in the same location in both the S. ratti and S. hyicus genomes. The genes for predicted virulence factors, surface and extracellular proteins found in S. ratti CCM 9025T genome are shown in Table 3.
The accessory genome is often associated with virulence and antimicrobial resistance and has an important role in the ability of species to colonize particular hosts or persist in the environment. To determine the genomic diversity within all Hyicus-Intermedius species group representatives, the pangenome was analyzed. Type strains within the group have shown an extensive accessory genome whose profile correlates well with the phylogenetic relationship of individual species (Figure 5). The comparison indicates a distinct accessory genes group which is shared with S. hyicus, S. chromogenes and S. agnetis species (Figure 5).
The genome harbors various types of insertion sequences (ISs) from the IS200/IS605, IS3 and IS6 families. In addition to the above ISs, a 9.1-kb composite transposon with 28.3% G+C content is integrated into the genome and flanked by two elements from the IS3 family (Figure 4). The transposon encodes a gene for the radical S-adenosyl-L-methionine (SAM) enzyme (LN051_08115), which is involved in a number of metabolic processes, including post-transcriptional and post-translational modifications, and a gene encoding a YcaO–like protein (LN051_08100) which is responsible for the synthesis of thiazole/oxazole-modified microcin antibiotics [39]. The SagB/ThcOx family dehydrogenase gene (LN051_08090) for a membrane-associated N-acetylglucosaminidase that cleaves polymerized glycan strands to their physiological length and a major facilitator superfamily (MFS) transporter gene (LN051_08095), are also part of the transposon. Since these genes are linked to modulating antibiotic resistance in methicillin-resistant S. aureus [40], we hypothesize that these genes may be responsible for the penicillin resistance of this strain.
One prophage designated vB_SraS_LR1 with a typical siphoviral modular structure is integrated at the 13-bp-long putative att site 5′ AAAATCAACYTTT 3′ adjacent to the tRNAArg gene (locus tag LN051_09185) and exhibits 75.9% identity and 43% coverage with the S. hyicus phage EW (= RG = NCTC 9856) [41], which was previously misclassified according to its genomic sequence [42] as an S. aureus bacteriophage belonging to the Phietavirus genus [43]. Moderate similarity was found (77.3% identity and 30% coverage) to the S. hyicus phage PMBT9 from the Siphoviridae family [44].
A 15-kb long phage-inducible chromosomal island designated SrRICCM9025 was identified in the genome of CCM 9025T (Figure 4). It has 30.1% G+C content and harbours a site-specific integrase (LN051_02785), IS6 family transposase (LN051_02790), the gene virE encoding virulence-associated E family protein (LN051_02795), and new putative phosphotransferase genes which may be related to antimicrobial resistance. SrRICCM9025 exhibits partial sequence similarity to the Macrococcus island McRImsr [45] in the proteins for integration and replication but otherwise demonstrates a distinct gene structure (Figure 6). It is integrated adjacent to the lacA, lacB, lacD, lacG genes and genes for lactose/cellobiose-specific phosphotransferase system genes, which are required for galactose 6-phosphate isomerase activity, described as part of the S. aureus lactose operon lacABCDFEG [46].
Despite the fact that mobile elements were found in the genome, the identification of the CRISPR/Cas system in the S. ratti CCM 9025T genome is consistent with the need to limit the uptake of foreign DNA. Strain CCM 9025T carries a 431-bp-long CRISPR loci with 6 spacers flanked by 36-bp-long direct repeats (DR: GTTTTAGTACTCTGTAATTTTAGTATAAGGTATTC) and putative cas genes encoding CRISPR-associated endonuclease Cas9 and the proteins Cas1, Cas2 and Csn1 typical for CRISPR system type II-A. Cas9 exhibited 79% amino acid identity to the type II CRISPR-associated Cas9 of S. agnetis (WP_107390356) [47] and S. chromogenes (WP_145399953), and 74% amino acid identity to Cas9 of S. hyicus ATCC 11249T (WP_167696241) [48]. Spacers target bacteriophage-related sequences, but no significant similarity to staphylococcal phage genomes was found, indicating a gap in the knowledge of phages infecting this host.

2.4. Description of Staphylococcus ratti sp. nov.

Staphylococcus ratti (rat’ti. L. gen. n. ratti of the rat) cells are Gram-stain-positive spherical cocci occurring predominantly in pairs and clusters. They are non-spore-forming and non-motile. Colonies on tryptic soy agar (TSA) are circular with a whole margin and are flat, smooth, shiny, 1–2 mm in diameter, aerobic and white. They demonstrate positive hemolytic activity on blood agar. Growth occurs in the presence of 15% NaCl at 20 °C and 45 °C, but not at 15 °C or 48 °C. Cells contain catalase, arginine dihydrolase, urease, β-glucuronidase, and hyaluronidase, and demonstrate nitrate reduction and hydrolysis of gelatin and DNA. Hydrolysis of esculin and Tween 80 is negative. Cells demonstrate weak anaerobic growth in thioglycolate medium and are coagulase-, thermostable nuclease-, oxidase-, pyrrolidonyl arylamidase-, Voges–Proskauer test (acetoin)- and ornithine decarboxylase-negative. Cells are susceptible to furazolidone (100 µg) and resistant to bacitracin (10 IU). Cells are acid phosphatase-, alkaline phosphatase-, esterase (C4)-, esterase lipase (C8)- and α-chymotrypsin (weak)-positive in the API ZYM kit. Cells are lipase (C14)-, leucine arylamidase-, valine arylamidase-, cystine arylamidase-, trypsin-, naphthol-AS-Bi-phosphohydrolase-, α-galactosidase-, β-galactosidase-, α-glucosidase-, β-glucosidase-, N-acetyl-β-glucosaminidase-, α-mannosidase- and α-fucosidase-negative with the API ZYM kit. Acid is produced from glycerol, ribose, galactose (weak), D-glucose, D-fructose, mannose, N-acetyl glucosamine, lactose, sucrose, and trehalose. Acid is not produced from erythritol, D-arabinose, L-arabinose, D-xylose, L-xylose, adonitol, β-methyl-D-xyloside, sorbose, rhamnose, dulcitol, inositol, mannitol, sorbitol, α-methyl-D-mannoside, α-methyl-D-glucoside, amygdaline, arbutine, salicin, cellobiose, maltose, melibiose, inulin, melezitose, D-raffinose, starch, glycogen, xylitol, β-gentiobiose, D-turanose, D-lyxose, D-tagatose, D-fucose, L-fucose, D-arabitol, L-arabitol, gluconate, 2 keto-gluconate, and 5 keto-gluconate. Cells are susceptible to ceftazidin, cephalothin, erythromycin, gentamicin, chloramphenicol, imipenem, kanamycin, neomycin, novobiocin, oxacillin, trimethoprim, sulphamethoxazole/trimethoprim (cotrimoxazole), tetracycline rifampicin, vancomycin and fusidic acid. Cells are resistant to ampicillin, clindamycin, penicillin G, and polymyxin B. Ciprofloxacin resistance is intermediate. Cells are susceptible to lysostaphin but resistant to lysozyme.
Utilisation of D-trehalose, sucrose, β-methyl-D-glucoside, N-acetyl-D-glucosamine, N-acetyl neuraminic acid, α-D-glucose, D-mannose, D-fructose, glycerol, D-glucose-6-PO4, L-histidine, L-serine, pectin and acetic acid is positive, and utilisation of D-turanose is borderline with Biolog MicroPlate GEN III, protocol A. There is negative utilisation of dextrin, D-maltose, D-cellobiose, gentiobiose, stachyose, D-raffinose, α-D-lactose, D-melibiose, D-salicin, N-acetyl-β-D-mannosamine, N-acetyl-D-galactosamine, D-galactose, 3-methyl glucose, D-fucose, L-fucose, L-rhamnose, inosine, D-sorbitol, D-mannitol, D-arabitol, myo-inositol, D-fructose-6-PO4, D-aspartic acid, D-serine, gelatin, glycyl-L-proline, L-alanine, L-arginine, L-aspartic acid, L-glutamic acid, L-pyroglutamic acid, D-galacturonic acid, D-galactonic acid lactone, D-gluconic acid, D-glucuronic acid, glucuronamide, mucic acid, quinic acid, D-saccharic acid, p-hydroxy phenylacetic acid, methyl pyruvate, D-lactic acid methyl ester, L-lactic acid, citric acid, α-keto glutaric acid, D-malic acid, L-malic acid, bromo-succinic acid, tween 40, γ-amino-butyric acid, α-hydroxy-butyric acid, β-hydroxy-D,L-butyric acid, α-keto butyric acid, acetoacetic acid, propionic acid and formic acid.
The peptidoglycan is of the type A3α L-Lys–Gly3–4, MK-7 is the major menaquinone, and the major fatty acids are iso-C15:0, anteiso-C15:0, iso-C17:0 and anteiso-C17:0.
The type strain is NRL/St 03/464T (= CCM 9025T = LMG 31873T = DSM 111348T), isolated in 2003 from a swab from the nostrils of a healthy laboratory rat. The DNA G+C content of the type strain is 36.1 mol%, calculated from the whole genomic sequence. The GenBank/ENA/DDBJ accession number for the 16S rRNA gene is OL352091. The complete genome sequence is available under GenBank accession number CP086654.

3. Materials and Methods

3.1. Bacterial Strains and Cultivation

The isolate NRL St 03/464T was isolated from a single animal and maintained as glycerol stock at −70 °C until analyzed. Reference strains of the phylogenetic relatives S. agnetis CCM 8869T, S. hyicus CCM 2368T, and S. chromogenes CCM 3387T were obtained from the Czech Collection of Microorganisms (Brno, Czech Republic). All cultivations were performed at 30 °C for 24 h unless stated otherwise in the test specifications. Type strain Staphylococcus ratti NRL/St 03/464T has been deposited as publicly accessible in the Czech Collection of Microorganisms (CCM), the German Collection of Microorganisms and Cell Cultures (DSMZ), and the BCCM/LMG Bacteria collection.

3.2. Phenotypic Characterization

Extensive phenotypic characterization using the commercial kits STAPHYtest 24 (Erba Lachema, Brno, Czechia) and API ZYM (bioMérieux, Marcy l’Etoile, France), phenotypic fingerprinting using the Biolog system with the identification test panel GEN III MicroPlate (Biolog, Hayward, CA, USA), and conventional biochemical, physiological, and growth tests relevant for the genus Staphylococcus were done as described previously [49,50,51]. The antibiotic resistance pattern was tested by the disc diffusion method on Mueller–Hinton agar (Oxoid, Basingstoke, UK). A set of discs (Oxoid) generally used for Gram-positive cocci were applied: ampicillin (10 µg), oxacillin (1 µg), ceftazidime (30 µg), cephalothin (30 µg), ciprofloxacin (5 µg), clindamycin (2 µg), erythromycin (15 µg), gentamicin (10 µg), chloramphenicol (30 µg), imipenem (10 µg), kanamycin (30 µg), neomycin (10 µg), novobiocin (5 µg), penicillin G (1 IU), rifampicin (5 µg), trimethoprim (5 µg), cotrimoxazole (25 µg), tetracycline (30 µg), vancomycin (30 µg), fusidic acid (10 µg) and polymyxin B (300 U). EUCAST standards and manufacturer’s recommendations (Oxoid) were strictly followed for cultivation and inhibition zone diameter measurement [52].

3.3. Transmission Electron Microscopy

A 200-mesh carbon/formvar-coated grid was placed on a drop of suspension of bacteria in water for 20 min. Bacterial cells located on the grid were negatively stained with 2% ammonium molybdate and treated with UV light. A Morgagni 268D Philips (ThermoFisher Scientific, Amsterdam, The Netherlands) transmission electron microscope was used to visualize bacterial cells.

3.4. Matrix-Assisted Laser Desorption/Ionization Time-of-Flight Mass Spectrometry (MALDI-TOF MS) Analysis

MALDI-TOF MS was performed in automatic acquisition mode as described previously [53] in a Microflex LT MALDI-TOF spectrometer (Bruker Daltonics, Bremen, Germany) by using MBT Compass 4.1 software (Bruker Daltonics). Identification was made using the Bruker’s database MBT Compass Library Revision L 2020 (9607 MSP).

3.5. Chemotaxonomic Analyses

Cellular fatty acids (FAs) were extracted from all compared Staphylococcus strains grown on the Trypticase soy broth agar (TSBA) plates under the same cultivation conditions with a cultivation temperature of 37 °C for 24 h to reach the late-exponential stage of growth according to the four-quadrant streak method. The extraction of cellular FAs was performed according to the standard protocol recommended by the MIDI Microbial Identification System [54]. Extracted FAs were identified using an Agilent 7890B gas chromatograph (Agilent Technologies, Santa Clara, CA, USA) according to the Standard Protocol of the Sherlock Identification System (MIDI Sherlock version 6.2, MIDI database RTSBA version 6.21).
Isolation and structure analyses of the peptidoglycan were performed according to published protocols and some modifications [35,55,56]. In brief, the amino acid composition of total hydrolysate (4 N HCl at 100 °C for 16 h) of the peptidoglycan was analyzed by gas chromatography/mass spectrometry (protocol 10 [35]). The partial hydrolysate (4 N HCl, 100 °C, 45 min) of the peptidoglycan was analyzed by high-resolution liquid chromatography mass spectrometry (LC-MS) as described in [35,55]. Enantiomeric analysis was performed by liquid chromatography as described recently [56].
Respiratory quinones were extracted and analyzed as described previously [57]. Their identity was confirmed by mass spectrometry as described previously [55].

3.6. Genotypic Analysis by (GTG)5-PCR Fingerprinting

Rep-PCR fingerprinting using the (GTG)5 primer was carried out according to Švec et al. [58]. Numerical analysis of the (GTG)5-PCR fingerprints was done using BioNumerics version 7.6 (Applied Maths, Kortrijk, Belgium).

3.7. Phylogenetic Analysis Based on 16S rRNA, Housekeeping Genes and Core Genome

The 16S rRNA gene sequences were amplified from crude boiled cell extracts and sequenced by Sanger sequencing in the Eurofins MWG Operon sequencing facility (Ebersberg, Germany) with previously described primers [59]. The partial rpoB gene was sequenced as described previously [60]. Initial identification of the strain to the genus level was based on pairwise sequence alignment and calculation of similarity values with the algorithm used in the EzBioCloud database v.2021.07.07 [61]. 16S rRNA gene sequences obtained from PCR products were aligned with those extracted from WGS data using RNAmmer version 1.2 [62]. The multilocus sequence data of six housekeeping genes that are commonly used in phylogenetic studies of the Staphylococcus genus were extracted from whole-genome sequence assemblies of type strains available in the NCBI Assembly resource [63], including the NCTC 3000 project (https://www.sanger.ac.uk/resources/downloads/bacteria/nctc/ accessed on 27 October 2021) and FDA-ARGOS project [64]. The partial gene sequences used correspond to the following gene coordinates of S. aureus: 1420..1974 for rpoB, 270..826 for groEL, 23..911 for dnaJ, 49..929 for gap, 383..1032 for tufA, and 50..480 for the sodA gene. The phylogenetic analyses were performed with the software MEGA X [65]. Genetic distances were corrected using the Tamura–Nei model [66], and the evolutionary history was inferred using the maximum likelihood (ML) and neighbor-joining (NJ) methods using a bootstrap test based on 500 replications [67]. The up-to-date bacterial core gene (UBCG) pipeline version 3.0 was used for whole-genome phylogenetic analysis based on core gene sequences [68].
The ANI and dDDH values were calculated using the web-based genome-to-genome distance calculator (GGDC) version 3.0 [69] and FastANI [70], respectively.

3.8. Genome Sequencing and Bioinformatics Analyses

Total genomic DNA was extracted using a GenElute Bacterial Genomic DNA kit (Sigma-Aldrich, St. Louis, MO, USA) from pure culture colonies cultivated on Colombia sheep blood agar (Oxoid). The preparation of DNA libraries with a Nextera XT DNA Library Preparation Kit (Illumina, San Diego, CA, USA) and whole-genome sequencing on the Illumina platform were conducted externally (LGC Genomics, Berlin, Germany) using 2 × 150 bp paired-end reads on the NextSeq sequencing platform (Illumina).
For sequencing using the Oxford Nanopore platform, bacterial DNA was isolated as described previously [71]. The library was prepared using the SQK-RBK004 rapid barcoding kit (Oxford Nanopore Technologies, Oxford, UK) according to the manufacturer’s instructions. Libraries were sequenced with FLO-MIN106 flow cells (R9.4.1) in a MinION device (Oxford Nanopore Technologies, Oxford, UK). The device was controlled with the software MinKNOW version 4.1.2 (Oxford Nanopore Technologies, Oxford, UK).
Basecalling, demultiplexing and barcode trimming were performed using stand-alone ONT Guppy software version 5.0.11 using the config file dna_r9.4.1_450bps_sup.cfg with the default minimum q-score threshold, i.e., 10. The MinION reads were subsequently filtered by mapping to Illumina reads using Filtlong version 0.2.1 (https://github.com/rrwick/Filtlong accessed on 20 September 2021) with a minimum length of 1500 bp and quality threshold set to 95%. Only data that exceeded these thresholds was used in the assembly. The quality of reads was assessed with FastQC version 0.11.9 (http://www.bioinformatics.babraham.ac.uk/projects/fastqc accessed on 20 September 2021) and NanoStat [72]. Complete bacterial genome sequences were obtained using a hybrid assembly with Unicycler version 0.4.9 [73] using SPAdes version 3.12.0 [74], and the parameters chosen were bold mode and k-mers 21,55,77,99,127. The resulting assembly was polished with Pilon version 1.24 [75].
For pangenome analysis, the complete genomes were initially annotated with Prokka version 1.14.6 [76] and clustered with Roary [77]; the results were then visualized with the script roary_plots.py, which is provided in the Roary package. Further, the genome was annotated using the NCBI Prokaryotic Genome Annotation Pipeline [78]. Sequences were manipulated and inspected in the cross-platform bioinformatics software Ugene version 38.1 [79]. The multiple sequence alignment was visualised using EasyFig version 2.2.5 [80]. Gene content was further examined manually with NCBI BLAST (https://blast.ncbi.nlm.nih.gov accessed on 4 November 2021), and VGEs were identified with PHASTER [81], PhiSpy version 3.4 [82], IslandViewer 4 [83], and ISFinder [84]. The CRISPR/Cas system was characterized with CRISPRCasTyper [85]. Virulence factors were predicted using the VFanalyzer tool available at the Virulence Factors Database [86].

Author Contributions

Conceptualization, R.P., I.S. and V.K.; data analysis, V.K., I.S., I.M. and P.P.; investigation, P.P., S.K., P.Š., T.B., M.N.-S., T.G. and E.S.; writing—original draft preparation, V.K., R.P., I.S., P.P., S.K. and I.M.; writing—review and editing, R.P., I.S., P.Š. and J.D.; visualization V.K. and R.P.; supervision, R.P. and I.M. All authors have read and agreed to the published version of the manuscript.

Funding

The students involved in this study were supported by the Grant Agency of Masaryk University (MUNI/A/1522/2020). T.B. gratefully acknowledges funding by the Ministry of Health of the Czech Republic (NU21J-05-00035). P.P. was supported by the Ministry of Health, Czech Republic—conceptual development of research organization (The National Institute of Public Health—NIPH, 75010330). Computational resources were supplied by the project e-INFRA, supported by the Ministry of Education, Youth and Sports of the Czech Republic (LM2018140).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The complete genome sequence of strain NRL/St 03/464T (= CCM 9025T) has been deposited in GenBank/ENA/DDBJ under accession number CP086654 (chromosome) and CP086655 and CP086656 (plasmids). The accession number for the 16S rRNA gene is OL352091. The associated BioProject and BioSample accession numbers are PRJNA779216 and SAMN22456156, respectively. The raw reads were deposited into the SRA database with accession numbers SRR16917662 (Illumina NextSeq 500) and SRR16917661 (Oxford Nanopore MinION).

Acknowledgments

Daniel Krsek (NRL for Diagnostic Electron Microscopy of Infectious Agents, National Institute of Public Health, Prague, Czech Republic) is gratefully acknowledged for transmission electron microscopy, Anika Wasner for technical assistance and Peter Schumann for helpful discussions (both Leibniz Institute DSMZ).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Götz, F.; Bannerman, T.; Schleifer, K.-H. The Genera Staphylococcus and Macrococcus. In The Prokaryotes; Dworkin, M., Falkow, S., Rosenberg, E., Schleifer, K.-H., Stackebrandt, E., Eds.; Springer: New York, NY, USA, 2006; pp. 5–75. [Google Scholar] [CrossRef]
  2. Nováková, D.; Pantůček, R.; Hubálek, Z.; Falsen, E.; Busse, H.-J.; Schumann, P.; Sedláček, I. Staphylococcus microti sp. nov., isolated from the common vole (Microtus arvalis). Int. J. Syst. Evol. Microbiol. 2010, 60, 566–573. [Google Scholar] [CrossRef] [Green Version]
  3. Hauschild, T. Phenotypic and genotypic identification of staphylococci isolated from wild small mammals. Syst. Appl. Microbiol. 2001, 24, 411–416. [Google Scholar] [CrossRef]
  4. Hauschild, T.; Slizewski, P.; Masiewicz, P. Species distribution of staphylococci from small wild mammals. Syst. Appl. Microbiol. 2010, 33, 457–460. [Google Scholar] [CrossRef] [PubMed]
  5. Shimizu, A.; Ozaki, J.; Kawano, J.; Saitoh, Y.; Kimura, S. Distribution of Staphylococcus species on animal skin. J. Vet. Med. Sci. 1992, 54, 355–357. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Raafat, D.; Mrochen, D.M.; Al’Sholui, F.; Heuser, E.; Ryll, R.; Pritchett-Corning, K.R.; Jacob, J.; Walther, B.; Matuschka, F.R.; Richter, D.; et al. Molecular epidemiology of methicillin-susceptible and methicillin-resistant Staphylococcus aureus in wild, captive and laboratory rats: Effect of habitat on the nasal S. aureus population. Toxins 2020, 12, 80. [Google Scholar] [CrossRef] [Green Version]
  7. Mrochen, D.M.; Grumann, D.; Schulz, D.; Gumz, J.; Trube, P.; Pritchett-Corning, K.; Johnson, S.; Nicklas, W.; Kirsch, P.; Martelet, K.; et al. Global spread of mouse-adapted Staphylococcus aureus lineages CC1, CC15, and CC88 among mouse breeding facilities. Int. J. Med. Microbiol. 2018, 308, 598–606. [Google Scholar] [CrossRef] [PubMed]
  8. Hanses, F.; Roux, C.; Dunman, P.M.; Salzberger, B.; Lee, J.C. Staphylococcus aureus gene expression in a rat model of infective endocarditis. Genome Med. 2014, 6, 93. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Power, M.E.; Olson, M.E.; Domingue, P.A.; Costerton, J.W. A rat model of Staphylococcus aureus chronic osteomyelitis that provides a suitable system for studying the human infection. J. Med. Microbiol. 1990, 33, 189–198. [Google Scholar] [CrossRef]
  10. Donnelly, T.M.; Stark, D.M. Susceptibility of laboratory rats, hamsters, and mice to wound infection with Staphylococcus aureus. Am. J. Vet. Res. 1985, 46, 2634–2638. [Google Scholar]
  11. Devriese, L.A. Staphylococci in healthy and diseased animals. J. Appl. Microbiol. 1990, 69, 71S–80S. [Google Scholar] [CrossRef]
  12. Bannoehr, J.; Guardabassi, L. Staphylococcus pseudintermedius in the dog: Taxonomy, diagnostics, ecology, epidemiology and pathogenicity. Vet. Dermatol. 2012, 23, 253-e52. [Google Scholar] [CrossRef] [PubMed]
  13. Devriese, L.A.; Vancanneyt, M.; Baele, M.; Vaneechoutte, M.; De Graef, E.; Snauwaert, C.; Cleenwerck, I.; Dawyndt, P.; Swings, J.; Decostere, A.; et al. Staphylococcus pseudintermedius sp. nov., a coagulase-positive species from animals. Int. J. Syst. Evol. Microbiol. 2005, 55, 1569–1573. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Igimi, S.; Takahashi, E.; Mitsuoka, T. Staphylococcus schleiferi subsp. coagulans subsp. nov., isolated from the external auditory meatus of dogs with external ear otitis. Int. J. Syst. Bacteriol. 1990, 40, 409–411. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Newstead, L.L.; Harris, J.; Goodbrand, S.; Varjonen, K.; Nuttall, T.; Paterson, G.K. Staphylococcus caledonicus sp. nov. and Staphylococcus canis sp. nov. isolated from healthy domestic dogs. Int. J. Syst. Evol. Microbiol. 2021, 71, 004878. [Google Scholar] [CrossRef] [PubMed]
  16. Vrbovská, V.; Sedláček, I.; Zeman, M.; Švec, P.; Kovařovic, V.; Šedo, O.; Laichmanová, M.; Doskař, J.; Pantůček, R. Characterization of Staphylococcus intermedius group isolates associated with animals from Antarctica and emended description of Staphylococcus delphini. Microorganisms 2020, 8, 204. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Perreten, V.; Kania, S.A.; Bemis, D. Staphylococcus ursi sp. nov., a new member of the ‘Staphylococcus intermedius group’ isolated from healthy black bears. Int. J. Syst. Evol. Microbiol. 2020, 70, 4637–4645. [Google Scholar] [CrossRef]
  18. Worthing, K.; Pang, S.; Trott, D.J.; Abraham, S.; Coombs, G.W.; Jordan, D.; McIntyre, L.; Davies, M.R.; Norris, J. Characterisation of Staphylococcus felis isolated from cats using whole genome sequencing. Vet. Microbiol. 2018, 222, 98–104. [Google Scholar] [CrossRef] [Green Version]
  19. Guardabassi, L.; Loeber, M.E.; Jacobson, A. Transmission of multiple antimicrobial-resistant Staphylococcus intermedius between dogs affected by deep pyoderma and their owners. Vet. Microbiol. 2004, 98, 23–27. [Google Scholar] [CrossRef]
  20. Boerlin, P.; Eugster, S.; Gaschen, F.; Straub, R.; Schawalder, P. Transmission of opportunistic pathogens in a veterinary teaching hospital. Vet. Microbiol. 2001, 82, 347–359. [Google Scholar] [CrossRef]
  21. Fudaba, Y.; Nishifuji, K.; Andresen, L.O.; Yamaguchi, T.; Komatsuzawa, H.; Amagai, M.; Sugai, M. Staphylococcus hyicus exfoliative toxins selectively digest porcine desmoglein 1. Microb. Pathog. 2005, 39, 171–176. [Google Scholar] [CrossRef]
  22. Casanova, C.; Iselin, L.; von Steiger, N.; Droz, S.; Sendi, P. Staphylococcus hyicus bacteremia in a farmer. J. Clin. Microbiol. 2011, 49, 4377–4378. [Google Scholar] [CrossRef] [Green Version]
  23. Foissac, M.; Lekaditi, M.; Loutfi, B.; Ehrhart, A.; Dauchy, F.A. Spondylodiscitis and bacteremia due to Staphylococcus hyicus in an immunocompetent man. Germs 2016, 6, 106–110. [Google Scholar] [CrossRef] [Green Version]
  24. Taponen, S.; Supre, K.; Piessens, V.; Van Coillie, E.; De Vliegher, S.; Koort, J.M.K. Staphylococcus agnetis sp. nov., a coagulase-variable species from bovine subclinical and mild clinical mastitis. Int. J. Syst. Evol. Microbiol. 2012, 62, 61–65. [Google Scholar] [CrossRef] [PubMed]
  25. Poulsen, L.L.; Thofner, I.; Bisgaard, M.; Olsen, R.H.; Christensen, J.P.; Christensen, H. Staphylococcus agnetis, a potential pathogen in broiler breeders. Vet. Microbiol. 2017, 212, 1–6. [Google Scholar] [CrossRef] [PubMed]
  26. Al-Rubaye, A.A.; Couger, M.B.; Ojha, S.; Pummill, J.F.; Koon, J.A., 2nd; Wideman, R.F., Jr.; Rhoads, D.D. Genome analysis of Staphylococcus agnetis, an agent of lameness in broiler chickens. PLoS ONE 2015, 10, e0143336. [Google Scholar] [CrossRef] [PubMed]
  27. Devriese, L.A.; Baele, M.; Vaneechoutte, M.; Martel, A.; Haesebrouck, F. Identification and antimicrobial susceptibility of Staphylococcus chromogenes isolates from intramammary infections of dairy cows. Vet. Microbiol. 2002, 87, 175–182. [Google Scholar] [CrossRef]
  28. Andresen, L.O.; Ahrens, P.; Daugaard, L.; Bille-Hansen, V. Exudative epidermitis in pigs caused by toxigenic Staphylococcus chromogenes. Vet. Microbiol. 2005, 105, 291–300. [Google Scholar] [CrossRef] [PubMed]
  29. Andrews, A.H.; Lamport, A. Isolation of Staphylococcus chromogenes from an unusual case of impetigo in a goat. Vet. Rec. 1997, 140, 584. [Google Scholar] [CrossRef]
  30. Schmidt, T.; Kock, M.M.; Ehlers, M.M. Diversity and antimicrobial susceptibility profiling of staphylococci isolated from bovine mastitis cases and close human contacts. J. Dairy Sci. 2015, 98, 6256–6269. [Google Scholar] [CrossRef]
  31. Lamers, R.P.; Muthukrishnan, G.; Castoe, T.A.; Tafur, S.; Cole, A.M.; Parkinson, C.L. Phylogenetic relationships among Staphylococcus species and refinement of cluster groups based on multilocus data. BMC Evol. Biol. 2012, 12, 171. [Google Scholar] [CrossRef] [Green Version]
  32. Goris, J.; Konstantinidis, K.T.; Klappenbach, J.A.; Coenye, T.; Vandamme, P.; Tiedje, J.M. DNA-DNA hybridization values and their relationship to whole-genome sequence similarities. Int. J. Syst. Evol. Microbiol. 2007, 57, 81–91. [Google Scholar] [CrossRef] [Green Version]
  33. Schleifer, K.-H.; Bell, J.A. Staphylococcus. In Bergey’s Manual of Systematics of Archaea and Bacteria, Online ed.; Whitman, W.B., Ed.; Wiley & Sons and Bergey’s Manual Trust: Hoboken, NJ, USA, 2015. [Google Scholar] [CrossRef]
  34. Nahaie, M.R.; Goodfellow, M.; Minnikin, D.E.; Hájek, V. Polar lipid and isoprenoid quinone composition in the classification of Staphylococcus. J. Gen. Microbiol. 1984, 130, 2427–2437. [Google Scholar] [CrossRef] [Green Version]
  35. Schumann, P. Peptidoglycan Structure. Methods Microbiol. 2011, 38, 101–129. [Google Scholar] [CrossRef]
  36. Devriese, L.A.; Hajek, V.; Oeding, P.; Meyer, S.A.; Schleifer, K.H. Staphylococcus hyicus (Sompolinsky 1953) comb. nov. and Staphylococcus hyicus subsp. chromogenes subsp. nov. Int. J. Syst. Bacteriol. 1978, 28, 482–490. [Google Scholar] [CrossRef] [Green Version]
  37. Shwani, A.; Adkins, P.R.F.; Ekesi, N.S.; Alrubaye, A.; Calcutt, M.J.; Middleton, J.R.; Rhoads, D.D. Whole-genome comparisons of Staphylococcus agnetis isolates from cattle and chickens. Appl. Environ. Microbiol. 2020, 86, e00484-20. [Google Scholar] [CrossRef] [PubMed]
  38. Luthje, P.; von Kockritz-Blickwede, M.; Schwarz, S. Identification and characterization of nine novel types of small staphylococcal plasmids carrying the lincosamide nucleotidyltransferase gene lnu(A). J. Antimicrob. Chemother. 2007, 59, 600–606. [Google Scholar] [CrossRef] [Green Version]
  39. Burkhart, B.J.; Schwalen, C.J.; Mann, G.; Naismith, J.H.; Mitchell, D.A. YcaO-dependent posttranslational amide activation: Biosynthesis, structure, and function. Chem. Rev. 2017, 117, 5389–5456. [Google Scholar] [CrossRef] [PubMed]
  40. Chan, Y.G.; Frankel, M.B.; Missiakas, D.; Schneewind, O. SagB glucosaminidase is a determinant of Staphylococcus aureus glycan chain length, antibiotic susceptibility, and protein secretion. J. Bacteriol. 2016, 198, 1123–1136. [Google Scholar] [CrossRef] [Green Version]
  41. Dean, B.A.; Williams, R.E.O.; Hall, F.; Corse, J. Phage typing of coagulase-negative staphylococci and micrococci. Epidemiol. Infect. (J. Hyg.) 1973, 71, 261–270. [Google Scholar] [CrossRef] [Green Version]
  42. Kwan, T.; Liu, J.; DuBow, M.; Gros, P.; Pelletier, J. The complete genomes and proteomes of 27 Staphylococcus aureus bacteriophages. Proc. Natl. Acad. Sci. USA 2005, 102, 5174–5179. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Gutierrez, D.; Adriaenssens, E.M.; Martinez, B.; Rodriguez, A.; Lavigne, R.; Kropinski, A.M.; Garcia, P. Three proposed new bacteriophage genera of staphylococcal phages: “3alikevirus”, “77likevirus” and “Phietalikevirus”. Arch. Virol. 2014, 159, 389–398. [Google Scholar] [CrossRef] [Green Version]
  44. Tetens, J.; Sprotte, S.; Thimm, G.; Wagner, N.; Brinks, E.; Neve, H.; Holzel, C.S.; Franz, C.M.A.P. First molecular characterization of Siphoviridae-like bacteriophages infecting Staphylococcus hyicus in a case of exudative epidermitis. Front. Microbiol. 2021, 12, 653501. [Google Scholar] [CrossRef]
  45. Schwendener, S.; Dona, V.; Perreten, V. The novel macrolide resistance genes mef(D), msr(F), and msr(H) are present on resistance islands in Macrococcus canis, Macrococcus caseolyticus, and Staphylococcus aureus. Antimicrob. Agents Chemother. 2020, 64, e00160-20. [Google Scholar] [CrossRef]
  46. Rosey, E.L.; Oskouian, B.; Stewart, G.C. Lactose metabolism by Staphylococcus aureus: Characterization of lacABCD, the structural genes of the tagatose 6-phosphate pathway. J. Bacteriol. 1991, 173, 5992–5998. [Google Scholar] [CrossRef] [Green Version]
  47. Naushad, S.; Barkema, H.W.; Luby, C.; Condas, L.A.; Nobrega, D.B.; Carson, D.A.; De Buck, J. Comprehensive phylogenetic analysis of bovine non-aureus staphylococci species based on whole-genome sequencing. Front. Microbiol. 2016, 7, 1990. [Google Scholar] [CrossRef] [Green Version]
  48. Haft, D.H.; Selengut, J.; Mongodin, E.F.; Nelson, K.E. A guild of 45 CRISPR-associated (Cas) protein families and multiple CRISPR/Cas subtypes exist in prokaryotic genomes. PLoS Comput. Biol. 2005, 1, e60. [Google Scholar] [CrossRef]
  49. Pantůček, R.; Švec, P.; Dajcs, J.J.; Machová, I.; Černohlavková, J.; Šedo, O.; Gelbíčová, T.; Mašlaňová, I.; Doškař, J.; Zdráhal, Z.; et al. Staphylococcus petrasii sp. nov. including S. petrasii subsp. petrasii subsp. nov. and S. petrasii subsp. croceilyticus subsp. nov., isolated from human clinical specimens and human ear infections. Syst. Appl. Microbiol. 2013, 36, 90–95. [Google Scholar] [CrossRef] [PubMed]
  50. Mannerová, S.; Pantůček, R.; Doškař, J.; Švec, P.; Snauwaert, C.; Vancanneyt, M.; Swings, J.; Sedláček, I. Macrococcus brunensis sp. nov., Macrococcus hajekii sp. nov. and Macrococcus lamae sp. nov., from the skin of llamas. Int. J. Syst. Evol. Microbiol. 2003, 53, 1647–1654. [Google Scholar] [CrossRef] [PubMed]
  51. Mašlaňová, I.; Wertheimer, Z.; Sedláček, I.; Švec, P.; Indráková, A.; Kovařovic, V.; Schumann, P.; Spröer, C.; Králová, S.; Šedo, O.; et al. Description and comparative genomics of Macrococcus caseolyticus subsp. hominis subsp. nov., Macrococcus goetzii sp. nov., Macrococcus epidermidis sp. nov., and Macrococcus bohemicus sp. nov., novel macrococci from human clinical material with virulence potential and suspected uptake of foreign DNA by natural transformation. Front. Microbiol. 2018, 9, 1178. [Google Scholar] [CrossRef] [PubMed]
  52. EUCAST. Breakpoint Tables for Interpretation of MICs and Zone Diameters. The European Committee on Antimicrobial Susceptibility Testing: Version 6.0. 2016. Available online: https://www.eucast.org (accessed on 10 June 2019).
  53. Pantůček, R.; Sedláček, I.; Indraková, A.; Vrbovská, V.; Mašlaňová, I.; Kovařovic, V.; Švec, P.; Králová, S.; Krištofová, L.; Kekláková, J.; et al. Staphylococcus edaphicus sp. nov., isolated in Antarctica, harbors the mecC gene and genomic islands with a suspected role in adaptation to extreme environments. Appl. Env. Microbiol. 2018, 84, e01746-17. [Google Scholar] [CrossRef] [Green Version]
  54. Sasser, M. Identification of Bacteria by Gas Chromatography of Cellular Fatty Acids, MIDI Technical Note 101, Revision July 2006 ed.; MIDI Inc.: Newark, DE, USA, 1990. [Google Scholar]
  55. Schumann, P.; Kalensee, F.; Cao, J.; Criscuolo, A.; Clermont, D.; Kohler, J.M.; Meier-Kolthoff, J.P.; Neumann-Schaal, M.; Tindall, B.J.; Pukall, R. Reclassification of Haloactinobacterium glacieicola as Occultella glacieicola gen. nov., comb. nov., of Haloactinobacterium album as Ruania alba comb. nov, with an emended description of the genus Ruania, recognition that the genus names Haloactinobacterium and Ruania are heterotypic synonyms and description of Occultella aeris sp. nov., a halotolerant isolate from surface soil sampled at an ancient copper smelter. Int. J. Syst. Evol. Microbiol. 2021, 71, 004769. [Google Scholar] [CrossRef]
  56. Kämpfer, P.; McInroy, J.A.; Clermont, D.; Neumann-Schaal, M.; Criscuolo, A.; Busse, H.-J.; Glaeser, S.P. Leucobacter soli sp. nov., from soil amended with humic acid. Int. J. Syst. Evol. Microbiol. 2021, 71, 005156. [Google Scholar] [CrossRef] [PubMed]
  57. Vieira, S.; Huber, K.J.; Neumann-Schaal, M.; Geppert, A.; Luckner, M.; Wanner, G.; Overmann, J. Usitatibacter rugosus gen. nov., sp. nov. and Usitatibacter palustris sp. nov., novel members of Usitatibacteraceae fam. nov. within the order Nitrosomonadales isolated from soil. Int. J. Syst. Evol. Microbiol. 2021, 71, 004631. [Google Scholar] [CrossRef] [PubMed]
  58. Švec, P.; Pantůček, R.; Petráš, P.; Sedláček, I.; Nováková, D. Identification of Staphylococcus spp. using (GTG)5-PCR fingerprinting. Syst. Appl. Microbiol. 2010, 33, 451–456. [Google Scholar] [CrossRef] [PubMed]
  59. Pantůček, R.; Sedláček, I.; Petráš, P.; Koukalová, D.; Švec, P.; Štětina, V.; Vancanneyt, M.; Chrastinová, L.; Vokurková, J.; Růžičková, V.; et al. Staphylococcus simiae sp. nov., isolated from South American squirrel monkeys. Int. J. Syst. Evol. Microbiol. 2005, 55, 1953–1958. [Google Scholar] [CrossRef]
  60. Mellmann, A.; Becker, K.; von Eiff, C.; Keckevoet, U.; Schumann, P.; Harmsen, D. Sequencing and staphylococci identification. Emerg. Infect. Dis. 2006, 12, 333–336. [Google Scholar] [CrossRef]
  61. Yoon, S.H.; Ha, S.M.; Kwon, S.; Lim, J.; Kim, Y.; Seo, H.; Chun, J. Introducing EzBioCloud: A taxonomically united database of 16S rRNA gene sequences and whole-genome assemblies. Int. J. Syst. Evol. Microbiol. 2017, 67, 1613–1617. [Google Scholar] [CrossRef]
  62. Lagesen, K.; Hallin, P.; Rodland, E.A.; Staerfeldt, H.H.; Rognes, T.; Ussery, D.W. RNAmmer: Consistent and rapid annotation of ribosomal RNA genes. Nucleic Acids Res. 2007, 35, 3100–3108. [Google Scholar] [CrossRef]
  63. Kitts, P.A.; Church, D.M.; Thibaud-Nissen, F.; Choi, J.; Hem, V.; Sapojnikov, V.; Smith, R.G.; Tatusova, T.; Xiang, C.; Zherikov, A.; et al. Assembly: A resource for assembled genomes at NCBI. Nucleic Acids Res. 2016, 44, D73–D80. [Google Scholar] [CrossRef]
  64. Sichtig, H.; Minogue, T.; Yan, Y.; Stefan, C.; Hall, A.; Tallon, L.; Sadzewicz, L.; Nadendla, S.; Klimke, W.; Hatcher, E.; et al. FDA-ARGOS is a database with public quality-controlled reference genomes for diagnostic use and regulatory science. Nat. Commun. 2019, 10, 3313. [Google Scholar] [CrossRef] [Green Version]
  65. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K. MEGA X: Molecular Evolutionary Genetics Analysis across computing platforms. Mol. Biol. Evol. 2018, 35, 1547–1549. [Google Scholar] [CrossRef] [PubMed]
  66. Tamura, K.; Nei, M. Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Mol. Biol. Evol. 1993, 10, 512–526. [Google Scholar] [CrossRef]
  67. Felsenstein, J. Confidence limits on phylogenies: An approach using the bootstrap. Evolution 1985, 39, 783–791. [Google Scholar] [CrossRef]
  68. Na, S.I.; Kim, Y.O.; Yoon, S.H.; Ha, S.M.; Baek, I.; Chun, J. UBCG: Up-to-date bacterial core gene set and pipeline for phylogenomic tree reconstruction. J. Microbiol. 2018, 56, 280–285. [Google Scholar] [CrossRef]
  69. Meier-Kolthoff, J.P.; Carbasse, J.S.; Peinado-Olarte, R.L.; Göker, M. TYGS and LPSN: A database tandem for fast and reliable genome-based classification and nomenclature of prokaryotes. Nucleic Acids Res. 2021, 49, gkab902. [Google Scholar] [CrossRef]
  70. Jain, C.; Rodriguez, R.L.; Phillippy, A.M.; Konstantinidis, K.T.; Aluru, S. High throughput ANI analysis of 90K prokaryotic genomes reveals clear species boundaries. Nat. Commun. 2018, 9, 5114. [Google Scholar] [CrossRef] [Green Version]
  71. Fišarová, L.; Botka, T.; Du, X.; Mašlaňová, I.; Bárdy, P.; Pantůček, R.; Benešík, M.; Roudnický, P.; Winstel, V.; Larsen, J.; et al. Staphylococcus epidermidis phages transduce antimicrobial resistance plasmids and mobilize chromosomal islands. mSphere 2021, 6, e00223-21. [Google Scholar] [CrossRef]
  72. De Coster, W.; D’Hert, S.; Schultz, D.T.; Cruts, M.; Van Broeckhoven, C. NanoPack: Visualizing and processing long-read sequencing data. Bioinformatics 2018, 34, 2666–2669. [Google Scholar] [CrossRef]
  73. Wick, R.R.; Judd, L.M.; Gorrie, C.L.; Holt, K.E. Unicycler: Resolving bacterial genome assemblies from short and long sequencing reads. PLoS Comput. Biol. 2017, 13, e1005595. [Google Scholar] [CrossRef] [Green Version]
  74. 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] [Green Version]
  75. Walker, B.J.; Abeel, T.; Shea, T.; Priest, M.; Abouelliel, A.; Sakthikumar, S.; Cuomo, C.A.; Zeng, Q.; Wortman, J.; Young, S.K.; et al. Pilon: An integrated tool for comprehensive microbial variant detection and genome assembly improvement. PLoS ONE 2014, 9, e112963. [Google Scholar] [CrossRef]
  76. Seemann, T. Prokka: Rapid prokaryotic genome annotation. Bioinformatics 2014, 30, 2068–2069. [Google Scholar] [CrossRef] [PubMed]
  77. Page, A.J.; Cummins, C.A.; Hunt, M.; Wong, V.K.; Reuter, S.; Holden, M.T.; Fookes, M.; Falush, D.; Keane, J.A.; Parkhill, J. Roary: Rapid large-scale prokaryote pan genome analysis. Bioinformatics 2015, 31, 3691–3693. [Google Scholar] [CrossRef] [PubMed]
  78. Li, W.; O’Neill, K.R.; Haft, D.H.; DiCuccio, M.; Chetvernin, V.; Badretdin, A.; Coulouris, G.; Chitsaz, F.; Derbyshire, M.K.; Durkin, A.S.; et al. RefSeq: Expanding the Prokaryotic Genome Annotation Pipeline reach with protein family model curation. Nucleic Acids Res. 2021, 49, D1020–D1028. [Google Scholar] [CrossRef]
  79. Okonechnikov, K.; Golosova, O.; Fursov, M.; The UGENE Team. Unipro UGENE: A unified bioinformatics toolkit. Bioinformatics 2012, 28, 1166–1167. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  80. Sullivan, M.J.; Petty, N.K.; Beatson, S.A. Easyfig: A genome comparison visualizer. Bioinformatics 2011, 27, 1009–1010. [Google Scholar] [CrossRef]
  81. Arndt, D.; Grant, J.R.; Marcu, A.; Sajed, T.; Pon, A.; Liang, Y.; Wishart, D.S. PHASTER: A better, faster version of the PHAST phage search tool. Nucleic Acids Res. 2016, 44, W16–W21. [Google Scholar] [CrossRef] [Green Version]
  82. Edwards, R.; Decewicz, P.; Katelyn; Daniel, S.; Laurasisk. linsalrob/PhiSpy: Dropped prophages (v.4.2.19). In Zenodo; CERN: Geneva, Switzerland, 2021. [Google Scholar] [CrossRef]
  83. Bertelli, C.; Laird, M.R.; Williams, K.P.; Simon Fraser University Research Computing Group; Lau, B.Y.; Hoad, G.; Winsor, G.L.; Brinkman, F.S.L. IslandViewer 4: Expanded prediction of genomic islands for larger-scale datasets. Nucleic Acids Res. 2017, 45, W30–W35. [Google Scholar] [CrossRef]
  84. Siguier, P.; Perochon, J.; Lestrade, L.; Mahillon, J.; Chandler, M. ISfinder: The reference centre for bacterial insertion sequences. Nucleic Acids Res. 2006, 34, D32–D36. [Google Scholar] [CrossRef] [Green Version]
  85. Russel, J.; Pinilla-Redondo, R.; Mayo-Munoz, D.; Shah, S.A.; Sorensen, S.J. CRISPRCasTyper: Automated identification, annotation, and classification of CRISPR-Cas loci. CRISPR J. 2020, 3, 462–469. [Google Scholar] [CrossRef]
  86. Liu, B.; Zheng, D.; Jin, Q.; Chen, L.; Yang, J. VFDB 2019: A comparative pathogenomic platform with an interactive web interface. Nucleic Acids Res. 2019, 47, D687–D692. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Evolutionary analyses of S. hyicus-S. intermedius phylogenetic complex including S. ratti sp. nov. (a) Unrooted phylogenetic tree based on complete 16S rRNA gene sequences extracted from whole-genomic sequencing data (GenBank accession numbers are in parentheses). The evolutionary history was inferred by using the maximum likelihood method and Tamura–Nei model. Filled circles indicate that the corresponding nodes were also obtained in the tree constructed by the neighbor-joining method. The percentage of 500 tree replications above 50% in which the associated taxa clustered together is shown next to the branches. The tree is drawn to scale, with branch lengths measured in the number of substitutions per site. There were a total of 1551 positions in the final dataset. (b) Unrooted maximum likelihood tree based on multilocus sequence analysis of concatenated nucleotide sequences from six loci, rpoB, hsp60, dnaJ, tufA, gap and sodA, extracted from whole genome assemblies (accession numbers are in parentheses). There were a total of 3972 positions in the final dataset. Bootstrap probability values (percentages of 500 tree replications) greater than 50% are shown at branch points. The evolutionary distances are given as the number of substitutions per site. (c) Nucleotide core gene set phylogenetic tree of S. ratti sp. nov. and phylogenetically related species. (d) Protein sequence-based phylogenetic tree of the core gene set of S. ratti sp. nov. and phylogenetically related species. Trees (c) and (d) were constructed using up-to-date bacterial core gene set (UBCG; concatenated alignment of 92 core genes). The maximum likelihood tree was inferred using RAxML software and set to 100 replicates. Gene support indices are given at branching points (maximal possible value is 92). Bar, 0.1 substitution per position.
Figure 1. Evolutionary analyses of S. hyicus-S. intermedius phylogenetic complex including S. ratti sp. nov. (a) Unrooted phylogenetic tree based on complete 16S rRNA gene sequences extracted from whole-genomic sequencing data (GenBank accession numbers are in parentheses). The evolutionary history was inferred by using the maximum likelihood method and Tamura–Nei model. Filled circles indicate that the corresponding nodes were also obtained in the tree constructed by the neighbor-joining method. The percentage of 500 tree replications above 50% in which the associated taxa clustered together is shown next to the branches. The tree is drawn to scale, with branch lengths measured in the number of substitutions per site. There were a total of 1551 positions in the final dataset. (b) Unrooted maximum likelihood tree based on multilocus sequence analysis of concatenated nucleotide sequences from six loci, rpoB, hsp60, dnaJ, tufA, gap and sodA, extracted from whole genome assemblies (accession numbers are in parentheses). There were a total of 3972 positions in the final dataset. Bootstrap probability values (percentages of 500 tree replications) greater than 50% are shown at branch points. The evolutionary distances are given as the number of substitutions per site. (c) Nucleotide core gene set phylogenetic tree of S. ratti sp. nov. and phylogenetically related species. (d) Protein sequence-based phylogenetic tree of the core gene set of S. ratti sp. nov. and phylogenetically related species. Trees (c) and (d) were constructed using up-to-date bacterial core gene set (UBCG; concatenated alignment of 92 core genes). The maximum likelihood tree was inferred using RAxML software and set to 100 replicates. Gene support indices are given at branching points (maximal possible value is 92). Bar, 0.1 substitution per position.
Pathogens 11 00051 g001
Figure 2. Transmission electron microscopy image of type strain Staphylococcus ratti sp. nov. CCM 9025T negatively stained with ammonium molybdate. Bar represents 500 nm (original magnification × 10,000).
Figure 2. Transmission electron microscopy image of type strain Staphylococcus ratti sp. nov. CCM 9025T negatively stained with ammonium molybdate. Bar represents 500 nm (original magnification × 10,000).
Pathogens 11 00051 g002
Figure 3. Distinct DNA banding patterns obtained in identification of Staphylococcus ratti sp. nov. and the type strains of related species based on repetitive PCR fingerprinting with (GTG)5 primer. The dendrogram based on cluster analysis of rep-PCR fingerprints was calculated with Pearson’s correlation coefficients with an unweighted pair group method with arithmetic average (UPGMA) clustering method (r, expressed as percentage similarity values).
Figure 3. Distinct DNA banding patterns obtained in identification of Staphylococcus ratti sp. nov. and the type strains of related species based on repetitive PCR fingerprinting with (GTG)5 primer. The dendrogram based on cluster analysis of rep-PCR fingerprints was calculated with Pearson’s correlation coefficients with an unweighted pair group method with arithmetic average (UPGMA) clustering method (r, expressed as percentage similarity values).
Pathogens 11 00051 g003
Figure 4. Whole-genome alignment of Staphylococcus chromogenes NCTC 10530T (assembly accession no. GCA_900458195.1), Staphylococcus ratti sp. nov. CCM 9025T (GCA_020883535.1) and Staphylococcus hyicus ATCC 11249T (GCA_000816085.1) chromosomes. Blast map shows nucleotide sequence identity above 60%. The location of variable genetic elements is color coded as in the legend.
Figure 4. Whole-genome alignment of Staphylococcus chromogenes NCTC 10530T (assembly accession no. GCA_900458195.1), Staphylococcus ratti sp. nov. CCM 9025T (GCA_020883535.1) and Staphylococcus hyicus ATCC 11249T (GCA_000816085.1) chromosomes. Blast map shows nucleotide sequence identity above 60%. The location of variable genetic elements is color coded as in the legend.
Pathogens 11 00051 g004
Figure 5. Pangenome analysis of the Hyicus-Intermedius species group showing genomic diversity within the type strains. Gene clusters (n = 6056) distinguished at 70% blastp identity were grouped by Roary. The pangenome matrix including 992 core genes present in <16 genomes, 2360 shell genes present in 2–15 genomes and 2704 cloud genes present in <2 genomes where the genes were either present or absent is visualized by the Roary plot on the right. The evolutionary insights between species based on the pangenome is shown by the tree on the left.
Figure 5. Pangenome analysis of the Hyicus-Intermedius species group showing genomic diversity within the type strains. Gene clusters (n = 6056) distinguished at 70% blastp identity were grouped by Roary. The pangenome matrix including 992 core genes present in <16 genomes, 2360 shell genes present in 2–15 genomes and 2704 cloud genes present in <2 genomes where the genes were either present or absent is visualized by the Roary plot on the right. The evolutionary insights between species based on the pangenome is shown by the tree on the left.
Pathogens 11 00051 g005
Figure 6. Structure comparison of the phage-inducible chromosomal island SrRICCM9025 from Staphylococcus ratti CCM 9025T, the partial sequence of McRImsr-like island from Macrococcus canis KM0218 and the homologous genomic region from Macrococcus caseolyticus 19Msa0687. GenBank accession numbers and depicted regions are indicated below the species designation. The Blast map shows protein identities above 30%. The genes are color coded by their function as in the legend.
Figure 6. Structure comparison of the phage-inducible chromosomal island SrRICCM9025 from Staphylococcus ratti CCM 9025T, the partial sequence of McRImsr-like island from Macrococcus canis KM0218 and the homologous genomic region from Macrococcus caseolyticus 19Msa0687. GenBank accession numbers and depicted regions are indicated below the species designation. The Blast map shows protein identities above 30%. The genes are color coded by their function as in the legend.
Pathogens 11 00051 g006
Table 1. Phenotypic differentiation of Staphylococcus ratti sp. nov. from closely related Staphylococcus spp. type strains.
Table 1. Phenotypic differentiation of Staphylococcus ratti sp. nov. from closely related Staphylococcus spp. type strains.
Feature S. ratti sp. nov.
CCM 9025T
S. agnetis
CCM 8869T
S. hyicus
CCM 2368T
S. chromogenes
CCM 3387T
Pigmentwhitewhitewhiteorange
Growth at 15 °C++
Growth in 15% NaCl+++
Hemolysis+
Voges–Proskauer test 1+
Acid from trehalose+++
Tween 80 hydrolysis++
Alkaline phosphatase+++
β-glucuronidase 2+w+
Urease++
+, positive; w, weakly positive; −, negative; all data were taken from this study. 1 Tube test; 2 STAPHYtest 24 kit.
Table 2. Cellular fatty acid composition (as a percentage of the total) of S. ratti sp. nov. CCM 9025T, S. chromogenes CCM 3387T, S. agnetis CCM 8869T and S. hyicus CCM 2368T. Values of less than 1% are not shown.
Table 2. Cellular fatty acid composition (as a percentage of the total) of S. ratti sp. nov. CCM 9025T, S. chromogenes CCM 3387T, S. agnetis CCM 8869T and S. hyicus CCM 2368T. Values of less than 1% are not shown.
Fatty AcidS. ratti sp. nov.
CCM 9025T
S. chromogenes
CCM 3387T
S. agnetis
CCM 8869T
S. hyicus
CCM 2368T
iso-C14:0TR1.0TRTR
C14:0TRTR1.01.0
iso-C15:035.016.543.933.7
anteiso-C15:024.136.425.132.5
iso-C16:01.51.5TR1.2
C16:01.51.52.42.7
iso-C17:010.86.09.27.4
anteiso-C17:010.213.44.86.9
C18:02.65.52.94.0
iso-C19:04.43.22.12.1
anteiso-C19:02.32.9TR1.0
C20:05.810.46.16.2
TR, trace amounts < 1%. Cultivation on TSBA medium for 24 h at 37 °C.
Table 3. Candidate virulence factors predicted in Staphylococcus ratti sp. nov. Locus tags in the genome representing homologs with known and previously predicted virulence factors are shown.
Table 3. Candidate virulence factors predicted in Staphylococcus ratti sp. nov. Locus tags in the genome representing homologs with known and previously predicted virulence factors are shown.
Function and RoleVirulence FactorsRelated GenesPrediction in CCM 9025T Genome
AdherenceClumping factor BclfBLN051_01230
Fibronectin binding proteinsfnbALN051_04265
Ser-Asp rich fibrinogen-binding proteinssdrCLN051_00305
EnzymesCysteine proteasesspBLN051_01195
Hyaluronate lyasehysALN051_02175
LipasegehLN051_01425
ThermonucleasenucLN051_06665
Secretion systemType VII secretion systemesaALN051_10885; LN051_10890
esaBLN051_10875
essCLN051_10865
esxALN051_10895
Surface protein anchoringLipoprotein diacylglyceryl transferaselgtLN051_09220
Lipoprotein-specific signal peptidase IIlspALN051_07315
Immune evasionCapsuleUndeterminedLN051_02355; LN051_03230; LN051_07970
Toxinβ-hemolysinhlbLN051_01075
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kovařovic, V.; Sedláček, I.; Petráš, P.; Králová, S.; Mašlaňová, I.; Švec, P.; Neumann-Schaal, M.; Botka, T.; Gelbíčová, T.; Staňková, E.; et al. Staphylococcus ratti sp. nov. Isolated from a Lab Rat. Pathogens 2022, 11, 51. https://doi.org/10.3390/pathogens11010051

AMA Style

Kovařovic V, Sedláček I, Petráš P, Králová S, Mašlaňová I, Švec P, Neumann-Schaal M, Botka T, Gelbíčová T, Staňková E, et al. Staphylococcus ratti sp. nov. Isolated from a Lab Rat. Pathogens. 2022; 11(1):51. https://doi.org/10.3390/pathogens11010051

Chicago/Turabian Style

Kovařovic, Vojtěch, Ivo Sedláček, Petr Petráš, Stanislava Králová, Ivana Mašlaňová, Pavel Švec, Meina Neumann-Schaal, Tibor Botka, Tereza Gelbíčová, Eva Staňková, and et al. 2022. "Staphylococcus ratti sp. nov. Isolated from a Lab Rat" Pathogens 11, no. 1: 51. https://doi.org/10.3390/pathogens11010051

APA Style

Kovařovic, V., Sedláček, I., Petráš, P., Králová, S., Mašlaňová, I., Švec, P., Neumann-Schaal, M., Botka, T., Gelbíčová, T., Staňková, E., Doškař, J., & Pantůček, R. (2022). Staphylococcus ratti sp. nov. Isolated from a Lab Rat. Pathogens, 11(1), 51. https://doi.org/10.3390/pathogens11010051

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