Since its discovery in the early 2000s, methicillin-resistant Staphylococcus aureus (MRSA) clonal complex 398 (CC398) has become a rapidly emerging cause of human infections, most often associated with livestock exposure. We applied whole-genome sequence typing to characterize a diverse collection of CC398 isolates ( n = 89), including MRSA and methicillin-susceptible S. aureus (MSSA) from animals and humans spanning 19 countries and four continents. We identified 4,238 single nucleotide polymorphisms (SNPs) among the 89 core genomes. Minimal homoplasy (consistency index = 0.9591) was detected among parsimony-informative SNPs, allowing for the generation of a highly accurate phylogenetic reconstruction of the CC398 clonal lineage. Phylogenetic analyses revealed that MSSA from humans formed the most ancestral clades. The most derived lineages were composed predominantly of livestock-associated MRSA possessing three different staphylococcal cassette chromosome mec element (SCC mec ) types (IV, V, and VII-like) including nine subtypes. The human-associated isolates from the basal clades carried phages encoding human innate immune modulators that were largely missing among the livestock-associated isolates. Our results strongly suggest that livestock-associated MRSA CC398 originated in humans as MSSA. The lineage appears to have undergone a rapid radiation in conjunction with the jump from humans to livestock, where it subsequently acquired tetracycline and methicillin resistance. Further analyses are required to estimate the number of independent genetic events leading to the methicillin-resistant sublineages, but the diversity of SCC mec subtypes is suggestive of strong and diverse antimicrobial selection associated with food animal production. IMPORTANCE Modern food animal production is characterized by densely concentrated animals and routine antibiotic use, which may facilitate the emergence of novel antibiotic-resistant zoonotic pathogens. Our findings strongly support the idea that livestock-associated MRSA CC398 originated as MSSA in humans. The jump of CC398 from humans to livestock was accompanied by the loss of phage-carried human virulence genes, which likely attenuated its zoonotic potential, but it was also accompanied by the acquisition of tetracycline and methicillin resistance. Our findings exemplify a bidirectional zoonotic exchange and underscore the potential public health risks of widespread antibiotic use in food animal production.

The Escherichia coli sequence type 131 (ST131) clone is notorious for extraintestinal infections, fluoroquinolone resistance, and extended-spectrum beta-lactamase (ESBL) production, attributable to a CTX-M-15-encoding mobile element. Here, we applied pulsed-field gel electrophoresis (PFGE) and whole-genome sequencing to reconstruct the evolutionary history of the ST131 clone. PFGE-based cluster analyses suggested that both fluoroquinolone resistance and ESBL production had been acquired by multiple ST131 sublineages through independent genetic events. In contrast, the more robust whole-genome-sequence-based phylogenomic analysis revealed that fluoroquinolone resistance was confined almost entirely to a single, rapidly expanding ST131 subclone, designated H 30-R. Strikingly, 91% of the CTX-M-15-producing isolates also belonged to a single, well-defined clade nested within H 30-R, which was named H 30-Rx due to its more extensive resistance. Despite its tight clonal relationship with H 30Rx, the CTX-M-15 mobile element was inserted variably in plasmid and chromosomal locations within the H 30-Rx genome. Screening of a large collection of recent clinical E. coli isolates both confirmed the global clonal expansion of H 30-Rx and revealed its disproportionate association with sepsis (relative risk, 7.5; P < 0.001). Together, these results suggest that the high prevalence of CTX-M-15 production among ST131 isolates is due primarily to the expansion of a single, highly virulent subclone, H 30-Rx. IMPORTANCE We applied an advanced genomic approach to study the recent evolutionary history of one of the most important Escherichia coli strains in circulation today. This strain, called sequence type 131 (ST131), causes multidrug-resistant bladder, kidney, and bloodstream infections around the world. The rising prevalence of antibiotic resistance in E. coli is making these infections more difficult to treat and is leading to increased mortality. Past studies suggested that many different ST131 strains gained resistance to extended-spectrum cephalosporins independently. In contrast, our research indicates that most extended-spectrum-cephalosporin-resistant ST131 strains belong to a single highly pathogenic subclone, called H 30-Rx. The clonal nature of H 30-Rx may provide opportunities for vaccine or transmission prevention-based control strategies, which could gain importance as H 30-Rx and other extraintestinal pathogenic E. coli subclones become resistant to our best antibiotics.

Anthrax, caused by the bacterium Bacillus anthracis, is a disease of historical and current importance that is found throughout the world. The basis of its historical transmission is anecdotal and its true global population structure has remained largely cryptic. Seven diverse B. anthracis strains were whole-genome sequenced to identify rare single nucleotide polymorphisms (SNPs), followed by phylogenetic reconstruction of these characters onto an evolutionary model. This analysis identified SNPs that define the major clonal lineages within the species. These SNPs, in concert with 15 variable number tandem repeat (VNTR) markers, were used to subtype a collection of 1,033 B. anthracis isolates from 42 countries to create an extensive genotype data set. These analyses subdivided the isolates into three previously recognized major lineages (A, B, and C), with further subdivision into 12 clonal sub-lineages or sub-groups and, finally, 221 unique MLVA15 genotypes. This rare genomic variation was used to document the evolutionary progression of B. anthracis and to establish global patterns of diversity. Isolates in the A lineage are widely dispersed globally, whereas the B and C lineages occur on more restricted spatial scales. Molecular clock models based upon genome-wide synonymous substitutions indicate there was a massive radiation of the A lineage that occurred in the mid-Holocene (3,064–6,127 ybp). On more recent temporal scales, the global population structure of B. anthracis reflects colonial-era importation of specific genotypes from the Old World into the New World, as well as the repeated industrial importation of diverse genotypes into developed countries via spore-contaminated animal products. These findings indicate humans have played an important role in the evolution of anthrax by increasing the proliferation and dispersal of this now global disease. Finally, the value of global genotypic analysis for investigating bioterrorist-mediated outbreaks of anthrax is demonstrated.