kSNP3.0: SNP detection and phylogenetic analysis of genomes without genome alignment or reference genome

Transient and resident pathogens: Intra-facility genetic diversity of Listeria monocytogenes and Salmonella from food production environments

Food production facilities are often routinely tested over time for the presence of foodborne pathogens (e.g., Listeria monocytogenes or Salmonella enterica subsp. enterica). Strains detected in a single sampling event can be classified as transient; positive findings of the same strain across multiple sampling events can be classified as resident pathogens. We analyzed whole-genome sequence (WGS) data from 4,758 isolates (L. monocytogenes = 3,685; Salmonella = 1,073) from environmental samples taken by FDA from 536 U.S. facilities. Our primary objective was to determine the frequency of transient or resident pathogens within food production facilities. Strains were defined as isolates from the same facility that are less than 50 SNP (single-nucleotide polymorphisms) different from one another. Resident pathogens were defined as strains that had more than one isolate collected >59 days apart and from the same facility. We found 1,076 strains (median = 1 and maximum = 21 strains per facility); 180 were resident pathogens, 659 were transient, and 237 came from facilities that had only been sampled once. As a result, 21% of strains (180/ 839) from facilities with positive findings and that were sampled multiple times were found to be resident pathogens; nearly 1 in 4 (23%) of L. monocytogenes strains were found to be resident pathogens compared to 1 in 6 (16%) of Salmonella strains. Our results emphasize the critical importance of preventing the colonization of food production environments by foodborne pathogens, since when colonization does occur, there is an appreciable chance it will become a resident pathogen that presents an ongoing potential to contaminate product.

Spread of Carbapenem-Resistant Klebsiella pneumoniae Clinical Isolates Producing NDM-Type Metallo-β-Lactamase in Myanmar

A total of 38 isolates of carbapenem-resistant Klebsiella pneumoniae harboring blaNDM were obtained during surveillance of 10 hospitals in Myanmar. Of these 38 isolates, 19 (50%) harbored genes encoding 16S rRNA methylases, such as armA or rmtB. The K. pneumoniae strains tested belonged to 17 sequence types (STs), including the high-risk clonal lineages ST101 and ST147. The ST101 and ST147 isolates carried IncFII plasmids harboring blaNDM-5 and IncFIB(pQil) plasmids harboring blaNDM-1, respectively. These results indicate that IncFII plasmids harboring blaNDM-5 and IncFIB(pQil) plasmids harboring blaNDM-1 have been spreading in K. pneumoniae ST101 and ST147 isolates, respectively, in Myanmar. IMPORTANCE The emergence of carbapenem-resistant K. pneumoniae has become a serious problem in medical settings worldwide. The present study demonstrated that carbapenem-resistant K. pneumoniae strains have been spreading in medical settings in Myanmar. In particular, plasmid genes encoding NDMs and 16S rRNA methylases have been spreading in K. pneumoniae high-risk clones.

Invasive Pneumococcal Disease Clusters Disproportionally Impact Persons Experiencing Homelessness, Injecting Drug Users, and the Western United States

BACKGROUND

Invasive pneumococcal disease (IPD) isolates forming genomic clusters can reflect rapid disease transmission between vulnerable individuals.

METHODS

We performed whole genome sequencing of 2820 IPD isolates recovered during 2019 through Centers for Disease Control and Prevention's Active Bacterial Core surveillance to provide strain information (serotypes, resistance, genotypes), and 2778 of these genomes were analyzed to detect highly related genomic clusters.

RESULTS

Isolates from persons experiencing homelessness (PEH) were more often within genomic clusters than those from persons not experiencing homelessness (PNEH) (105/198 [53.0%] vs 592/2551 [23.2%]; P < .001). The 4 western sites accounted for 33.4% (929/2778) of isolates subjected to cluster analysis yet accounted for 48.7% (343/705) of clustering isolates (P < .001) and 75.8% (150/198) of isolates recovered from PEH (P < .001). Serotypes most frequent among PEH were (in rank order) 12F, 4, 3, 9N, 8, 20, and 22F, all of which were among the 10 serotypes exhibiting the highest proportions of clustering isolates among all cases. These serotypes accounted for 44.9% (1265/2820) of all IPD cases and are included within available vaccines.

CONCLUSIONS

We identified serotype-specific and geographic differences in IPD transmission. We show the vulnerability of PEH within different regions to rapidly spreading IPD transmission networks representing several pneumococcal serotypes included in available vaccines.

De Novo Long-Read Whole-Genome Assemblies and the Comparative Pan-Genome Analysis of Ascochyta Blight Pathogens Affecting Field Pea

Ascochyta Blight (AB) is a major disease of many cool-season legumes globally. In field pea, three fungal pathogens have been identified to be responsible for this disease in Australia, namely Peyronellaea pinodes, Peyronellaea pinodella and Phoma koolunga. Limited genomic resources for these pathogens have been generated, which has hampered the implementation of effective management strategies and breeding for resistant cultivars. Using Oxford Nanopore long-read sequencing, we report the first high-quality, fully annotated, near-chromosome-level nuclear and mitochondrial genome assemblies for 18 isolates from the Australian AB complex. Comparative genome analysis was performed to elucidate the differences and similarities between species and isolates using phylogenetic relationships and functional diversity. Our data indicated that P. pinodella and P. koolunga are heterothallic, while P. pinodes is homothallic. More homology and orthologous gene clusters are shared between P. pinodes and P. pinodella compared to P. koolunga. The analysis of the repetitive DNA content showed differences in the transposable repeat composition in the genomes and their expression in the transcriptomes. Significant repeat expansion in P. koolunga's genome was seen, with strong repeat-induced point mutation (RIP) activity being evident. Phylogenetic analysis revealed that genetic diversity can be exploited for species marker development. This study provided the much-needed genetic resources and characterization of the AB species to further drive research in key areas such as disease epidemiology and host-pathogen interactions.

Characterization of Adherent-Invasive Escherichia coli (AIEC) Outer Membrane Proteins Provides Potential Molecular Markers to Screen Putative AIEC Strains

Adherent-invasive E. coli (AIEC) is a pathotype associated with the etiopathogenesis of Crohn's disease (CD), albeit with an as-yet unclear role. The main pathogenic mechanisms described for AIEC are adherence to epithelial cells, invasion of epithelial cells, and survival and replication within macrophages. A few virulence factors have been described as participating directly in these phenotypes, most of which have been evaluated only in AIEC reference strains. To date, no molecular markers have been identified that can differentiate AIEC from other E. coli pathotypes, so these strains are currently identified based on the phenotypic characterization of their pathogenic mechanisms. The identification of putative AIEC molecular markers could be beneficial not only from the diagnostic point of view but could also help in better understanding the determinants of AIEC pathogenicity. The objective of this study was to identify molecular markers that contribute to the screening of AIEC strains. For this, we characterized outer membrane protein (OMP) profiles in a group of AIEC strains and compared them with the commensal E. coli HS strain. Notably, we found a set of OMPs that were present in the AIEC strains but absent in the HS strain. Moreover, we developed a PCR assay and performed phylogenomic analyses to determine the frequency and distribution of the genes coding for these OMPs in a larger collection of AIEC and other E. coli strains. As result, it was found that three genes (chuA, eefC, and fitA) are widely distributed and significantly correlated with AIEC strains, whereas they are infrequent in commensal and diarrheagenic E. coli strains (DEC). Additional studies are needed to validate these markers in diverse strain collections from different geographical regions, as well as investigate their possible role in AIEC pathogenicity.

Anthroponotic-Based Transfer of Staphylococcus to Dog: A Case Study

Although usually harmless, Staphylococcus spp. can cause nosocomial and community-onset skin and soft tissue infections in both humans and animals; thus, it is considered a significant burden for healthcare systems worldwide. Companion animals have been identified as potential reservoirs of pathogenic Staphylococcus with specific reference to Methicillin Resistant Staphylococcus aureus (MRSA). In this study, we investigated the circulation and the genetic relationships of a collection of Staphylococcus spp. isolates in a family composed of four adults (a mother, father, grandmother, and grandfather), one child, and a dog, which were sampled over three years. The routes of transmission among humans and between humans and the dog werelyzed. The results displayed the circulation of many Staphylococcus lineages, belonging to different species and sequence types (ST) and being related to both human and pet origins. However, among the observed host-switch events, one of them clearly underpinnthroponotic route from a human to a dog. This suggests that companion animals can potentially have a role as a carrier of Staphylococcus, thus posing a serious concern about MRSA spreading within human and animal microbial communities.

Analysis of the Taxonomy, Synteny, and Virulence Factors for Soft Rot Pathogen Pectobacterium aroidearum in Amorphophallus konjac Using Comparative Genomics

Bacterial soft rot is a devastating disease for a wide range of crops, vegetables, and ornamental plants including konjac (Amorphophallus konjac). However, the pangenome and genomic plasticity of the konjac soft rot pathogens is little explored. In this study, we reported the complete genome sequences of 11 bacterial isolates that can cause typical soft rot symptoms in konjac by in vitro and in vivo pathogenicity tests. Based on in silico DNA-DNA hybridization, average nucleotide identity and phylogenomic analysis, all 11 isolates were determined to be Pectobacterium aroidearum. In addition, synteny analysis of these genomes revealed considerable chromosomal inversions, one of which is triggered by homologous recombination of ribose operon. Pangenome analysis and COG enrichment analysis showed that the pangenome of P. aroidearum is open and that accessory genes are enriched in replication, recombination, and repair. Variations in type IV secretion system and type VI secretion system were found, while plant cell wall degrading enzymes were conserved. Furthermore, sequence analyses also provided evidence for the presence of a type V secretion system in Pectobacterium. These findings advance our understanding of the pathogenicity determinants, genomic plasticity, and evolution of P. aroidearum.

