Genomic Epidemiology and Global Population Structure of Exfoliative Toxin A-Producing Staphylococcus aureus Strains Associated With Staphylococcal Scalded Skin Syndrome

ST913-IVa-t991 Methicillin-Resistant Staphylococcus aureus among Pediatric Patients, Israel

In Israel, prevalence of sequence type 913, staphylococcal cassette chromosome mecIVa, spa type t991 methicillin-resistant Staphylococcus aureus lineage has surged among pediatric populations, predominantly in Arab and Orthodox Jewish communities. Antimicrobial resistance patterns vary by demographics. This lineage's spread and microevolution in the Middle East underscore the need for ongoing surveillance.

Average nucleotide identity-based Staphylococcus aureus strain grouping allows identification of strain-specific genes in the pangenome

Staphylococcus aureus causes both hospital- and community-acquired infections in humans worldwide. Due to the high incidence of infection, S. aureus is also one of the most sampled and sequenced pathogens today, providing an outstanding resource to understand variation at the bacterial subspecies level. We processed and downsampled 83,383 public S. aureus Illumina whole-genome shotgun sequences and 1,263 complete genomes to produce 7,954 representative substrains. Pairwise comparison of average nucleotide identity revealed a natural boundary of 99.5% that could be used to define 145 distinct strains within the species. We found that intermediate frequency genes in the pangenome (present in 10%-95% of genomes) could be divided into those closely linked to strain background ("strain-concentrated") and those highly variable within strains ("strain-diffuse"). Non-core genes had different patterns of chromosome location. Notably, strain-diffuse genes were associated with prophages; strain-concentrated genes were associated with the vSaβ genome island and rare genes (<10% frequency) concentrated near the origin of replication. Antibiotic resistance genes were enriched in the strain-diffuse class, while virulence genes were distributed between strain-diffuse, strain-concentrated, core, and rare classes. This study shows how different patterns of gene movement help create strains as distinct subspecies entities and provide insight into the diverse histories of important S. aureus functions.

IMPORTANCE

We analyzed the genomic diversity of Staphylococcus aureus, a globally prevalent bacterial species that causes serious infections in humans. Our goal was to build a genetic picture of the different strains of S. aureus and which genes may be associated with them. We reprocessed >84,000 genomes and subsampled to remove redundancy. We found that individual samples sharing >99.5% of their genome could be grouped into strains. We also showed that a portion of genes that are present in intermediate frequency in the species are strongly associated with some strains but completely absent from others, suggesting a role in strain specificity. This work lays the foundation for understanding individual gene histories of the S. aureus species and also outlines strategies for processing large bacterial genomic data sets.

A look at staphylococci from the one health perspective

Staphylococcus aureus and other staphylococcal species are resident and transient multihost colonizers as well as conditional pathogens. Especially S. aureus represents an excellent model bacterium for the "One Health" concept because of its dynamics at the human-animal interface and versatility with respect to host adaptation. The development of antimicrobial resistance plays another integral part. This overview will focus on studies at the human-animal interface with respect to livestock farming and to companion animals, as well as on staphylococci in wildlife. In this context transmissions of staphylococci and of antimicrobial resistance genes between animals and humans are of particular significance.

Average Nucleotide Identity based Staphylococcus aureus strain grouping allows identification of strain-specific genes in the pangenome

Staphylococcus aureus causes both hospital and community acquired infections in humans worldwide. Due to the high incidence of infection S. aureus is also one of the most sampled and sequenced pathogens today, providing an outstanding resource to understand variation at the bacterial subspecies level. We processed and downsampled 83,383 public S. aureus Illumina whole genome shotgun sequences and 1,263 complete genomes to produce 7,954 representative substrains. Pairwise comparison of core gene Average Nucleotide Identity (ANI) revealed a natural boundary of 99.5% that could be used to define 145 distinct strains within the species. We found that intermediate frequency genes in the pangenome (present in 10-95% of genomes) could be divided into those closely linked to strain background ("strain-concentrated") and those highly variable within strains ("strain-diffuse"). Non-core genes had different patterns of chromosome location; notably, strain-diffuse associated with prophages, strain-concentrated with the vSaβ genome island and rare genes (<10% frequency) concentrated near the origin of replication. Antibiotic genes were enriched in the strain-diffuse class, while virulence genes were distributed between strain-diffuse, strain-concentrated, core and rare classes. This study shows how different patterns of gene movement help create strains as distinct subspecies entities and provide insight into the diverse histories of important S. aureus functions.

Multiclonal human origin and global expansion of an endemic bacterial pathogen of livestock

