Mobile elements drive recombination hotspots in the core genome of Staphylococcus aureus

Within-host genomic evolution of methicillin-resistant Staphylococcus aureus in long-term carriers

Assessing the genomic evolution of Staphylococcus aureus can help us understand how the bacteria adapt to its environment. In this study, we aimed to assess the mutation rate within 144 methicillin-resistant Staphylococcus aureus (MRSA) carriers with a carriage time from 4 to 11 years, including some carriers who belonged to the same households. We found that 23 of the 144 individuals had completely different MRSA types over time and were therefore not long-term carriers of the same MRSA. From the remaining 121 individuals, we performed whole-genome sequencing (WGS) on 424 isolates and then compared these pairwise using core genome multilocus sequence typing (cgMLST) and single-nucleotide polymorphism (SNP) analyses. We found a median within-host mutation rate in long-term MRSA carriers of 4.9 (3.4-6.9) SNPs/genome/year and 2.7 (1.8-4.2) allelic differences/genome/year, when excluding presumed recombination. Furthermore, we stratified the cohort into subgroups and found no significant difference between the median mutation rate of members of households, individuals with presumed continued exposure, e.g., from travel and persons without known continued exposure. Finally, we found that SNPs occurred at random within the genes in our cohort. KEY POINTS: • Median mutation rate within long-term MRSA carriers of 4.9 (3.4-6.9) SNPs/genome/year • Similar median mutation rates in subgroups (households, travelers) • No hotspots for SNPs within the genome.

Homologous Recombination Shapes the Architecture and Evolution of Bacterial Genomes

Homologous recombination is a key evolutionary force that varies considerably across bacterial species. However, how the landscape of homologous recombination varies across genes and within individual genomes has only been studied in a few species. Here, we used Approximate Bayesian Computation to estimate the recombination rate along the genomes of 145 bacterial species. Our results show that homologous recombination varies greatly along bacterial genomes and shapes many aspects of genome architecture and evolution. The genomic landscape of recombination presents several key signatures: rates are highest near the origin of replication in most species, patterns of recombination generally appear symmetrical in both replichores (i.e. replicational halves of circular chromosomes) and most species have genomic hotpots of recombination. Furthermore, many closely related species share conserved landscapes of recombination across orthologs indicating that recombination landscapes are conserved over significant evolutionary distances. We show evidence that recombination drives the evolution of GC-content through increasing the effectiveness of selection and not through biased gene conversion, thereby contributing to an ongoing debate. Finally, we demonstrate that the rate of recombination varies across gene function and that many hotspots of recombination are associated with adaptive and mobile regions often encoding genes involved in pathogenicity.

Mechanisms of acquisition of the vanA operon among vancomycin-resistant Staphylococcus aureus genomes: The tip of the iceberg?

BACKGROUND

Vancomycin is frequently used as a last line of defence against infections due to multidrug-resistant Staphylococcus aureus (S. aureus). A recent finding described the acquisition of vancomycin-resistant S. aureus strains by the integration of an enterococcal plasmid containing the vanA operon into the S. aureus chromosome via homologous recombination involving a specific integration site called locus L2.

METHODS

To characterise all mechanisms of acquisition of vanA, this study analysed the 15 706 S. aureus genomes to look for vanA and described its genetic environment.

RESULTS

A complete vanA operon was found in 25 S. aureus strains isolated from 12 patients, including nine co-isolated with vancomycin-resistant Enterococcus strains. VanA was found within transposon Tn1546-like elements on 17 plasmids and eight chromosomes. VanA might be acquired through conjugation of enterococcal and staphylococcal plasmids, transposition of Tn1546 carrying vanA and plasmid integration into the chromosome. Further, L2 was detected in 2087 genomes (13.3%) of S. aureus strains across different continents. Six potential chromosomal hotspots for integration of the entire vanA-containing enterococcal plasmid were identified by homologous recombination via L2.

CONCLUSIONS

These findings suggest that the recently described scenario in a New York patient could be reproduced anywhere. Surveillance of this possibility is mandatory, especially in patients with vancomycin-resistant Enterococcus infection or colonisation.

Isolation and characterization of novel Staphylococcus aureus bacteriophage Hesat from dairy origin

A novel temperate phage, named Hesat, was isolated by the incubation of a dairy strain of Staphylococcus aureus belonging to spa-type t127 with either bovine or ovine milk. Hesat represents a new species of temperate phage within the Phietavirus genus of the Azeredovirinae subfamily. Its genome has a length of 43,129 bp and a GC content of 35.11% and contains 75 predicted ORFs, some of which linked to virulence. This includes (i) a pathogenicity island (SaPln2), homologous to the type II toxin-antitoxin system PemK/MazF family toxin; (ii) a DUF3113 protein (gp30) that is putatively involved in the derepression of the global repressor Stl; and (iii) a cluster coding for a PVL. Genomic analysis of the host strain indicates Hesat is a resident prophage. Interestingly, its induction was obtained by exposing the bacterium to milk, while the conventional mitomycin C-based approach failed. The host range of phage Hesat appears to be broad, as it was able to lyse 24 out of 30 tested S. aureus isolates. Furthermore, when tested at high titer (108 PFU/ml), Hesat phage was also able to lyse a Staphylococcus muscae isolate, a coagulase-negative staphylococcal strain. KEY POINTS: • A new phage species was isolated from a Staphylococcus aureus bovine strain. • Pathogenicity island and PVL genes are encoded within phage genome. • The phage is active against most of S. aureus strains from both animal and human origins.

CRISPR-COPIES: an in silico platform for discovery of neutral integration sites for CRISPR/Cas-facilitated gene integration