Genomic Degeneration and Reduction in the Fish Pathogen Mycobacterium shottsii

Mycobacterium shottsii is a dysgonic, nonpigmented mycobacterium originally isolated from diseased striped bass (Morone saxatilis) in the Chesapeake Bay, USA. Genomic analysis reveals that M. shottsii is a Mycobacterium ulcerans/Mycobacterium marinum clade (MuMC) member, but unlike the superficially similar M. pseudoshottsii, also isolated from striped bass, it is not an M. ulcerans ecovar, instead belonging to a transitional group of strains basal to proposed "Aronson" and "M" lineages. Although phylogenetically distinct from the human pathogen M. ulcerans, the M. shottsii genome shows parallel but nonhomologous genomic degeneration, including massive accumulation of pseudogenes accompanied by proliferation of unique insertion sequences (ISMysh01, ISMysh03), large-scale deletions, and genomic reorganization relative to typical M. marinum strains. Coupled with its observed ecological characteristics and loss of chromogenicity, the genomic structure of M. shottsii is suggestive of evolution toward a state of obligate pathogenicity, as observed for other Mycobacterium spp., including M. ulcerans, M. tuberculosis, and M. leprae. IMPORTANCE Morone saxatilis (striped bass) is an ecologically and economically important finfish species on the United States east coast. Mycobacterium shottsii and Mycobacterium pseudoshottsii were originally described in the early 2000s as novel species from outbreaks of visceral and dermal mycobacteriosis in this species. Biochemical and genetic characterization place these species within the Mycobacterium ulcerans/M. marinum clade (MuMC), and M. pseudoshottsii has been proposed as an ecovar of M. ulcerans. Here, we describe the complete genome of M. shottsii, demonstrating that it is clearly not an M. ulcerans ecovar; however, it has undergone parallel genomic modification suggestive of a transition to obligate pathogenicity. As in M. ulcerans, the M. shottsii genome demonstrates widespread pseudogene formation driven by proliferation of insertion sequences, as well as genomic reorganization. This work clarifies the phylogenetic position of M. shottsii relative to other MuMC members and provides insight into processes shaping its genomic structure.

The Number and Type of Chaperone-Usher Fimbriae Reflect Phylogenetic Clade Rather than Host Range in Salmonella

Salmonella is one of the most successful foodborne pathogens worldwide, owing in part to its ability to colonize or infect a wide range of hosts. Salmonella serovars are known to encode a variety of different fimbriae (hairlike organelles that facilitate binding to surfaces); however, the distribution, number, and sequence diversity of fimbriae encoded across different lineages of Salmonella were unknown. We queried whole-genome sequence (WGS) data for 242 Salmonella enterica subsp. enterica (subspecies enterica) isolates from the top 217 serovars associated with isolation from humans and agricultural animals; this effort identified 2,894 chaperone-usher (CU)-type fimbrial usher sequences, representing the most conserved component of CU fimbriae. On average, isolates encoded 12 different CU fimbrial ushers (6 to 18 per genome), although the distribution varied significantly (P = 1.328E-08) by phylogenetic clade, with isolates in section Typhi having significantly fewer fimbrial ushers than isolates in clade A2 (medians = 10 and 12 ushers, respectively). Characterization of fimbriae in additional non-enterica subspecies genomes suggested that 8 fimbrial ushers were classified as being unique to subspecies enterica isolates, suggesting that the majority of fimbriae were most likely acquired prior to the divergence of subspecies enterica. Characterization of mobile elements suggested that plasmids represent an important vehicle facilitating the acquisition of a wide range of fimbrial ushers, particularly for the acquisition of fimbriae from other Gram-negative genera. Overall, our results suggest that differences in the number and type of fimbriae encoded most likely reflect differences in phylogenetic clade rather than differences in host range. IMPORTANCE Fimbriae of the CU assembly pathway represent important organelles that mediate Salmonella's interactions with host tissues and abiotic surfaces. Our analyses provide a comprehensive overview of the diversity of CU fimbriae in Salmonella spp., highlighting that the majority of CU fimbriae are distributed broadly across multiple subspecies and suggesting that acquisition most likely occurred prior to the divergence of subspecies enterica. Our data also suggest that plasmids represent the primary vehicles facilitating the horizontal transfer of diverse CU fimbriae in Salmonella. Finally, the observed high sequence similarity between some ushers suggests that different names may have been assigned to closely related fimbrial ushers that likely should be represented by a single designation. This highlights the need to establish standard criteria for fimbria classification and nomenclature, which will also facilitate future studies seeking to associate virulence factors with adaptation to or differences in the likelihood of causing disease in a given host.

Whole genome sequencing analyses revealed that Salmonella enterica serovar Dublin strains from Brazil belonged to two predominant clades

Salmonella Dublin is a cattle-associated serovar sporadically causing disease in humans. S. Dublin strains isolated in Brazil and in other countries were analyzed to determine their phylogenetic relationships, the presence of genes, plasmids, genomic regions related to virulence and antimicrobial resistance genes repertoire, using WGS analyses. Illumina was used to sequence the genome of 112 S. Dublin strains isolated in Brazil from humans (n = 82) and animals (n = 30) between 1983 and 2016. Furthermore, 87 strains from other countries were analyzed. WGSNP analysis revealed three different clades, in which the strains from Brazil belonged to two clades, A and C. Most of the genes and genomic regions searched varied among the strains studied. The siderophore genes iroB and iroC were exclusively found in strains from Brazil and pegD gene, related to fimbrial adherence determinants, were positive in 124 strains from clades A and B but absent in all the strains from clade C (n = 71). Eleven plasmid replicons were found in the strains from Brazil, and nine were exclusively found in strains from other countries. The antimicrobial resistance genes mdsA and mdsB, that encode an efflux pump, were found in all the strains studied. The strains from Brazil carried other resistance genes, such as tet(A) (n = 11), tet(B) (n = 4) and tet(C) (n = 4), blaTEM-1 (n = 4), catA1 (n = 1), aadA1 (n = 1), and sul1 (n = 1). In conclusion, S. Dublin strains isolated in Brazil presented some few unique genes not found in strains from other countries and were allocated into two distinct clades with strains of human and animal origin epidemiologically related. This fact stresses the zoonotic potential of S. Dublin circulating in Brazil for more than 30 years.

Comparative Genomics Provides Insights Into Genetic Diversity of Clostridium tyrobutyricum and Potential Implications for Late Blowing Defects in Cheese

Clostridium tyrobutyricum has been recognized as the main cause of late blowing defects (LBD) in cheese leading to considerable economic losses for the dairy industry. Although differences in spoilage ability among strains of this species have been acknowledged, potential links to the genetic diversity and functional traits remain unknown. In the present study, we aimed to investigate and characterize genomic variation, pan-genomic diversity and key traits of C. tyrobutyricum by comparing the genomes of 28 strains. A comparative genomics analysis revealed an “open” pangenome comprising 9,748 genes and a core genome of 1,179 genes shared by all test strains. Among those core genes, the majority of genes encode proteins related to translation, ribosomal structure and biogenesis, energy production and conversion, and amino acid metabolism. A large part of the accessory genome is composed of sets of unique, strain-specific genes ranging from about 5 to more than 980 genes. Furthermore, functional analysis revealed several strain-specific genes related to replication, recombination and repair, cell wall, membrane and envelope biogenesis, and defense mechanisms that might facilitate survival under stressful environmental conditions. Phylogenomic analysis divided strains into two clades: clade I contained human, mud, and silage isolates, whereas clade II comprised cheese and milk isolates. Notably, these two groups of isolates showed differences in certain hypothetical proteins, transcriptional regulators and ABC transporters involved in resistance to oxidative stress. To the best of our knowledge, this is the first study to provide comparative genomics of C. tyrobutyricum strains related to LBD. Importantly, the findings presented in this study highlight the broad genetic diversity of C. tyrobutyricum, which might help us understand the diversity in spoilage potential of C. tyrobutyricum in cheese and provide some clues for further exploring the gene modules responsible for the spoilage ability of this species.

Genomic Epidemiology of MBL-Producing Pseudomonas putida Group Isolates in Poland

Introduction

Pseudomonas putida group are described as low-incidence opportunistic pathogens, but also as a significant reservoir of antimicrobial resistance (AMR) genes, including those of metallo-β-lactamases (MBLs). Our objective was the molecular and genomic characterization of MBL-producing P. putida (MPPP) group isolates from Poland, focusing on population structures, successful genotypes and MBL-encoding integrons.

Methods

During a country-wide MBL surveillance in Pseudomonas spp., 59 non-duplicate MPPP isolates were collected from 36 hospitals in 23 towns from 2003 to 2016. All of the isolates were subjected to whole-genome sequencing (WGS), followed by species identification, multi-locus sequence typing (MLST), single-nucleotide polymorphism (SNP)-based phylogenetic/clonality analysis, resistome determination, and susceptibility testing.

Results

The study collection comprised 12 species, of which P. alloputida (n = 19), P. monteilii (n = 15), and P. asiatica (n = 11) prevailed, while the others were P. kurunegalensis, P. putida, P. soli, P. mosselii, P. juntendi, and four potentially new species. MLST classified the isolates into 23 sequence types (STs) of which 21 were new, with three main clones, namely P. alloputida ST69, P.monteilii ST95 and P. asiatica ST15. The isolates produced VIM-like MBLs only, largely VIM-2 (n = 40), encoded by 24 different class 1 integrons (ten new), a number of which occurred also in P. aeruginosa and/or Enterobacterales in Poland. The plasmid pool was dominated by IncP-9, IncP-2, and pMOS94-like types. Multiple isolates were extensively drug-resistant.