Most new pathogens of humans and animals arise via switching events from distinct host species. However, our understanding of the evolutionary and ecological drivers of successful host adaptation, expansion, and dissemination are limited. Staphylococcus aureus is a major bacterial pathogen of humans and a leading cause of mastitis in dairy cows worldwide. Here we trace the evolutionary history of bovine S. aureus using a global dataset of 10,254 S. aureus genomes including 1,896 bovine isolates from 32 countries in 6 continents. We identified 7 major contemporary endemic clones of S. aureus causing bovine mastitis around the world and traced them back to 4 independent host-jump events from humans that occurred up to 2,500 y ago. Individual clones emerged and underwent clonal expansion from the mid-19th to late 20th century coinciding with the commercialization and industrialization of dairy farming, and older lineages have become globally distributed via established cattle trade links. Importantly, we identified lineage-dependent differences in the frequency of host transmission events between humans and cows in both directions revealing high risk clones threatening veterinary and human health. Finally, pangenome network analysis revealed that some bovine S. aureus lineages contained distinct sets of bovine-associated genes, consistent with multiple trajectories to host adaptation via gene acquisition. Taken together, we have dissected the evolutionary history of a major endemic pathogen of livestock providing a comprehensive temporal, geographic, and gene-level perspective of its remarkable success.

Case Report: Necrotizing fasciitis caused by Staphylococcus aureus positive for a new sequence variant of exfoliative toxin E

Objectives: Necrotizing fasciitis (NF) caused by S. aureus is a rare, aggressive and rapidly progressing superficial fascia infection with a high mortality rate. The aim of this study was to identify virulence-related genes from a complete genome sequence of a methicillin-susceptible S. aureus (MSSA) isolate recovered from a monomicrobial case of NF.

Materials and methods: The MSSA isolate UMCG579 was cultured from a pus collection from the subcutis of a patient with NF. The genome of isolate UMCG579 was sequenced using MinION (Oxford Nanopore) and MiSeq (illumina) platforms.

Results: The genome of the UMCG579 isolate was composed of a 2,741,379 bp chromosome and did not harbor any plasmids. Virulence factor profiling identified multiple pore-forming toxin genes in the UMCG579 chromosome, including the Panton-Valentine leukocidin (PVL) genes, and none of the superantigen genes. The UMCG579 isolate harbored a new sequence variant of the recently described ete gene encoding exfoliative toxin (type E). A search in the GenBank database revealed that the new sequence variant (ete2) was exclusively found among isolates (n = 115) belonging to MLST CC152. While the majority of S. aureus ete-positive isolates were recovered from animal sources, S. aureus ete2-positive isolates originated from human carriers and human infections. Comparative genome analysis revealed that the ete2 gene was located on a 8777 bp genomic island.

Conclusion: The combination of two heterogeneously distributed potent toxins, ETE2 and PVL, is likely to enhance the pathogenic ability of S. aureus isolates. Since anti-virulence therapies for the treatment of S. aureus infections continue to be explored, the understanding of specific pathogenetic mechanisms may have an important prophylactic and therapeutic value. Nevertheless, the exact contribution of ETE sequence variants to S. aureus virulence in NF infections must be determined.

Genome Evolution of Invasive Methicillin-Resistant Staphylococcus aureus in the Americas

Staphylococcus aureus causes a variety of debilitating and life-threatening diseases, and thus remains a challenging global health threat. S. aureus is remarkably diverse, yet only a minority of methicillin-resistant S. aureus (MRSA) clones have caused pandemic proportions of diseases. The genetic drivers of the successful dissemination of some clones across wide geographical expanses remain poorly understood. We analyzed 386 recently published MRSA genomes from bloodstream infections sampled in North, Central, and South America from 2011 to 2018. Here, we show that MRSA-associated bloodstream infections were attributable to two genetically distinct lineages. One lineage consisted almost exclusively of sequence type (ST) 8, which emerged in 1964. A second lineage emerged in 1986 and consisted of STs 5, 105, and 231. The two lineages have simultaneously disseminated across geographically distant sites. Sublineages rapidly diverged within locations in the early 2000s. Their diversification was associated with independent acquisitions of unique variants of the mobile mecA-carrying chromosomal cassette and distinct repertoires of antimicrobial resistance genes. We show that the evolution and spread of invasive multidrug-resistant MRSA in the Americas was driven by transcontinental dissemination, followed by more recent establishment and divergence of local pathogen populations. Our study highlights the need for continued international surveillance of high-risk clones to control the global health threat of multidrug resistance.

IMPORTANCE Bloodstream infections due to S. aureus cause significant patient morbidity and mortality worldwide, exacerbated by the emergence and spread of methicillin resistant S. aureus (MRSA). This study provides important insights on the evolution and long-distance geographic expansion of two distinct MRSA lineages that predominate in bloodstream infections in the past 5 decades. The success of these two lineages partly lies on their acquisition of a diverse set of antimicrobial resistance genes and of unique variants of the mobile genetic element SCCmec that carries the gene mecA conferring resistance to beta-lactams. High-risk antimicrobial resistant clones can therefore rapidly disseminate across long distances and establish within local communities within a short period of time. These results have important implications for global initiatives and local epidemiological efforts to monitor and control invasive MRSA infections and transcontinental spread of multidrug resistance.

Meaningful Use of Pathogen Genomic Data

