1. Introduction
Staphylococcus aureus is a Gram-positive bacterium inhabiting healthy individuals’ nostrils and skin. However, it has become an important opportunistic pathogen in communities and hospitals [
1]. It causes severe skin infections, pneumonia, endocarditis, and bloodstream infections (BSIs) [
2]. BSIs caused by
S. aureus infections have high morbidity and mortality if not treated timeously [
3]. The most significant risk factors for
S. aureus BSIs are intravascular devices, surgical procedures, and a debilitated immune system [
4].
Methicillin-resistant
S. aureus (MRSA) has become a significant cause of BSIs. MRSA poses a major public health threat because of multidrug resistance to different antibiotic classes that limit treatment options [
5]. Resistance is mediated by the
mecA gene, found on a mobile genetic element (MGE) known as the staphylococcal cassette chromosome
mec (SCC
mec) [
6]. Methicillin-susceptible
S. aureus (MSSA) is also emerging as a causative agent of BSIs [
7] and has been reported to display high virulence and multidrug resistance [
8].
The pathogenicity of
S. aureus depends on its ability to produce a wide array of virulence factors involved in adhesion, invasion of host tissues, immune system evasion, and biofilm formation [
9,
10]. Virulence factors and multiple resistance genes can be transmitted by horizontal gene transfer (HGT) [
11] on diverse MGEs, amongst which plasmids are reported as the primary sources for dissemination [
4].
The epidemiology of
S. aureus strains indicates that its molecular characteristics continually change over time, resulting in new clones, which vary by region. In a study in the United States, ST5 and ST8 were the most prevalent sequence types [
12]. In South Africa, ST612 is dominant in the hospital environment [
13]. The ST612-IV [2B], belonging to spa type t1257, was identified as a typical clone in clinical settings [
14] and sporadically in poultry settings [
15]. The ST5 and ST8 clones are commonly associated with BSIs and the pandemic lineages of
S. aureus, such as the clonal complex CC8 and CC5 [
16]. Notably, the sequence types ST612, ST5, ST8, and ST72 have displayed high resistance to most antibiotic drug classes and are challenging to treat [
16].
Multidrug-resistant (MDR)
S. aureus infections pose a serious clinical concern. A high incidence of pathogenic MDR MRSA has been reported, and the data suggest that its prevalence is increasing in Africa [
17]. A recent South African study investigating the genetic relatedness of hospital-acquired-associated MRSA isolates in two hospitals revealed that all isolates were resistant to aminoglycosides and β-lactams. All the isolates carried the
aacA-aphD and
mecA-resistant genes and clusters of virulence genes [
18]. We elaborate on the genomic characteristics, antibiotic resistance, virulence, pathogenicity, phylogenomics, and clonal diversity of six
S. aureus clinical strains implicated in bloodstream infections at a regional hospital in the KwaZulu-Natal province in South Africa
3. Discussion
We studied the genomic characteristics of six MDR S aureus isolates implicated in BSIs. The study analysed the resistome, virulome, mobilome, phylogeny, and genetic environment of the resistance genes using WGS and bioinformatics. The genomes analysed herein were predominantly recovered from patients ≤1-year-old.
There was a diversity of ARGs encoding resistance to different antibiotics and good concordance between the observed phenotypic and genotypic resistance. The incidence of ARG’s encoding resistance to β-lactams, aminoglycosides, macrolides, fosfomycin, trimethoprim, tetracycline, and genes coding multidrug resistance (MDR) efflux pumps (
norA,
mepR,
arlR,
mgrA, and
lmrS) was not dependent on the clonal type. Resistance genes found in this study included
erm genes which mediate macrolide–lincosamide–streptogramin B (MLSB) resistance, the aminoglycoside-modifying enzymes, MFS antibiotic efflux pump (
msrA,
norA,
norC), which also confers resistance to fluoroquinolones [
25]. The
erm(C) and
erm(A) genes that are commonly found in MLS
B-resistant
S. aureus were found in erythromycin and clindamycin resistant isolates (
Table 2), which was expected since resistance to erythromycin co-selects resistance to other antibiotics, such as streptogramin B (MLS
B) and lincosamides [
26]. The
ermC gene is among the primary
erm type which facilitates ribosome methylation of the 23S rRNA, triggering conformational changes, resulting in drug binding inhibition [
27], and have been reported in clinical
S. aureus isolates from South Africa [
28]. In this study, the
ermC encoding macrolide resistance was carried on a plasmid, on a contig that had the closest nucleotide homology to plasmids from
S. epidermidis strain TMDU-137 plasmid p5, complete sequence (CP093178.1), implying the likelihood of horizontal transfer of
ermC genes in clinical
S. aureus isolates. The
ermC are often plasmid-mediated, resulting in high resistance to macrolides in
S. aureus [
29].
The
blaZ gene which inactivates penicillin through hydrolysis of the beta-lactam ring was observed in all six isolates that were phenotypically resistant to penicillin. The
blaZ genes have also been isolated in clinical isolates of Staphylococci in South Africa [
30]. In this study, the
blaZ genes were found on contigs with closest homology to either chromosomes or plasmids. This agrees with a study done in Spain that analysed ARGs presence in chromosomes and plasmids from the genomes of
S. aureus. WGS analysis of
S. aureus revealed that
blaZ (n=2) were located on chromosomic contigs, while
blaZ was found in plasmid contigs in three isolates [
31]. It is important to note that most
blaZ and associated MGEs from isolates belonging to ST5 (S13, S34) isolated from the intensive care unit (ICU), and paediatric ward (S29) were located on contigs that had the closest homology to plasmids, implying that plasmids play a crucial role in mobilizing the
blaZ gene in clinical
S. aureus isolates. The S29 isolate belonging to t045-CC5 lineage, carried assortment of ARGs encoding resistance to different antibiotics (Table 4). Similar ARGs in MRSA lineage t045-CC5-MRSA were also reported in a study conducted in South Africa, where t045-CC5 MRSA lineages obtained from different clinical samples from South Africa and Nigeria reported that t045 lineages were MDR, suggesting that this lineage is hospital-associated, and their multidrug resistance nature may compromise treatment [
32].
Also, the
blaZ genes, heavy metal genes and associated MGEs were carried on either plasmid or chromosome. The
blaZ and
cadAC genes were found on the genetic element recombinase:
blaI:
blaR1:
blaZ:
cadC:
cadA for isolates S24 (MSSA) that was from the ICU, and S29 (MRSA) from the paediatric ward, suggesting co-selection of heavy metal resistance dissemination and adaptation in different wards. The
cadA gene confers a high resistance to cadmium and other heavy metals like zinc and lead in
S. aureus isolates [
33]. The
cadA was associated with a plasmid, similar to the findings of a study that was done by Al-Trat et al. (2023) in Malaysia who used WGS to analyse the plasmid content of clinical MRSA isolates, and reported that heavy metal resistance plasmids harboured cadmium resistance genes with the majority being
cadAC [
34]. The HMRGs have been reported to trigger co-selection mechanism with antibiotics, which may complicate treatment [
35]. This may pose a challenge especially among patients in the ICU where broad-spectrum antibiotics are often used.
Tetracycline resistance genes (
tetK, tetM) were observed in two isolates. Isolate S11 carried
tet(K) associated with a genetic context: plasmid recombination:
tet(K) that had a high similarity to
S. epidermidis BPH0662, plasmid: 1 (LT614820.1), which could be significant in mobilising TET resistant genes. Also, the
tet(M) was bracketed by integrase and IS256 in isolate S31. The IS
256 is a retrotransposon that can mobilize the resistance genes through a copy-and-paste mechanism and has been shown to confer a robust genomic plasticity in MRSA strains [
36]. We found mutations in
gyrA, gyrB parC,
parE, and
rpoB, implicated in fluoroquinolone and rifampicin resistance. However, the detection of double mutations of gyrA/parC could be associated with high levels of resistance to fluoroquinolones [
37].
We found that ARGs, and virulence genes were associated with MGEs, which may enable their transfer within and between plasmids and chromosomes [
38]. In this study, the
mecA gene was located on IS1182 in two isolates, surrounded by recombinase in genetic context
mecA:mecR1::IS1182::recombinase. The insertion sequence IS1182 was present in 2/3 MRSA strains that contained
mecA. IS1182 has been shown to occur close to the SCC
mec element and increase resistance through inactivating
lytH gene encoding a putative lytic enzyme in pathogenic MRSA isolates [
39].
MLST typing, clonal complex,
spa typing, and SCC
mec typing were used to analyse the molecular characteristics of the
S. aureus isolates. Four ST types and two clonal clusters (CCs) were found among the six clinical isolates in this study, with ST5, the most predominant complex clonal CC5 and CC8. Generally, clonal lineages ST5, ST8, ST152, and ST612 are among the most commonly reported in hospital environments, along with other sequence types of
S. aureus [
40].
S. aureus ST5, belonging to CC5, was predominant in this study and was previously reported among patients with bloodstream infections at Ruijin Hospital in Shanghai [
3]. The detection of clonal complexes CC5 and CC8 agrees with a study by Smith et al. [
16], which also found CC8, and CC5 were predominant in a study that analysed the genomic epidemiology of MRSA and MSSA from bloodstream infections in the USA. Their results revealed that the MDR phenotype observed in strains belonging to CC5 and CC8 was responsible for the occurrence of multidrug and methicillin resistance in the
S. aureus population. MRSA strains belonging to CC8 and CC5 are frequently associated with global outbreaks and have been identified in Africa [
41].
The
spa typing revealed six different spa types, suggesting a non-clonal MRSA and MSSA distribution. The detection of
spa types t1257, t045 and t355 agrees with a study conducted in South Africa, which analysed the diversity of SCC
mec elements and spa types in
S. aureus isolates from blood culture in the Gauteng, KwaZulu-Natal, Free State, and Western Cape provinces [
14], in which t037 and t1257 were the most common and predominated throughout the seven-year study period. In this study, some antibiotic resistance genes were associated with specific MRSA clones belonging to
spa types t1257 and t045, t9475. Shittu et al. (2021) found the
spa types t045 and t1257 to be the most prevalent and associated with genes conferring resistance to aminoglycosides, trimethoprim, macrolides and tetracycline in clinical isolates of
S. aureus from South Africa and Nigeria.
The analysis of SCC
mec types revealed the presence of SCC
mec type IVd (2B) and SCC
mec type I (B) carrying the
mecA gene, which occurred in tandem with
mecR1 in both isolates. However, one MRSA (S11) isolate had a non-typeable SCC
mec element cassette due to the missing cassette chromosome recombinase (
ccr) gene complex [
42]. The
ccr gene complex is an essential component required to facilitate the integration or excision of the SCC
mec element in the staphylococcal chromosome, and their loss has also been reported [
43]. The SCC
mec IV detected in our study is associated with the spa type t1257, previously reported in South Africa in
S. aureus obtained from poultry isolates [
15], implying its possible transfer between humans and animals.
We found different MRSA genotypes ST612-t1257-CC8, ST8-t9475-CC8, and ST5-t045-CC5, suggesting that MRSA isolates were not clonally and epidemiologically related. The ST612-t1257-CC8 identified in this study, is an endemic MRSA clone which have been reported in animal and clinical settings [
14,
15]. The ST5-I-MRSA, known as the pandemic British EMRSA-3 clone, was detected in the paediatric ward. This is similar to a study conducted in South Africa, where the t045-MRSA strain occurred in paediatric patients [
18]. The isolation of t045-ST5-MRSA strain could confirm its successful persistence in the hospital and its capacity to cause infections in neonatal and paediatric wards [
44].
Several virulence factors, including adherence, immune invasion, toxins, and exoenzymes associated with invasive infections, were detected in our isolates. The virulence genes encoding clumping factor proteins (
clfA and
clfB) are involved in the pathogenesis of
S. aureus, including bacteremia [
9]. Consistent with pathogenic
S. aureus strains isolated in various environments globally, our isolates were characterized by
icaADBC operon and (
sdrC,
sdrD,
sdrE) involved in biofilm-forming genes [
45]. Most strains haboured genes, including the alpha and gamma-hemolysin genes (
hlgA,
hlgB,
hlgC,
hly/hla,
hlb), and the
ica operon associated with pathogenicity and adhesion. Additionally, our isolates were characterised by various toxins, including
lukE/D genes, and panton-valentine leukocidin (PVL)
lukS-PV/
lukF-PV genes in one MSSA, and MRSA strains. The expression of these PVL toxin genes in
S. aureus isolates, lyses host cells and promotes virulence of the bacteria [
46], which might worsen the outcomes
S. aureus infection. Consistent with clinical
S. aureus strains, our isolates were characterised by a capsular polysaccharide (CP) serotype 8, which shields the bacterial pathogen from host immune defence mechanisms associated with increased virulence in BSIs [
47].
Most virulence genes including those encoding SEs,
sak,
hlg,
luk,
scn clfA,
sbi, and associated MGEs were carried on chromosomes in the majority of isolates. The
ica gene operon and
vra genes were found to be associated with ISs (IS30, IS6) and recombinase for S11 (ST8) isolate from the surgical ward. The
ica genes
vraDEH genes have been shown to play an important role in biofilm formation [
48], daptomycin resistance in
S. aureus [
49]. which could enhance antibiotic resistance traits and chronic infection. The occurrence of ST8-t9475 MRSA strains co-harbouring
ica genes and genes encoding daptomycin resistance in ST8 MRSA could be advantageous to the ST8-t9475 colonization, invasion, and survival in the surgical ward. The virulence genes encoding SEs,
eap, scn, sak, sph,
lukH, and
cadA, were found on a contig that had high sequence similarity to
S. aureus strain ER10678.3 plasmid pER10678.3A.1 (CP051928.1), implying that they are mobilized by plasmids. Virulence genes, including those encoding hla/hld, toxin production, and biofilm formation, are plasmid-mediated [
50], thus could easily facilitate their transfer resulting in highly pathogenic strains that may be difficult to treat.
Phylogenomic analyses revealed that the clinical isolates in this study clustered mainly with clinical isolates from hospital patients (
Figure 2). ST5 study isolates were closely related to clinical isolates from South Africa suggesting possible dissemination of ST5 strains and adaptation in hospital environments. Furthermore, ST152 isolate was closely related to ST152 strains from Egypt and Ghana, implying a possible spread and epidemiological linkage between these isolates. ST152-PVL-producing
S. aureus isolates is particularly frequent and widespread in West and Central Africa [
51], and livestock [
52]. The ST152- PVL-positive MSSA, has also been reported from cutaneous abscesses among mine workers at a gold mine in Gauteng, South Africa [
53]. Identifying ST152 in livestock and humans suggests animal–human transmission, which requires further investigation. ST8 and ST612 isolates were closely related to ST8 isolated from Tanzania, indicating that ST612 is a double-locus variant of ST8. ST8 and ST612 isolates are potentially multidrug-resistant and highly virulent strains associated with hospital outbreaks [
54].