Conclusions

This study, being one of the most comprehensive analyses of MPPP so far, has shown high diversity of the isolates in general, with three apparently international lineages, each internally diversified by MBL-encoding structures.

Supplementary Information

The online version contains supplementary material available at 10.1007/s40121-022-00659-z.

Profiling Fusobacterium infection at high taxonomic resolution reveals lineage-specific correlations in colorectal cancer

The bacterial genus Fusobacterium promotes colorectal cancer (CRC) development, but an understanding of its precise composition at the species level in the human gut and the relevant association with CRC is lacking. Herein, we devise a Fusobacterium rpoB amplicon sequencing (FrpoB-seq) method that enables the differentiation of Fusobacterium species and certain subspecies in the microbiota. By applying this method to clinical tissue and faecal samples from CRC patients, we detect 62 Fusobacterium species, including 45 that were previously undescribed. We additionally reveal that Fusobacterium species may display different lineage-dependent functions in CRC. Specifically, a lineage (designated L1) including F. nucleatum, F. hwasookii, F. periodonticum and their relatives (rather than any particular species alone) is overabundant in tumour samples and faeces from CRC patients, whereas a non-enriched lineage (designated L5) represented by F. varium and F. ulcerans in tumours has a positive association with lymphovascular invasion.

Prevalence and characteristics of the mcr-1 gene in retail meat samples in Zhejiang Province, China

Considering the serious threat to food safety and public health posed by pathogens with colistin resistance, colistin was banned as a growth promoter in 2017 in China. In recent years, the resistance rate of Escherichia coli isolated from animal intestines or feces to colistin has decreased. However, the prevalence and characteristics of the mcr-1 gene in retail meat have not been well explored. Herein, 106 mcr-1-negative and 16 mcr-1-positive E. coli isolates were randomly recovered from 120 retail meat samples and screened using colistin. The 106 E. coli isolates showed maximum resistance to sulfafurazole (73.58%) and tetracycline (62.26%) but susceptibility to colistin (0.00%). All 16 mcr-1-positive E. coli isolates showed resistance to colistin, were extended spectrum beta-lactamase (ESBL)-positive and exhibited complex multidrug resistance (MDR). For these 16 isolates, 17 plasmid replicons and 42 antibiotic resistance genes were identified, and at least 7 antibiotic resistance genes were found in each isolate. Acquired disinfectant resistance genes were identified in 75.00% (12/16) of the isolates. Furthermore, comparative genomic and phylogenetic analysis results indicated that these 16 mcr-1-positive E. coli isolates and the most prevalent mcr-1-harboring IncI2 plasmid in this study were closely related to other previously reported mcr-1-positive E. coli isolates and the IncI2 plasmid, respectively, showing their wide distribution. Taken together, our findings showed that retail meat products were a crucial reservoir of mcr-1 during the colistin ban period and should be continuously monitored.

The c-di-GMP Phosphodiesterase PipA (PA0285) Regulates Autoaggregation and Pf4 Bacteriophage Production in Pseudomonas aeruginosa PAO1

In Pseudomonas aeruginosa PAO1, 41 genes encode proteins predicted to be involved in the production or degradation of c-di-GMP, a ubiquitous secondary messenger that regulates a variety of physiological behaviors closely related to biofilm and aggregate formation. Despite extensive effort, the entire picture of this important signaling network is still unclear, with one-third of these proteins remaining uncharacterized. Here, we show that the deletion of pipA, which produces a protein containing two PAS domains upstream of a GGDEF-EAL tandem, significantly increased the intracellular c-di-GMP level and promoted the formation of aggregates both on surfaces and in planktonic cultures. However, this regulatory effect was not contributed by either of the two classic pathways modulating biofilm formation, exopolysaccharide (EPS) overproduction or motility inhibition. Transcriptome sequencing (RNA-Seq) data revealed that the expression levels of 361 genes were significantly altered in a ΔpipA mutant strain compared to the wild type (WT), indicating the critical role of PipA in PAO1. The most remarkably downregulated genes were located on the Pf4 bacteriophage gene cluster, which corresponded to a 2-log reduction in the Pf4 phage production in the ΔpipA mutant. The sizes of aggregates in ΔpipA cultures were affected by exogenously added Pf4 phage in a concentration-dependent manner, suggesting the quantity of phage plays a part in regulating the formation of aggregates. Further analysis demonstrated that PipA is highly conserved across 83 P. aeruginosa strains. Our work therefore for the first time showed that a c-di-GMP phosphodiesterase can regulate bacteriophage production and provided new insights into the relationship between bacteriophage and bacterial aggregation. IMPORTANCE The c-di-GMP signaling pathways in P. aeruginosa are highly organized and well coordinated, with different diguanylate cyclases and phosphodiesterases playing distinct roles in a complex network. Understanding the function of each enzyme and the underlying regulatory mechanisms not only is crucial for revealing how bacteria decide the transition between motile and sessile lifestyles, but also greatly facilitates the development of new antibiofilm strategies. This work identified bacteriophage production as a novel phenotypic output controlled transcriptionally by a phosphodiesterase, PipA. Further analysis suggested that the quantity of phage may be important in regulating autoaggregation, as either a lack of phage or overproduction was associated with higher levels of aggregation. Our study therefore extended the scope of c-di-GMP-controlled phenotypes and discovered a potential signaling circuit that can be target for biofilm treatment.

Whole-Genome Analysis Reveals That Bacteriophages Promote Environmental Adaptation of Staphylococcus aureus via Gene Exchange, Acquisition, and Loss

The study of bacteriophages is experiencing a resurgence owing to their antibacterial efficacy, lack of side effects, and low production cost. Nonetheless, the interactions between Staphylococcus aureus bacteriophages and their hosts remain unexplored. In this study, whole-genome sequences of 188 S. aureus bacteriophages—20 Podoviridae, 56 Herelleviridae, and 112 Siphoviridae—were obtained from the National Center for Biotechnology Information (NCBI, USA) genome database. A phylogenetic tree was constructed to estimate their genetic relatedness using single-nucleotide polymorphism analysis. Comparative analysis was performed to investigate the structural diversity and ortholog groups in the subdividing clusters. Mosaic structures and gene content were compared in relation to phylogeny. Phylogenetic analysis revealed that the bacteriophages could be distinguished into three lineages (I–III), including nine subdividing clusters and seven singletons. The subdividing clusters shared similar mosaic structures and core ortholog clusters, including the genes involved in bacteriophage morphogenesis and DNA packaging. Notably, several functional modules of bacteriophages 187 and 2368A shared more than 95% nucleotide sequence identity with prophages in the S. aureus strain RJ1267 and the Staphylococcus pseudintermedius strain SP_11306_4, whereas other modules exhibited little nucleotide sequence similarity. Moreover, the cluster phages shared similar types of holins, lysins, and DNA packaging genes and harbored diverse genes associated with DNA replication and virulence. The data suggested that the genetic diversity of S. aureus bacteriophages was likely due to gene replacement, acquisition, and loss among staphylococcal phages, which may have crossed species barriers. Moreover, frequent module exchanges likely occurred exclusively among the subdividing cluster phages. We hypothesize that during evolution, the S. aureus phages enhanced their DNA replication in host cells and the adaptive environment of their host.

Colistin Resistance and Molecular Characterization of the Genomes of mcr-1-Positive Escherichia coli Clinical Isolates

Research on resistance against polymyxins induced by the mcr-1 gene is gaining interest. In this study, using agar dilution method, polymerase chain reaction, and comparative genomic analysis, we investigated the colistin resistance mechanism of clinical E. coli isolates. The minimum inhibitory concentration (MIC) analysis results revealed that of the 515 isolates tested, bacteria with significantly increased MIC levels against colistin were isolated in 2019. Approximately one-fifth (17.14% to 19.65%) of the isolates showed MIC values ≥1 mg/L against colistin in 2015, 2016, and 2017. However, in 2019, up to three-quarters (74.11%, 146/197) of the isolates showed MIC values ≥1 mg/L against colistin indicating an increase in colistin resistance. Six isolates (EC7518, EC4968, EC3769, EC16, EC117, EC195, 1.13%, 6/515) were found to carry the mcr-1 gene and a novel mcr-1 variant with Met2Ile mutation was identified in EC3769. All six strains showed higher MIC levels (MIC=4 mg/L) than any mcr-1-negative strains (MIC ≤ 2 mg/L). Whole-genome sequencing of the six mcr-1-positive isolates revealed that EC195 carried the highest number of resistance genes (n = 28), nearly a half more than those of the following EC117 (n = 19). Thus, EC195 showed a wider resistance spectrum and higher MIC levels against the antimicrobials tested than the other five isolates. Multi-locus sequence typing demonstrated that these mcr-1-positive strains belonged to six different sequence types. The six mcr-1 genes were located in three different incompatibility group plasmids (IncI2, IncHI2 and IncX4). The genetic context of mcr-1 was related to a sequence derived from Tn6330 (ISApl1-mcr-1-pap2-ISApl1). Investigations into the colistin resistance mechanism and characterization of the molecular background of the mcr genes may help trace the development and spread of colistin resistance in clinical settings.

Comparison of Two Distinct Subpopulations of Klebsiella pneumoniae ST16 Co-Occurring in a Single Patient