The CRISPR/Cas system has emerged as a powerful tool for genome editing in metabolic engineering and human gene therapy. However, locating the optimal site on the chromosome to integrate heterologous genes using the CRISPR/Cas system remains an open question. Selecting a suitable site for gene integration involves considering multiple complex criteria, including factors related to CRISPR/Cas-mediated integration, genetic stability, and gene expression. Consequently, identifying such sites on specific or different chromosomal locations typically requires extensive characterization efforts. To address these challenges, we have developed CRISPR-COPIES, a COmputational Pipeline for the Identification of CRISPR/Cas-facilitated intEgration Sites. This tool leverages ScaNN, a state-of-the-art model on the embedding-based nearest neighbor search for fast and accurate off-target search, and can identify genome-wide intergenic sites for most bacterial and fungal genomes within minutes. As a proof of concept, we utilized CRISPR-COPIES to characterize neutral integration sites in three diverse species: Saccharomyces cerevisiae, Cupriavidus necator, and HEK293T cells. In addition, we developed a user-friendly web interface for CRISPR-COPIES (https://biofoundry.web.illinois.edu/copies/). We anticipate that CRISPR-COPIES will serve as a valuable tool for targeted DNA integration and aid in the characterization of synthetic biology toolkits, enable rapid strain construction to produce valuable biochemicals, and support human gene and cell therapy applications.

Impact of Homologous Recombination on Core Genome Evolution and Host Adaptation of Pectobacterium parmentieri

Homologous recombination is a major force mechanism driving bacterial evolution, host adaptability, and acquisition of novel virulence traits. Pectobacterium parmentieri is a plant bacterial pathogen distributed worldwide, primarily affecting potatoes, by causing soft rot and blackleg diseases. The goal of this investigation was to understand the impact of homologous recombination on the genomic evolution of P. parmentieri. Analysis of P. parmentieri genomes using Roary revealed a dynamic pan-genome with 3,742 core genes and over 55% accessory genome variability. Bayesian population structure analysis identified 7 lineages, indicating species heterogeneity. ClonalFrameML analysis displayed 5,125 recombination events, with the lineage 4 exhibiting the highest events. fastGEAR analysis identified 486 ancestral and 941 recent recombination events ranging from 43 bp to 119 kb and 36 bp to 13.96 kb, respectively, suggesting ongoing adaptation. Notably, 11% (412 genes) of the core genome underwent recent recombination, with lineage 1 as the main donor. The prevalence of recent recombination (double compared to ancient) events implies continuous adaptation, possibly driven by global potato trade. Recombination events were found in genes involved in vital cellular processes (DNA replication, DNA repair, RNA processing, homeostasis, and metabolism), pathogenicity determinants (type secretion systems, cell-wall degrading enzymes, iron scavengers, lipopolysaccharides (LPS), flagellum, etc.), antimicrobial compounds (phenazine and colicin) and even CRISPR-Cas genes. Overall, these results emphasize the potential role of homologous recombination in P. parmentieri's evolutionary dynamics, influencing host colonization, pathogenicity, adaptive immunity, and ecological fitness.

Genomic analysis of Ralstonia pickettii reveals the genetic features for potential pathogenicity and adaptive evolution in drinking water

Ralstonia pickettii, the most critical clinical pathogen of the genus Ralstonia, has been identified as a causative agent of numerous harmful infections. Additionally, Ralstonia pickettii demonstrates adaptability to extreme environmental conditions, such as those found in drinking water. In this study, we conducted a comprehensive genomic analysis to investigate the genomic characteristics related to potential pathogenicity and adaptive evolution in drinking water environments of Ralstonia pickettii. Through phylogenetic analysis and population genetic analysis, we divided Ralstonia pickettii into five Groups, two of which were associated with drinking water environments. The open pan-genome with a large and flexible gene repertoire indicated a high genetic plasticity. Significant differences in functional enrichment were observed between the core- and pan-genome of different groups. Diverse mobile genetic elements (MGEs), extensive genomic rearrangements, and horizontal gene transfer (HGT) events played a crucial role in generating genetic diversity. In drinking water environments, Ralstonia pickettii exhibited strong adaptability, and the acquisition of specific adaptive genes was potentially facilitated by genomic islands (GIs) and HGT. Furthermore, environmental pressures drove the adaptive evolution of Ralstonia pickettii, leading to the accumulation of unique mutations in key genes. These mutations may have a significant impact on various physiological functions, particularly carbon metabolism and energy metabolism. The presence of virulence-related elements associated with macromolecular secretion systems, virulence factors, and antimicrobial resistance indicated the potential pathogenicity of Ralstonia pickettii, making it capable of causing multiple nosocomial infections. This study provides comprehensive insights into the potential pathogenicity and adaptive evolution of Ralstonia pickettii in drinking water environments from a genomic perspective.

Epidemiological characterization of clinical isolates of meticillin resistant Staphylococcus aureus through multilocus sequence typing and staphylococcal cassette chromosome mec typing in Northwest Iran

BACKGROUND

Methicillin-resistant Staphylococcus aureus (MRSA), is considered a potential and aggressive nosocomial pathogen. It accounts for 50% of S. aureus isolates in tertiary hospitals in Iran, however, there is no sufficient evolutionary and epidemiological investigation about this medically important bacterium. We aimed to study the lineage and evolution of MRSA in Northwest Iran during 2021-2022 based on the obtained phenotypic and genotypic characteristics.

MATERIALS AND METHODS

Seventy-two non-duplicate MRSA isolates were collected from 3 referral hospitals in Tabriz, Ardebil, and Urmia cities. The antimicrobial susceptibility patterns were determined by disk diffusion test and micro broth dilution methods. Thereafter 4 virulence genes (eta, etb, pvl, tst) and 5 types of staphylococcal cassette chromosome mec (SCCmec) were detected by PCR. In the final step, representative isolates were selected to be studied by Multilocus sequence typing (MLST).

RESULTS

The highest resistance was observed to erythromycin and clindamycin at a rate of 76.4%, followed by ciprofloxacin (61.1%), gentamicin (54.2%), rifampin (38.9%), and co-trimoxazole (27.8%). All isolates were susceptible to vancomycin. The virulence genes of etb, pvl, tst, and eta were detected in 50%, 29.2%, 21.8%, and 13.9% of isolates, respectively. SCCmec types III and I were the most prevalent types, followed by types IV, II, and V. MLST analysis revealed 6 sequence types: ST6854, ST5282, ST127, ST7804, ST1607, and ST7784. Two MLST-based clonal complexes (CC8, and CC97) were identified as well.

CONCLUSION

The ST numbers were non-repetitive. CC8 as a pandemic clone and an individual lineage and clinically significant clade was reported as the most prevalent clonal complex. It is essential periodic evaluations of antibiotic susceptibility patterns and study the evolutionary characteristics of medical-challenging microorganisms in particular MRSA to effectively treat and restrict the outbreaks.

Staphylococcus brunensis sp. nov. isolated from human clinical specimens with a staphylococcal cassette chromosome-related genomic island outside of the rlmH gene bearing the ccrDE recombinase gene complex

Novel species of coagulase-negative staphylococci, which could serve as reservoirs of virulence and antimicrobial resistance factors for opportunistic pathogens from the genus Staphylococcus, are recognized in human and animal specimens due to advances in diagnostic techniques. Here, we used whole-genome sequencing, extensive biotyping, MALDI-TOF mass spectrometry, and chemotaxonomy to characterize five coagulase-negative strains from the Staphylococcus haemolyticus phylogenetic clade obtained from human ear swabs, wounds, and bile. Based on the results of polyphasic taxonomy, we propose the species Staphylococcus brunensis sp. nov. (type strain NRL/St 16/872T = CCM 9024T = LMG 31872T = DSM 111349T). The genomic analysis revealed numerous variable genomic elements, including staphylococcal cassette chromosome (SCC), prophages, plasmids, and a unique 18.8 kb-long genomic island SbCIccrDE integrated into the ribosomal protein L7 serine acetyltransferase gene rimL. SbCIccrDE has a cassette chromosome recombinase (ccr) gene complex with a typical structure found in SCCs. Based on nucleotide and amino acid identity to other known ccr genes and the distinct integration site that differs from the canonical methyltransferase gene rlmH exploited by SCCs, we classified the ccr genes as novel variants, ccrDE. The comparative genomic analysis of SbCIccrDE with related islands shows that they can accumulate virulence and antimicrobial resistance factors creating novel resistance elements, which reflects the evolution of SCC. The spread of these resistance islands into established pathogens such as Staphylococcus aureus would pose a great threat to the healthcare system. IMPORTANCE The coagulase-negative staphylococci are important opportunistic human pathogens, which cause bloodstream and foreign body infections, mainly in immunocompromised patients. The mobile elements, primarily the staphylococcal cassette chromosome mec, which confers resistance to methicillin, are the key to the successful dissemination of staphylococci into healthcare and community settings. Here, we present a novel species of the Staphylococcus genus isolated from human clinical material. The detailed analysis of its genome revealed a previously undescribed genomic island, which is closely related to the staphylococcal cassette chromosome and has the potential to accumulate and spread virulence and resistance determinants. The island harbors a set of conserved genes required for its mobilization, which we recognized as novel cassette chromosome recombinase genes ccrDE. Similar islands were revealed not only in the genomes of coagulase-negative staphylococci but also in S. aureus. The comparative genomic study contributes substantially to the understanding of the evolution and pathogenesis of staphylococci.

Atopic dermatitis-derived Staphylococcus aureus strains: what makes them special in the interplay with the host

Background

Atopic dermatitis (AD) is a chronic inflammatory skin condition whose pathogenesis involves genetic predisposition, epidermal barrier dysfunction, alterations in the immune responses and microbial dysbiosis. Clinical studies have shown a link between Staphylococcus aureus and the pathogenesis of AD, although the origins and genetic diversity of S. aureus colonizing patients with AD is poorly understood. The aim of the study was to investigate if specific clones might be associated with the disease.

Methods

WGS analyses were performed on 38 S. aureus strains, deriving from AD patients and healthy carriers. Genotypes (i.e. MLST, spa-, agr- and SCCmec-typing), genomic content (e.g. virulome and resistome), and the pan-genome structure of strains have been investigated. Phenotypic analyses were performed to determine the antibiotic susceptibility, the biofilm production and the invasiveness within the investigated S. aureus population.

Results

Strains isolated from AD patients revealed a high degree of genetic heterogeneity and a shared set of virulence factors and antimicrobial resistance genes, suggesting that no genotype and genomic content are uniquely associated with AD. The same strains were characterized by a lower variability in terms of gene content, indicating that the inflammatory conditions could exert a selective pressure leading to the optimization of the gene repertoire. Furthermore, genes related to specific mechanisms, like post-translational modification, protein turnover and chaperones as well as intracellular trafficking, secretion and vesicular transport, were significantly more enriched in AD strains. Phenotypic analysis revealed that all of our AD strains were strong or moderate biofilm producers, while less than half showed invasive capabilities.

Conclusions

We conclude that in AD skin, the functional role played by S. aureus may depend on differential gene expression patterns and/or on post-translational modification mechanisms rather than being associated with peculiar genetic features.

Leaky barriers to gene sharing between locally co-existing coagulase-negative Staphylococcus species

Coagulase-negative Staphylococcus (CoNS) are opportunistic pathogens implicated in many human and animal infections. The evolutionary history of CoNS remains obscure because of the historical lack of recognition for their clinical importance and poor taxonomic sampling. Here, we sequenced the genomes of 191 CoNS isolates representing 15 species sampled from diseased animals diagnosed in a veterinary diagnostic laboratory. We found that CoNS are important reservoirs of diverse phages, plasmids and mobilizable genes encoding antimicrobial resistance, heavy metal resistance, and virulence. Frequent exchange of DNA between certain donor-recipient partners suggests that specific lineages act as hubs of gene sharing. We also detected frequent recombination between CoNS regardless of their animal host species, indicating that ecological barriers to horizontal gene transfer can be surmounted in co-circulating lineages. Our findings reveal frequent but structured patterns of transfer that exist within and between CoNS species, which are driven by their overlapping ecology and geographical proximity.

Staphylococcus aureus host interactions and adaptation

Invasive Staphylococcus aureus infections are common, causing high mortality, compounded by the propensity of the bacterium to develop drug resistance. S. aureus is an excellent case study of the potential for a bacterium to be commensal, colonizing, latent or disease-causing; these states defined by the interplay between S. aureus and host. This interplay is multidimensional and evolving, exemplified by the spread of S. aureus between humans and other animal reservoirs and the lack of success in vaccine development. In this Review, we examine recent advances in understanding the S. aureus-host interactions that lead to infections. We revisit the primary role of neutrophils in controlling infection, summarizing the discovery of new immune evasion molecules and the discovery of new functions ascribed to well-known virulence factors. We explore the intriguing intersection of bacterial and host metabolism, where crosstalk in both directions can influence immune responses and infection outcomes. This Review also assesses the surprising genomic plasticity of S. aureus, its dualism as a multi-mammalian species commensal and opportunistic pathogen and our developing understanding of the roles of other bacteria in shaping S. aureus colonization.

Integrated conjugative plasmid drives high frequency chromosomal gene transfer in Sulfolobus islandicus

Gene transfer in crenarchaea has been observed within natural and experimental populations of Sulfolobus. However, the molecular factors that govern how gene transfer and recombination manifest themselves in these populations is still unknown. In this study, we examine a plasmid-mediated mechanism of gene transfer in S. islandicus that results in localized high frequency recombination within the chromosome. Through chromosomal marker exchange assays with defined donors and recipients, we find that while bidirectional exchange occurs among all cells, those possessing the integrated conjugative plasmid, pM164, mobilize a nearby locus at a significantly higher frequency when compared to a more distal marker. We establish that traG is essential for this phenotype and that high frequency recombination can be replicated in transconjugants after plasmid transfer. Mapping recombinants through genomic analysis, we establish the distribution of recombinant tracts with decreasing frequency at increasing distance from pM164. We suggest the bias in transfer is a result of an Hfr (high frequency recombination)-like conjugation mechanism in this strain. In addition, we find recombinants containing distal non-selected recombination events, potentially mediated by a different host-encoded marker exchange (ME) mechanism.

Comparative Genomics of Streptococcus oralis Identifies Large Scale Homologous Recombination and a Genetic Variant Associated with Infection

The viridans group streptococci (VGS) are a large consortium of commensal streptococci that colonize the human body. Many species within this group are opportunistic pathogens causing bacteremia and infective endocarditis (IE), yet little is known about why some strains cause invasive disease. Identification of virulence determinants is complicated by the difficulty of distinguishing between the closely related species of this group. Here, we analyzed genomic data from VGS that were isolated from blood cultures in patients with invasive infections and oral swabs of healthy volunteers and then determined the best-performing methods for species identification. Using whole-genome sequence data, we characterized the population structure of a diverse sample of Streptococcus oralis isolates and found evidence of frequent recombination. We used multiple genome-wide association study tools to identify candidate determinants of invasiveness. These tools gave consistent results, leading to the discovery of a single synonymous single nucleotide polymorphism (SNP) that was significantly associated with invasiveness. This SNP was within a previously undescribed gene that was conserved across the majority of VGS species. Using the growth in the presence of human serum and a simulated infective endocarditis vegetation model, we were unable to identify a phenotype for the enriched allele in laboratory assays, suggesting a phenotype may be specific to natural infection. These data highlighted the power of analyzing natural populations for gaining insight into pathogenicity, particularly for organisms with complex population structures like the VGS. IMPORTANCE The viridians group streptococci (VGS) are a large collection of closely related commensal streptococci, with many being opportunistic pathogens causing invasive diseases, such as bacteremia and infective endocarditis. Little is known about virulence determinants in these species, and there is a distinct lack of genomic information available for the VGS. In this study, we collected VGS isolates from invasive infections and healthy volunteers and performed whole-genome sequencing for a suite of downstream analyses. We focused on a diverse sample of Streptococcus oralis genomes and identified high rates of recombination in the population as well as a single genome variant highly enriched in invasive isolates. The variant lies within a previously uncharacterized gene, nrdM, which shared homology with the anaerobic ribonucleoside triphosphate reductase, nrdD, and was highly conserved among VGS. This work increased our knowledge of VGS genomics and indicated that differences in virulence potential among S. oralis isolates were, at least in part, genetically determined.

ggMOB: Elucidation of genomic conjugative features and associated cargo genes across bacterial genera using genus-genus mobilization networks

Horizontal gene transfer mediated by conjugation is considered an important evolutionary mechanism of bacteria. It allows organisms to quickly evolve new phenotypic properties including antimicrobial resistance (AMR) and virulence. The frequency of conjugation-mediated cargo gene exchange has not yet been comprehensively studied within and between bacterial taxa. We developed a frequency-based network of genus-genus conjugation features and candidate cargo genes from whole-genome sequence data of over 180,000 bacterial genomes, representing 1,345 genera. Using our method, which we refer to as ggMOB, we revealed that over half of the bacterial genomes contained one or more known conjugation features that matched exactly to at least one other genome. Moreover, the proportion of genomes containing these conjugation features varied substantially by genus and conjugation feature. These results and the genus-level network structure can be viewed interactively in the ggMOB interface, which allows for user-defined filtering of conjugation features and candidate cargo genes. Using the network data, we observed that the ratio of AMR gene representation in conjugative versus non-conjugative genomes exceeded 5:1, confirming that conjugation is a critical force for AMR spread across genera. Finally, we demonstrated that clustering genomes by conjugation profile sometimes correlated well with classical phylogenetic structuring; but that in some cases the clustering was highly discordant, suggesting that the importance of the accessory genome in driving bacterial evolution may be highly variable across both time and taxonomy. These results can advance scientific understanding of bacterial evolution, and can be used as a starting point for probing genus-genus gene exchange within complex microbial communities that include unculturable bacteria. ggMOB is publicly available under the GNU licence at https://ruiz-hci-lab.github.io/ggMOB/.

How to survive pig farming: Mechanism of SCCmec element deletion and metabolic stress adaptation in livestock-associated MRSA

Previous research on methicillin susceptible Staphylococcus aureus (MSSA) belonging to livestock-associated (LA-) sequence type (ST) 398, isolated from pigs and their local surroundings, indicated that differences between these MSSA and their methicillin resistant predecessors (MRSA) are often limited to the absence of the staphylococcal cassette chromosome mec (SCCmec) and few single nucleotide polymorphisms. So far, our understanding on how LA-MRSA endure the environmental conditions associated with pig-farming as well as the putative impact of this particular environment on the mobilisation of SCCmec elements is limited. Thus, we performed in-depth genomic and transcriptomic analyses using the LA-MRSA ST398 strain IMT38951 and its methicillin susceptible descendant. We identified a mosaic-structured SCCmec region including a putative replicative SCCmecVc which is absent from the MSSA chromosome through homologous recombination. Based on our data, such events occur between short repetitive sequences identified within and adjacent to two distinct alleles of the large cassette recombinase genes C (ccrC). We further evaluated the global transcriptomic response of MRSA ST398 to particular pig-farm associated conditions, i.e., contact with host proteins (porcine serum) and a high ammonia concentration. Differential expression of global regulators involved in stress response control were identified, i.e., ammonia-induced alternative sigma factor B-depending activation of genes for the alkaline shock protein 23, the heat shock response and the accessory gene regulator (agr)-controlled transcription of virulence factors. Exposure to serum transiently induced the transcription of distinct virulence factor encoding genes. Transcription of genes reported for mediating the loss of methicillin resistance, especially ccrC, was not significantly different compared to the unchallenged controls. We concluded that, from an evolutionary perspective, bacteria may save energy by incidentally dismissing a fully replicative SCCmec element in contrast to the induction of ccr genes on a population scale. Since the genomic SCCmec integration site is a hot-spot of recombination, occasional losses of elements of 16 kb size may restore capacities for the uptake of foreign genetic material. Subsequent spread of resistance, on the other hand, might depend on the autonomous replication machinery of the deleted SCCmec elements that probably enhance chances for reintegration of SCCmec into susceptible genomes by mere multiplication.

Molecular Epidemiology of Staphylococcus aureus in China Reveals the Key Gene Features Involved in Epidemic Transmission and Adaptive Evolution

Staphylococcus aureus is a Gram-positive pathogen that causes various infections in humans and domestic animals. In China, S. aureus is the most common Gram-positive pathogen that causes clinical infections. However, there are few comprehensive genome-based molecular epidemiology studies to investigate the genotypic background of the major S. aureus clones that are epidemic in China. Here, four S. aureus isolates that were recovered from hospital personnel were sequenced. In combination with whole-genome sequencing (WGS) data of 328 S. aureus strains as references, we performed a comprehensive molecular epidemiology study to reveal the molecular epidemic characterization of S. aureus that is epidemic in China. It was found that 332 S. aureus isolates were phylogenetically categorized into 4 major epidemic groups with different epidemiology phenotypes. Each group has exclusive features in virulence genotypic profiles, antimicrobial resistance genotypic profiles, core and pangenome features representing the differences involved in genetic features, evolutionary processes, and potential future evolutionary directions. Moreover, a comparative core genome analysis of 332 S. aureus isolates indicated several key genes that contributed to differences in molecular epidemic characterization and promoted the adaptive evolutionary process of each group. This study provides a comprehensive understanding of molecular epidemiological characteristics and adaptive evolutionary directions of major S. aureus clones that are epidemic in China. IMPORTANCE Staphylococcus aureus is an important Gram-positive pathogen that is epidemic worldwide and causes various infections in humans and domestic animals. However, there has been relatively little research on comprehensive molecular epidemiology in China. In this research, we reconstructed the phylogenetic relationship based on whole-genome data of strains almost all over China, screened for resistance and virulence genes, and took core and pan genome analysis to perform a comprehensive molecular epidemiology study of S. aureus that is epidemic in China. Our results highlight that there are 4 major epidemic groups with different epidemiology phenotypes after phylogenetic categorization with exclusive genetic features in virulence genotypic profiles, antimicrobial-resistance genotypic profiles, and core and pangenome features, and we found key gene features involved in epidemic transmission and adaptive evolution. Our findings are critical in describing molecular characteristic profiles of S. aureus infection, which could update existing preventive measures and take appropriate strategies.

A scalable analytical approach from bacterial genomes to epidemiology

Recent years have seen a remarkable increase in the practicality of sequencing whole genomes from large numbers of bacterial isolates. The availability of this data has huge potential to deliver new insights into the evolution and epidemiology of bacterial pathogens, but the scalability of the analytical methodology has been lagging behind that of the sequencing technology. Here we present a step-by-step approach for such large-scale genomic epidemiology analyses, from bacterial genomes to epidemiological interpretations. A central component of this approach is the dated phylogeny, which is a phylogenetic tree with branch lengths measured in units of time. The construction of dated phylogenies from bacterial genomic data needs to account for the disruptive effect of recombination on phylogenetic relationships, and we describe how this can be achieved. Dated phylogenies can then be used to perform fine-scale or large-scale epidemiological analyses, depending on the proportion of cases for which genomes are available. A key feature of this approach is computational scalability and in particular the ability to process hundreds or thousands of genomes within a matter of hours. This is a clear advantage of the step-by-step approach described here. We discuss other advantages and disadvantages of the approach, as well as potential improvements and avenues for future research. This article is part of a discussion meeting issue 'Genomic population structures of microbial pathogens'.

Metagenomic methylation patterns resolve bacterial genomes of unusual size and structural complexity

The plasticity of bacterial and archaeal genomes makes examining their ecological and evolutionary dynamics both exciting and challenging. The same mechanisms that enable rapid genomic change and adaptation confound current approaches for recovering complete genomes from metagenomes. Here, we use strain-specific patterns of DNA methylation to resolve complex bacterial genomes from long-read metagenomic data of a marine microbial consortium, the “pink berries” of the Sippewissett Marsh (USA). Unique combinations of restriction-modification (RM) systems encoded by the bacteria produced distinctive methylation profiles that were used to accurately bin and classify metagenomic sequences. Using this approach, we finished the largest and most complex circularized bacterial genome ever recovered from a metagenome (7.9 Mb with >600 transposons), the finished genome of Thiohalocapsa sp. PB-PSB1 the dominant bacteria in the consortia. From genomes binned by methylation patterns, we identified instances of horizontal gene transfer between sulfur-cycling symbionts (Thiohalocapsa sp. PB-PSB1 and Desulfofustis sp. PB-SRB1), phage infection, and strain-level structural variation. We also linked the methylation patterns of each metagenome-assembled genome with encoded DNA methyltransferases and discovered new RM defense systems, including novel associations of RM systems with RNase toxins.

Core genes can have higher recombination rates than accessory genes within global microbial populations

Recombination is essential to microbial evolution, and is involved in the spread of antibiotic resistance, antigenic variation, and adaptation to the host niche. However, assessing the impact of homologous recombination on accessory genes which are only present in a subset of strains of a given species remains challenging due to their complex phylogenetic relationships. Quantifying homologous recombination for accessory genes (which are important for niche-specific adaptations) in comparison to core genes (which are present in all strains and have essential functions) is critical to understanding how selection acts on variation to shape species diversity and genome structures of bacteria. Here, we apply a computationally efficient, non-phylogenetic approach to measure homologous recombination rates in the core and accessory genome using >100,000 whole genome sequences from Streptococcus pneumoniae and several additional species. By analyzing diverse sets of sequence clusters, we show that core genes often have higher recombination rates than accessory genes, and for some bacterial species the associated effect sizes for these differences are pronounced. In a subset of species, we find that gene frequency and homologous recombination rate are positively correlated. For S. pneumoniae and several additional species, we find that while the recombination rate is higher for the core genome, the mutational divergence is lower, indicating that divergence-based homologous recombination barriers could contribute to differences in recombination rates between the core and accessory genome. Homologous recombination may therefore play a key role in increasing the efficiency of selection in the most conserved parts of the genome.

From a large-scale genomic analysis of insertion sequences to insights into their regulatory roles in prokaryotes

Background

Insertion sequences (ISs) are mobile repeat sequences and most of them can copy themselves to new host genome locations, leading to genome plasticity and gene regulation in prokaryotes. In this study, we present functional and evolutionary relationships between IS and neighboring genes in a large-scale comparative genomic analysis.

Results

IS families were located in all prokaryotic phyla, with preferential occurrence of IS3, IS4, IS481, and IS5 families in Alpha-, Beta-, and Gammaproteobacteria, Actinobacteria and Firmicutes as well as in eukaryote host-associated organisms and autotrophic opportunistic pathogens. We defined the concept of the IS-Gene couple (IG), which allowed to highlight the functional and regulatory impacts of an IS on the closest gene. Genes involved in transcriptional regulation and transport activities were found overrepresented in IG. In particular, major facilitator superfamily (MFS) transporters, ATP-binding proteins and transposases raised as favorite neighboring gene functions of IS hotspots. Then, evolutionary conserved IS-Gene sets across taxonomic lineages enabled the classification of IS-gene couples into phylum, class-to-genus, and species syntenic IS-Gene couples. The IS5, IS21, IS4, IS607, IS91, ISL3 and IS200 families displayed two to four times more ISs in the phylum and/or class-to-genus syntenic IGs compared to other IS families. This indicates that those families were probably inserted earlier than others and then subjected to horizontal transfer, transposition and deletion events over time. In phylum syntenic IG category, Betaproteobacteria, Crenarchaeota, Calditrichae, Planctomycetes, Acidithiobacillia and Cyanobacteria phyla act as IS reservoirs for other phyla, and neighboring gene functions are mostly related to transcriptional regulators. Comparison of IS occurrences with predicted regulatory motifs led to ~ 26.5% of motif-containing ISs with 2 motifs per IS in average. These results, concomitantly with short IS-Gene distances, suggest that those ISs would interfere with the expression of neighboring genes and thus form strong candidates for an adaptive pairing.

Conclusions

All together, our large-scale study provide new insights into the IS genetic context and strongly suggest their regulatory roles.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08678-3.

Landscape of mobile genetic elements and their antibiotic resistance cargo in prokaryotic genomes

Prokaryotic Mobile Genetic Elements (MGEs) such as transposons, integrons, phages and plasmids, play important roles in prokaryotic evolution and in the dispersal of cargo functions like antibiotic resistance. However, each of these MGE types is usually annotated and analysed individually, hampering a global understanding of phylogenetic and environmental patterns of MGE dispersal. We thus developed a computational framework that captures diverse MGE types, their cargos and MGE-mediated horizontal transfer events, using recombinases as ubiquitous MGE marker genes and pangenome information for MGE boundary estimation. Applied to ∼84k genomes with habitat annotation, we mapped 2.8 million MGE-specific recombinases to six operational MGE types, which together contain on average 13% of all the genes in a genome. Transposable elements (TEs) dominated across all taxa (∼1.7 million occurrences), outnumbering phages and phage-like elements (<0.4 million). We recorded numerous MGE-mediated horizontal transfer events across diverse phyla and habitats involving all MGE types, disentangled and quantified the extent of hitchhiking of TEs (17%) and integrons (63%) with other MGE categories, and established TEs as dominant carriers of antibiotic resistance genes. We integrated all these findings into a resource (proMGE.embl.de), which should facilitate future studies on the large mobile part of genomes and its horizontal dispersal.

Cryptic niche differentiation of novel sediment ecotypes of Rugeria pomeroyi correlates with nitrate respiration

Marine intertidal sediments fluctuate in redox conditions and nutrient availability, and they are also known as an important sink of nitrogen mainly through denitrification, yet how denitrifying bacteria adapt to this dynamic habitat remains largely untapped. Here, we investigated novel intertidal benthic ecotypes of the model pelagic marine bacterium Ruegeria pomeroyi DSS-3 with a population genomic approach. While differing by only 1.3% at the 16S rRNA gene level, members of the intertidal benthic ecotypes are complete denitrifiers whereas the pelagic ecotype representative (DSS-3) is a partial denitrifier lacking a nitrate reductase. The intertidal benthic ecotypes are further differentiated by using non-homologous nitrate reductases and a different set of genes that allow alleviating oxidative stress and acquiring organic substrates. In the presence of nitrate, the two ecotypes showed contrasting growth patterns under initial oxygen concentrations at 1 vol% versus 7 vol% and supplemented with different carbon sources abundant in intertidal sediments. Collectively, this combination of evidence indicates that there are cryptic niches in coastal intertidal sediments that support divergent evolution of denitrifying bacteria. This knowledge will in turn help understand how these benthic environments operate to effectively remove nitrogen.

Evolutionary Processes Driving the Rise and Fall of Staphylococcus aureus ST239, a Dominant Hybrid Pathogen

Selection plays a key role in the spread of antibiotic resistance, but the evolutionary drivers of clinically important resistant strains remain poorly understood. Here, we use genomic analyses and competition experiments to study Staphylococcus aureus ST239, a prominent MRSA strain that is thought to have been formed by large-scale recombination between ST8 and ST30. Genomic analyses allowed us to refine the hybrid model for the origin of ST239 and to date the origin of ST239 to 1920 to 1945, which predates the clinical introduction of methicillin in 1959. Although purifying selection has dominated the evolution of ST239, parallel evolution has occurred in genes involved in antibiotic resistance and virulence, suggesting that ST239 has evolved toward an increasingly pathogenic lifestyle. Crucially, ST239 isolates have low competitive fitness relative to both ST8 and ST30 isolates, supporting the idea that fitness costs have driven the demise of this once-dominant pathogen strain.

Quantifying bacterial evolution in the wild: A birthday problem for Campylobacter lineages

Measuring molecular evolution in bacteria typically requires estimation of the rate at which nucleotide changes accumulate in strains sampled at different times that share a common ancestor. This approach has been useful for dating ecological and evolutionary events that coincide with the emergence of important lineages, such as outbreak strains and obligate human pathogens. However, in multi-host (niche) transmission scenarios, where the pathogen is essentially an opportunistic environmental organism, sampling is often sporadic and rarely reflects the overall population, particularly when concentrated on clinical isolates. This means that approaches that assume recent common ancestry are not applicable. Here we present a new approach to estimate the molecular clock rate in Campylobacter that draws on the popular probability conundrum known as the 'birthday problem'. Using large genomic datasets and comparative genomic approaches, we use isolate pairs that share recent common ancestry to estimate the rate of nucleotide change for the population. Identifying synonymous and non-synonymous nucleotide changes, both within and outside of recombined regions of the genome, we quantify clock-like diversification to estimate synonymous rates of nucleotide change for the common pathogenic bacteria Campylobacter coli (2.4 x 10-6 s/s/y) and Campylobacter jejuni (3.4 x 10-6 s/s/y). Finally, using estimated total rates of nucleotide change, we infer the number of effective lineages within the sample time frame-analogous to a shared birthday-and assess the rate of turnover of lineages in our sample set over short evolutionary timescales. This provides a generalizable approach to calibrating rates in populations of environmental bacteria and shows that multiple lineages are maintained, implying that large-scale clonal sweeps may take hundreds of years or more in these species.

Multilocus sequence analysis reveals genetic diversity in Staphylococcus aureus isolate of goat with mastitis persistent after treatment with enrofloxacin

Staphylococcus aureus is one of the main bacterial agents responsible for cases of mastitis in ruminants, playing an important role in the persistence and chronicity of diseases treated with antimicrobials. Using the multilocus sequence typing technique, network approaches and study of the population diversity of microorganisms, we performed analyzes of S. aureus (ES-GPM) isolated from goats with persistent mastitis (GPM). The most strains of ES-GPM were categorically different phylogenetically from the others and could be divided into two lineages: one with a majority belonging to ES-GPM and the other to varied strains. These two lineages were separated by 27 nuclear polymorphisms. The 43 strains comprised 22 clonal complexes (CCs), of which the ES-GPM strains were present in CC133, CC5 and a new complex formed by the sequence type 4966. The genetic diversity of some alleles showed be greater diversity and polymorphism than others, such as of the aroE and yqiL genes less than glpF gene. In addition, the sequences ES-GPM to the arc gene and glpF alleles showed the greatest number of mutations for ES-GPM in relation to non-ES-GPM. Therefore, this study identified genetic polymorphisms characteristic of S. aureus isolated from milk of goats diagnosed with persistent mastitis after the failed treatment with the antibiotic enrofloxacin. This study may help in the future to identify and discriminate this agent in cases of mastitis, and with that, the most appropriate antibiotic treatment can be performed in advance of the appearance of persistent mastitis caused by the agent, reducing the chances of premature culling and animal suffering.

The global population structure and beta-lactamase repertoire of the opportunistic pathogen Serratia marcescens

Serratia marcescens is a global spread nosocomial pathogen. This rod-shaped bacterium displays a broad host range and worldwide geographical distribution. Here we analyze an international collection of this multidrug-resistant, opportunistic pathogen from 35 countries to infer its population structure. We show that S. marcescens comprises 12 lineages; Sm1, Sm4, and Sm10 harbor 78.3% of the known environmental strains. Sm5, Sm6, and Sm7 comprise only human-associated strains which harbor smallest pangenomes, genomic fluidity and lowest levels of core recombination, indicating niche specialization. Sm7 and Sm9 lineages exhibit the most concerning resistome; bla_KPC-2 plasmid is widespread in Sm7, whereas Sm9, also an anthropogenic-exclusive lineage, presents highest plasmid/lineage size ratio and plasmid-diversity encoding metallo-beta-lactamases comprising bla_NDM-1. The heterogeneity of resistance patterns of S. marcescens lineages elucidated herein highlights the relevance of surveillance programs, using whole-genome sequencing, to provide insights into the molecular epidemiology of carbapenemase producing strains of this species.

Characterizing a novel SCCmec with a composite structure from a clinical strain of Staphylococcus hominis, C34847

Staphylococcal cassette chromosome mec (SCCmec) has predominantly been described in methicillin-resistant Staphylococcus aureus. However, studies have indicated that coagulase-negative staphylococci (CoNS) carry a larger diversity of SCC elements. We characterized a composite SCCmec element carrying an uncharacterized ccr1 and type A mec gene combination, in conjunction with a secondary element bearing ccr4, from a clinical strain of S. hominis. The element's complex structure points to a high degree of recombination occurring in SCCmec in CoNS.

Overview of the risks of Staphylococcus aureus infections and their control by bacteriophages and bacteriophage-encoded products

Staphylococcus aureus is the leading cause of secondary infections in hospitals and a challenging pathogen in food industries. Decades after it was first reported, β-lactam-resistant S. aureus remains a subject of intense research owing to the ever-increasing issue of drug resistance. S. aureus bacteriophages (phages) or their encoded products are considered an alternative to antibiotics as they have been shown to be effective in treating some S. aureus-associated infections. In this review, we present a concise collection of the literature on the pathogenic potential of S. aureus and examine the prospects of using S. aureus phages and their encoded products as antimicrobials.

The Phosphoarginine Phosphatase PtpB from Staphylococcus aureus Is Involved in Bacterial Stress Adaptation during Infection

Staphylococcus aureus continues to be a public health threat, especially in hospital settings. Studies aimed at deciphering the molecular and cellular mechanisms that underlie pathogenesis, host adaptation, and virulence are required to develop effective treatment strategies. Numerous host-pathogen interactions were found to be dependent on phosphatases-mediated regulation. This study focused on the analysis of the role of the low-molecular weight phosphatase PtpB, in particular, during infection. Deletion of ptpB in S. aureus strain SA564 significantly reduced the capacity of the mutant to withstand intracellular killing by THP-1 macrophages. When injected into normoglycemic C57BL/6 mice, the SA564 ΔptpB mutant displayed markedly reduced bacterial loads in liver and kidney tissues in a murine S. aureus abscess model when compared to the wild type. We also observed that PtpB phosphatase-activity was sensitive to oxidative stress. Our quantitative transcript analyses revealed that PtpB affects the transcription of various genes involved in oxidative stress adaptation and infectivity. Thus, this study disclosed first insights into the physiological role of PtpB during host interaction allowing us to link phosphatase-dependent regulation to oxidative bacterial stress adaptation during infection.

Prevalence and distribution of multilocus sequence types of Staphylococcus aureus isolated from bulk tank milk and cows with mastitis in Pennsylvania

A total of 163 S. aureus isolates; 113 from mastitic milk (MM) and 50 from bulk tank milk (BTM) (2008, 2013–2015) submitted for bacteriologic analysis at the Penn State Animal Diagnostic Laboratory were examined for their phenotypic and genotypic characteristics. Multi-locus sequence typing (MLST) analysis identified 16 unique sequence types (STs) which belonged to eight clonal complexes (CCs). Majority of the isolates were variants of CC97 (68.7%) and CC151 (25.1%). CC97 comprised of seven STs, of which two were new STs (ST3273, ST3274), while CC151 comprised of three STs of which ST3272 was identified for the first time. Several farms had more than one ST type that were either members of the same clonal complex or unrelated STs. On one farm, six different STs of both categories were seen over the years within the farm. It was observed that ST352 and ST151 were the two main clonal populations in cattle not only in Pennsylvania but also globally. Most isolates were susceptible to all the antibiotics evaluated. 6.7% of isolates showed resistance to vancomycin and penicillin. Two isolates of ST398 showed multidrug resistance (>3 antibiotics) against clindamycin, erythromycin, tetracycline, and penicillin. It was noted that 59 of 163 (36.2%) isolates encoded for enterotoxigenic genes. Enterotoxin genes seg/sei accounted for ~85% of enterotoxin positive isolates. Toxic shock syndrome gene tsst-1 alone was positive in two isolates (ST352, ST 2187). 97.5% of CC151 isolates were enterotoxin seg/sei positive. Most isolates were positive for lukED (95%) and lukAB (96.3%) leukotoxin genes. Bovine specific bi-component leucocidin lukMF’ was present in 54% of isolates. A prominent observation of this study was the explicit association of lukMF’ with lineages ST151 and ST352. In conclusion, the findings of the study, suggest that small number of S. aureus STs types (ST352, ST2187, ST3028, and ST151) are associated with majority of cases of bovine mastitis in Pennsylvania dairy farms. It was observed that one ST of S. aureus predominated in the herd and this ST can coexist with several other ST types of S. aureus strains. When STs were interpreted along with virulence, leucocidin genes and antimicrobial resistance, ST-variants allowed better interpretation of the S. aureus molecular epidemiologic findings specifically for tracing recurrence or persistence of infections in cow over time, among cows in the herd, and between herds in Pennsylvania.

Genes Influencing Phage Host Range in Staphylococcus aureus on a Species-Wide Scale

Staphylococcus aureus is a human pathogen that causes serious diseases, ranging from skin infections to septic shock. Bacteriophages (phages) are both natural killers of S. aureus, offering therapeutic possibilities, and important vectors of horizontal gene transfer (HGT) in the species. Here, we used high-throughput approaches to understand the genetic basis of strain-to-strain variation in sensitivity to phages, which defines the host range. We screened 259 diverse S. aureus strains covering more than 40 sequence types for sensitivity to eight phages, which were representatives of the three phage classes that infect the species. The phages were variable in host range, each infecting between 73 and 257 strains. Using genome-wide association approaches, we identified putative loci that affect host range and validated their function using USA300 transposon knockouts. In addition to rediscovering known host range determinants, we found several previously unreported genes affecting bacterial growth during phage infection, including trpA, phoR, isdB, sodM, fmtC, and relA We used the data from our host range matrix to develop predictive models that achieved between 40% and 95% accuracy. This work illustrates the complexity of the genetic basis for phage susceptibility in S. aureus but also shows that with more data, we may be able to understand much of the variation. With a knowledge of host range determination, we can rationally design phage therapy cocktails that target the broadest host range of S. aureus strains and address basic questions regarding phage-host interactions, such as the impact of phage on S. aureus evolution.IMPORTANCEStaphylococcus aureus is a widespread, hospital- and community-acquired pathogen, many strains of which are antibiotic resistant. It causes diverse diseases, ranging from local to systemic infection, and affects both the skin and many internal organs, including the heart, lungs, bones, and brain. Its ubiquity, antibiotic resistance, and disease burden make new therapies urgent. One alternative therapy to antibiotics is phage therapy, in which viruses specific to infecting bacteria clear infection. In this work, we identified and validated S. aureus genes that influence phage host range-the number of strains a phage can infect and kill-by testing strains representative of the diversity of the S. aureus species for phage host range and associating the genome sequences of strains with host range. These findings together improved our understanding of how phage therapy works in the bacterium and improve prediction of phage therapy efficacy based on the predicted host range of the infecting strain.

Global Epidemiology and Evolutionary History of Staphylococcus aureus ST45

Staphylococcus aureus ST45 is a major global MRSA lineage with huge strain diversity and a high clinical impact. It is one of the most prevalent carrier lineages but also frequently causes severe invasive disease, such as bacteremia. Little is known about its evolutionary history. In this study, we used whole-genome sequencing to analyze a large collection of 451 diverse ST45 isolates from 6 continents and 26 countries. De novo-assembled genomes were used to understand genomic plasticity and to perform coalescent analyses. The ST45 population contained two distinct sublineages, which correlated with the isolates' geographical origins. One sublineage primarily consisted of European/North American isolates, while the second sublineage primarily consisted of African and Australian isolates. Bayesian analysis predicted ST45 originated in northwestern Europe about 500 years ago. Isolation time, host, and clinical symptoms did not correlate with phylogenetic groups. Our phylogenetic analyses suggest multiple acquisitions of the SCCmec element and key virulence factors throughout the evolution of the ST45 lineage.

Impact of homologous recombination on core genome phylogenies

BACKGROUND

Core genome phylogenies are widely used to build the evolutionary history of individual prokaryote species. By using hundreds or thousands of shared genes, these approaches are the gold standard to reconstruct the relationships of large sets of strains. However, there is growing evidence that bacterial strains exchange DNA through homologous recombination at rates that vary widely across prokaryote species, indicating that core genome phylogenies might not be able to reconstruct true phylogenies when recombination rate is high. Few attempts have been made to evaluate the robustness of core genome phylogenies to recombination, but some analyses suggest that reconstructed trees are not always accurate.

RESULTS

In this study, we tested the robustness of core genome phylogenies to various levels of recombination rates. By analyzing simulated and empirical data, we observed that core genome phylogenies are relatively robust to recombination rates; nevertheless, our results suggest that many reconstructed trees are not completely accurate even when bootstrap supports are high. We found that some core genome phylogenies are highly robust to recombination whereas others are strongly impacted by it, and we identified that the robustness of core genome phylogenies to recombination is highly linked to the levels of selective pressures acting on a species. Stronger selective pressures lead to less accurate tree reconstructions, presumably because selective pressures more strongly bias the routes of DNA transfers, thereby causing phylogenetic artifacts.

CONCLUSIONS

Overall, these results have important implications for the application of core genome phylogenies in prokaryotes.

Staphylococcal Protein A (spa) Locus Is a Hot Spot for Recombination and Horizontal Gene Transfer in Staphylococcus pseudintermedius

Staphylococcus pseudintermedius is a major canine pathogen but can also occasionally infect humans. Identification of genetic factors contributing to the virulence and clonal success of multidrug-resistant S. pseudintermedius clones is critical for the development of therapeutics against this pathogen. Here, we characterized the genome sequences of a global collection of 622 S. pseudintermedius isolates. We show that all major clones, besides carrying core virulence genes, which are present in all strains, carry one or more lineage-specific genes. Many of these genes have been acquired from other bacterial species through a horizontal gene transfer mechanism. Importantly, we have discovered that the staphylococcal protein A gene (spa), a widely used marker for molecular typing of S. pseudintermedius strains and a potential vaccine candidate antigen, is deleted in 62% of strains. Furthermore, the spa locus in S. pseudintermedius acts as a reservoir to accumulate lineage-associated genes with adaptive functions.

Staphylococcus pseudintermedius is a major canine pathogen but also occasionally colonizes and infects humans. Multidrug-resistant methicillin-resistant S. pseudintermedius (MDR MRSP) strains have emerged globally, making treatment and control of this pathogen challenging. Sequence type 71 (ST71), ST68, and ST45 are the most widespread and successful MDR MRSP clones. The potential genetic factors underlying the clonal success of these and other predominant clones remain unknown. Characterization of the pangenome, lineage-associated accessory genes, and genes acquired through horizontal gene transfer from other bacteria is important for identifying such factors. Here, we analyzed genome sequence data from 622 S. pseudintermedius isolates to investigate the evolution of pathogenicity across lineages. We show that the predominant clones carry one or more lineage-associated virulence genes. The gene encoding staphylococcal protein A (SpA), a key virulence factor involved in immune evasion and a potential vaccine antigen, is deleted in 62% of isolates. Most importantly, we have discovered that the spa locus is a hot spot for recombination and horizontal gene transfer in S. pseudintermedius, where genes related to restriction modification, prophage immunity, mercury resistance, and nucleotide and carbohydrate metabolism have been acquired in different lineages. Our study also establishes that ST45 is composed of two distinct sublineages that differ in their accessory gene content and virulence potential. Collectively, this study reports several previously undetected lineage-associated genetic factors that may have a role in the clonal success of the major MDR MRSP clones. These data provide a framework for future experimental studies on S. pseudintermedius pathogenesis and for developing novel therapeutics against this pathogen.

IMPORTANCEStaphylococcus pseudintermedius is a major canine pathogen but can also occasionally infect humans. Identification of genetic factors contributing to the virulence and clonal success of multidrug-resistant S. pseudintermedius clones is critical for the development of therapeutics against this pathogen. Here, we characterized the genome sequences of a global collection of 622 S. pseudintermedius isolates. We show that all major clones, besides carrying core virulence genes, which are present in all strains, carry one or more lineage-specific genes. Many of these genes have been acquired from other bacterial species through a horizontal gene transfer mechanism. Importantly, we have discovered that the staphylococcal protein A gene (spa), a widely used marker for molecular typing of S. pseudintermedius strains and a potential vaccine candidate antigen, is deleted in 62% of strains. Furthermore, the spa locus in S. pseudintermedius acts as a reservoir to accumulate lineage-associated genes with adaptive functions.

Molecular investigations on a chimeric strain of Staphylococcus aureus sequence type 80

A PVL-positive, methicillin-susceptible Staphylococcus aureus was cultured from pus from cervical lymphadenitis of a patient of East-African origin. Microarray hybridisation assigned the isolate to clonal complex (CC) 80 but revealed unusual features, including the presence of the ORF-CM14 enterotoxin homologue and of an ACME-III element as well as the absence of etD and edinB. The isolate was subjected to both, Illumina and Nanopore sequencing allowing characterisation of deviating regions within the strain´s genome. Atypical features of this strain were attributable to the presence of two genomic regions that originated from other S. aureus lineages and that comprised, respectively, 3% and 1.4% of the genome. One deviating region extended from walJ to sirB. It comprised ORF-CM14 and the ACME-III element. A homologous but larger fragment was also found in an atypical S. aureus CC1/ST567 strain whose lineage might have served as donor of this genomic region. This region itself is a chimera comprising fragments from CC1 as well as fragments of unknown origin. The other deviating region comprised the region from htsB to ecfA2, i.e., another 3% of the genome. It was very similar to CC1 sequences. Either this suggests an incorporation of CC1 DNA into the study strain, or alternatively a recombination event affecting “canonical” CC80. Thus, the study strain bears witness of several recombination events affecting supposedly core genomic genes. Although the exact mechanism is not yet clear, such chimerism seems to be an additional pathway in the evolution of S. aureus. This could facilitate also a transmission of virulence and resistance factors and therefore offer an additional evolutionary advantage.

Mobile genetic elements responsible for discordant Staphylococcus aureus phenotypes and genotypes in the same blood culture bottle

Approximately 15-20% of the S. aureus genome contains mobile genetic elements that can cause discrepancies between phenotypic and genotypic identification methods. Three blood culture bottles (each from a different patient) that showed discordant results, were shown to contain 2 S. aureus isolates after additional subcultures. One bottle had MRSA and MSSA that by DNA sequence analysis differed only by 31 kb; however, the deletions encompassed parts of SCCmec including mecA and SCC_M1. The second bottle contained MRSA and MSSA that differed by 124 kb; the MSSA was missing the entire SCCmec and spa regions. The last bottle contained 2 MRSA, one with ACME II disrupting SCCmec and a 24 bp spa deletion. The deletions in SCCmec and the other elements gave rise to the discrepancies between molecular and the original culture results. Such discrepancies should prompt a search for additional strains in the blood culture bottle.

Efflux pump activity potentiates the evolution of antibiotic resistance across S. aureus isolates

The rise of antibiotic resistance in many bacterial pathogens has been driven by the spread of a few successful strains, suggesting that some bacteria are genetically pre-disposed to evolving resistance. Here, we test this hypothesis by challenging a diverse set of 222 isolates of Staphylococcus aureus with the antibiotic ciprofloxacin in a large-scale evolution experiment. We find that a single efflux pump, norA, causes widespread variation in evolvability across isolates. Elevated norA expression potentiates evolution by increasing the fitness benefit provided by DNA topoisomerase mutations under ciprofloxacin treatment. Amplification of norA provides a further mechanism of rapid evolution in isolates from the CC398 lineage. Crucially, chemical inhibition of NorA effectively prevents the evolution of resistance in all isolates. Our study shows that pre-existing genetic diversity plays a key role in shaping resistance evolution, and it may be possible to predict which strains are likely to evolve resistance and to optimize inhibitor use to prevent this outcome.

The Bacillus subtilis dCMP deaminase ComEB acts as a dynamic polar localization factor for ComGA within the competence machinery

Bacillus subtilis can import DNA from the environment by an uptake machinery that localizes to a single cell pole. We investigated the roles of ComEB and of the ATPase ComGA during the state of competence. We show that ComEB plays an important role during competence, possibly because it is necessary for the recruitment of GomGA to the cell pole. ComEB localizes to the cell poles even upon expression during exponential phase, indicating that it can serve as polar marker. ComEB is also a deoxycytidylate monophosphate (dCMP) deaminase, for the function of which a conserved cysteine residue is important. However, cysteine-mutant ComEB is still capable of natural transformation, while a comEB deletion strain is highly impaired in competence, indicating that ComEB confers two independent functions. Single-molecule tracking (SMT) reveals that both proteins exchange at the cell poles between bound and unbound in a time scale of a few milliseconds, but turnover of ComGA increases during DNA uptake, whereas the mobility of ComEB is not affected. Our data reveal a highly dynamic role of ComGA during DNA uptake and an unusual role for ComEB as a mediator of polar localization, localizing by diffusion-capture on an extremely rapid time scale and functioning as a moonlighting enzyme.

Sequential Monte Carlo with transformations

This paper examines methodology for performing Bayesian inference sequentially on a sequence of posteriors on spaces of different dimensions. For this, we use sequential Monte Carlo samplers, introducing the innovation of using deterministic transformations to move particles effectively between target distributions with different dimensions. This approach, combined with adaptive methods, yields an extremely flexible and general algorithm for Bayesian model comparison that is suitable for use in applications where the acceptance rate in reversible jump Markov chain Monte Carlo is low. We use this approach on model comparison for mixture models, and for inferring coalescent trees sequentially, as data arrives.

Comprehensive Analysis Reveals the Evolution and Pathogenicity of Aeromonas, Viewed from Both Single Isolated Species and Microbial Communities

The genus Aeromonas is a common gastrointestinal pathogen associated with human and animal infections. Due to the high level of cross-species similarity, their evolutionary dynamics and genetic diversity are still fragmented. Hereby, we investigated the pan-genomes of 29 Aeromonas species, as well as Aeromonas species in microbial communities, to clarify their evolutionary dynamics and genetic diversity, with special focus on virulence factors and horizontal gene transfer events. Our study revealed an open pan-genome of Aeromonas containing 10,144 gene families. These Aeromonas species exhibited different functional constraints, with the single-copy core genes and most accessory genes experiencing purifying selection. The significant congruence between core genome and pan-genome trees revealed that core genes mainly affected evolutionary divergences of Aeromonas species. Gene gains and losses revealed a high level of genome plasticity, exhibited by hundreds of gene expansions and contractions, horizontally transferred genes, and mobile genetic elements. The selective constraints shaped virulence gene pools of these Aeromonas strains, where genes encoding hemolysin were ubiquitous. Of these strains, Aeromonas aquatica MX16A seemed to be more resistant, as it harbored most resistance genes. Finally, the virulence factors of Aeromonas in microbial communities were quite dynamic in response to environment changes. For example, the virulence diversity of Aeromonas in microbial communities could reach levels that match some of the most virulent Aeromonas species (such as A. hydrophila) in penetrated-air and modified-air packaging. Our work shed some light onto genetic diversity, evolutionary history, and functional features of Aeromonas, which could facilitate the detection and prevention of infections.IMPORTANCE Aeromonas has long been known as a gastrointestinal pathogen, yet it has many species whose evolutionary dynamics and genetic diversity had been unclear until now. We have conducted pan-genome analysis for 29 Aeromonas species and revealed a high level of genome plasticity exhibited by hundreds of gene expansions and contractions, horizontally transferred genes, and mobile genetic elements. These species also contained many virulence factors both identified from single isolated species and microbial community. This pan-genome study could elevate the level for detection and prevention of Aeromonas infections.

The recombination-cold region as an epidemiological marker of recombinogenic opportunistic pathogen Mycobacterium avium

Background

The rapid identification of lineage remains a challenge in the genotyping of clinical isolates of recombinogenic pathogens. The chromosome of Mycobacterium avium subsp. hominissuis (MAH), an agent of Mycobacterium avium complex (MAC) lung disease, is often mosaic and is composed of chromosomal segments originating from different lineages. This makes it difficult to infer the MAH lineage in a simple experimental set-up. To overcome this difficulty, we sought to identify chromosomal marker genes containing lineage-specific alleles by genome data mining.

Results

We conducted genetic population structure analysis, phylogenetic analysis, and a survey of historical recombination using data from 125 global MAH isolates. Six MAH lineages (EA1, EA2, SC1, SC2, SC3, and SC4) were identified in the current dataset. One P-450 gene (locus_tag MAH_0788/MAV_0940) in the recombination-cold region was found to have multiple alleles that could discriminate five lineages. By combining the information about allele type from one additional gene, the six MAH lineages as well as other M. avium subspecies were distinguishable. A recombination-cold region of 116 kb contains an insertion hotspot and is flanked by a mammalian cell-entry protein operon where allelic variants have previously been reported to occur. Hence, we speculate that the acquisition of lineage- or strain-specific insertions has introduced homology breaks in the chromosome, thereby reducing the chance of interlineage recombination.

Conclusions

The allele types of the newly identified marker genes can be used to predict major lineages of M. avium. The single nucleotide polymorphism typing approach targeting multiallelic loci in recombination-cold regions will facilitate the epidemiological study of MAC, and may also be useful for equivalent studies of other nontuberculous mycobacteria potentially carrying mosaic genomes.

panX: pan-genome analysis and exploration

Horizontal transfer, gene loss, and duplication result in dynamic bacterial genomes shaped by a complex mixture of different modes of evolution. Closely related strains can differ in the presence or absence of many genes, and the total number of distinct genes found in a set of related isolates—the pan-genome—is often many times larger than the genome of individual isolates. We have developed a pipeline that efficiently identifies orthologous gene clusters in the pan-genome. This pipeline is coupled to a powerful yet easy-to-use web-based visualization for interactive exploration of the pan-genome. The visualization consists of connected components that allow rapid filtering and searching of genes and inspection of their evolutionary history. For each gene cluster, panX displays an alignment, a phylogenetic tree, maps mutations within that cluster to the branches of the tree and infers gain and loss of genes on the core-genome phylogeny. PanX is available at pangenome.de. Custom pan-genomes can be visualized either using a web server or by serving panX locally as a browser-based application.

Long-Term Intrahost Evolution of Methicillin Resistant Staphylococcus aureus Among Cystic Fibrosis Patients With Respiratory Carriage

Staphylococcus aureus is the most commonly identified airway colonizer of cystic fibrosis (CF) patients, and infections with methicillin-resistant S. aureus (MRSA) are associated with poor outcomes. Yet, little is known about the intrahost evolution of S. aureus among CF patients. We investigated convergent evolution and adaptation of MRSA among four CF patients with long-term respiratory carriage. For each patient, we performed whole-genome sequencing on an average of 21 isolates (range: 19–23) carried for a mean of 1,403 days (range: 903–1,679), including 25 pairs of isolates collected on the same day. We assessed intrahost diversity, population structure, evolutionary history, evidence of switched intergenic regions (IGRs), and signatures of adaptation in the context of patient age, antibiotic treatment, and co-colonizing microbes. Phylogenetic analysis delineated distinct multilocus sequence type ST5 (n = 3) and ST72 (n = 1) clonal populations in addition to sporadic, non-clonal isolates, and uncovered a putative transmission event. Variation in antibiotic resistance was observed within clonal populations, even among isolates collected on the same day. Rates of molecular evolution ranged from 2.21 to 8.64 nucleotide polymorphisms per year, and lineage ages were consistent with acquisition of colonization in early childhood followed by subsequent persistence of multiple sub-populations. Selection analysis of 1,622 core genes present in all four clonal populations (n = 79) found 11 genes variable in three subjects – most notably, ATP-dependent protease clpX, 2-oxoglutarate dehydrogenase odhA, fmtC, and transcription-repair coupling factor mfd. Only one gene, staphylococcal protein A (spa), was found to have evidence of gene-wide diversifying selection. We identified three instances of intrahost IGR switching events, two of which flanked genes related to quorum sensing. The complex microbial ecology of the CF airway poses challenges for management. We illustrate appreciable intrahost diversity as well as persistence of a dominant lineage. We also show that intrahost adaptation is a continual process, despite purifying selective pressure, and provide targets that should be investigated further for their function in CF adaptation.

Genome hypermobility by lateral transduction

Genetic transduction is a major evolutionary force that underlies bacterial adaptation. Here we report that the temperate bacteriophages of Staphylococcus aureus engage in a distinct form of transduction we term lateral transduction. Staphylococcal prophages do not follow the previously described excision-replication-packaging pathway but instead excise late in their lytic program. Here, DNA packaging initiates in situ from integrated prophages, and large metameric spans including several hundred kilobases of the S. aureus genome are packaged in phage heads at very high frequency. In situ replication before DNA packaging creates multiple prophage genomes so that lateral-transducing particles form during normal phage maturation, transforming parts of the S. aureus chromosome into hypermobile regions of gene transfer.

Vaccination with the Staphylococcus aureus secreted proteins EapH1 and EapH2 impacts both S. aureus carriage and invasive disease

INTRODUCTION

There is a need for an efficacious vaccine reducing infections due to Staphylococcus aureus, a common cause of community and hospital infection. Infecting organisms originate from S. aureus populations colonising the nares and bowel. Antimicrobials are widely used to transiently reduce S. aureus colonisation prior to surgery, a practice which is selecting for resistant S. aureus isolates. S. aureus secretes multiple proteins, including the protease inhibitors extracellular adhesion protein homologue 1 and 2 (EapH1 and EapH2).

METHODS

Mice were vaccinated intramuscularly or intranasally with Adenovirus serotype 5 and Modified Vaccinia Ankara viral vectors expressing EapH1 and EapH2 proteins, or with control viruses. Using murine S. aureus colonisation models, we monitored S. aureus colonisation by sequential stool sampling. Monitoring of S. aureus invasive disease after intravenous challenge was performed using bacterial load and abscess numbers in the kidney.

RESULTS

Intramuscular vaccination with Adenovirus serotype 5 and Modified Vaccinia Ankara viral vectors expressing EapH1 and EapH2 proteins significantly reduces bacterial recovery in the murine renal abscess model of infection, but the magnitude of the effect is small. A single intranasal vaccination with an adenoviral vaccine expressing these proteins reduced S. aureus gastrointestinal (GI) tract colonisation.

CONCLUSION

Vaccination against EapH1 / EapH2 proteins may offer an antibiotic independent way to reduce S. aureus colonisation, as well as contributing to protection against S. aureus invasive disease.

A statistical method to identify recombination in bacterial genomes based on SNP incompatibility

BACKGROUND

Phylogeny estimation for bacteria is likely to reflect their true evolutionary histories only if they are highly clonal. However, recombination events could occur during evolution for some species. The reconstruction of phylogenetic trees from an alignment without considering recombination could be misleading, since the relationships among strains in some parts of the genome might be different than in others. Using a single, global tree can create the appearance of homoplasy in recombined regions. Hence, the identification of recombination breakpoints is essential to better understand the evolutionary relationships of isolates among a bacterial population.

RESULTS

Previously, we have developed a method (called ACR) to detect potential breakpoints in an alignment by evaluating compatibility of polymorphic sites in a sliding window. To assess the statistical significance of candidate breakpoints, we propose an extension of the algorithm (ptACR) that applies a permutation test to generate a null distribution for comparing the average local compatibility. The performance of ptACR is evaluated on both simulated and empirical datasets. ptACR is shown to have similar sensitivity (true positive rate) but a lower false positive rate and higher F1 score compared to basic ACR. When used to analyze a collection of clinical isolates of Staphylococcus aureus, ptACR finds clear evidence of recombination events in this bacterial pathogen, and is able to identify statistically significant boundaries of chromosomal regions with distinct phylogenies.

CONCLUSIONS

ptACR is an accurate and efficient method for identifying genomic regions affected by recombination in bacterial genomes.

A genetic linkage map for the salmon louse (Lepeophtheirus salmonis): evidence for high male:female and inter-familial recombination rate differences

A salmon louse (Lepeophtheirus salmonis salmonis) genetic linkage map was constructed to serve as a genomic resource for future investigations into the biology of this important marine parasitic copepod species, and to provide insights into the inheritance patterns of genetic markers in this species. SNP genotyping of 8 families confirmed the presence of 15 linkage groups based upon the assignment of 93,773 markers. Progeny sample size weight adjusted map sizes in males (with the exception of SL12 and SL15) ranged in size from 96.50 cM (SL11) to 134.61 cM (SL06), and total combined map steps or bins ranged from 143 (SL09) to 203 (SL13). The SL12 male map was the smallest linkage group with a weight-averaged size of 3.05 cM with 6 recombination bins. Male:female specific recombination rate differences are 10.49:1 and represent one of the largest reported sex-specific differences for any animal species. Recombination ratio differences (M:F) ranged from 1.0 (SL12) to 29:1 (SL15). The number of markers exhibiting normal Mendelian segregation within the sex linkage group SL15 was extremely low (N = 80) in comparison to other linkage groups genotyped [range: 1459 (SL12)-10206 markers (SL05)]. Re-evaluation of Mendelian inheritance patterns of markers unassigned to any mapping parent according to hemizygous segregation patterns (models presented) identified matches for many of these markers to hemizygous patterns. The greatest proportion of these markers assigned to SL15 (N increased to 574). Inclusion of the hemizygous markers revised SL15 sex-specific recombination rate differences to 28:1. Recombination hot- and coldspots were identified across all linkage groups with all linkage groups possessing multiple peaks. Nine of 13 linkage groups evaluated possessed adjacent domains with hot-coldspot transitional zones. The most common pattern was for one end of the linkage to show elevated recombination in addition to internal regions. For SL01 and SL06, however, a terminal region with high recombination was not evident while a central domain possessing extremely high-recombination levels was present. High levels of recombination were weakly coupled to higher levels of SNP variation within domains, but this association was very strong for the central domains of SL01 and SL06. From the pooled paternal half-sib lots (several virgin females placed with 1 male), only 1 or two surviving family lots were obtained. Surviving families possessed parents where both the male and female possessed either inherently low or high recombination rates. This study provides insight into the organization of the sea louse genome, and describes large differences in recombination rate that exist among individuals of the same sex, and between the sexes. These differences in recombination rate may be coupled to the capabilities of this species to adapt to environmental and pharmaceutical treatments, given that family survivorship appears to be enhanced when parents have similar recombination levels.

The Enterococcus Cassette Chromosome, a Genomic Variation Enabler in Enterococci

Enterococcus faecium has a highly variable genome prone to recombination and horizontal gene transfer. Here, we have identified a novel genetic island with an insertion locus and mobilization genes similar to those of staphylococcus cassette chromosome elements SCCmec This novel element termed the enterococcus cassette chromosome (ECC) element was located in the 3' region of rlmH and encoded large serine recombinases ccrAB similar to SCCmec Horizontal transfer of an ECC element termed ECC::cat containing a knock-in cat chloramphenicol resistance determinant occurred in the presence of a conjugative rep _pLG1 plasmid. We determined the ECC::cat insertion site in the 3' region of rlmH in the E. faecium recipient by long-read sequencing. ECC::cat also mobilized by homologous recombination through sequence identity between flanking insertion sequence (IS) elements in ECC::cat and the conjugative plasmid. The ccrAB_Ent genes were found in 69 of 516 E. faecium genomes in GenBank. Full-length ECC elements were retrieved from 32 of these genomes. ECCs were flanked by attR and attL sites of approximately 50 bp. The attECC sequences were found by PCR and sequencing of circularized ECCs in three strains. The genes in ECCs contained an amalgam of common and rare E. faecium genes. Taken together, our data imply that ECC elements act as hot spots for genetic exchange and contribute to the large variation of accessory genes found in E. faecium IMPORTANCE Enterococcus faecium is a bacterium found in a great variety of environments, ranging from the clinic as a nosocomial pathogen to natural habitats such as mammalian intestines, water, and soil. They are known to exchange genetic material through horizontal gene transfer and recombination, leading to great variability of accessory genes and aiding environmental adaptation. Identifying mobile genetic elements causing sequence variation is important to understand how genetic content variation occurs. Here, a novel genetic island, the enterococcus cassette chromosome, is shown to contain a wealth of genes, which may aid E. faecium in adapting to new environments. The transmission mechanism involves the only two conserved genes within ECC, ccrAB _Ent, large serine recombinases that insert ECC into the host genome similarly to SCC elements found in staphylococci.

Equine Methicillin-Resistant Sequence Type 398 Staphylococcus aureus (MRSA) Harbor Mobile Genetic Elements Promoting Host Adaptation

Continuing introduction of multi-drug resistant, zoonotic pathogens such as methicillin-resistant Staphylococcus aureus (MRSA) in horse clinics challenges the biosafety of employees and animal patients. This study was aimed to determine the occurrence of mobile genetic elements facilitating survival in the early stages of invasive infection in different host species, including humans and horses, in MRSA carried by equine patients admitted to a large horse clinic. A total of 341 equine patients were investigated for carriage of MRSA by hygiene screening directly at hospital admission. MRSA were further investigated by antimicrobial susceptibility testing, whole-genome sequencing and genomic composition, including virulence factors involved in immune evasion and host adaption. From a total of 340 validated specimens from equine nostrils, 3.5% yielded positive results for MRSA. All MRSA were found to be closely related belonging to sequence type (ST) 398_t011 with up to four additional antimicrobial resistances. All MRSA harbored a specific Staphylococcal Pathogenicity Island (SaPIbov5) involved in facilitating survival in ruminant and equine plasma. Moreover, a β-hemolysin (hlb) converting ΦSa3 phage encoding the human-specific Immune Evasion Cluster (IEC) was present in 72% of the isolates. An equid-specific leukotoxin encoded by a further temperate phage (Saeq1) was only rarely detected (22%). Despite the absence of β-hemolysin production for all IEC-positive ST398, a prominent hemolysis zone was demonstrable on sheep blood agar. Thus, IEC might remain undetected among the ST398 lineage, since the presence of IEC is commonly associated with reduction of hemolysis in S. aureus belonging to other genetic backgrounds. Here we describe MRSA-ST398 harboring different mobile genetic elements encoding variants of immune evasion factors and toxins previously shown to contribute to S. aureus invasive diseases in specific host species or ecologic niches. We suggest these combinations contribute to the adaptation of MRSA belonging to ST398 with respect to epidemic spread across different habitats and hosts, and may therefore confer a host “generalist” phenotype.