Next Article in Journal
The Odor Release Regularity of Livestock and Poultry Manure and the Screening of Deodorizing Strains
Next Article in Special Issue
Hepatitis B Virus Integration into Transcriptionally Active Loci and HBV-Associated Hepatocellular Carcinoma
Previous Article in Journal
2-Methylisoborneol (2-MIB) Excretion by Pseudanabaena yagii under Low Temperature
Previous Article in Special Issue
Feline Coronaviruses Identified in Feline Effusions in Suspected Cases of Feline Infectious Peritonitis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Antibiotic Resistance and Pathogenomics of Staphylococci Circulating in Novosibirsk, Russia

1
Institute of Chemical Biology and Fundamental Medicine SB RAS, Lavrentiev Avenue 8, 630090 Novosibirsk, Russia
2
Joint-Stock Company Vector-Best, P.O. Box 121, 630117 Novosibirsk, Russia
3
Novosibirsk Research Institute of Traumatology and Orthopedics n.a. Ya. L. Tsivyan, Frunze 17, 630091 Novosibirsk, Russia
*
Authors to whom correspondence should be addressed.
Microorganisms 2021, 9(12), 2487; https://doi.org/10.3390/microorganisms9122487
Submission received: 6 November 2021 / Revised: 26 November 2021 / Accepted: 29 November 2021 / Published: 30 November 2021
(This article belongs to the Special Issue Microorganisms Associated with Infectious Disease)

Abstract

:
A total of 394 strains of staphylococci found in humans and pets in Novosibirsk, Siberian Russia, were characterized in terms of antibiotic resistance and corresponding genes. Two coagulase-positive and 17 coagulase-negative species were identified. The majority of isolates, with the exception of S. haemolyticus and hospital S. epidermidis isolates, were sensitive to most of the tested antibiotics, and isolates from pets displayed the lowest level of resistance. Nevertheless, methicillin-resistant (MRS) and/or multidrug-resistant (MDR) isolates were found in all prevailed species, including coagulase-negative. A set of genes corresponding to the detected resistance was identified: mecA (beta-lactam resistance), aac(6′)-Ie-aph(2″)-Ia, aph(3′)-IIIa, ant(4′)-Ia (aminoglycoside-modifying enzymes), ermA/ermC, and msrA (macrolide resistance). Complete genome analysis for ten MDR S. epidermidis and five MDR S. haemolyticus isolates revealed additional antibiotic resistance genes mphC, qacA/qacB, norA, dfrC/dfrG, lnuA, BseSR, and fosB. NorA, dfrC, and fosB were present in all S. epidermidis genomes, whereas mphC and msrA were identified in all S. haemolyticus ones. All investigated MDR S. epidermidis and four of five S. haemolyticus strains were moderate or strong biofilm producers, whereas multiple genes responsible for this function and for virulence and pathogenicity were identified mostly in S. epidermidis, but were less frequently represented in S. haemolyticus.

1. Introduction