The higher resistance rate to ceftazidime-avibactam (CZA) is mainly related to carbapenem resistance, especially New Delhi metallo-β-lactamase (NDM). The CZA-susceptible Klebsiella pneumoniae (K191663) and the later CZA-resistant isolates (K191724, K191725, K191773) co-producing NDM-4 and OXA-181 were obtained from the same hospitalized patient returning from Vietnam. Our study aims to elucidate the diversity of K. pneumoniae ST16 through comparative analysis of whole-genome sequencing (WGS) data and identify the potential evolution of plasmids by sequencing longitudinal clinical isolates during antibiotic treatment. Firstly, multilocus sequence typing analysis and phylogenic analysis suggested that these strains belong to the two lineages of K. pneumoniae ST16. Surprisingly, the CZA-resistant strains were closely related to K. pneumoniae ST16 described in South Korea, instead of the bla_NDM-4- or bla_OXA-181-carrying ST16 reported in Vietnam. Secondly, bla_NDM-4, bla_TEM-1B, and rmtB co-existed on a self-conjugative IncFII(Yp)-like plasmid, which played a significant role in CZA resistance. It could transfer into the recipient Escherichia coli J53 at high frequency, indicating the risk of mobile carbapenemases. In addition, the loss of 12-kbp fragment occurred in bla_NDM-4-positive isolate (K191773), which was likely caused by insertion sequence-mediated homologous recombination. Last but not least, as a repressor of acrAB operon system, acrR was truncated by a frameshift mutation in K191663. Thus, our study provided baseline information for monitoring the occurrence and development of bacterial resistance.

IMPORTANCE As a leading health care-acquired infection pathogen, Klebsiella pneumoniae is threatening a large number of inpatients due to its diverse antibiotic resistance and virulence factors. Heretofore, with a growing number of reports about the coexistence of several carbapenemases in carbapenem-resistant K. pneumoniae (CRKP), epidemiologic surveillance has been strengthened. Nevertheless, the nosocomial outbreaks by CRKP ST16 are gradually increasing worldwide. Our study provides a deeper insight into the diversification of clinical isolates of CRKP ST16 in China. In addition, the comparison analysis of resistant plasmids may reveal the transmission of carbapenemase-encoding genes. Furthermore, our study also highlights the importance of longitudinal specimen collection and continuous monitoring during the treatment, which play a crucial role in understanding the development of antibiotic resistance and the evolution of resistance plasmids.

Rapid Identification and Source Tracing of a Salmonella Typhimurium Outbreak in China by Metagenomic and Whole-Genome Sequencing

Salmonella spp. are among the most prevalent foodborne pathogens. Rapid identification of etiologic agents during foodborne outbreaks is of great importance. In this study, we report a traceback investigation of a Salmonella outbreak in China. Metagenomic sequencing of suspected food samples was performed on MinION and MiSeq platforms. Real-time nanopore sequencing analysis identified reads belonging to the Enterobacteriaceae family. MiSeq sequencing identified 63 reads specifically mapped to Salmonella. Conventional methods including quantitative-PCR and culture-based isolation confirmed as Salmonella enterica serovar Typhimurium. The foodborne outbreak of Salmonella Typhimurium was further recognized by whole-genome sequencing and pulsed-field gel electrophoresis analysis. Our study demonstrates the ability of metagenomic sequencing to rapidly identify enteric pathogens directly from food samples. These results highlight the capacity of metagenomic sequencing to deliver actionable information rapidly and to expedite the tracing and identification of etiologic agents during foodborne outbreaks.

Whole-Genome Sequencing-Based Re-Identification of Pseudomonas putida/fluorescens Clinical Isolates Identified by Biochemical Bacterial Identification Systems

The genus Pseudomonas, a complex Gram-negative genus, includes species isolated from various environments, plants, animals, and humans. We compared whole-genome sequencing (WGS) with clinical bacteriological methods and evaluated matrix-assisted laser desorption ionization-time of flight mass spectrometry (MALDI-TOF MS) to identify Pseudomonas species. Clinical isolates (N = 42) identified as P. putida or P. fluorescens by a bacterial identification system based on biochemical properties were reexamined by another identification system based on biochemical properties, two systems based on MALDI-TOF MS, and WGS. WGS revealed that 30 of the 42 isolates belonged to one of 14 known Pseudomonas species, respectively. The remaining 12 belonged to one of 9 proposed novel Pseudomonas species, respectively. MALDI-TOF MS analysis showed that the 9 novel species had unique major peaks. These results suggest that WGS is the optimal method to identify Pseudomonas species and that MALDI-TOF MS may complement WGS in identification. Based on their morphologic, physiologic, and biochemical properties, we propose nine novel Pseudomonas species.

IMPORTANCE Most of the clinical isolates, identified as P. putida or P. fluorescens, were misidentified in clinical laboratories. Whole-genome sequencing (WGS) revealed that these isolates belonged to different Pseudomonas species, including novel species. WGS is a gold-standard method to identify Pseudomonas species, and MALDI-TOF MS analysis has the potential to complement WGS to reliably identify them.

Purification and characterization of antifungal lipopeptide produced by Bacillus velezensis isolated from raw honey

Raw honey contains a diverse microbiota originating from honeybees, plants, and soil. Some gram-positive bacteria isolated from raw honey are known for their ability to produce secondary metabolites that have the potential to be exploited as antimicrobial agents. Currently, there is a high demand for natural, broad-spectrum, and eco-friendly bio-fungicides in the food industry. Naturally occurring antifungal products from food-isolated bacteria are ideal candidates for agricultural applications. To obtain novel antifungals from natural sources, we isolated bacteria from raw clover and orange blossom honey to evaluate their antifungal-producing potential. Two Bacillus velezensis isolates showed strong antifungal activity against food-isolated fungal strains. Antifungal compound production was optimized by adjusting the growth conditions of these bacterial isolates. Extracellular proteinaceous compounds were purified via ammonium sulfate precipitation, solid phase extraction, and RP-HPLC. Antifungal activity of purified products was confirmed by deferred overlay inhibition assay. Mass spectrometry (MS) was performed to determine the molecular weight of the isolated compounds. Whole genome sequencing (WGS) was conducted to predict secondary metabolite gene clusters encoded by the two antifungal-producing strains. Using MS and WGS data, we determined that the main antifungal compound produced by these two Bacillus velezensis isolates was iturin A, a lipopeptide exhibiting broad spectrum antifungal activity.

Evaluating a Cluster and the Overall Trend of Invasive Haemophilus influenzae Serotype b in Alaska 2005-2019

BACKGROUND

In 2019, 5 cases of invasive Haemophilus influenzae serotype b (Hib) occurred in the Anchorage region of Alaska over a period of 16 days. No cases had occurred in Alaska in the preceding 26 months.

METHODS

Alaska Hib isolates from 2005 through 2019 were analyzed using whole-genome sequencing (WGS). Rates were compared with the CDC's Active Bacterial Core surveillance (ABCs) data.

RESULTS

A total of 33 cases of invasive Hib occurred in Alaska from 2005 through 2019. Of the 5 cases associated with the cluster, 2 (40%) occurred in adults and all occurred in the Anchorage region. In contrast, only 14% (4/28) of the noncluster cases occurred in this region (P < 0.01). Two cluster cases were linked epidemiologically and the bacteria were nearly identical. The other 3 cluster cases were caused by 3 genetically distinct bacteria. When the full period was evaluated, the unadjusted rate of invasive Hib disease in Alaska was 15.5 times higher in Alaska Native (AN) people than non-AN people [1.3/100,000 vs. 0.07/100,000, 95% confidence intervals (CI): 10.2-22.5). The age-adjusted rate of invasive Hib disease in Alaska was 9.4 times higher than the ABCs rate (95% CI: 6.3-14.1).

CONCLUSIONS

While clustered in time and space, the 5 cases in 2019 were not due to a single bacterial strain. AN people continue to have elevated rates of invasive Hib infection compared with both non-AN people in Alaska and the ABCs population.

Dataset of genome sequence, de novo assembly, and functional annotation of Ruegeria sp. (PBVC088), a marine bacterium associated with the toxin-producing harmful dinoflagellate, Pyrodinium bahamense var. compressum

The dataset comprises a whole-genome sequence of Ruegeria sp. PBVC088, a symbiotic (Gram-negative) bacterium associated with Pyrodinium bahamense var. compressum, which has been associated with harmful algal blooms in the coastal waters of west Sabah, Malaysia. Harmful algal blooms contribute to economic losses for the aquaculture industry, as well as human illnesses and fatalities due to paralytic shellfish poisoning. Bacteria-algae dynamics have posited that the interaction is potentially responsible for the toxin production during a toxic harmful algal bloom event. Despite the expanding body of literature on the capabilities of these bacteria to metabolize, produce, and modify toxins autonomously, it has yet to be confirmed that these toxin-producing bacteria are capable of autonomous toxin synthesis. Saxitoxin, a paralytic shellfish poisoning toxin, is produced by a unique biosynthetic pathway, where the genetic basis for the saxitoxin production was first reported in the saxitoxin-producing cyanobacteria strain Cylindrospermopsis raciborskii T3 (NCBI accession no. DQ787200). The genes responsible for saxitoxin biosynthesis in dinoflagellates, have yet to be fully elucidated. The identification of cyanobacteria saxitoxin biosynthesis genes (sxt) may eventually lead to the identification of homologous genes within the dinoflagellates. Previous studies on the diversity of the bacterial communities associated with the same toxic P. bahamense harmful alga has been carried out by using both the culture-dependent 16S ribosomal RNA gene sequence analysis and culture-independent 16S metagenomic sequence analysis. This study extends the knowledge pertaining to the genomic aspect of an associated bacterium isolated from P. bahamense alga by adopting a whole genome sequencing approach. Here, we report the genome sequencing, de novo assembly, and annotation data of a bacterium, Ruegeria sp. PBVC088, associated with harmful alga P. bahamense, which can be referenced by researchers to identify the genes and pathways related to toxin biosynthesis from a much larger data set. The genome of Ruegeria sp. PBVC088 was sequenced using the Illumina MiSeq platform with 250 bp paired-end reads. The number of reads generated from the MiSeq sequencer was 1,135,484, with an estimated coverage of 100X. The estimated genome size for the marine bacterium was computed to be 5.78 Mb. Annotation of the genome predicted 5,689 gene sequences, which were assigned putative functions based on homology to existing protein sequences in public databases. In addition, annotation of genes related to saxitoxin biosynthesis pathway was also performed. Raw fastq reads and the final version of the genome assembly have been deposited in the National Center for Biotechnology Information (NCBI) (BioProject: PRJNA324753, WGS: LZNT00000000, SRA: SRR3646181). The genome data provided here are expected to better understand the genetic processes involved in saxitoxin biosynthesis in marine bacteria associated with dinoflagellates.