Population genomic analysis is a powerful tool to understand the evolutionary history of pathogens and the factors contributing to the success or failure of lineages. These studies have significant implications for human health, as evident from our ongoing tracking of SARS-CoV-2. In their article, Gill et al. (J. L. Gill, J. Hedge, D. J. Wilson, and R. C. MacLean, mBio 12:e02168-21, 2021, https://doi.org/10.1128/mBio.02168-21) demonstrate the utility of pathogen genomic data by comprehensively elucidating the origin of methicillin-resistant Staphylococcus aureus ST239. To accomplish this, they leveraged newly developed tools for querying large genomic data sets. Overall, these analyses rely on the availability of representative genomic data along with their associated metadata-information about where and when samples were collected, clinical and epidemiological characteristics, and phenotypic properties. However, in many instances, these data are missing. Here, I borrow the term "meaningful use" from the Health IT field to describe the need to maximize the utility of genomic data and make suggestions for how to address the current limitations.

Predominance of community-associated, methicillin-susceptible Staphylococcus aureus infections among hospitalized children and adolescents

Introduction. Staphylococcus aureus infections cause significant morbidity and mortality in children and adolescents.Gap Statement. There is limited data on the characteristics of S. aureus infections requiring hospitalization in childhood.Aim.To investigate the molecular epidemiology and antibiotic resistance of S. aureus clinical isolates from children and adolescents.Methodology.All S. aureus isolates recovered from patients aged <18 years, admitted to a referral hospital, with culture-proven invasive or non-invasive infections during the 4 year period 2015 to 2018 were analysed for antimicrobial resistance, virulence genes, PFGE and multilocus sequence typing (MLST). Cases were assigned to community-associated, community-onset healthcare-associated or hospital-associated infections based on epidemiological case definitions.Results.Among 139 S. aureus infections, 88.5 % (123/139) were caused by methicillin-susceptible isolates (MSSA) and 73.4 % (102/139) were classified as community-associated infections. tst and lukS/lukF-PV genes were more common among MRSA as compared to MSSA isolates (tst, p 0.04; lukS/lukF-PV, p 0.007). Invasive disease was noted in 22/139 patients (15.8 %). Staphylococcal scalded skin syndrome caused by fusidic-resistant MSSA increased over time (22.8 % in 2017-2018 vs 8.3 % in 2015-2016, OR 3.24; 95 % CI 1.10-8.36; P 0.03). By PFGE genotyping, 22 pulsotypes were identified. A total of five sequence types (STs) were identified among 58 isolates analysed by MLST. More than one third of MSSA isolates (40/123, 32.5 %) and 13/23 (56.5 %) of SSSS isolates belonged to pulsotype 1, classified as sequence type 121 (ST121). MRSA isolates were equally distributed to pulsotypes A (ST30), B (ST239), C (ST80), H (ST225). ST121 isolates carried fnbA (40/40), eta/etb genes (29/40), exhibited high resistance to fusidic acid and were increasingly resistant to mupirocin.Conclusion.In our population, community-associated MSSA was the predominant cause of S. aureus infections characterized by polyclonality, increasing resistance to fusidic acid and mupirocin. PFGE type 1 ST121 clone, harboured exfoliative toxin genes and was associated with rising trends of SSSS.

Three-Dimensional Structure Characterization and Inhibition Study of Exfoliative Toxin D From Staphylococcus aureus

The Staphylococcus aureus exfoliative toxins (ETs) are the main toxins that produce staphylococcal scalded skin syndrome (SSSS), an abscess skin disorder. The victims of the disease are usually newborns and kids, as well as grown-up people. Five ETs namely, exfoliative toxins A, B, C, D, and E have been identified in S. aureus. The three-dimensional (3D) structure of exfoliative toxins A, B, C and E is known, while that of exfoliative toxin D (ETD) is still unknown. In this work, we have predicted the 3D structure of ETD using protein modeling techniques (software used for 3D structure modeling comprising the MODELLER 9v19 program, SWISS-Model, and I-TESSER). The validation of the build model was done using PROCHECK (Ramachandran plot), ERRAT2, and Verify 3D programs. The results from 3D modeling show that the build model was of good quality as indicated by a GMQE score of 0.88 and by 91.1% amino acid residues in the most favored region of the Ramachandran plot, the ERRAT2 quality factor of 90.1%, and a verify3D score of >0.2 for 99.59% of amino acid residues. The 3D structure analysis indicates that the overall structure of ETD is similar to the chymotrypsin-like serine protease fold. The structure is composed of 13 β-strands and seven α-helices that fold into two well-defined six-strand β-barrels whose axes are roughly perpendicular to each other. The active site residues include histidine-97, aspartic acid-147, and serine-221. This represents the first structure report of ETD. Structural comparison with the other ETs shows some differences, particularly in the loop region, which also change the overall surface charge of these toxins. This may convey variable substrate specificity to these toxins. The inhibition of these toxins by natural (2S albumin and flocculating proteins from Moringa oleifera seeds) and synthetic inhibitors (suramin) was also carried out in this study. The results from docking indicate that the inhibitors bind near the C-terminal domain which may restrict the movement of this domain and may halt the access of the substrate to the active site of this enzyme. Molecular dynamic simulation was performed to see the effect of inhibitor binding to the enzyme. This work will further elucidate the structure-function relationship of this enzyme. The inhibition of this enzyme will lead to a new treatment for SSSS.