Staphylococci are Gram-positive facultative anaerobic bacteria that belong to the Staphylococcaceae family and Staphylococcus genus. To date, approximately 60 validated members of the genus have been described (https://lpsn.dsmz.de/genus/staphylococcus, accessed on 10 September 2021). Most bacteria from this genus normally inhabit the skin and mucosae of humans and animals and are also a part of soil microbial communities. At the same time, staphylococci can affect almost any organs and tissues of the human body, causing superficial and deep purulent abscesses, respiratory and urinary tract infections, purulent–necrotic processes in postoperative wounds, and food poisoning or intoxication [1].
Staphylococci are divided into coagulase-positive (able to produce coagulase) and coagulase-negative species. Among coagulase-positive species, Staphylococcus aureus is the most clinically important. A number of other coagulase-positive staphylococci (Staphylococcus delphini, Staphylococcus intermedius, and Staphylococcus pseudintermedius, belonging to the Staphylococcus intermedius group) are causative agents of the most common staphylococcal infections in veterinary medicine and can also be transmitted to humans through close contact with animals [2,3,4]. Coagulase-negative members of the genus are considered less virulent and are usually identified as commensals [1]. Nevertheless, cases of infections have become more frequent in people, the causative agents of which are coagulase-negative species Staphylococcus epidermidis, Staphylococcus saprophyticus, Staphylococcus haemolyticus, and Staphylococcus hominis [5]. The most common and studied of the coagulase-negative staphylococci is S. epidermidis, which is normally a commensal of the skin and mucous surfaces. However, the ability of S. epidermidis to form biofilms makes it a serious problem in surgery associated with implanted structures [6].
Notably, S. epidermidis and other coagulase-negative staphylococci are an important reservoir of antibiotic resistance genes. Diseases caused by methicillin-resistant staphylococci (MRS) that are resistant to beta-lactams are particularly difficult to treat, especially when methicillin resistance is associated with reduced sensitivity to antibiotics from other classes. For many years, coagulase-negative MRS strains were considered exclusively as hospital pathogens; however, the situation has changed for the worse as these pathogens are increasingly causing community-acquired infections [7,8]. It has been shown that staphylococci with multiple drug resistances (MDR) can be present in wastewater and in other places associated with human activity [9,10].
Most studies of staphylococci in Russia focused on the monitoring of antibiotic resistant isolates in hospitals and were often limited to characterizing only S. aureus [11,12,13,14,15]. Only a few studies of environmental, veterinary, and community-associated staphylococci in Russia have been published [16,17]. Data on genomes, virulence factors (VF), and antibiotic resistance genes (ARG) of staphylococci found in Russia are still limited, especially in Siberian Russia, which is located between the Far Eastern and European regions. A few S. aureus complete genome studies have been published [17,18,19]. No data regarding the genome characteristics of coagulase-negative staphylococci from Russia have been reported.
In this study, staphylococci found in humans and pets in Novosibirsk, Siberian Russia, were characterized in terms of antibiotic resistance and respective ARGs. In addition, five MDR S. haemolyticus and ten S. epidermidis isolates were examined for the ability to form biofilms. Their complete genomes were sequenced and analyzed, focusing on the genes responsible for pathogenicity factors and biofilm formation.

2. Materials and Methods

2.1. Bacterial Strain Isolation and Identification

Clinical specimens and pure cultures were obtained from a number of Novosibirsk hospitals and medical centers, including Railway Clinical Hospital, Department of Surgery of Purulent Wounds; Novosibirsk Research Institute of Traumatology and Orthopedics, Department of Endoscopic Joint Surgery; Scientific Institute of Clinical and Experimental Lymphology, Department of Diabetic Foot Therapy; Center of New Medical Technologies, Department of Gynecology; City Infectious Diseases Clinical Hospital No. 1, Department of Gastroenterology; and the Federal Center of Neurosurgery. The specimens included biopsy material, wound exudate, synovial liquid, cerebrospinal fluid, swabs, urine and fecal samples, skin and mucosal scrapings, sputum, and pure bacterial cultures; the last was obtained from the hospitals. Veterinary isolates were obtained from specimens taken from pets (cats and dogs) in veterinary clinics of Novosibirsk. The origin and amount of clinical and veterinary isolates are listed in Table 1. Ten-fold dilutions of specimens were prepared and the resulting cell suspensions were plated on mannitol salt agar (OXOID, Basingstoke, UK). Cells were grown at 35 °C in an aerobic atmosphere for 18–24 h, and individual colonies were passaged three times under the same conditions. Pure bacterial cultures were obtained, and bacteria of the genus Staphylococcus were determined by cultural and morphological characteristics. In each of the specimens, from one to three isolates were obtained, differing in colony morphology, growth rate, and biochemical properties. Pure cultures of staphylococci were deposited in the Collection of Extremophilic Microorganisms and Type Cultures of the Institute of Chemical Biology and Fundamental Medicine Siberian Branch of Russian Academy of Science (CEMTC ICBFM SB RAS).
Strains, conserved at −80 °C in a lysogeny broth medium (LB) containing 25% glycerol, were subcultured in LB medium and on LB agar plates. All cultures were grown at 37 °C for 48 h, and the colonies grown were used for PCR amplification of the 16S rRNA gene. Subsequent identification of the strains was performed by sequencing the 1308-bp PCR fragment of the 16S rRNA gene, as described previously [20]. All PCR amplicons were gel purified (0.6% SeaKem® GTG-agarose, Lonza, ME, USA) and sequenced (BigDye™ Terminator v.3.1 Cycle Sequencing Kit and ABI 3500 Genetic Analyzer, Applied Biosystems, Foster City, CA, USA). The obtained nucleotide sequences of 16S rRNA genes were compared with the corresponding nucleotide sequences extracted from the NCBI GenBank database (http://www.ncbi.nlm.nih.gov, accessed on 1 April 2021). The closest reference sequences of the genus Staphylococcus were used with a sequence similarity level of at least 98%, and species names were determined according to a recently published updated classification [21]. Additional identification was done for a number of isolates using GEN III OmniLog Plus ID System (Biolog, Inc., Hayward, CA, USA).
All patient data were collected anonymously. Isolates were considered as hospital specimens if they were obtained from patients who had been in the clinic for at least two weeks and were treated with antibiotics during that time. The isolates were considered as outpatient if they were obtained from patients who were not hospitalized or from patients who had been hospitalized only a few hours. The study was approved by a local ethics Committee of the Center for New Medical Technologies in Novosibirsk; protocol #2, 12 February 2014.

2.2. Susceptibility Testing

Susceptibility to five antibiotics was determined using a disk diffusion test according to the guidelines of EUCAST 10.0 (https://eucast.org, accessed on 10 September 2020). Cefoxitin (FOX, 30 µg), amikacin (AMK, 30 µg), gentamicin (CN, 10 µg), erythromycin (E, 15 µg), and clindamycin (DA, 2 µg) were examined. Disks with antibiotics (OXOID) were applied to the lawns of the investigated cultures on Mueller–Hinton agar (OXOID). Cefoxitin was used as a marker for the detection of methicillin resistance and strain S. aureus ATCC 25923 was used as a control of susceptibility.

2.3. Antibiotic Resistance Genes Detection

Bacterial DNA was extracted from bacterial cultures grown at 35 °C overnight in 2 mL of brain heart infusion broth (BHI Broth, BioMerieux, Marcy-l’Étoile, France) using a DNA extraction kit (Biolabmix, Novosibirsk, Russia) according to the manufacturer’s instructions. Subsequently, all strains were tested for the presence of genes that contribute to the resistance to various classes of antibiotics: mecA and blaZ genes (beta-lactam and penicillin resistance), aac(6′)-Ie-aph(2″)-Ia, ant(4′)-Ia, and aph(3′)-IIIa genes (aminoglycoside resistance), and erm A, erm C, and msrA genes (macrolide and lincosamide resistance) using a set of primers (Table S1).

2.4. Study of Biofilm Formation

The study of staphylococcus biofilm formation with the crystal violet staining method was done as described previously [22]. Briefly, an overnight growth culture in tryptic soy broth (TSB) medium (Merch, Darmstadt, Germany) was adjusted to a final OD600 of 1.00 ± 0.05 by adding sterile TSB. This OD-adjusted suspension was then diluted 100-fold with TSB containing 1% glucose to obtain the initial bacterial suspension. The wells of a 96-well microplate (TPP, Trasadingen, Switzerland) were inoculated with 200 μL of the suspension and incubated at 35 °C for 24 h. A control well was inoculated with a sterile TSB medium containing 1% of glucose. Each experiment was done in triplicate in three technical repeats. The medium was removed and the wells were washed three times with sterile saline solution. The wells were air dried for 45 min, and the adherent cells were stained with a 0.1% water solution of crystal violet at room temperature for 30 min. Next, the excess crystal violet was removed and the wells were washed five times with 300 μL sterile saline solution. The dye was dissolved with 150 μL of 96% ethanol and the absorbance of each well was read at 595 nm in a microplate reader (Uniplan, Moscow, Russia). The results were processed and interpreted as described previously [22]. In brief, the average OD values (ODav) were calculated for all tested strains and for a negative control (non-inoculated medium). The cutoff value (ODc) was defined as three standard deviations above the mean OD of the negative control (ODc = ODav-neg + 3 × SD). The strains were divided into three categories, based on the calculated ODc: ODc < ODav ≤ 2 × ODc was defined as a weak biofilm producer; 2 × ODc < ODav ≤ 4 × ODc was defined as a moderate biofilm producer; 4ODc ≤ ODav was defined as a strong biofilm producer.

2.5. Complete Genome Sequencing and Analysis

Bacterial DNA was extracted from bacterial cultures grown at 35 °C overnight in 10 mL of brain heart infusion broth (BioMerieux, Marcy-l’Étoile, France) using a DNA extraction kit (Biolabmix, Novosibirsk, Russia) according to the manufacturer’s instructions.
The extracted DNA was further used for a paired-end library construction, which was performed using a NEB Next DNA Ultra library prep kit (New England Biolabs, Ipswich, MA, USA). Sequencing was carried out using a MiSeq Benchtop Sequencer (Illumina, Inc., San Diego, CA, USA) and a MiSeq Reagent Kit 2 × 250 v. 2.0 (Illumina, San Diego, CA, USA). The obtained sequences were assembled de novo using the SPAdes genome assembler v.3.15.2 (http://cab.spbu.ru/software/spades, accessed on 1 August 2021). Annotation was performed using Rapid Annotation Subsystem Technology (RAST) v.2.0 (https://rast.nmpdr.org, accessed on 2 September 2021). The genomes and their RAST annotations are located in Supplementary Data S1 and S2, respectively. The obtained sequences were analyzed for the presence of antibiotic resistance genes (ARG) and virulence factors (VF) using a Resistance Genes Identifier (https://card.mcmaster.ca/analyze/rgi, accessed on 2 September 2021) and the Virulence Factor Database (http://www.mgc.ac.cn/VFs, (accessed on 3 September 2021), respectively. In addition, the genomes were checked for the presence of ARGs and VFs manually, using RAST annotation and BLAST search. The sequence types for MDR Staphylococcus strains were determined using the PubMLST database (https://pubmlst.org, accessed on 15 September 2021), and the clonal complexes were identified using a BURST analysis [23].
All Staphylococcus sequences/genomes are publicly available through NSBI resources under BioProject PRJNA774949. Unassembled sequencing reads have been deposited in the NCBI Short Read Archive (SRA) under accessions SRR16591584-SRR16591598. Whole Genome Shotgun projects have been deposited at DDBJ/ENA/GenBank under accessions JAJFNL000000000, JAJFNM000000000, JAJFNN000000000, JAJFNO000000000, JAJFNP000000000, JAJFNQ000000000, JAJFNR000000000, JAJFNS000000000, JAJFNT000000000, JAJFNU000000000, JAJFNV000000000, JAJFNW000000000, JAJFNX000000000, JAJFNY000000000, and JAJFNZ000000000. Versions described in this paper are JAJFNL010000000, JAJFNM010000000, JAJFNN010000000, JAJFNO010000000, JAJFNP010000000, JAJFNQ010000000, JAJFNR010000000, JAJFNS010000000, JAJFNT010000000, JAJFNU010000000, JAJFNV010000000, JAJFNW010000000, JAJFNX010000000, JAJFNY010000000, and JAJFNZ010000000.

3. Results

3.1. Staphylococcus Strain Isolation and Identification

In 2014–2020, a collection of 394 strains of staphylococci from humans and pets (dogs and cats) was isolated in Novosibirsk, Russia (Table 1). Nineteen Staphylococcus species were revealed (Table 2), and five strains were identified as belonging to the newly established Mammaliicoccus genus [21].
S. aureus, S. epidermidis, S. haemolyticus, and S. hominis were the major species in clinical samples, whereas S. pseudintermedius was the dominant species in veterinary samples. The distribution of coagulase-positive and coagulase-negative staphylococci varied in different panels (Figure 1). Among the strains obtained from hospitalized patients, S. aureus strains prevailed (~70%); the share of S. epidermidis strains was approximately 19%, and S. haemolyticus, S. warneri, and other staphylococci accounted for 11%. Outpatient isolates had a different ratio of species: S. aureus, S. epidermidis, and S. haemolyticus were at approximately 35%, 37%, and 9%, respectively. Veterinary isolates contained ~39%, 20%, and 13% of the S. pseudintermedius, S. aureus, and S coagulans strains, respectively (Figure 1, Table 2). The 16S rRNA gene sequences of all investigated Staphylococcus strains (except S. aureus strains) were deposited in the NCBI GenBank database. A list of the accession numbers for 16S rRNA gene sequences is given in Table 2.

3.2. Antibiotic Resistance and Resistance-Encoding Genes

All Staphylococcus isolates were tested for antibiotic resistance (Figure 2). The majority of them (~85%, 330/389) were methicillin-sensitive (MSS), including 91% of S. aureus, 70% of S. epidermidis, 85% of S. hominis, and all isolates from other species. The exception was S. haemolyticus isolates: ~60% (13/22) of them were MRS.
Nevertheless, MRS and MDR isolates resistant to three or more classes of antibiotics were found in all predominant species (Figure 2). All identified MDR isolates were methicillin-resistant, with the exception of S. pseudintermedius MDR isolates, which were sensitive to cefoxitin, but resistant to the tested non-beta-lactam drugs. Notably, five of seven S. aureus MDR hospital isolates were obtained from purulent diabetic ulcers and six S. epidermidis MDR hospital isolates were associated with prosthetic joint and postoperative wound infections. Isolates from the other studied Staphylococcus species (N = 63) were sensitive to most of the tested antibiotics; however, 10 of 63 commensal isolates were resistant to erythromycin, including all S. devriesei isolates, two S. simulans isolates, and one isolate each of M. lentus, S. borealis, S. warneri, and S. pasteuri. The S. pasteuri isolate was also resistant to gentamicin.
In addition, all isolates were assayed for antibiotic resistance genes (ARGs). The presence of the ARGs among the prevailing Staphylococcus species is shown in Table 3. The mecA gene was identified in all MRS isolates (Table 3); approximately half of all S. aureus and half of outpatient S. epidermidis isolates contained the blaZ gene, encoding penicillin resistance. The majority of aminoglycoside-resistant Staphylococcus isolates (59/63) contained a single aac(6′)-Ie-aph(2″)-Ia gene (n = 34), or its combinations with aph(3′)-IIIa (n = 14) or ant(4′)-Ia (n = 11) genes. The erythromycin resistance of S. aureus isolates was encoded only by the ermA/ermC genes, not msrA. In contrast, the outpatient erythromycin-resistant S. epidermidis isolates mostly contained msrA genes (n = 29). A combination of ermA and msrA was revealed in only one S. hominis isolate (Table 3).
In general, there was a clear correlation between the resistance and the presence of the corresponding ARGs. The exception was seven S. pseudintermedius erythromycin-resistant isolates and two S. epidermidis clindamycin-resistant isolates; in all these isolates the ermA, ermC, and msrA genes were not detected. Probably, other ARGs are responsible for resistance in these isolates. Conversely, genes encoding aminoglycoside resistance were determined in four S. haemolyticus and two S.hominis isolates; all of them were sensitive to aminoglycosides. The same was observed in erythromycin-sensitive strain of S. haemolyticus, in which the msrA gene was detected (Table 3). The lack of resistance in these cases is likely due to the mutant variants of the detected gene and/or the lack of its activity.

3.3. MDR Isolate Sequence Type Identification

Complete genomes were determined for ten and five MDR isolates of S. epidermidis and S. haemolyticus, respectively. The sequence types for staphylococcus isolates (Table 4) were determined using the PubMLST database (https://pubmlst.org, accessed on 15 September 2021) and the BURST analysis was used for grouping the isolates into clonal complexes [23]. Eight sequence types were found among ten S. epidermidis isolates; six of them were members of the CC5, one (ST23) was determined to be a member of CC23, and two ST 20 isolates were evaluated as singletons. All S. haemolyticus isolates were grouped into the clonal complex 3 (Table 4).

3.4. Analysis of the Investigated Genomes for the Presence of Antibiotic Resistance Genes

An ARG search was performed using Resistance Gene Identifier (https://card.mcmaster.ca/analyze/rgi, accessed on 2 September 2021). The default selection criteria, which identified genes based on the strict or perfect mode, were used. In addition to those previously identified (Table 3), multiple ARGs have been found in the investigated genomes (Table 5). The main differences between S. epidermidis and S. haemolyticus isolates were as follows: all S. epidermidis genomes contained norA (quinolone resistance), dfrC (diaminopyrimidine resistance), and fosB (fosfomycin resistance), and nine of ten genomes possessed an aminocoumarin-resistant variant of gyrB. None of these genes were revealed in S. haemolyticus genomes. On the contrary, msrA and mphC (macrolide resistance) were found in all S. haemolyticus genomes but not in S. epidermidis genomes. Both the S. haemolyticus and S. epidermidis genomes had one to three different genes encoding aminoglycoside-modifying enzymes, and the aac(6′)-Ie-aph(2″)-Ia gene was identified in all isolates. FusC (fusidic acid resistance) and cat8 (chloramphenicol resistance) were rare and found only in S. epidermidis CEMTC 3750 and in S. haemolyticus CEMTC1553 isolates (Table 5). Genes encoding vancomycin resistance were not identified; the cluster of bacitracin resistance genes bce was found in all studied genomes.
In addition, a number of genes encoding different MDR pumps mediating the efflux were found in the genomes. Among them, genes responsible for biocide resistance were identified, including the gene encoding acriflavin resistance protein (found in all genomes) and qacA encoding antiseptic resistance protein (identified in eight S. epidermidis and four S. haemolyticus genomes, except the 3107, 3117, and 2688 isolates).

3.5. Virulence Factor Identification

A number of VFs were identified in the investigated genomes, including genes encoding toxins, exoenzymes, and immune evasion factors (Table 6). All investigated Staphylococcus genomes contained the genes lip and nuc, encoding lipase and thermonuclease, respectively. In addition, two hemolysin-encoding genes (hemolysin III and hemolysin, containing the CBS domain) were found in all studied genomes. The pgs ABCDE cluster, responsible for the synthesis of the surface-attached poly-gamma-glutamate (PGA), was also identified in all S. epidermidis and S. haemolyticus genomes (Table 6). The PGA capsule is produced by a number of coagulase-negative staphylococci and is absent in S. aureus. It is a key factor of pathogen survival during infection and efficiently shelters it from components of innate host defense, including antimicrobial peptides and neutrophil phagocytosis [24,25]. Only S. epidermidis genomes contained the gene hlb, encoding β-hemolysin, and a set of the genes sspA, sspB, and geh, encoding exoenzymes.
The number of virulence factors in the studied S. epidermidis genomes was lower compared to the known pathogenic S. aureus genomes. This is typical for coagulase-negative staphylococci, which are usually not as pathogenic as S. aureus isolates [5]. Notably, the number of VFs identified in the genomes of S. haemolyticus was lower than in the genomes of studied S. epidermidis, which may partially be due to insufficient information about the virulence factors of S. haemolyticus [5].
At the same time, a capsular of investigated S. haemolyticus isolates probably has a more complex structure compared to S. epidermidis, as additional genes encoding enzymes of capsular synthesis were revealed in the genomes of S. haemolyticus strains CEMTC 2119 and CEMTC 3413 (Table 6). These genes possess more than 50% similarity with the closely related cap5 and cap8 gene clusters of S. aureus, which are responsible for CP5 and CP8 (the two main CP serotypes usual for clinical S. aureus strains) [26,27].

3.6. In Vitro Biofilm Formation by MDR S. epidermidis and S. haemolyticus Isolates and Factors Responsible for Adhesion and Biofilm Formation

The studied isolates have demonstrated a different ability to form biofilms (Figure 3). All isolates were divided into three categories (weak, moderate, and strong biofilm producers), as described previously [22]. Based on the cutoff value (ODc), set as three standard deviations above the mean OD of negative control and calculated as 0.15, S. haemolyticus CEMTC 3413 was classified as a weak biofilm producer (0.15 < ODav ≤ 0.3), S. epidermidis CEMTC 1833, CEMTC 3750, and S. haemolyticus CEMTC 2119, CEMTC 3601 were defined as moderate biofilm producers (0.3 < ODav ≤ 0.6), and the remaining isolates were defined as strong biofilm producers (0.6 ≤ ODav). The S. epidermidis CEMTC 3824 and CEMTC 3918 isolates showed the highest biofilm formation.
Two distinct mechanisms of biofilm formation have been identified previously. One involves a number of cell-wall-associated proteins (CWA), responsible for adhesion to the host cells or abiotic surfaces and subsequent biofilm production, and the other requires the synthesis and secretion of a polysaccharide intercellular adhesin (PIA) [28,29]. The genome analysis revealed a number of genes that promote biofilm formation (Table 7). All S. epidermidis and S. haemolyticus genomes contained atlE (autolysin), ebpS (elastin binding protein), and sas family genes, all of them encoding CWA proteins responsible for adhesion. A number of sdr family genes, encoding Ser-Asp-rich fibrinogen-binding proteins, were found in all S. epidermidis genomes; however, two of five S. haemolyticus genomes contained a single sdrC gene. The aap gene (which encodes an accumulation-associated protein, Aap) was found in the genomes of all S. epidermidis strains capable of strong biofilm production, except strain 1827.
No clear correlation between the presence of ica-operon and the ability to form biofilms was found (Table 7, Figure 3). Six S. epidermidis genomes contained ica-operon, associated with a high level of biofilm formation. Most of them, except strain 1833, were strong biofilm producers. Notably, S. epidermidis strains 3824 and 3918, which differed considerably from the others in terms of their ability to produce biofilms (Figure 3), had a combination of ica-operon, aap, pls, and four sdr-genes. At the same time, genes encoding phenol-soluble modulins (PSMs) were absent from both strains (Table 7). PSMs disrupt non-covalent bonds between cells inside the biofilm, promoting channel formation and cell release from the mature biofilm; mutant S. epidermidis strains, lacking PSM, were unable to form channels and produced more substantial biofilms [31,32,33].
In general, the investigated S. epidermidis strains were strong or moderate biofilm producers and their genomes contained multiple genes responsible for this function. At the same time, only two of five S. haemolyticus strains were capable of intense biofilm formation, and the well-known genes responsible for the formation of biofilms were poorly represented in their genomes. Presumably, their genomes contain some unidentified genes responsible for biofilm formation. Notably, both S. haemolyticus strains CEMTC 3413 and CEMTC 2119, possessing weak or moderate ability for biofilm formation, contained a putative capsule synthesis operon (Table 7). This may be one of the reasons for the weak biofilm formation—only non-encapsulated cells are able to adhere to the extracellular matrix or to endothelial cells [34].

4. Discussion

The widespread, sometimes unjustified use of antibiotics in recent decades has led to a significant increase in the proportion of antibiotic-resistant and MDR strains among nosocomial agents [35,36,37,38]. In this study, a wide spectrum of staphylococci species found in hospitalized patients and outpatients in Novosibirsk (with more than 1.5 million inhabitants) were investigated for their antibiotic resistance. In addition, staphylococci isolated from pets in close contact with humans were involved in the study, because such staphylococci can be a reservoir of pathogenic MRS and MDR strains [39,40].
Only two coagulase-positive species were identified: S. aureus (in humans and pets) and S. pseudintermedius (only in pets); such staphylococci dominated in samples from hospitalized patients and pets (~70% and ~60%, respectively). Seventeen coagulase-negative staphylococci were revealed, with the most variety in outpatients (15 species), where coagulase-negative staphylococci were in the majority (~75%). Two coagulase-negative species, S. coagulans and S. felis, were found only in veterinary samples.
Our investigation showed that the majority of isolates of all found species, with the exception of S. haemolyticus and hospital S. epidermidis isolates, were sensitive to most of the antibiotics tested and the number of MRS and MDR strains was not high. Most S. haemolyticus isolates were MRS, which corresponded to the known data on a high level of methicillin resistance in S. haemolyticus [5]. The animal isolates had the lowest antibiotic resistance to tested antibiotics, and only one veterinary S. haemolyticus MRS and four S. pseudintermedius MDR isolates (all methicillin-sensitive) were found among 61 isolates from domestic animals. Nevertheless, MRS and/or MDR isolates were found in all prevailed species, including hospital, outpatient, and veterinary isolates.
The data on MDR and MRS coagulase-negative staphylococci in different countries and regions are heterogeneous and depend on the country, region, type of hospital, etc. It is also known that hospital, community-associated, and environmental coagulase-negative isolates differ significantly in their antibiotic resistance, with the highest level of resistance a characteristic of hospital isolates. Unfortunately, the situation is changing for the worse in recent years. For instance, 17% of community-associated isolates, 30% of healthcare personnel isolates, and ~87% of hospital S. epidermidis isolates were MRS in 2013 in Shanghai, China [41]. In 2018, a higher percentage of MRSE isolates (76.5%) was recovered from healthcare personnel in two public hospitals in Tianjin, China [42]. Another study revealed that 20.5% of environmental CoNS staphylococci detected in the environment of a university in Thailand were MRS and 61.0% of these MRS isolates were MDR [43].
The beta-lactam resistance was mediated by the mecA gene, which was found in all MRS isolates. Aminoglycoside-modifying enzymes were represented in all species and encoded mainly by the aac(6′)-Ie-aph(2″)-Ia gene and its combination with aph(3′)-IIIa or ant(4′)-Ia genes. At the same time, the macrolide resistance genes distribution was variable: only ermA/ermC genes were found in S. aureus isolates; in contrast, macrolide-resistant S. epidermidis, S. haemolyticus, and S. hominis isolates were mostly found to contain the msrA gene. Perhaps this indicates separate ways of the resistance genes transfer between different species of coagulase-negative staphylococci. None of the macrolide-resistant veterinary S. pseudintermedius isolates contained ermA/ermC or msrA genes. Other mechanisms (efflux pump or antibiotic-modifying enzymes) likely mediated macrolide resistance in these isolates.
Complete genome sequencing and analysis, which was performed for ten MDR S. epidermidis and five MDR S. haemolyticus isolates, identified other ARGs and genes associated with biofilm formation and virulence. The species-specific distribution of a number of ARGs was revealed in the genomes; norA, dfrC, and fosB genes were found in all S. epidermidis genomes, whereas mphC and msrA were identified in all S. haemolyticus ones.
The distribution of the macrolide-resistance-encoding genes msrA and ermA/ermC, which were found in complete genomes, was in contrast to the usual distribution of these genes among S. epidermidis isolates. Notably, almost all macrolide-resistant S. epidermidis isolates for which genome sequences were determined in this study were isolated from hospitalized patients; possibly, ermA/ermC genes were transmitted to S. epidermidis isolates from the nosocomial strains of S. aureus.
Six S. epidermidis MDR hospital isolates were associated with prosthetic joint and post-operative wound infections. The last association may be mediated by one of the major factors of S. epidermidis pathogenicity, their biofilm formation ability [5,6]. Indeed, four of the six S. epidermidis MDR hospital isolates were strong biofilm producers and the remaining two isolates demonstrated moderate biofilm-forming activity. In addition, all these strains contained the antiseptic resistance gene qacA and multiple ARGs. All of these factors were found to correlate with poor patient outcome, as was shown previously [6].
In addition, the MLST profiling revealed the genetic diversity of isolates: eight different ST were found among ten S. epidermidis isolates. The observed genetic diversity is in accordance with previously published data [44,45]. The majority of identified S. epidermidis sequence types belonged to clonal complex CC5, which has been previously reported as one of the leading causative agents for bloodstream and prosthetic joint infections [45,46,47]. A recent study [48] has shown that ST5 isolates can be detected as matched commensal/invasive pairings in the same human organism, unlike ST2, which is probably purely a hospital pathogen. Presumably, ST5 isolates are able to evolve in one organism from commensal to infectious variants.
A number of S. epidermidis isolates have been previously reported in Russia (https://pubmlst.org, accessed on 15 September 2021). Among the sequence types identified in this study, the ST5, ST20, and ST23 isolates were found in 2008–2009 in Moscow and Nizhny Novgorod, European Russia, whereas ST786 was obtained from a hospitalized patient in 2018 in Moscow [49]. Others, namely ST17, ST152, and ST210, were found for the first time in Russia. All eight STs were revealed previously as both colonizing and infectious agents. Isolate ST210 was reported at first to be associated with samples taken from healthy people [41]; however, a number of reports of ST210 isolates causing infection were subsequently published [50,51]. The ST17, ST20, and ST23 isolates were previously identified as both colonizing and infectious agents [44]. A set of ST152 isolates were previously revealed as human-colonizing staphylococci (isolated in the USA) and as mastitis-associated agents in cows (isolated in Greece), whereas ST786 was identified in China as a community-associated isolate (https://pubmlst.org, accessed on 15 September 2021).
S. haemolyticus isolates were identified as belonging to four different ST; all of them were grouped into CC3. STs of S. haemolyticus isolates from Russia have not been previously reported (https://pubmlst.org, accessed on 15 September 2021).
In conclusion, the majority of isolates of 19 staphylococci species found in Novosibirsk in 2014–2020 were sensitive to most of the tested antibiotics. The exception was S. haemolyticus and hospital S. epidermidis isolates. Nevertheless, MRS and/or MDR isolates were found in all prevailed (S. aureus, S. epidermidis, S. haemolyticus, S. hominis, and S. pseudintermedius) species. ARGs (mecA, aac(6′)-Ie-aph(2″)-Ia, aph(3′)-IIIa, ant(4′)-Ia, ermA/ermC, and msrA) were identified in staphylococcus isolates, and a clear correlation with the corresponding resistance was revealed. Complete genomes were sequenced and analyzed for ten MDR S. epidermidis and five MDR S. haemolyticus isolates. Antibiotic resistance genes mphC, qacA/qacB, norA, dfrC/dfrG, lnuA, BseSR, fosB, and multiple genes responsible for virulence and pathogenicity were found. All investigated MDR S. epidermidis and four of five S. haemolyticus strains were moderate or strong biofilm producers, whereas multiple genes responsible for this function were identified mostly in S. epidermidis genomes and were less represented in S. haemolyticus genomes.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/microorganisms9122487/s1, Table S1: Primers, specific to antibiotic resistance encoding genes; Data S1: Staphylococcus genome sequences; Data S2: RAST annotation of the genomes.

Author Contributions

Conceptualization, V.M. and N.T.; methodology, A.B., V.P. and N.F.; validation, S.C., A.T. and V.M.; formal analysis, A.B. and V.M.; investigation, A.B., V.F., T.K. and Y.K.; resources, V.M.; data curation, A.T., E.Z. and A.B.; writing—original draft preparation, A.B. and V.M.; writing—review and editing, N.T.; visualization, A.B. and V.M.; supervision, N.T.; project administration, V.M.; funding acquisition, V.M. All authors have read and agreed to the published version of the manuscript.

Funding

The research was funded by Russian Foundation of Basic Research and government of Novosibirsk region, project number 20-44-540012; Collection of Extremophilic Microorganisms and Type Cultures of the ICBFM SB RAS was funded with Russian State Budget Project 121031300043-8.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and was approved by the local Ethics Committee of the Center for New Medical Technologies in Novosibirsk; protocol # 2, 12 February 2014.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available in supplementary, Data S1: Staphylococcus genome sequences; Data S2: RAST annotation of the genomes.

Acknowledgments

The authors thank Kaverina Galina for excellent technical support. Authors are grateful to Ganichev D.A. (Railway Clinical Hospital); Nimaev V.V. (Scientific Institute of Clinical and Experimental Lymphology); Yakovets E.A. (Center of New Medical Technologies); Pozdnyakova L.L. (City Infectious Diseases Clinical Hospital No. 1) for specimens obtained.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Götz, F.; Bannerman, T.; Schleifer, K.H. The Prokaryotes; Springer: New York, NY, USA, 2006; pp. 5–75. [Google Scholar]
  2. Somayaji, R.; Priyantha, M.A.; Rubin, J.E.; Church, D. Human infections due to Staphylococcus pseudintermedius, an emerging zoonosis of canine origin: Report of 24 cases. Diagn. Microbiol. Infect. Dis. 2016, 85, 471–476. [Google Scholar] [CrossRef]
  3. Darlow, C.A.; Paidakakos, N.; Sikander, M.; Atkins, B. A spinal infection with Staphylococcus pseudintermedius. BMJ Case Rep. 2017, 2017, bcr2017221260. [Google Scholar] [CrossRef] [PubMed]
  4. Magleby, R.; Bemis, D.A.; Kim, D.; Carroll, K.C.; Castanheira, M.; Kania, S.A.; Jenkins, S.G.; Westblade, L.F. First reported human isolation of Staphylococcus delphini. Diagn. Microbiol. Infect. Dis. 2019, 94, 274–276. [Google Scholar] [CrossRef] [PubMed]
  5. Argemi, X.; Hansmann, Y.; Prola, K.; Prévost, G. Coagulase-Negative Staphylococci Pathogenomics. Int. J. Mol. Sci. 2019, 20, 1215. [Google Scholar] [CrossRef] [Green Version]
  6. Post, V.; Harris, L.G.; Morgenstern, M.; Mageiros, L.; Hitchings, M.D.; Méric, G.; Pascoe, B.; Sheppard, S.K.; Richards, R.G.; Moriarty, T.F. A comparative genomics study of Staphylococcus epidermidis from orthopedic device-related infections correlated with patient outcome. J. Clin. Microbiol. 2017, 55, 3089–3103. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Otto, M. Staphylococcus epidermidis—The ‘accidental’ pathogen. Nat. Rev. Microbiol. 2009, 7, 555–567. [Google Scholar] [CrossRef] [Green Version]
  8. Namvar, A.E.; Bastarahang, S.; Abbasi, N.; Ghehi, G.S.; Farhadbakhtiarian, S.; Arezi, P.; Hosseini, M.; Baravati, S.Z.; Jokar, Z.; Chermahin, S.G. Clinical characteristics of Staphylococcus epidermidis: A systematic review. GMS Hyg. Infect. Control 2014, 9, Doc23. [Google Scholar] [CrossRef]
  9. Gómez, P.; Lozano, C.; Benito, D.; Estepa, V.; Tenorio, C.; Zarazaga, M.; Torres, C. Characterization of staphylococci in urban wastewater treatment plants in Spain, with detection of methicillin resistant Staphylococcus aureus ST398. Environ. Pollut. 2016, 212, 71–76. [Google Scholar] [CrossRef]
  10. Xu, Z.; Misra, R.; Jamrozy, D.; Paterson, G.K.; Cutler, R.R.; Holmes, M.A.; Gharbia, S.; Mkrtchyan, H.V. Whole genome sequence and comparative genomics analysis of multi-drug resistant environmental Staphylococcus epidermidis ST59. G3 Genes Genomes Genet. 2018, 8, 2225–2230. [Google Scholar] [CrossRef] [Green Version]
  11. Stratchounski, L.S.; Dekhnich, A.V.; Kretchikov, V.A.; Edelstain, I.A.; Narezkina, A.D.; Afinogenov, G.E.; Akhmetova, L.I.; Boronina, L.G.; Gugutcidze, E.N.; Gudkova, L.V.; et al. Antimicrobial resistance of nosocomial strains of Staphylococcus aureus in Russia: Results of a prospective study. J. Chemother. 2005, 17, 54–60. [Google Scholar] [CrossRef]
  12. Dekhnich, A.V.; Nikulin, A.A.; Ryabkova, E.L.; Krechikova, O.I.; Sukhorukova, M.V.; Kozlov, R.S.; ROSNET Study Group. Epidemiology of resistance of S. aureus strains isolated from patients in the ICU of Russian hospitals: Results of a multicenter study. Clin. Microbiol. Antimicrob. Chemother. 2008, 10, 333–344. (In Russian) [Google Scholar]
  13. Sukhorukova, M.V.; Skleenova, E.Y.; Ivanchik, N.V.; Timokhova, A.V.; Edelstein, M.V.; Dekhnich, A.V.; Kozlov, R.S.; Shek, E.A.; «MARATHON» Study Group. Antimicrobial resistance of nosocomial Staphylococcus aureus isolates in Russia: Results of national multicenter surveillance study «MARATHON» 2011–2012. Clin. Microbiol. Antimicrob. Chemother. 2014, 16, 280–286. (In Russian) [Google Scholar]
  14. Romanov, A.V.; Dekhnich, A.V.; Sukhorukova, M.V.; Skleenova, E.Y.; Ivanchik, N.V.; Edelstein, M.V.; Kozlov, R.S.; the «MARATHON» Study Group. Antimicrobial resistance of nosocomial Staphylococcus aureus isolates in Russia: Results of multicenter epidemiological study «MARATHON» 2013–2014. Clin. Microbiol. Antimicrob. Chemother. 2017, 19, 57–62. (In Russian) [Google Scholar]
  15. Barantsevich, N.E.; Barantsevich, E.P. Species diversity and methicillin resistance in Staphylococcus spp. in nosocomial infections. Clin. Microbiol. Antimicrob. Chemother. 2019, 21, 207–211. (In Russian) [Google Scholar] [CrossRef] [Green Version]
  16. Kozlova, Y.N.; Fomenko, N.V.; Morozova, V.V.; Saranina, I.; Tikunov, A.; Ganichev, D.A.; Samokhin, A.; Pavlov, V.; Rozhnova, O.; Bondar, I.A.; et al. Genetic and biochemical characterization of staphylococci occurring in Novosibirsk, Russia. Vavilovskii Zhurnal Genet. Sel. 2017, 21, 952–958. [Google Scholar] [CrossRef] [Green Version]
  17. Fursova, K.; Sorokin, A.; Sokolov, S.; Dzhelyadin, T.; Shulcheva, I.; Shchannikova, M.; Nikanova, D.; Artem’eva, O.; Zinovieva, N.; Brovko, F. Virulence factors and phylogeny of Staphylococcus aureus associated with bovine mastitis in Russia based on genome sequences. Front. Vet. Sci. 2020, 7, 135. [Google Scholar] [CrossRef]
  18. Khokhlova, O.E.; Hung, W.C.; Wan, T.W.; Iwao, Y.; Takano, T.; Higuchi, W.; Yachenko, S.V.; Teplyakova, O.V.; Kamshilova, V.V.; Kotlovsky, Y.V.; et al. Healthcare- and community-associated methicillin-resistant Staphylococcus aureus (MRSA) and fatal pneumonia with pediatric deaths in Krasnoyarsk, Siberian Russia: Unique MRSA’s multiple virulence factors, genome, and stepwise evolution. PLoS ONE 2015, 10, e0128017. [Google Scholar] [CrossRef]
  19. Wan, T.W.; Khokhlova, O.E.; Iwao, Y.; Higuchi, W.; Hung, W.C.; Reva, I.V.; Singur, O.A.; Gostev, V.V.; Sidorenko, S.V.; Peryanova, O.V.; et al. Complete circular genome sequence of successful ST8/SCCmecIV community-associated methicillin-resistant Staphylococcus aureus (OC8) in Russia: One-megabase genomic inversion, IS256’s spread, and evolution of Russia ST8-IV. PLoS ONE 2016, 11, e0164168. [Google Scholar] [CrossRef]
  20. Morozova, V.; Kozlova, Y.; Shedko, E.; Kurilshikov, A.; Babkin, I.; Tupikin, A.; Yunusova, A.; Chernonosov, A.; Baykov, I.; Kondratov, I. Lytic Bacteriophage PM16 Specific for Proteus mirabilis: A Novel Member of the genus phiKMVvirus. Arch. Virol. 2016, 161, 2457–2472. [Google Scholar] [CrossRef]
  21. Madhaiyan, M.; Wirth, J.S.; Saravanan, V.S. Phylogenomic analyses of the Staphylococcaceae family suggest the reclassification of five species within the genus Staphylococcus as heterotypic synonyms, the promotion of five subspecies to novel species, the taxonomic reassignment of five Staphylococcus species to Mammaliicoccus gen. nov., and the formal assignment of Nosocomiicoccus to the family Staphylococcaceae. Int. J. Syst. Evol. Microbiol. 2020, 70, 5926–5936. [Google Scholar] [CrossRef]
  22. Stepanović, S.; Vuković, D.; Hola, V.; Di Bonaventura, G.; Djukić, S.; Cirković, I.; Ruzicka, F. Quantification of biofilm in microtiter plates: Overview of testing conditions and practical recommendations for assessment of biofilm production by staphylococci. APMIS 2007, 115, 891–899. [Google Scholar] [CrossRef]
  23. Jolley, K.A.; Bray, J.E.; Maiden, M.C.J. Open-access bacterial population genomics: BIGSdb software, the PubMLST.org website and their applications. Wellcome Open Res. 2018, 3, 124. [Google Scholar] [CrossRef] [PubMed]
  24. Kocianova, S.; Vuong, C.; Yao, Y.; Voyich, J.M.; Fischer, E.R.; DeLeo, F.R.; Otto, M. Key role of poly-gamma-DL-glutamic acid in immune evasion and virulence of Staphylococcus epidermidis. J. Clin. Investig. 2005, 115, 688–694. [Google Scholar] [CrossRef] [PubMed]
  25. Takeuchi, F.; Watanabe, S.; Baba, T.; Yuzawa, H.; Ito, T.; Morimoto, Y.; Kuroda, M.; Cui, L.; Takahashi, M.; Ankai, A.; et al. Whole-genome sequencing of Staphylococcus haemolyticus uncovers the extreme plasticity of its genome and the evolution of human-colonizing staphylococcal species. J. Bacteriol. 2005, 187, 7292–7308. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Flahaut, S.; Vinogradov, E.; Kelley, K.A.; Brennan, S.; Hiramatsu, K.; Lee, J.C. Structural and biological characterization of a capsular polysaccharide produced by Staphylococcus haemolyticus. J. Bacteriol. 2008, 190, 1649–1657. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Keinhörster, D.; Salzer, A.; Duque-Jaramillo, A.; George, S.E.; Marincola, G.; Lee, J.C.; Weidenmaier, C.; Wolz, C. Revisiting the regulation of the capsular polysaccharide biosynthesis gene cluster in Staphylococcus aureus. Mol. Microbiol. 2019, 112, 1083–1099. [Google Scholar] [CrossRef] [Green Version]
  28. Arciola, C.R.; Campoccia, D.; Ravaioli, S.; Montanaro, L. Polysaccharide intercellular adhesin in biofilm: Structural and regulatory aspects. Front. Cell. Infect. Microbiol. 2015, 5, 7. [Google Scholar] [CrossRef] [Green Version]
  29. Foster, T.J. Surface Proteins of Staphylococcus epidermidis. Front. Microbiol. 2020, 11, 1829. [Google Scholar] [CrossRef]
  30. Clarke, S.R.; Harris, L.G.; Richards, R.G.; Foster, S.J. Analysis of Ebh, a 1.1-megadalton cell wall-associated fibronectin-binding protein of Staphylococcus aureus. Infect. Immun. 2002, 70, 6680–6687. [Google Scholar] [CrossRef] [Green Version]
  31. Wang, R.; Khan, B.A.; Cheung, G.Y.; Bach, T.H.; Jameson-Lee, M.; Kong, K.F.; Queck, S.Y.; Otto, M. Staphylococcus epidermidis surfactant peptides promote biofilm maturation and dissemination of biofilm-associated infection in mice. J. Clin. Investig. 2011, 121, 238–248. [Google Scholar] [CrossRef]
  32. Dastgheyb, S.S.; Villaruz, A.E.; Le, K.Y.; Tan, V.Y.; Duong, A.C.; Chatterjee, S.S.; Cheung, G.Y.C.; Joo, H.-S.; Hickok, N.J.; Otto, M. Role of phenol-soluble modulins in formation of Staphylococcus aureus biofilms in synovial fluid. Infect. Immun. 2015, 83, 2966–2975. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Le, K.Y.; Park, M.D.; Otto, M. Immune evasion mechanisms of Staphylococcus epidermidis biofilm infection. Front. Microbiol. 2018, 9, 359. [Google Scholar] [CrossRef] [PubMed]
  34. Pöhlmann-Dietze, P.; Ulrich, M.; Kiser, K.B.; Döring, G.; Lee, J.C.; Fournier, J.M.; Botzenhart, K.; Wolz, C. Adherence of Staphylococcus aureus to endothelial cells: Influence of capsular polysaccharide, global regulator agr, and bacterial growth phase. Infect. Immun. 2000, 68, 4865–4871. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Santajit, S.; Indrawattana, N. Mechanisms of antimicrobial resistance in ESKAPE pathogens. BioMed Res. Int. 2016, 2016, 2475067. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Suwantarat, N.; Carroll, K.C. Epidemiology and molecular characterization of multidrug-resistant Gram-negative bacteria in Southeast Asia. Antimicrob. Resist. Infect. Control 2016, 5, 15. [Google Scholar] [CrossRef] [Green Version]
  37. Pedroso, S.; Sandes, S.; Filho, R.; Nunes, A.C.; Serufo, J.C.; Farias, L.M.; Carvalho, M.; Bomfim, M.; Santos, S.G. Coagulase-negative staphylococci isolated from human bloodstream infections showed multidrug resistance profile. Microb. Drug Resist. 2018, 24, 635–647. [Google Scholar] [CrossRef]
  38. Asante, J.; Amoako, D.G.; Abia, A.; Somboro, A.M.; Govinden, U.; Bester, L.A.; Essack, S.Y. Review of clinically and epidemiologically relevant coagulase-negative staphylococci in Africa. Microbial. Drug Resist. 2020, 26, 951–970. [Google Scholar] [CrossRef] [Green Version]
  39. Ben Zakour, N.L.; Beatson, S.A.; van den Broek, A.H.; Thoday, K.L.; Fitzgerald, J.R. Comparative genomics of the Staphylococcus intermedius group of animal pathogens. Front. Cell. Infect. Microbiol. 2012, 2, 44. [Google Scholar] [CrossRef] [Green Version]
  40. Jantorn, P.; Heemmamad, H.; Soimala, T.; Indoung, S.; Saising, J.; Chokpaisarn, J.; Wanna, W.; Tipmanee, V.; Saeloh, D. Antibiotic Resistance Profile and Biofilm Production of Staphylococcus pseudintermedius Isolated from Dogs in Thailand. Pharmaceuticals 2021, 14, 592. [Google Scholar] [CrossRef] [PubMed]
  41. Du, X.; Zhu, Y.; Song, Y.; Li, T.; Luo, T.; Sun, G.; Yang, C.; Cao, C.; Lu, Y.; Li, M. Molecular analysis of Staphylococcus epidermidis strains isolated from community and hospital environments in China. PLoS ONE 2013, 8, e62742. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Xu, Z.; Cave, R.; Chen, L.; Yangkyi, T.; Liu, Y.; Li, K.; Meng, G.; Niu, K.; Zhang, W.; Tang, N.; et al. Antibiotic resistance and molecular characteristics of methicillin-resistant Staphylococcus epidermidis recovered from hospital personnel in China. J. Glob. Antimicrob. Resist. 2020, 22, 195–201. [Google Scholar] [CrossRef] [PubMed]
  43. Seng, R.; Leungtongkam, U.; Thummeepak, R.; Chatdumrong, W.; Sitthisak, S. High prevalence of methicillin-resistant coagulase-negative staphylococci isolated from a university environment in Thailand. Int. Microbiol. 2017, 20, 65–73. [Google Scholar] [CrossRef] [PubMed]
  44. Miragaia, M.; Thomas, J.C.; Couto, I.; Enright, M.C.; de Lencastre, H. Inferring a population structure for Staphylococcus epidermidis from multilocus sequence typing data. J. Bacteriol. 2007, 189, 2540–2552. [Google Scholar] [CrossRef] [Green Version]
  45. Mendes, R.E.; Deshpande, L.M.; Costello, A.J.; Farrell, D.J. Molecular epidemiology of Staphylococcus epidermidis clinical isolates from U.S. hospitals. Antimicrob. Agents Chemother. 2012, 56, 4656–4661. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Deplano, A.; Vandendriessche, S.; Nonhoff, C.; Dodemont, M.; Roisin, S.; Denis, O. National surveillance of Staphylococcus epidermidis recovered from bloodstream infections in Belgian hospitals. J. Antimicrob. Chemother. 2016, 71, 1815–1819. [Google Scholar] [CrossRef] [Green Version]
  47. Cremniter, J.; Sivadon-Tardy, V.; Caulliez, C.; Bauer, T.; Porcher, R.; Lortat-Jacob, A.; Piriou, P.; Judet, T.; Aegerter, P.; Herrmann, J.-L.; et al. Genetic analysis of glycopeptide-resistant Staphylococcus epidermidis strains from bone and joint infections. J. Clin. Microbiol. 2013, 51, 1014–1019. [Google Scholar] [CrossRef] [Green Version]
  48. Both, A.; Huang, J.; Qi, M.; Lausmann, C.; Weißelberg, S.; Büttner, H.; Lezius, S.; AFailla, V.; Christner, M.; Stegger, M.; et al. Distinct clonal lineages and within-host diversification shape invasive Staphylococcus epidermidis populations. PLoS Pathog. 2021, 17, e1009304. [Google Scholar] [CrossRef]
  49. Skachkova, T.S.; Zamyatin, M.N.; Orlova, O.A.; Yumtsunova, N.A.; Lashenkova, N.N.; Fomina, V.S.; Gusarov, V.G.; Mikhaylova, Y.V.; Shelenkov, A.A.; Goloveshkina, E.N.; et al. Monitoring methicillin-resistant staphylococcus strains in the Moscow medical and surgical center using molecular-biological methods. Epidemiol. Vaccinal Prev. 2021, 20, 44–50. [Google Scholar] [CrossRef]
  50. Lee, J.Y.H.; Monk, I.R.; Gonçalves da Silva, A.; Seemann, T.; Chua, K.Y.L.; Kearns, A.; Hill, R.; Woodford, N.; Bartels, M.D.; Strommenger, B.; et al. Global spread of three multidrug-resistant lineages of Staphylococcus epidermidis. Nat. Microbiol. 2018, 3, 1175–1185. [Google Scholar] [CrossRef]
  51. Arora, S.; Li, X.; Hillhouse, A.; Konganti, K.; Little, S.V.; Lawhon, S.D.; Threadgill, D.; Shelburne, S.; Hook, M. Staphylococcus epidermidis MSCRAMM SesJ is encoded in composite islands. mBio 2020, 11, e02911-19. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Species distribution of staphylococci isolates from hospitalized patients (A), outpatients (B), and pets (C). Colors: blue—S. aureus, light red—S. epidermidis, green—S. haemolyticus, violet—S. hominis, yellow—S. pseudintermedius, pink—S.coagulans, grey—other staphylococci species. The number of isolates for each species is represented in the circle segments.
Figure 1. Species distribution of staphylococci isolates from hospitalized patients (A), outpatients (B), and pets (C). Colors: blue—S. aureus, light red—S. epidermidis, green—S. haemolyticus, violet—S. hominis, yellow—S. pseudintermedius, pink—S.coagulans, grey—other staphylococci species. The number of isolates for each species is represented in the circle segments.
Microorganisms 09 02487 g001
Figure 2. Resistance of predominant Staphylococcus species. MSS isolates are marked with green, MDR isolates with rose, and MRS isolates with violet. The number of corresponding isolates is shown on chart columns. 1 Chart column marked with diagonal stripes correspond to four MDR S. pseudintermedius isolates sensitive to cefoxitin.
Figure 2. Resistance of predominant Staphylococcus species. MSS isolates are marked with green, MDR isolates with rose, and MRS isolates with violet. The number of corresponding isolates is shown on chart columns. 1 Chart column marked with diagonal stripes correspond to four MDR S. pseudintermedius isolates sensitive to cefoxitin.
Microorganisms 09 02487 g002
Figure 3. Biofilm formation by MDR S. epidermidis and S. haemolyticus isolates.
Figure 3. Biofilm formation by MDR S. epidermidis and S. haemolyticus isolates.
Microorganisms 09 02487 g003
Table 1. The origin of Staphylococcus isolates.
Table 1. The origin of Staphylococcus isolates.
PanelsOrigin of StrainsNumber of Isolates
Strains isolated from hospitalized patients (n = 105)Prosthetic joint and post-operative wound infections56
Purulent diabetic ulcers43
Decubitus ulcers3
Sputum from a patient with pneumonia2
Cerebrospinal fluid1
Strains from outpatients (n = 228)Respiratory tract infections83
Skin and mucosal infections59
Urinary tract infections21
Bacterial vaginosis25
Fecal samples from patients with diarrhea36
Purulent diabetic ulcers4
Strains from pets (39 from dogs and 22 from cats) (n = 61)Skin and mucosal infections57
Fecal samples from animals with diarrhea4
Table 2. Species identification of clinical and veterinary staphylococci isolates.
Table 2. Species identification of clinical and veterinary staphylococci isolates.
NoSpecies (Number of Isolates)Isolation Source (Number of Isolates)GenBank Identifier for 16S rRNA Gene
1Staphylococcus aureus (167)Hospital (74)No ID
Outpatient (81)No ID
Veterinary (12)No ID
2Staphylococcus auricularis (4)Outpatient (4)MZ014399–MZ014402
3Staphylococcus borealis (2)Outpatient (2)MZ014403, MZ014404
4Staphylococcus capitis (4)Outpatient (4)MZ014405–MZ014407, MZ014411
5Staphylococcus caprae (5)Hospital (1)MZ014409
Outpatient (4)MZ014408, MZ014410, MZ014412, MZ014413
6Staphylococcus carnosus (2)Outpatient (2)MZ014414, MZ014415
7Staphylococcus casei (2)Outpatient (1)MZ014416
8Staphylococcus coagulans/S. schleiferi subsp. coagulans (8)Veterinary (8)MW979964–MW979970
9Staphylococcus cohnii (2)Outpatient (2)MZ014417, MZ014418
10Staphylococcus devriesei (4)Outpatient (4)MZ014419–MZ014422
11Staphylococcus equorum (3)Outpatient (1)MZ014423
Veterinary (2)MW979955, MW979971
12Staphylococcus epidermidis (104)Hospital (20)MZ027385–MZ027397, MZ040881, MZ040882, MZ041685
Outpatient (84)MZ027349–MZ027358, MZ027360–MZ027362, MZ027364–MZ027366, MZ040893–MZ040916, MZ047203–MZ047210, MZ049531–MZ049537,MZ723062
13Staphylococcus felis (5)Veterinary (5)MW979972–MW979976
14Staphylococcus haemolyticus (22)Hospital (2)MZ027398, MZ027399
Outpatient (19)MZ027367–MZ027375, MZ040917–MZ040921, MZ047213–MZ047216, MZ723059, MZ723059
Veterinary (1)MW979956
15Staphylococcus hominis (14)Hospital (2)MZ014434, MZ014435
Outpatient (11)MZ014424–MZ014433
Veterinary (1)MW979977
16Staphylococcus pasteuri (1)Outpatient (1)MZ047217
17Staphylococcus pseudintermedius (24)Veterinary (24)MW965646, MW979953, MW979978–MW979995, MW979997, MW979998
18Staphylococcus simulans (8)Hospital (1)MZ014437
Outpatient (5)MZ014436, MZ014438–MZ014441
Veterinary (2)MW979999, MW980000
19Staphylococcus warneri (8)Hospital (2)MZ014445, MZ014447
Outpatient (5)MZ014442–MZ014444, MZ014446, MZ014448
Veterinary (1)MW979957
20Mammaliicoccus sciuri/Staphylococcus sciuri (3) 1Hospital (1)MW965541
Veterinary (2)MW965545, MW965551
21Mammaliicoccus vitulinus/Staphylococcus vitulinus (1) 1Veterinary (1)MW965552
22Mammaliicoccus lentus/Staphylococcus lentus (1) 1Veterinary (1)MW965543
1 The species were moved to the newly established Mammaliicoccus genus [21], previous species names are separated by slash.
Table 3. Antimicrobial resistance profiles and resistance-encoding genes of prevailing Staphylococcus isolates.
Table 3. Antimicrobial resistance profiles and resistance-encoding genes of prevailing Staphylococcus isolates.
AntimicrobialResistance-Encoding
Gene(s) 1
Number of Isolates with Resistance Profiles and ARGs That Were Identified in These Isolates
S. aureus
Hospital (n = 74)Outpatient (n = 81)Veterinary (n = 12)
Cefoxitin 1320
mecA1320
Penicillin N.d 2N.d.N.d.
blaZ32383
Gentamicin and/or Amikacin 1220
aac(6′)-Ie-aph(2″)-Ia810
aac(6′)-Ie-aph(2″)-Ia + ant(4′)-Ia110
aac(6′)-Ie-aph(2″)-Ia + aph(3′)-IIIa300
Erythromycin 050
ermC050
Erythromycin and Clindamycin 750
ermA500
ermC250
S. epidermidis
Hospital (n = 20)Outpatient (n = 84)Veterinary (n = 0)
Cefoxitin 13180
mecA13180
Penicillin N.d.N.d.N.d.
blaZ6390
Gentamicin and/or Amikacin 9190
aac(6′)-Ie-aph(2″)-Ia4110
ant(4′)-Ia030
aac(6′)-Ie-aph(2″)-Ia + ant(4′)-Ia520
aac(6′)-Ie-aph(2″)-Ia + aph(3′)-IIIa030
Erythromycin 6340
ermC250
msrA4290
Clindamycin 1 31 30
Erythromycin and Clindamycin 420
ermA300
ermC100
msrA020
S. haemolyticus
Hospital (n = 2)Outpatient (n = 19)Veterinary (n = 1)
Cefoxitin 2101
mecA2101
Penicillin N.d.N.d.N.d.
blaZ130
Gentamicin and/or Amikacin 2101
aac(6′)-Ie-aph(2″)-Ia171
ant(4′)-Ia04 40
aac(6′)-Ie-aph(2″)-Ia + aph(3′)-IIIa130
Erythromycin 2130
ermC030
msrA2101 4
S. hominis
Hospital (n = 2)Outpatient (n = 19)Veterinary (n = 1)
Cefoxitin 020
mecA020
Penicillin N.d.N.d.N.d.
blaZ130
Gentamicin and/or Amikacin 020
ant(4′)-Ia01 40
aph(3′)-IIIa001 4
aac(6′)-Ie-aph(2)-Ia + ant(4′)-Ia020
Erythromycin 050
msrA040
ermA + msrA010
S. pseudintermedius
Hospital (n = 0)Outpatient (n = 0)Veterinary (n = 24)
Penicillin N.d.N.d.N.d.
blaZ004
Gentamicin and/or Amikacin 006
aac(6′)-Ie-aph(2″)-Ia001
aph(3′)-IIIa001
aac(6′)-Ie-aph(2″)-Ia + aph(3′)-IIIa004
Erythromycin and Clindamycin 007 3
1mecA, aac(6′)-Ie-aph(2″)-Ia, aph(3′)-IIIa, ant(4′)-Ia, ermA/ermC, and msrA genes were checked in all isolates; 2 not determined; 3 isolates had clindamycin or erythromycin resistant profiles, but investigated genes were not found; 4 isolates contained genes, encoding aminoglycoside resistance, but had no aminoglycoside-resistant phenotype.
Table 4. Source of isolation, resistance characteristics, and sequence types of MDR S. epidermidis and MDR S. haemolyticus isolates.
Table 4. Source of isolation, resistance characteristics, and sequence types of MDR S. epidermidis and MDR S. haemolyticus isolates.
NoSpeciesCEMTC No of Isolate/16S rRNA GenBank IDSource of IsolationResistanceSequence TypeClonal Complex
1S. epidermidis1827/MZ027386Hospital (swab from purulent diabetic ulcer)FOX, AK, CN, E, DAST 23CC23
2S. epidermidis1833/MZ027390Hospital (pure culture, post-operative wound infection)FOX, AK, CN, E, DAST 23CC23
3S. epidermidis2300/MZ027350Outpatient (nasal swab, rhinitis)FOX, AK, CNST 20S 1
4S. epidermidis3054/MZ027389Hospital (pure culture, post-operative wound infection)FOX, CN, E, DAST 5CC5
5S. epidermidis3107/MZ027352Outpatient (nasal swab, rhinitis)FOX, CN, EST 152CC5
6S. epidermidis3117/MZ027364Outpatient (faeces, diarrhea)FOX, CN, EST 152CC5
7S. epidermidis3750/MZ027394Hospital (biopsy material, prosthetic joint infection)FOX, AK, CNST 210CC5
8S. epidermidis3824/MZ027395Hospital (biopsy material, prosthetic joint infection)FOX, AK, CN, EST 786CC5
9S. epidermidis3918/MZ027396Hospital (pure culture, post-operative wound infection)FOX, AK, CN, EST 20S 1
10S. epidermidis4066/MZ027397Hospital (pure culture, post-operative wound infection)FOX, AK, CN, E, DAST 17CC5
11S. haemolyticus1553/MZ027371Outpatient (urine sample, pyelonephritis)FOX, AK, CN, EST 1CC3
12S. haemolyticus2119/MZ027368Outpatient (purulent diabetic ulcer)FOX, CN, EST 3CC3
13S. haemolyticus2688/MZ723059Outpatient (faeces, diarrhea)FOX, CN, EST 3CC3
14S. haemolyticus3413/MZ027399Hospital (sputum, ventilator-associated pneumonia)FOX, AK, CN, EST 42CC3
15S. haemolyticus3601/MZ027370Outpatient (skin infection, skin scraping)FOX, AK, CN, EST 8CC3
1 Singleton.
Table 5. Antibiotic resistance genes revealed in the genomes of MDR S. epidermidis and S. haemolyticus isolates.
Table 5. Antibiotic resistance genes revealed in the genomes of MDR S. epidermidis and S. haemolyticus isolates.
Antibiotic ClassesARGS. epidermidis Isolates 1S. haemolyticus Isolates 1
182718332300305431073117375038243918406615532119268834133601
PenicillinsblaZ 2
CephalosporinsmecA
mecI
mecR1
Macrolides and LincosamidesermA
ermC
msrA
mphC
Aminoglycosidesaph(3′)-IIIa
ant(4′)-Ia
aac(6’)-Ie-aph(2″)-Ia
aad(6)
ant(9)
QuinolonesnorA
qacA
qacB
Chlorampheni-
col
catA8
catA7
Fusidic acidfusC
Diaminopyrimi-
dines
dfrC
dfrG
AminocoumarinsAminocoumarin resistant gyrB
LincosamideslnuA
Tetracyclinestet(K)
tet(45)
FosfomycinfosB
BacitracinBceSR
1S. epidermidis isolates are marked with blue, S. haemolyticus isolates are marked with green; 2 the “tick” symbol means that corresponding antibiotic resistance gene was revealed in the isolate.
Table 6. Virulence factors identified in the genomes of S. epidermidis and S. haemolyticus isolates.
Table 6. Virulence factors identified in the genomes of S. epidermidis and S. haemolyticus isolates.
Species and No in CEMTCExoenzymesToxinsImmune Evasion
Cysteine ProteaseSerine V8 ProteaseLipasesThermonucleaseHemo-
lysin III
Hemolysin, Containing CBS Domainβ-HemolysinPGA, Capsule
S. epidermidis 1827sspBsspAgeh, lipnuc+ 1+ 1hlbcapA, pgsABCDE 2
S. epidermidis 1833sspBsspAgeh, lipnuc++hlbcapA, pgsABCDE
S. epidermidis 2300sspBsspAgeh, lipnuc++hlbcapA, pgsABCDE
S. epidermidis 3054sspBsspAgeh, lipnuc++hlbcapA, pgsABCDE
S. epidermidis 3107sspBsspAgeh, lipnuc++hlbcapA, pgsABCDE
S. epidermidis 3117sspBsspAgeh, lipnuc++hlbcapA, pgsABCDE
S. epidermidis 3750sspBsspAgeh, lipnuc++hlbcapA, pgsABCDE
S. epidermidis 3824sspBsspAgeh, lipnuc++hlbcapA, pgsABCDE
S. epidermidis 3918sspBsspAgeh, lipnuc++hlbcapA, pgsABCDE
S. epidermidis 4066sspBsspAgeh, lipnuc++hlbcapA, pgsABCDE
S. haemolyticus 1553 lipnuc++ capA, pgsABCDE
S. haemolyticus 2119sspBsspAgeh, lipnuc++ capA, pgsABCDE, cap5L, cap5F, cap8C, cap8E, cap8I, cap8H, cap 8M, cap8N2
S. haemolyticus 2688 lipnuc++ capA, pgsABCDE, cap5A
S. haemolyticus 3413 lipnuc++ capA, pgsABCDE, cap5L, cap5F, cap8C, cap8E, cap8I, cap8H, cap 8M, cap8N2
S. haemolyticus 3601 lipnuc++ capA, pgsABCDE, cap5A, cap8C3
1 the “plus” symbol means that gene, encoding the appropriate hemolysin, was revealed in the isolate; 2 pgsABCDE cluster responsible for synthesis of the surface-attached poly-gamma-glutamate (PGA); 3 putative capsule operon, possessing ≥50% identity to S. aureus cap5 or cap8 locus.
Table 7. Adherence and biofilm formation factors, identified in the genomes of MDR S. epidermidis and S. haemolyticus isolates.
Table 7. Adherence and biofilm formation factors, identified in the genomes of MDR S. epidermidis and S. haemolyticus isolates.
Species and
No in CEMTC
Ability to Form BiofilmAdherence and Biofilm Formation Genes
PIAPSMAapPlsSas-Family ProteinsAutoly-
sin
Fibronectin Binding Protein 1Elastin Binding ProteinSdr-Family Proteins
S. epidermidis 1827strongicaADBC, icaRpsmβ1, psmβ2 plsSasA, SasC, SasFatlEebhebpsdrG, sdrH
S. epidermidis
1833
moderateicaADBC, icaR plsSasA, SasC, SasFatlEebhebpsdrG
S. epidermidis
2300
strongicaADBC, icaRpsmβ1, psmβ2aapplsSasA, SasC, SasFatlEebhebpsdrC, sdrF, sdrH
S. epidermidis
3054
strong psmβ1, psmβ2aap SasA, SasC, SasFatlEebhebpsdrF, sdrG, sdrH
S. epidermidis
3107
strong psmβ1, psmβ2aap SasA, SasC, SasFatlEebhebpsdrF, sdrG, sdrH
S. epidermidis
3117
strong psmβ1, psmβ2aap SasA, SasC, SasFatlEebhebpsdrF, sdrG, sdrH
S. epidermidis
3750
moderate psmβ1, psmβ2 plsSasA, SasC, SasFatlEebhebpsdrE, sdrF, sdrG, sdrH
S. epidermidis
3824
strongicaADBC, icaR aapplsSasA, SasC, SasFatlE ebpsdrE, sdrF, sdrG, sdrH
S. epidermidis
3918
strongicaADBC, icaR aapplsSasA, SasC, SasF, SasGatlEebhebpsdrE, sdrF, sdrG, sdrH
S. epidermidis
4066
strongicaADBC, icaRpsmβ1, psmβ2aap SasA, SasC, SasFatlE ebpsdrF, sdrG, sdrH
S. haemolyticus
1553
strong aap SasA, SasC, SasF, SasGatlE ebp
S. haemolyticus
2119
moderate SasA, SasC, SasF, SasGatlE ebp
S. haemolyticus
2688
strong SasA, SasC, SasF, SasGatlE ebp
S. haemolyticus
3413
weak aap SasA, SasC, SasF, SasGatlE ebpsdrC
S. haemolyticus
3601
moderate SasA, SasC, SasF, SasGatlE ebpsdrC
1 Extracellular matrix-binding protein (Embp) homologue [30].
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bardasheva, A.; Tikunov, A.; Kozlova, Y.; Zhirakovskaia, E.; Fedorets, V.; Fomenko, N.; Kalymbetova, T.; Chretien, S.; Pavlov, V.; Tikunova, N.; et al. Antibiotic Resistance and Pathogenomics of Staphylococci Circulating in Novosibirsk, Russia. Microorganisms 2021, 9, 2487. https://doi.org/10.3390/microorganisms9122487

AMA Style

Bardasheva A, Tikunov A, Kozlova Y, Zhirakovskaia E, Fedorets V, Fomenko N, Kalymbetova T, Chretien S, Pavlov V, Tikunova N, et al. Antibiotic Resistance and Pathogenomics of Staphylococci Circulating in Novosibirsk, Russia. Microorganisms. 2021; 9(12):2487. https://doi.org/10.3390/microorganisms9122487

Chicago/Turabian Style

Bardasheva, Alevtina, Artem Tikunov, Yuliya Kozlova, Elena Zhirakovskaia, Valeriya Fedorets, Natalya Fomenko, Tatyana Kalymbetova, Svetlana Chretien, Vitaliy Pavlov, Nina Tikunova, and et al. 2021. "Antibiotic Resistance and Pathogenomics of Staphylococci Circulating in Novosibirsk, Russia" Microorganisms 9, no. 12: 2487. https://doi.org/10.3390/microorganisms9122487

APA Style

Bardasheva, A., Tikunov, A., Kozlova, Y., Zhirakovskaia, E., Fedorets, V., Fomenko, N., Kalymbetova, T., Chretien, S., Pavlov, V., Tikunova, N., & Morozova, V. (2021). Antibiotic Resistance and Pathogenomics of Staphylococci Circulating in Novosibirsk, Russia. Microorganisms, 9(12), 2487. https://doi.org/10.3390/microorganisms9122487

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