Draft genome sequence data of a clinical Enterococcus faecalis isolate SHH039 from a patient with cholecystitis from a tertiary care hospital in Sabah, Malaysia

An Enterococcus faecalis strain SHH039 was isolated from a 68 year old man who was hospitalised with cholecystitis. The genomic sequence of this isolate which had a size of 2,990,081 bp and 2,663 proteins with functional assignments is presented here. Analysis of the genome revealed Enterococcus faecalis with multiple antibiotic resistance genes which may be associated with acute cholecystitis. It may be not clear if the infection symptoms are the consequence of enterococci manifestation. However, this opportunistic organism may play a minor role in the disease.

Integration of bioinformatic and chemoproteomic tools for the study of enzyme conservation in closely related bacterial species

Activity-based protein profiling (ABPP) is a commonly utilized technique to globally characterize the endogenous activity of multiple enzymes within a related family. While it has been used extensively to identify enzymes that are differentially active across various mammalian tissues, recent efforts have expanded this technique to studying bacteria. As ABPP is applied to diverse sets of bacterial strains found in microbial communities, there is also an increasing need for robust tools for assessing the conservation of enzymes across closely related bacterial species and strains. In this chapter, we detail the integration of gel-based ABPP with basic bioinformatic tools to enable the analysis of enzyme activity, distribution, and homology. We use as an example the family of serine hydrolases identified in the skin commensal bacterium Staphylococcus epidermidis.

Phylogenetic Structure and Comparative Genomics of Multi-National Invasive Haemophilus influenzae Serotype a Isolates

Recent reports have indicated a rise of invasive disease caused by Haemophilus influenzae serotype a (Hia) in North America and some European countries. The whole-genome sequences for a total of 410 invasive Hia isolates were obtained from 12 countries spanning the years of 1998 to 2019 and underwent phylogenetic and comparative genomic analysis in order to characterize the major strains causing disease and the genetic variation present among factors contributing to virulence and antimicrobial resistance. Among 410 isolate sequences received, 408 passed our quality control and underwent genomic analysis. Phylogenetic analysis revealed that the Hia isolates formed four genetically distinct clades: clade 1 (n = 336), clade 2 (n = 13), clade 3 (n = 3) and clade 4 (n = 56). A low diversity subclade 1.1 was found in clade 1 and contained almost exclusively North American isolates. The predominant sequence types in the Hia collection were ST-56 (n = 125), ST-23 (n = 98) and ST-576 (n = 51), which belonged to clade 1, and ST-62 (n = 54), which belonged to clade 4. Clades 1 and 4 contained predominantly North American isolates, and clades 2 and 3 predominantly contained European isolates. Evidence of the presence of capsule duplication was detected in clade 1 and 2 isolates. Seven of the virulence genes involved in endotoxin biosynthesis were absent from all Hia isolates. In general, the presence of known factors contributing to β-lactam antibiotic resistance was low among Hia isolates. Further tests for virulence and antibiotic susceptibility would be required to determine the impact of these variations among the isolates.

Benchmarking the topological accuracy of bacterial phylogenomic workflows using in silico evolution

Phylogenetic analyses are widely used in microbiological research, for example to trace the progression of bacterial outbreaks based on whole-genome sequencing data. In practice, multiple analysis steps such as de novo assembly, alignment and phylogenetic inference are combined to form phylogenetic workflows. Comprehensive benchmarking of the accuracy of complete phylogenetic workflows is lacking. To benchmark different phylogenetic workflows, we simulated bacterial evolution under a wide range of evolutionary models, varying the relative rates of substitution, insertion, deletion, gene duplication, gene loss and lateral gene transfer events. The generated datasets corresponded to a genetic diversity usually observed within bacterial species (≥95 % average nucleotide identity). We replicated each simulation three times to assess replicability. In total, we benchmarked 19 distinct phylogenetic workflows using 8 different simulated datasets. We found that recently developed k-mer alignment methods such as kSNP and ska achieve similar accuracy as reference mapping. The high accuracy of k-mer alignment methods can be explained by the large fractions of genomes these methods can align, relative to other approaches. We also found that the choice of de novo assembly algorithm influences the accuracy of phylogenetic reconstruction, with workflows employing SPAdes or skesa outperforming those employing Velvet. Finally, we found that the results of phylogenetic benchmarking are highly variable between replicates. We conclude that for phylogenomic reconstruction, k-mer alignment methods are relevant alternatives to reference mapping at the species level, especially in the absence of suitable reference genomes. We show de novo genome assembly accuracy to be an underappreciated parameter required for accurate phylogenomic reconstruction.

The Association Between Onset of Staphylococcal Non-menstrual Toxic Shock Syndrome With Inducibility of Toxic Shock Syndrome Toxin-1 Production

Non-menstrual toxic shock syndrome (non-mTSS) is a life-threatening disease caused by Staphylococcus aureus strains producing superantigens, such as staphylococcal enterotoxins A, B, C, and toxic shock syndrome toxin-1 (TSST-1). However, little is known about why the TSS cases are rare, although S. aureus strains frequently carry a tst gene, which encodes TSST-1. To answer this question, the amount of TSST-1 produced by 541 clinical isolates was measured in both the presence and absence of serum supplementation to growth media. Then a set of S. aureus strains with similar genetic backgrounds isolated from patients presenting with non-mTSS and those with clinical manifestations other than non-mTSS was compared for their TSST-1 inducibility by human serum, and their whole-genome sequences were determined. Subsequently, the association of mutations identified in the tst promoter of non-mTSS strains with TSST-1 inducibility by human serum was evaluated by constructing promoter replacement mutants and green fluorescent protein (GFP) reporter recombinants. Results showed that 39 out of 541 clinical isolates (7.2%), including strains isolated from non-mTSS patients, had enhanced production of TSST-1 in the presence of serum. TSST-1 inducibility by human serum was more clearly seen in non-mTSS strains of clonal complex (CC)-5. Moreover, the whole-genome sequence analysis identified a set of sequence variations at a putative SarA-binding site of the tst promoter. This sequence variation was proven to be partially responsible for the induction of TSST-1 production by human serum. We conclude that the onset of staphylococcal toxic shock syndrome caused by TSST-1-producing CC-5 strains seem at least partially initiated by serum induction of TSST-1, which is regulated by the mutation of putative SarA-binding site at the tst promoter.

Escherichia fergusonii, an Underrated Repository for Antimicrobial Resistance in Food Animals

A total of 1,400 samples of food animals (pigs, chickens, and ducks) were collected between July and September 2019 in China to uncover the prevalence of E. fergusonii and its potential role in the evolution of antimicrobial resistance (AMR). An isolation of E. fergusonii was performed and pulsed-field gel electrophoresis (PFGE) was used to uncover the genetic relationship. The AMR of E. fergusonii isolates was comprehensively characterized using broth microdilution-based antimicrobial susceptibility testing, S1-PFGE, southern hybridization, whole-genome sequencing, and in-depth bioinformatics analysis. As a result, a total of 133 E. fergusonii isolates were obtained. These isolates could be grouped into 41 PFGE subclades, suggesting a diverse genetic relationship. The resistance phenotypes of sulfafurazole (97.74%) and tetracycline (94.74%) were the most frequently found. Of the E. fergusonii isolates, 51.88% were extended spectrum beta-lactamase (ESBL)-positive. Forty-three different AMR genes were revealed based on 25 genome sequences harboring mcr-1. Briefly, aph(6)-Id, aph(3'')-Ib and tet(A) genes were the most frequently observed, with the highest rate being 76.00% (19/25). Three mcr-1-harboring plasmids were identified after Nanopore sequencing, including pTB31P1 (IncHI2-IncHI2A, 184,652 bp), pTB44P3 (IncI2, 62,882 bp), and pTB91P1 (IncHI2-IncHI2A, 255,882 bp). Additionally, 25 E. fergusonii isolates harboring mcr-1 were clustered together with other E. fergusonii isolates from different regions and sources available in GenBank, suggesting a possible random process of mcr-1 transmission in E. fergusonii. In conclusion, E. fergusonii is widespread in food animals in China and might be an important reservoir of AMR genes, especially mcr-1, and facilitate the evolution of AMR. IMPORTANCEE. fergusonii, a member of the genus Escherichia, has been reported to transmit via the food chain and cause diseases in humans. However, the prevalence of multidrug-resistant E. fergusonii, especially mcr-1-positive E. fergusonii isolates, has rarely been reported. Here, we collected 1,400 samples from food animals in three provinces of China and obtained 133 E. fergusonii isolates (9.5%). We found that the prevalence of E. fergusonii isolates was diverse, with high levels of antimicrobial resistance. Among them, 18.8% E. fergusonii isolates carried the colistin resistance gene mcr-1. Thus, E. fergusonii may facilitate the evolution of colistin resistance as a reservoir of mcr-1. As far as we know, the prevalence and AMR of E. fergusonii in the food animals in this study was first reported in China. These findings increase our understanding of the role of E. fergusonii in public health and the evolution of antibiotic resistance.

