- Research article
- Open Access
Genome-wide association studies reveal candidate genes associated to bacteraemia caused by ST93-IV CA-MRSA
BMC Genomics volume 22, Article number: 418 (2021)
The global emergence of community-associated methicillin-resistant Staphylococcus aureus (CA-MRSA) has seen the dominance of specific clones in different regions around the world with the PVL-positive ST93-IV as the predominant CA-MRSA clone in Australia. In this study we applied a genome-wide association study (GWAS) approach on a collection of Australian ST93-IV MRSA genomes to screen for genetic traits that might have assisted the ongoing transmission of ST93-IV in Australia. We also compared the genomes of ST93-IV bacteraemia and non-bacteraemia isolates to search for potential virulence genes associated with bacteraemia.
Based on single nucleotide polymorphism phylogenetics we revealed two distinct ST93-IV clades circulating concurrently in Australia. One of the clades contained isolates primarily isolated in the northern regions of Australia whilst isolates in the second clade were distributed across the country. Analyses of the ST93-IV genome plasticity over a 15-year period (2002–2017) revealed an observed gain in accessory genes amongst the clone’s population. GWAS analysis on the bacteraemia isolates identified two gene candidates that have previously been associated to this kind of infection.
Although this hypothesis was not tested here, it is possible that the emergence of a ST93-IV clade containing additional virulence genes might be related to the high prevalence of ST93-IV infections amongst the indigenous population living in the northern regions of Australia. More importantly, our data also demonstrated that GWAS can reveal candidate genes for further investigations on the pathogenesis and evolution of MRSA strains within a same lineage.
Over the last three decades, community-associated methicillin-resistant Staphylococcus aureus (CA-MRSA) has emerged globally. Although polyclonal, a small number of CA-MRSA clones are dominant in different regions of the world such as multilocus sequence type (ST) 8-IV (USA300) in North America, ST80-IV in Europe and Northern Africa, ST59-IV/V in Asia, ST772-V and ST22-IV in the Indian subcontinent, and ST30-IV in the West Pacific region . Transmission of the dominant clones in other regions has occurred, and characteristically they harbour the lukS/F-PV genes that encode the Panton-Valentine leukocidin (PVL) toxin .
In Australia, the dominant CA-MRSA clone is PVL-positive ST93-IV. Colloquially known as the “Queensland CA-MRSA clone”, ST93-IV was first described in the early 2000 s. Although known to cause severe infections including necrotizing pneumonia, ST93-IV is typically associated with skin and soft tissue infections . Reported across Australia, the clone is frequently isolated in the indigenous Australian population where its dominance is believed to be linked to overcrowding , poor hygiene and healthcare . Using whole genome sequencing (WGS) and temporal and geographical analysis, ST93 has been shown to be an early diverging and recombinant lineage genetically related to ST59/ST121 and to an unknown S. aureus lineage that emerged in the 1970 s in the North Western region of Australia . Although earlier studies into the genetic diversity of ST93 showed multiple rearrangements of the spa sequence, the core regions of the genome were very stable . However in 2014, Stinear et al. suggested ST93 clone was under pressure for adaptive change due to a reduction in both exotoxin expression and oxacillin minimum inhibitory concentration .
To screen for potential association between gene content and disease, genome-wide association studies (GWAS) can be performed by analysing single nucleotide polymorphisms (SNPs), and the accessory genes provided by WGS data. For example, GWAS performed on isolates from children with acute S. aureus osteomyelitis selected a number of virulence gene candidates potentially associated to the severity of disease . In contrast, when applied to S. aureus bacteraemia isolates, no obvious associations in the number of virulence genes present in isolates from patients with and without S. aureus infective endocarditis were identified . GWAS can also be used to examine the evolution of a bacterial clone. For example recent GWAS performed on livestock-associated CC398 MRSA, showed the clone frequently lost antimicrobial resistance genes and acquired human specific virulence genes in relation to the origin of the host .
In this study, we performed GWAS on a collection of Australian ST93 MRSA bacteraemia isolates collected over a three-year period (2015–2017) and a collection of previously published ST93 MRSA genomes (2002–2012). Phylogenetic analysis of the genomes was performed by examining SNPs in the core genome and investigating the absence/presence of accessory genes. To screen for potential genetic traits that may have assisted the ongoing transmission of ST93-IV in Australia we correlated the absence and presence of accessory genes in the ST93-IV genomes to time, location and whether they originated from a bloodstream infection.
The 423 ST93-IV were isolated across Australia from the following states and mainland territories: Northern Territory (n = 141), Queensland (n = 98), New South Wales (n = 64), Western Australia (n = 54), Victoria (n = 43), South Australia (n = 19), Australia Capital Territory (n = 3) and Tasmania (n = 1). Overall, there were 302 bacteraemia and 121 non-bacteraemia isolates. The non-bacteraemia isolates were limited to four geographical regions: New South Wales, Victoria, Western Australia and Northern Territory.
Based on core genome SNPs, the rooted phylogeny based on 1383 SNPs depicted the ST93 population to cluster primarily in two main clades (Fig. 1). Clade 1 contained 111 bacteraemia isolates predominantly from northern Australia whilst clade 2 contained 185 bacteraemia and 119 non-bacteraemia isolates collected across Australia.
Comparison between Principal Component Analysis (PCA) and Phylogenetic Clustering
By examining the presence and absence of accessory genes, PCA identified two distinct clusters (Fig. 2). Isolates in the two PCA clusters correlated with isolates in the two SNP derived phylogenetic clades.
GWAS Comparison between Bacteraemia and Non-bacteraemia ST93 Isolates
GWAS revealed nine accessory genes correlated with the bacteremia isolates (p < 0.001 and odds ratio > 1) (Table 1). However, seven of these genes were clade 1 specific and were not considered bacteraemia factors (Supplementary Table 2).
Because the majority of clade 1 genomes were bacteremia isolates, GWAS was repeated without clade 1 genomes to remove a possible selection bias. The results for both GWAS showed that the two genes that correlated with bacteraemia were hsdM (type I restriction enzyme EcoKI M protein) and clfA (clumping factor A) (Supplementary Table 3). Overall, of the 302 bacteraemia isolates, 76 % (n = 230) carried both genes; 16 % (n = 49) carried one of the genes, and the remaining 7 % (n = 23) carried neither gene. Only 43 and 45 % of the non-bacteraemia isolates carried the clfA and hsdM genes respectively.
The seven clade 1 specific accessory genes were ohrR (organic hydroperoxide resistance transcriptional regulator), acul (putative acrylyl-Coa reductase), ypuA (hypothetical protein), hutl_2 (hypothetical protein), entE (enterotoxin E), soj (chromosome-partitioning ATPase) and entA_2 (enterotoxin A) (Fig. 1). Approximately 88 % (n = 98/111) of the clade 1 genomes harboured all seven genes, with seven isolates containing none of the seven genes. The seven genes were located on five different contigs, with entE and acuI co-located with soj and ohR respectively.
Genomic diversity of ST93 over Time and Location
No significant differences in the presence or absence of accessory genes over time or location were identified.
Recombination/rearrangement of the ST93 genome
When we analysed conserved gene neighbourhoods, we observed two genes affected by re-arrangements correlating to bacteraemia, sdrF (serine-aspartate repeat-containing protein F) and pls (surface protein) (Supplementary Table 4). Analysis of the genes show that inversions occurred in regions containing sdrF and pls (Supplementary Figure 1).
In the current study we have identified two distinct ST93-IV clades circulating concurrently in Australia. The identification of the two clades by SNP analysis of the core region was supported by the PCA based on the absence and presence of genes matrix. The clade 1 isolates were primarily isolated in the northern regions of Australia spread over three states/territories (Western Australia, Northern Territory and Queensland) whilst the clade 2 isolates were distributed across the country. Based on genomic data of the van Hal et al.  historic ST93-IV isolates that were located at the root of the phylogenetic tree, we believe the two clades recently diverged from a common ancestor.
Clade 1 isolates differed from the clade 2 isolates by having acquired up to seven additional accessory genes. The known biological significance of these accessory genes varies. The entA and entE genes, encode the superantigen enterotoxins A and E respectively and play an important role in serious staphylococcal infections by triggering an overexpression of inflammatory mediators . The ohrR gene, which has previously been identified in Pseudomonas aeruginosa  and Bacillus subtilis , is known to increase an organism’s resistance to oxidative stress. The ability to resist peroxide provides the organism a growth advantage and increases its survival in host cells . The soj gene, a parA homologue involved in chromosome segregation during DNA replication, is not normally found in S. aureus . Typically, chromosome segregation in S. aureus is performed by the parB homologue spo0J, which was identified in all ST93-IV genomes. In Bacillus subtilis, soj and spo0J are present and work together to prevent premature midcell Z ring assembly . By having acquired soj, clade 1 isolates might have an advantage over non-clade 1 isolates as represented by a more efficient DNA replication system. The roles of the three remaining accessory genes, acul (a putative protein), and ypuA and hutl (both hypothetical proteins) are not known. The acquisition of the seven accessory genes, which are likely to have originated on mobile genetic elements, may explain the high rates of ST93-IV skin infections amongst indigenous children living in the northern regions of Australia . Further studies are required to determine if clade 1 has become the predominant ST93-IV strain in the region’s indigenous communities and the role of these additional genes in the expansion and fitness of this pathogen.
Based on the variability of the ST93-IV accessory genes over time and location we attempted to identify clade 1 or 2 specific subclades. Despite minor accessory gene variations occurred in a small number of isolates (for example, four isolates contained qacA [antiseptic resistance protein], qacR [HTH-type transcriptional regulator] and tnsB [transposon] which were all located on the same contig), no important difference in the absence or presence of accessory genes related to specific subclades were observed.
GWAS for Bacteraemia vs. Non-bacteraemia MRSA
In 2017 a GWAS performed by Lilje and colleagues was not able to identify genetic differences between S. aureus bacteraemia and non-bacteraemia genomes . Their results however may have been influenced by studying a variety of S. aureus lineages and clonal complexes. To identify if specific genetic factors are harboured by S. aureus bacteraemia genomes our study was limited to a single S. aureus lineage. After accounting for a possible clade 1 selection bias GWAS identified two genes associated with the ST93-IV bacteraemia isolates. The hsdM gene has recently been shown to be a hotspot for chromosome rearrangements in staphylococcus which cause phenotype switching associated with persistent infections . The clfA gene, which mediates staphylococcal binding to fibrin-coated surfaces has previously shown to be highly expressed during rat models in infective endocarditis , while clfA mutants developed milder systemic inflammation in mice models .
Chromosome rearrangements of genes may lead to altered gene expression . The 23 bacteraemic genomes that did not harbor hsdM and clfA all carried rearrangements of the pls and sdrF, genes. The pls and sdrF genes encode surface proteins. Pls, which mediates bacterial aggregation and binding to glycolipids and human epithelial cells [22, 23], has been shown in mice models to be an important factor in causing sepsis . SdrF, which is a microbial surface components recognising adhesive matrix molecule (MSCRAMM), allows staphylococcus to attach to and colonise host cells . Among the sdr gene detected amongst the different S. aureus clones, the sdr gene in the ST93 strain JKD6159 is the most diverse suggesting sdr acquisition by horizontal gene transfer. In the Huping et al. study the sdr in the ST93 genome was classified as sdrC . However, an updated annotation database has identified the gene as sdrF which had previously only been reported in S. epidermidis. SdrF adheres to human keratinocytes and epithelial cells facilitating S. epidermidis colonisation of the skin .
GWAS is a powerful tool to screen for potential associations using large datasets. However, other factors related to bacterium-host evolution may also pressure for genetic diversification. For example, patient’s age and prior medical condition, which are factors associated with MRSA bacteraemia. In the current study we selected accessory genes and gene rearrangements that show significant statistical associations with ST93-IV bacteraemia. However, to validate the scientific impact of these findings future ex vivo and in vivo investigations using gene-knockouts and expression clones are required. Furthermore, to determine if these genes are bacteraemia determinants other S. aureus lineages should be examined. Finally, phylogenetic analysis has shown ST93-IV has recently gained accessory virulence genes which might be contributing to the clone’s persistence in the Australian indigenous communities.
Bacterial Strains and Genome Assembly
A total of 300 ST93 MRSA bacteraemia isolates were identified in the 2015 , 2016  and 2017  Australian Group on Antimicrobial Resistance (AGAR) Australian Staphylococcus aureus Sepsis Outcome Programs (ASSOPs). All isolates collected were from patients with systemic infections. As part of ASSOP, MRSA isolates were referred to a central reference laboratory where genomic libraries were prepared using the Illumina Nextera® XT DNA Library Prep Kit (Illumina, United States) according to the manufacturer’s protocol. WGS was performed on the Miseq or Nextseq platforms using the Miseq Reagent Kit V3 (600 cycle) and the Nextseq 500/550 Mid Output Kit V2.5 (300 cycles), respectively. The raw sequence reads were assembled de novo using SPAdes V3.12 . Sequencing quality control was determined based on average sequencing depth. Thirty-one had genomes less with than 40x coverage and therefore were excluded. The MLST profiles of the remaining 269 genomes were determined using the mlst tool described by Seeman et al. .
All sequence data obtained from this study were deposited to the NCBI Sequence Read Archive under BioProject ID PRJNA644215.
Using the chromosome of S. aureus CC398 reference strain SO395 (GenBank accession ID AM990992) as the reference genome, the bacterial variant calling tool snippy V4.1.0  was used to extract and align SNPs from the core genome. The 423 ST93 genomes were used to generate a rooted maximum parsimony phylogenetic tree using MEGA V10.1.7  with the following parameters; bootstrap value: 1000, nucleotide substitution model and the SPR model for the MP search method. Phylogenetic clades were defined as a cluster of isolates sharing multiple common SNP mutations. The iTOL V3 web service was used to visualise the phylogenetic tree and the corresponding metadata .
Genome-Wide Association Study (GWAS)
Genes from the 423 assembled S. aureus genome sequences were annotated with Prokka V1.13  using default parameters and the pan-genome was extracted by Roary V3.12.0  using the -s option of no paralog splitting. The pan-genome matrix from Roary containing of gene presence or absence for each genome was used as input for Scoary V1.6.16  with the following traits; SNP phylogeny clades, location (states and territories), year of isolation, clade specific genes and whether the isolate was from a bloodstream infection. Adapting the method described by Arnoud H. M. van Vliet , genes returning a Bonferroni corrected p value ≤ 10− 5 and odds ratio > 1 were further investigated. In addition to Scoary analysis, principal component analysis (PCA) of binomial variables on the pan-genome matrix was performed for determination of association with the statistical package R version 3.5.1  and ggplot2 V3.2.1 to confirm relationships between genes identified in Scoary and traits.
Detecting gene rearrangements
The pan-genome matrix was compared with and without the -s option. Scoary was used on both pan-genome matrices using bacteremia as phenotype. Identification of genes correlating to bacteremia were compared between both sets of data. Genes associated to bacteremia were extracted along with neighbouring genes and aligned against the ST93 genome JDK6159 using Artemis comparison tool  to visualise the rearrangement structure.
Availability of data and materials
The data that support the findings of this study are openly available on the SRA database under Bioproject: PRJNA644215, Accession: SRX8689588-SRX8689856. Collection of ST93 S. aureus used as reference in this study is available on ENA (Supplementary Table 1). An additional reference genome used is also available from GenBank (Accession: AM990992).
Australian Group on Antimicrobial Resistance
Australian Staphylococcus aureus Sepsis Outcome Programs
Community-Associated Methicillin-Resistant Staphylococcus aureus
Genome Wide Association Study
microbial surface components recognising adhesive matrix molecule
Multi Locus Sequence Typing
Principal Component Analysis
Single Nucleotide Polymorphisms
Whole Genome Sequencing
Lakhundi S, Zhang K. Methicillin-Resistant Staphylococcus aureus: Molecular Characterization, Evolution, and Epidemiology. Clin Microbiol Rev. 2018;31(4).
Coombs GW, Goering RV, Chua KY, Monecke S, Howden BP, Stinear TP, et al. The molecular epidemiology of the highly virulent ST93 Australian community Staphylococcus aureus strain. PLoS One. 2012;7(8):e43037.
Coombs GW, Daley DA, Mowlaboccus S, Lee YT, Pang S. Australian Group on Antimicrobial R. Australian Group on Antimicrobial Resistance (AGAR) Australian Staphylococcus aureus Sepsis Outcome Programme (ASSOP) Annual Report 2018. Commun Dis Intell (2018). 2020;44.
Peleg AY, Munckhof WJ. Fatal necrotising pneumonia due to community-acquired methicillin-resistant Staphylococcus aureus (MRSA). Med J Aust. 2004;181(4):228–9.
van Hal SJ, Steinig EJ, Andersson P, Holden MTG, Harris SR, Nimmo GR, et al. Global Scale Dissemination of ST93: A Divergent Staphylococcus aureus Epidemic Lineage That Has Recently Emerged From Remote Northern Australia. Front Microbiol. 2018;9:1453.
Uhlemann AC, Otto M, Lowy FD, DeLeo FR. Evolution of community- and healthcare-associated methicillin-resistant Staphylococcus aureus. Infect Genet Evol. 2014;21:563–74.
Stinear TP, Holt KE, Chua K, Stepnell J, Tuck KL, Coombs G, et al. Adaptive change inferred from genomic population analysis of the ST93 epidemic clone of community-associated methicillin-resistant Staphylococcus aureus. Genome Biol Evol. 2014;6(2):366–78.
Collins A, Wakeland EK, Raj P, Kim MS, Kim J, Tareen NG, et al. The impact of Staphylococcus aureus genomic variation on clinical phenotype of children with acute hematogenous osteomyelitis. Heliyon. 2018;4(6):e00674.
Lilje B, Rasmussen RV, Dahl A, Stegger M, Skov RL, Fowler VG, et al. Whole-genome sequencing of bloodstream Staphylococcus aureus isolates does not distinguish bacteraemia from endocarditis. Microb Genom. 2017;3(11).
Sieber RN, Larsen AR, Urth TR, Iversen S, Moller CH, Skov RL, et al. Genome investigations show host adaptation and transmission of LA-MRSA CC398 from pigs into Danish healthcare institutions. Sci Rep. 2019;9(1):18655.
Liu Q, Mazhar M, Miller LS. Immune and Inflammatory Reponses to Staphylococcus aureus Skin Infections. Curr Dermatol Rep. 2018;7(4):338–49.
Atichartpongkul S, Fuangthong M, Vattanaviboon P, Mongkolsuk S. Analyses of the regulatory mechanism and physiological roles of Pseudomonas aeruginosa OhrR, a transcription regulator and a sensor of organic hydroperoxides. J Bacteriol. 2010;192(8):2093–101.
Fuangthong M, Atichartpongkul S, Mongkolsuk S, Helmann JD. OhrR is a repressor of ohrA, a key organic hydroperoxide resistance determinant in Bacillus subtilis. J Bacteriol. 2001;183(14):4134–41.
Li X, Tao J, Han J, Hu X, Chen Y, Deng H, et al. The gain of hydrogen peroxide resistance benefits growth fitness in mycobacteria under stress. Protein Cell. 2014;5(3):182–5.
Murray H, Errington J. Dynamic control of the DNA replication initiation protein DnaA by Soj/ParA. Cell. 2008;135(1):74–84.
Hajduk IV, Mann R, Rodrigues CDA, Harry EJ. The ParB homologs, Spo0J and Noc, together prevent premature midcell Z ring assembly when the early stages of replication are blocked in Bacillus subtilis. Mol Microbiol. 2019;112(3):766–84.
Bowen AC, Mahe A, Hay RJ, Andrews RM, Steer AC, Tong SY, et al. The Global Epidemiology of Impetigo: A Systematic Review of the Population Prevalence of Impetigo and Pyoderma. PLoS One. 2015;10(8):e0136789.
Guerillot R, Kostoulias X, Donovan L, Li L, Carter GP, Hachani A, et al. Unstable chromosome rearrangements in Staphylococcus aureus cause phenotype switching associated with persistent infections. Proc Natl Acad Sci U S A. 2019;116(40):20135–40.
Hanses F, Roux C, Dunman PM, Salzberger B, Lee JC. Staphylococcus aureus gene expression in a rat model of infective endocarditis. Genome Med. 2014;6(10):93.
Josefsson E, Higgins J, Foster TJ, Tarkowski A. Fibrinogen binding sites P336 and Y338 of clumping factor A are crucial for Staphylococcus aureus virulence. PLoS One. 2008;3(5):e2206.
Raeside C, Gaffe J, Deatherage DE, Tenaillon O, Briska AM, Ptashkin RN, et al. Large chromosomal rearrangements during a long-term evolution experiment with Escherichia coli. mBio. 2014;5(5):e01377-14.
Huesca M, Peralta R, Sauder DN, Simor AE, McGavin MJ. Adhesion and virulence properties of epidemic Canadian methicillin-resistant Staphylococcus aureus strain 1: identification of novel adhesion functions associated with plasmin-sensitive surface protein. J Infect Dis. 2002;185(9):1285–96.
Roche FM, Meehan M, Foster TJ. The Staphylococcus aureus surface protein SasG and its homologues promote bacterial adherence to human desquamated nasal epithelial cells. Microbiology. 2003;149(Pt 10):2759–67.
Josefsson E, Juuti K, Bokarewa M, Kuusela P. The surface protein Pls of methicillin-resistant Staphylococcus aureus is a virulence factor in septic arthritis. Infect Immun. 2005;73(5):2812–7.
Foster TJ, Geoghegan JA, Ganesh VK, Hook M. Adhesion, invasion and evasion: the many functions of the surface proteins of Staphylococcus aureus. Nat Rev Microbiol. 2014;12(1):49–62.
Xue H, Lu H, Zhao X. Sequence diversities of serine-aspartate repeat genes among Staphylococcus aureus isolates from different hosts presumably by horizontal gene transfer. PLoS One. 2011;6(5):e20332.
Trivedi S, Uhlemann AC, Herman-Bausier P, Sullivan SB, Sowash MG, Flores EY, et al. The Surface Protein SdrF Mediates Staphylococcus epidermidis Adherence to Keratin. J Infect Dis. 2017;215(12):1846–54.
Coombs GW, Daley DA, Lee YT, Pang S, Bell JM, Turnidge JD, et al Australian Group on Antimicrobial Resistance (AGAR) Australian Staphylococcus aureus Sepsis Outcome Programme (ASSOP) Annual Report 2015. Commun Dis Intell (2018). 2018;42.
Coombs GW, Daley DA, Lee YT, Pang S. Australian Group on Antimicrobial R. Australian Group on Antimicrobial Resistance (AGAR) Australian Staphylococcus aureus Sepsis Outcome Programme (ASSOP) Annual Report 2016. Commun Dis Intell (2018). 2018;42.
Coombs GW, Daley DA, Lee YT, Pang S. Australian Group on Antimicrobial Resistance (AGAR) Australian Staphylococcus aureus Sepsis Outcome Programme (ASSOP) Annual Report 2017. Commun Dis Intell (2018). 2019;43.
Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, et al. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol. 2012;19(5):455–77.
Jolley KA, Maiden MC. BIGSdb. Scalable analysis of bacterial genome variation at the population level. BMC Bioinformatics. 2010;11:595.
Seemann T. Snippy: fast bacterial variant calling from NGS reads. 2015.
Kumar S, Stecher G, Li M, Knyaz C, Tamura K. MEGA X: Molecular Evolutionary Genetics Analysis across Computing Platforms. Mol Biol Evol. 2018;35(6):1547–9.
Letunic I, Bork P. Interactive tree of life (iTOL) v3: an online tool for the display and annotation of phylogenetic and other trees. Nucleic Acids Res. 2016;44(W1):W242-5.
Seemann T. Prokka: rapid prokaryotic genome annotation. Bioinformatics. 2014;30(14):2068–9.
Page AJ, Cummins CA, Hunt M, Wong VK, Reuter S, Holden MT, et al. Roary: rapid large-scale prokaryote pan genome analysis. Bioinformatics. 2015;31(22):3691–3.
Brynildsrud O, Bohlin J, Scheffer L, Eldholm V. Rapid scoring of genes in microbial pan-genome-wide association studies with Scoary. Genome Biol. 2016;17(1):238.
van Vliet AH. Use of pan-genome analysis for the identification of lineage-specific genes of Helicobacter pylori. FEMS Microbiol Lett. 2017;364(2).
R_core_Team R: A Language and Environment for Statistical Computing. 2019.
Carver TJ, Rutherford KM, Berriman M, Rajandream MA, Barrell BG, Parkhill J. ACT: the Artemis Comparison Tool. Bioinformatics. 2005;21(16):3422–3.
This work received no specific grant from any funding agency.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Table 1: ST93-IV isolates used in this study with bacteremia gene matrix
Supplementary Table 2: Statistics supporting clade specific genes
Supplementary Table 3: Statistics supporting bacteremia in clade 2
Supplementary Table 4: Gene rearrangements statistics associated to bacteremia
Supplementary Figure 1: Rearrangements of A) pls and B) sdrF gene neighbourhoods.
About this article
Cite this article
Pang, S., Daley, D.A., Sahibzada, S. et al. Genome-wide association studies reveal candidate genes associated to bacteraemia caused by ST93-IV CA-MRSA. BMC Genomics 22, 418 (2021). https://doi.org/10.1186/s12864-021-07738-4
- Staphylococcus aureus