Long-term spatiotemporal variation of antimicrobial resistance genes within the Serratia marcescens population and transmission of S. marcescens revealed by public whole-genome datasets

The development of antimicrobial resistance (AMR) is accelerated by the selective pressure exerted by the widespread use of antimicrobial drugs, posing an increasing danger to public health. However, long-term spatiotemporal variation in AMR genes in microorganisms, particularly in bacterial pathogens in response to antibiotic consumption, is not fully understood. Here, we used the NCBI RefSeq database to collect 478 whole-genome sequences for Serratia marcescens ranging from 1961 up to 2019, to document global long-term AMR trends in S. marcescens populations. In total, 100 AMR gene subtypes (16 AMR gene types) were detected in the genomes of S. marcescens populations. We identified 3 core resistance genes in S. marcescens genomes, and a high diversity of AMR genes was observed in S. marcescens genomes after corresponding antibiotics were discovered and introduced into clinical practice, suggesting the adaptation of S. marcescens populations to challenges with therapeutic antibiotics. Our findings indicate spatiotemporal variation of AMR genes in S. marcescens populations in relation to antibiotic consumption and suggest the potential transmission of S. marcescens isolates harboring AMR genes among countries and between the environment and the clinic, representing a public health threat that necessitates international solidarity to overcome.

Whole-Genome Sequencing of Streptomycin-Resistant Mycobacterium tuberculosis Strain SBH145 from Sabah, Malaysia

This paper reports on the whole-genome sequencing of a streptomycin-resistant Mycobacterium tuberculosis strain that was isolated from a patient with pulmonary tuberculosis in Sabah state of Malaysian Borneo. The strain belongs to the EAI2-Manila family of lineage 1 and is clustered with M. tuberculosis strains from the Philippines, India, and Taiwan.

Population genomics reveals historical and ongoing recombination in the Fusarium oxysporum species complex

The Fusarium oxysporum species complex (FOSC) is a group of closely related plant pathogens long-considered strictly clonal, as sexual stages have never been recorded. Several studies have questioned whether recombination occurs in FOSC, and if it occurs its nature and frequency are unknown. We analysed 410 assembled genomes to answer whether FOSC diversified by occasional sexual reproduction interspersed with numerous cycles of asexual reproduction akin to a model of predominant clonal evolution (PCE). We tested the hypothesis that sexual reproduction occurred in the evolutionary history of FOSC by examining the distribution of idiomorphs at the mating locus, phylogenetic conflict and independent measures of recombination from genome-wide SNPs and genes. A phylogenomic dataset of 40 single copy orthologs was used to define structure a priori within FOSC based on genealogical concordance. Recombination within FOSC was tested using the pairwise homoplasy index and divergence ages were estimated by molecular dating. We called SNPs from assembled genomes using a k-mer approach and tested for significant linkage disequilibrium as an indication of PCE. We clone-corrected and tested whether SNPs were randomly associated as an indication of recombination. Our analyses provide evidence for sexual or parasexual reproduction within, but not between, clades of FOSC that diversified from a most recent common ancestor about 500 000 years ago. There was no evidence of substructure based on geography or host that might indicate how clades diversified. Competing evolutionary hypotheses for FOSC are discussed in the context of our results.

Novel Putative Transposable Element Associated with the Subtype E5 Botulinum Toxin Gene Cluster of Neurotoxigenic Clostridium butyricum Type E Strains from China

Previously, a whole-genome comparison of three Clostridium butyricum type E strains from Italy and the United States with different C. botulinum type E strains indicated that the bont/e gene might be transferred between the two clostridia species through transposition. However, transposable elements (TEs) have never been identified close to the bont/e gene. Herein, we report the whole genome sequences for four neurotoxigenic C. butyricum type E strains that originated in China. An analysis of the obtained genome sequences revealed the presence of a novel putative TE upstream of the bont/e gene in the genome of all four strains. Two strains of environmental origin possessed an additional copy of the putative TE in their megaplasmid. Similar putative TEs were found in the megaplasmids and, less frequently, in the chromosomes of several C. butyricum strains, of which two were neurotoxigenic C. butyricum type E strains, and in the chromosome of a single C. botulinum type E strain. We speculate that the putative TE might potentially transpose the bont/e gene at the intracellular and inter-cellular levels. However, the occasional TE occurrence in the clostridia genomes might reflect rare transposition events.

Development of a Genomics-Based Approach To Identify Putative Hypervirulent Nontyphoidal Salmonella Isolates: Salmonella enterica Serovar Saintpaul as a Model

While differences in human virulence have been reported across nontyphoidal Salmonella (NTS) serovars and associated subtypes, a rational and scalable approach to identify Salmonella subtypes with differential ability to cause human diseases is not available. Here, we used NTS serovar Saintpaul (S. Saintpaul) as a model to determine if metadata and associated whole-genome sequence (WGS) data in the NCBI Pathogen Detection (PD) database can be used to identify (i) subtypes with differential likelihoods of causing human diseases and (ii) genes and single nucleotide polymorphisms (SNPs) potentially responsible for such differences. S. Saintpaul SNP clusters (n = 211) were assigned different epidemiology types (epi-types) based on statistically significant over- or underrepresentation of human clinical isolates, including human associated (HA; n = 29), non-human associated (NHA; n = 23), and other (n = 159). Comparative genomic analyses identified 384 and 619 genes overrepresented among isolates in 5 HA and 4 NHA SNP clusters most significantly associated with the respective isolation source. These genes included 5 HA-associated virulence genes previously reported to be present on Gifsy-1/Gifsy-2 prophages. Additionally, premature stop codons in 3 and 7 genes were overrepresented among the selected HA and NHA SNP clusters, respectively. Tissue culture experiments with strains representing 4 HA and 3 NHA SNP clusters did not reveal evidence for enhanced invasion or intracellular survival for HA strains. However, the presence of sodCI (encoding a superoxide dismutase), found in 4 HA and 1 NHA SNP clusters, was positively correlated with intracellular survival in macrophage-like cells. Post hoc analyses also suggested a possible difference in intracellular survival among S. Saintpaul lineages. IMPORTANCE Not all Salmonella isolates are equally likely to cause human disease, and Salmonella control strategies may unintentionally focus on serovars and subtypes with high prevalence in source populations but are rarely associated with human clinical illness. We describe a framework leveraging WGS data in the NCBI PD database to identify Salmonella subtypes over- and underrepresented among human clinical cases. While we identified genomic signatures associated with HA/NHA SNP clusters, tissue culture experiments failed to identify consistent phenotypic characteristics indicative of enhanced human virulence of HA strains. Our findings illustrate the challenges of defining hypo- and hypervirulent S. Saintpaul and potential limitations of phenotypic assays when evaluating human virulence, for which in vivo experiments are essential. Identification of sodCI, an HA-associated virulence gene associated with enhanced intracellular survival, however, illustrates the potential of the framework and is consistent with prior work identifying specific genomic features responsible for enhanced or reduced virulence of nontyphoidal Salmonella.

Comparision of biological and genomic characteristics of five virulent bacteriophages against Enterobacter hormaechei

Enterobacter hormaechei is a zoonotic bacteria that may cause respiratory diseases in animals and neonatal sepsis in humans. Bacteriophages are increasingly considered as potential biocontrol agents to control pathogens in the food industry. In this study, five E. hormaechei virulent phages, named as Ehp-YZU08, Ehp-YZU10, Ehp-YZU9-1, Ehp-YZU9-2 and Ehp-YZU9-3, were isolated from sewage in China and analyzed for their biological and whole-genome characteristics, and a comparative genomic analysis was performed to study the functional genes and phylogenetic evolution of phages. The results showed that four of the phage strains belong to the Podoviridae family and one belongs to the Myoviridae family. The burst sizes were 70-283 PFU/cell after a latent period of 5-40 min. Phages were able to survive in a pH range of 5-10 and resist temperatures up to 60 °C for 60 min. The sequencing results showed that the full length of the genomes of the five phages ranged from 39,502 to 173,418 bp. Each phage contained multiple genes related to phage replication, and genes related to bacterial virulence or drug resistance were not found. The five phages belonged to three different groups by a construction of a phylogenetic tree, and the significant genetic evolutionary distance from each E. hormaechei phage was observed. The inhibition assay showed that all five phages could completely inhibit the growth of E. hormaechei at 37 °C within 8 h, suggesting that the phages in this study have great potential for the development of biocontrol agents against E. hormaechei in the food industry.

Prevalence and Molecular Characterization of Fluoroquinolone-Resistant Escherichia coli in Healthy Children

Faecal E. coli can act as reservoirs for resistance genes. Here, we analyzed prevalence of drug resistance in faecal E. coli isolated from healthy children at a single kindergarten in Beijing, China, then used whole genome sequencing to characterize fluoroquinolone-non-susceptible strains. Our results revealed high resistance to ampicillin (54.0%), trimethoprim/sulphurmethoxazole (47.5%) and tetracycline (58.9%) among 576 faecal E. coli isolates, 49.2% of which exhibited multidrug resistance. A total of 113 E. coli isolates were not susceptible to ciprofloxacin, with four sequence types, namely ST1193 (25.7%), ST773 (13.3%), ST648 (8.8%) and ST131 (7.1%) found to be the most prevalent (54.9%). With regards to resistance to quinolones, we detected chromosomal mutations in gyrA, parC, and parE in 111 (98.2%), 105 (92.9%), and 67 (61.1%) isolates, respectively. bla _CTX-M (37.2%) was the major ESBL gene, whereas bla _CTX-M-14 (12.4%) and bla _CTX-M-27 (11.5%) were the most frequent subtypes. A total of 90 (79.6%) ExPEC and 65 (57.5%) UPEC isolates were classified. Overall, these findings revealed clonal spread of certain prevalent STs, namely ST1193, ST773, ST648 and ST131 E. coli isolates in healthy children within a single kindergarten in Beijing, China, affirming the seriousness of the multidrug resistance problem and potential pathogenicity of E. coli isolates in healthy children. Therefore, there is an urgent need for increased surveillance to enhance control of this problem.

Whole-Genome Analysis of Multidrug-Resistant Salmonella Enteritidis Strains Isolated from Poultry Sources in Korea

The Salmonella Enterica subsp. Enterica serovar Enteritidis is one of main serovars isolated from human patients with food poisoning and poultry without clinical signs. Consumption of poultry products contaminated with Salmonella Enteritidis is a common source of human salmonellosis; 82 Salmonella spp. were isolated from 291 samples of retail chicken meat, 201 one-day-old chicks, 30 internal organs of chickens, 156 chicken eggs, 100 duck eggs, 38 straw bedding samples, 18 samples of retail duck meat, and 19 swab samples from slaughterhouses in 2019 and 2020. An antibiotic susceptibility test was performed for all isolates, revealing 33 multidrug-resistant (MDR) strains. The whole genome of 33 MDR strains isolated in 2019 and 2020 and 10 strains isolated in 2011, 2012, and 2017 was sequenced using the MinION sequencing protocol. Within these 43 samples, 5 serovars were identified: S. Enteritidis, S. Agona, S. Virchow, S. Albany, and S. Bareilly. The most common serovar was S. Enteritidis (26/43), which showed the highest resistance to ampicillin (100%), followed by nalidixic acid (90%) and colistin (83%). Core genome multilocus sequence typing analysis showed that the S. Enteritidis strains isolated from different sources and in different years were clustered together. In addition, the S. Enteritidis strains isolated since 2011 consistently harbored the same antibiotic resistance patterns.

Comparative Genomic and Pan-Genomic Characterization of Staphylococcus epidermidis From Different Sources Unveils the Molecular Basis and Potential Biomarkers of Pathogenic Strains

Coagulase-negative Staphylococcus (CoNS) is the most common pathogen causing traumatic endophthalmitis. Among which, Staphylococcus epidermidis is the most common species that colonizes human skin, eye surfaces, and nasal cavity. It is also the main cause of nosocomial infection, specially foreign body-related bloodstream infections (FBR-BSIs). Although some studies have reported the genome characteristics of S. epidermidis, the genome of ocular trauma-sourced S. epidermidis strain and a comprehensive understanding of its pathogenicity are still lacking. Our study sequenced, analyzed, and reported the whole genomes of 11 ocular trauma-sourced samples of S. epidermidis that caused traumatic endophthalmitis. By integrating publicly available genomes, we obtained a total of 187 S. epidermidis samples from healthy and diseased eyes, skin, respiratory tract, and blood. Combined with pan-genome, phylogenetic, and comparative genomic analyses, our study showed that S. epidermidis, regardless of niche source, exhibits two founder lineages with different pathogenicity. Moreover, we identified several potential biomarkers associated with the virulence of S. epidermidis, including essD, uhpt, sdrF, sdrG, fbe, and icaABCDR. EssD and uhpt have high homology with esaD and hpt in Staphylococcus aureus, showing that the genomes of S. epidermidis and S. aureus may have communicated during evolution. SdrF, sdrG, fbe, and icaABCDR are related to biofilm formation. Compared to S. epidermidis from blood sources, ocular-sourced strains causing intraocular infection had no direct relationship with biofilm formation. In conclusion, this study provided additional data resources for studies on S. epidermidis and improved our understanding of the evolution and pathogenicity among strains of different sources.

Long-Read Sequencing Reveals Evolution and Acquisition of Antimicrobial Resistance and Virulence Genes in Salmonella enterica

Salmonella enterica is a significant and phylogenetically diverse zoonotic pathogen. To understand its genomic heterogeneity and antimicrobial resistance, we performed long-read sequencing on Salmonella isolated from retail meats and food animals. A collection of 134 multidrug-resistant isolates belonging to 33 serotypes were subjected to PacBio sequencing. One major locus of diversity among these isolates was the presence and orientation of Salmonella pathogenic islands (SPI), which varied across different serotypes but were largely conserved within individual serotypes. We also identified insertion of an IncQ resistance plasmid into the chromosome of fourteen strains of serotype I 4,[5],12:i:- and the Salmonella genomic island 1 (SGI-1) in five serotypes. The presence of various SPIs, SGI-1 and integrated plasmids contributed significantly to the genomic variability and resulted in chromosomal resistance in 55.2% (74/134) of the study isolates. A total of 93.3% (125/134) of isolates carried at least one plasmid, with isolates carrying up to seven plasmids. We closed 233 plasmid sequences of thirteen replicon types, along with twelve hybrid plasmids. Some associations between Salmonella isolate source, serotype, and plasmid type were seen. For instance, IncX plasmids were more common in serotype Kentucky from retail chicken. Plasmids IncC and IncHI had on average more than five antimicrobial resistance genes, whereas in IncX, it was less than one per plasmid. Overall, 60% of multidrug resistance (MDR) strains that carried >3 AMR genes also carried >3 heavy metal resistance genes, raising the possibility of co-selection of antimicrobial resistance in the presence of heavy metals. We also found nine isolates representing four serotypes that carried virulence plasmids with the spv operon. Together, these data demonstrate the power of long-read sequencing to reveal genomic arrangements and integrated plasmids with a high level of resolution for tracking and comparing resistant strains from different sources. Additionally, the findings from this study will help expand the reference set of closed Salmonella genomes that can be used to improve genome assembly from short-read data commonly used in One Health antimicrobial resistance surveillance.

Associations between Listeria monocytogenes genomic characteristics and adhesion to polystyrene at 8 °C

Listeria monocytogenes remains a threat to the food system and has led to numerous foodborne outbreaks worldwide. L. monocytogenes can establish itself in food production facilities by adhering to surfaces, resulting in increased resistance to environmental stressors. The aim of this study was to evaluate the adhesion ability of L. monocytogenes at 8 °C and to analyse associations between the observed phenotypes and genetic factors such as internalin A (inlA) genotypes, stress survival islet 1 (SSI-1) genotype, and clonal complex (CC). L. monocytogenes isolates (n = 184) were grown at 8 °C and 100% relative humidity for 15 days. The growth was measured by optical density at 600 nm every 24 h. Adherent cells were stained using crystal violet and quantified spectrophotometrically. Genotyping of inlA and SSI-1, multi-locus sequence typing, and a genome-wide association study (GWAS) were performed to elucidate the phenotype-genotype relationships in L. monocytogenes cold adhesion. Among all inlA genotypes, truncated inlA isolates had the highest mean adhered cells, ABS595nm = 0.30 ± 0.15 (Tukey HSD; P < 0.05), while three-codon deletion inlA isolates had the least mean adhered cells (Tukey HSD; P < 0.05). When SSI-1 was present, more cells adhered; less cells adhered when SSI-1 was absent (Welch's t-test; P < 0.05). Adhesion was associated with clonal complexes which have low clinical frequency, while reduced adhesion was associated with clonal complexes which have high frequency. The results of this study support that premature stop codons in the virulence gene inlA are associated with increased cold adhesion and that an invasion enhancing deletion in inlA is associated with decreased cold adhesion. This study also provides evidence to suggest that there is an evolutionary trade off between virulence and adhesion in L. monocytogenes. These results provide a greater understanding of L. monocytogenes adhesion which will aid in the development of strategies to reduce L. monocytogenes in the food system.

Extensive Genome Exploration of Clostridium botulinum Group III Field Strains

In animals, botulism is commonly sustained by botulinum neurotoxin C, D or their mosaic variants, which are produced by anaerobic bacteria included in Clostridium botulinum group III. In this study, a WGS has been applied to a large collection of C. botulinum group III field strains in order to expand the knowledge on these BoNT-producing Clostridia and to evaluate the potentiality of this method for epidemiological investigations. Sixty field strains were submitted to WGS, and the results were analyzed with respect to epidemiological information and compared to published sequences. The strains were isolated from biological or environmental samples collected in animal botulism outbreaks which occurred in Italy from 2007 to 2016. The new sequenced strains belonged to subspecific groups, some of which were already defined, while others were newly characterized, peculiar to Italian strains and contained genomic features not yet observed. This included, in particular, two new flicC types (VI and VII) and new plasmids which widen the known plasmidome of the species. The extensive genome exploration shown in this study improves the C. botulinum and related species classification scheme, enriching it with new strains of rare genotypes and permitting the highest grade of discrimination among strains for forensic and epidemiological applications.

In vivo and in silico Virulence Analysis of Leptospira Species Isolated From Environments and Rodents in Leptospirosis Outbreak Areas in Malaysia

The zoonotic disease leptospirosis is caused by pathogenic species of the genus Leptospira. With the advancement of studies in leptospirosis, several new species are being reported. It has always been a query, whether Leptospira species, serovars, and strains isolated from different geographical locations contribute to the difference in the disease presentations and severity. In an epidemiological surveillance study performed in Malaysia, we isolated seven novel intermediate and saprophytic species (Leptospira semungkisensis, Leptospira fletcheri, Leptospira langatensis, Leptospira selangorensis, Leptospira jelokensis, Leptospira perdikensis, Leptospira congkakensis) from environments and three pathogenic species from rodents (Leptospira borgpetersenii strain HP364, Leptospira weilii strain SC295, Leptospira interrogans strain HP358) trapped in human leptospirosis outbreak premises. To evaluate the pathogenic potential of these isolates, we performed an in vivo and in silico virulence analysis. Environmental isolates and strain HP364 did not induce any clinical manifestations in hamsters. Strain SC295 caused inactivity and weight loss with histopathological changes in kidneys, however, all hamsters survived until the end of the experiment. Strain HP358 showed a high virulent phenotype as all infected hamsters died or were moribund within 7 days postinfection. Lungs, liver, and kidneys showed pathological changes with hemorrhage as the main presentation. In silico analysis elucidated the genome size of strain HP358 to be larger than strains HP364 and SC295 and containing virulence genes reported in Leptospira species and a high number of specific putative virulence factors. In conclusion, L. interrogans strain HP358 was highly pathogenic with fatal outcome. The constituent of Leptospira genomes may determine the level of disease severity and that needs further investigations.

Phylogenetic Analysis Reveals Source Attribution Patterns for Campylobacter spp. in Tennessee and Pennsylvania

Campylobacteriosis is the most common bacterial foodborne illness in the United States and is frequently associated with foods of animal origin. The goals of this study were to compare clinical and non-clinical Campylobacter populations from Tennessee (TN) and Pennsylvania (PA), use phylogenetic relatedness to assess source attribution patterns, and identify potential outbreak clusters. Campylobacter isolates studied (n = 3080) included TN clinical isolates collected and sequenced for routine surveillance, PA clinical isolates collected from patients at the University of Pennsylvania Health System facilities, and non-clinical isolates from both states for which sequencing reads were available on NCBI. Phylogenetic analyses were conducted to categorize isolates into species groups and determine the population structure of each species. Most isolates were C. jejuni (n = 2132, 69.2%) and C. coli (n = 921, 29.9%), while the remaining were C. lari (0.4%), C. upsaliensis (0.3%), and C. fetus (0.1%). The C. jejuni group consisted of three clades; most non-clinical isolates were of poultry (62.7%) or cattle (35.8%) origin, and 59.7 and 16.5% of clinical isolates were in subclades associated with poultry or cattle, respectively. The C. coli isolates grouped into two clades; most non-clinical isolates were from poultry (61.2%) or swine (29.0%) sources, and 74.5, 9.2, and 6.1% of clinical isolates were in subclades associated with poultry, cattle, or swine, respectively. Based on genomic similarity, we identified 42 C. jejuni and one C. coli potential outbreak clusters. The C. jejuni clusters contained 188 clinical isolates, 19.6% of the total C. jejuni clinical isolates, suggesting that a larger proportion of campylobacteriosis may be associated with outbreaks than previously determined.

Characterization of Basal Transcriptomes Identifies Potential Metabolic and Virulence-Associated Adaptations Among Diverse Nontyphoidal Salmonella enterica Serovars

The zoonotic pathogen Salmonella enterica includes >2,600 serovars, which differ in the range of hosts they infect and the severity of disease they cause. To further elucidate the mechanisms behind these differences, we performed transcriptomic comparisons of nontyphoidal Salmonella (NTS) serovars with the model for NTS pathogenesis, S. Typhimurium. Specifically, we used RNA-seq to characterize the understudied NTS serovars S. Javiana and S. Cerro, representing a serovar frequently attributed to human infection via contact with amphibians and reptiles, and a serovar primarily associated with cattle, respectively. Whole-genome sequence (WGS) data were utilized to ensure that strains characterized with RNA-seq were representative of their respective serovars. RNA extracted from representative strains of each serovar grown to late exponential phase in Luria-Bertani (LB) broth showed that transcript abundances of core genes were significantly higher (p<0.001) than those of accessory genes for all three serovars. Inter-serovar comparisons identified that transcript abundances of genes in Salmonella Pathogenicity Island (SPI) 1 were significantly higher in both S. Javiana and S. Typhimurium compared to S. Cerro. Together, our data highlight potential transcriptional mechanisms that may facilitate S. Cerro and S. Javiana survival in and adaptation to their respective hosts and impact their ability to cause disease in others. Furthermore, our analyses demonstrate the utility of omics approaches in advancing our understanding of the diversity of metabolic and virulence mechanisms of different NTS serovars.

Investigating the Distribution of Strains of Erwinia amylovora and Streptomycin Resistance in Apple Orchards in New York Using Clustered Regularly Interspaced Short Palindromic Repeat Profiles: A 6-Year Follow-Up

Fire blight, caused by the bacterium Erwinia amylovora, is one of the most important diseases of apple. The antibiotic streptomycin is routinely used in the commercial apple industries of New York (NY) and New England to manage the disease. In 2002 and again, from 2011 to 2014, outbreaks of streptomycin resistance (SmR) were reported and investigated in NY. Motivated by new grower reports of control failures, we conducted a follow-up investigation of the distribution of SmR and E. amylovora strains for major apple production regions of NY over the last 6 years (2015 to 2020). Characterization of clustered regularly interspaced short palindromic repeat (CRISPR) profiles revealed that a few "cosmopolitan" strains were widely prevalent across regions, whereas many other "resident" strains were confined to one location. In addition, we uncovered novel CRISPR profile diversity in all investigated regions. SmR E. amylovora was detected only in a small area spanning two counties from 2017 to 2020 and was always associated with one CRISPR profile (41:23:38), which matched the profile of SmR E. amylovora, discovered in 2002. This suggests the original SmR E. amylovora was never fully eradicated and went undetected because of several seasons of low disease pressure in this region. Investigation of several representative isolates under controlled greenhouse conditions indicated significant differences in aggressiveness on 'Gala' apples. Potential implications of strain differences include the propensity of strains to become distributed across wide geographic regions and associated resistance management practices. Results from this work will directly influence sustainable fire blight management recommendations for commercial apple industries in NY state and other regions.

Population differentiation of Rhodobacteraceae along with coral compartments

Coral mucus, tissue, and skeleton harbor compositionally different microbiota, but how these coral compartments shape the microbial evolution remains unexplored. Here, we sampled bacteria inhabiting a prevalent coral species Platygyra acuta and sequenced genomes of 234 isolates comprising two populations in Rhodobacteraceae, an alphaproteobacterial lineage representing a significant but variable proportion (5-50%) of the coral microbiota. The Ruegeria population (20 genomes) contains three clades represented by eight, six, and six isolates predominantly sampled from the skeleton (outgroup), mucus (clade-M), and skeleton (clade-S), respectively. The clade-M possesses functions involved in the utilization of coral osmolytes abundant in the mucus (e.g., methylamines, DMSP, taurine, and L-proline), whereas the clade-S uniquely harbors traits that may promote adaptation to the low-energy and diurnally anoxic skeleton (e.g., sulfur oxidation and swimming motility). These between-clade genetic differences were largely supported by physiological assays. Expanded analyses by including genomes of 24 related isolates (including seven new genomes) from other marine environments suggest that clade-M and clade-S may have diversified in non-coral habitats, but they also consolidated a key role of distinct coral compartments in diversifying many of the above-mentioned traits. The unassigned Rhodobacteraceae population (214 genomes) varies only at a few dozen nucleotide sites across the whole genomes, but the number of between-compartment migration events predicted by the Slatkin-Maddison test supported that dispersal limitation between coral compartments is another key mechanism diversifying microbial populations. Collectively, our results suggest that different coral compartments represent ecologically distinct and microgeographically separate habitats that drive the evolution of the coral microbiota.

Phylogeographic Clustering Suggests that Distinct Clades of Salmonella enterica Serovar Mississippi Are Endemic in Australia, the United Kingdom, and the United States

Salmonella enterica serovar Mississippi is the 2nd and 14th leading cause of human clinical salmonellosis in the Australian island state of Tasmania and the United States, respectively. Despite its public health relevance, relatively little is known about this serovar. Comparison of whole-genome sequence (WGS) data of S. Mississippi isolates with WGS data for 317 additional S. enterica serovars placed one clade of S. Mississippi within S. enterica clade B (“clade B Mississippi”) and the other within section Typhi in S. enterica clade A (“clade A Mississippi”), suggesting that these clades evolved from different ancestors. Phylogenetic analysis of 364 S. Mississippi isolates from Australia, the United Kingdom, and the United States suggested that the isolates cluster geographically, with U.S. and Australian isolates representing different subclades (Ai and Aii, respectively) within clade A Mississippi and clade B isolates representing the predominant S. Mississippi isolates in the United Kingdom. Intraclade comparisons suggested that different mobile elements, some of which encode virulence factors, are responsible for the observed differences in gene content among isolates within these clades. Specifically, genetic differences among clade A isolates reflect differences in prophage contents, while differences among clade B isolates are due to the acquisition of a 47.1-kb integrative conjugative element (ICE). Phylogenies inferred from antigenic components (fliC, fljB, and O-antigen-processing genes) support that clade A and B Mississippi isolates acquired these loci from different ancestral serovars. Overall, these data support that different S. Mississippi phylogenetic clades are endemic in Australia, the United Kingdom, and the United States.

IMPORTANCE The number of known so-called “polyphyletic” serovars (i.e., phylogenetically distinct clades with the same O and H antigenic formulas) continues to increase as additional Salmonella isolates are sequenced. While serotyping remains a valuable tool for reporting and monitoring Salmonella, more discriminatory analyses for classifying polyphyletic serovars may improve surveillance efforts for these serovars, as we found that for S. Mississippi, distinct genotypes predominate at different geographic locations. Our results suggest that the acquisition of genes encoding O and H antigens from different ancestors led to the emergence of two Mississippi clades. Furthermore, our results suggest that different mobile elements contribute to the microevolution and diversification of isolates within these two clades, which has implications for the acquisition of novel adaptations, such as virulence factors.