Novel Temperate Phages of Salmonella enterica subsp. salamae and subsp. diarizonae and Their Activity against Pathogenic S. enterica subsp. enterica Isolates

Genetic characterization of four bacteriophages of Salmonella enterica derived from different geographic regions in China via genomic comparison

Based on the AT content > GC content in four Salmonella enterica lytic bacteriophage genomes, information entropy analysis revealed that overall nucleotide usage bias is shaped in the gene population. This genetic feature directly contributes to synonymous codons tending toward the A/T end rather than the C/G end. Furthermore, the interplay between the nucleotide composition constraint from the bacteriophage itself and the natural selection caused by outside environments forces our bacteriophages into similar evolutionary trends in terms of overall codon usage patterns. We identified the nucleotide composition constraint which plays an important role in shaping synonymous codon usage patterns including the keto skew at the first codon position, the pyrimidine skew at the second position and the AT skew at the third position. Although the four bacteriophages were isolated from different geographical regions in China, they display similar evolutionary trends in terms of genomic organization and synonymous codon usage, which are strongly influenced by the nucleotide composition constraint of the bacteriophage. The findings of the present study reveal important details of the evolutionary and host-pathogen interactions of Salmonella enterica, which will benefit the efficient utilization of phages for therapeutic and other applications.

Next-Generation Sequencing-Based Monitoring of Intestinal Bacteria and Bacteriophages Following Fecal Microbiota Transplantation in Inflammatory Bowel Diseases

Inflammatory bowel diseases (IBD) and acute graft-versus-host disease (GVHD) are associated with persistent intestinal dysfunction preceded by gut bacterial dysbiosis. There are limited data on intestinal bacteriophages in these conditions. The aim of the present work was to detect associations between dominant intestinal bacteria by means of 16S rRNA gene sequencing, and some clinically significant viruses detected with a customized primer panel for NGS-based study. The clinical group included patients with Crohn's disease (IBD, n = 9), or GVHD (n = 6) subjected to fecal microbiota transplantation (FMT) from healthy donors. The stool specimens were taken initially, and 5 times post-FMT until day 120. Using NGS approach, we have found a higher abundance of Proteobacterota phylum in GVHD, especially, at later terms post-FMT. Moreover, we found an early increase of Klebsiella and E. coli/Shigella abundance in GVHD, along with decreased relative content of Faecalibacterium. Upon evaluation of intestinal phageome, the relative amount of Caudoviricetes class was higher in GVHD. A significant correlation was found between Proteobacteria and Caudoviricetes, thus suggesting their association during the post-FMT period. Moreover, the relative amounts of five Caudoviricetes phage species showed distinct correlations with Klebsiella and Enterococcus ratios at different terms of FMT. In conclusion, parallel use of 16S rRNA gene sequencing and targeted NGS viral panel is a feasible and useful option for tracing specific viral strains in fecal microbiota. The developed array of viral primers may be extended to detect other phages infecting the clinically relevant bacteria.

Complete Genome Sequence of Bacteriophage EO1, Which Infects Both Escherichia coli O157:H7 and Shigella

The lytic bacteriophage EO1 has been newly isolated. This phage infects Escherichia coli O157:H7 and has a broad antibacterial spectrum, including against Shigella. The complete genome sequence of phage EO1 was determined; its full length is 166,941 bp, and it has a G+C content of 35.46%.

Tell me if you prefer bovine or poultry sectors and I'll tell you who you are: Characterization of Salmonella enterica subsp. enterica serovar Mbandaka in France

Introduction

In north-western France, Salmonella enterica susp. enterica serovar Mbandaka (S. Mbandaka) is most frequently isolated from bovine and dairy samples. While this serovar most often results in asymptomatic carriage, for a number of years it has caused episodes of abortions, which have serious economic consequences for the sector. Interestingly, this serovar is also isolated from Gallus gallus in the same geographic zone. Despite its prevalence in bovines in north-western France, S. Mbandaka has not been broadly studied at the genomic level, and its prevalence and host adaptation are still not fully understood.

Methods

In this study, we analyzed the genomic diversity of 304 strains of S. Mbandaka isolated from the bovine and poultry sectors in this area over a period of 5 years. A phylogenetic analysis was carried out and two approaches were followed to identify conserved genes and mutations related to host associations. The first approach targeted the genes compiled in the MEGARESv2, Resfinder, VFDB and SPI databases. Plasmid and phage contents were also investigated. The second approach refers to an in-house algorithm developed for this study that computes sensitivity, specificity, and accuracy of accessory genes and core variants according to predefined genomes groups.

Results and discussion

All the analyzed strains belong to the multi-locus sequence type profile ST413, and the phylogenomic analysis revealed main clustering by host (bovine and poultry), emphasizing the circulation of 12 different major clones, of which seven circulate in poultry and five in the bovine sector in France and a likely food production chain adaptation of these clones. All strains present resistance determinants including heavy metals and biocides that could explain the ability of this serovar to survive and persist in the environment, within herds, and in food processing plants. To explore the wild animal contribution to the spread of this serovar in north-western France, we retrieved S. Mbandaka genomes isolated from wild birds from EnteroBase and included them in the phylogenomic analysis together with our collection. Lastly, screening of accessory genes and major variants allowed us to identify conserved specific mutations characteristic of each major cluster. These mutations could be used to design useful probes for food safety surveillance.

Genomic characteristics and comparative genomics analysis of Salmonella enterica subsp. enterica serovar Thompson isolated from an outbreak in South Korea

Salmonella infections represent an important public health problem. In 2018, a multistate outbreak of S. enterica subsp. enterica serovar Thompson infection associated with contaminated chocolate cakes in schools was reported in South Korea. In this study, we sequenced the 37 S. Thompson strains isolated from chocolate cakes, egg whites, preserves, and cookware associated with the outbreak. In addition, we analyze the genomic sequences of 61 S. Thompson strains (37 chocolate cake-related outbreak strains, 4 strains isolated from outbreaks in South Korea and 20 strains available in the National Center for Biotechnology Information) to assess the genomic characteristics of outbreak-related strains by comparative genomics and phylogenetic analysis. The results showed that identically classified clusters divided strains into two clusters, sub-clusters A & I (with strains from 2018 in South Korea) and sub-clusters B & II (with strains from 2014 to 2015 in South Korea). S. Thompson isolated from South Korea were accurately distinguished from publicly-available strains. Unlike other S. Thompson genomes, those of chocolate cake outbreak-related strains had three Salmonella phages (SEN8, vB SosS Oslo, and SI7) integrated into their chromosome. Comparative genomics revealed several genes responsible for the specific genomic features of chocolate cake outbreak-related strains and three bacteriophages that may contribute to the pathogenicity of other S. Thompson strains.

The effect of a temperature-sensitive prophage on the evolution of virulence in an opportunistic bacterial pathogen

Viruses are key actors of ecosystems and have major impacts on global biogeochemical cycles. Prophages deserve particular attention as they are ubiquitous in bacterial genomes and can enter a lytic cycle when triggered by environmental conditions. We explored how temperature affects the interactions between prophages and other biological levels using an opportunistic pathogen, the bacterium Serratia marcescens, which harbours several prophages and that had undergone an evolution experiment under several temperature regimes. We found that the release of one of the prophages was temperature-sensitive and malleable to evolutionary changes. We further discovered that the virulence of the bacterium in an insect model also evolved and was positively correlated with phage release rates. We determined through analysis of genetic and epigenetic data that changes in the bacterial outer cell wall structure possibly explain this phenomenon. We hypothezise that the temperature-dependent phage release rate acted as a selection pressure on S. marcescens and that it resulted in modified bacterial virulence in the insect host. Our study system illustrates how viruses can mediate the influence of abiotic environmental changes to other biological levels and thus be involved in ecosystem feedback loops.

Genomic Analysis of a mcr-9.1-Harbouring IncHI2-ST1 Plasmid from Enterobacter ludwigii Isolated in Fish Farming

This study analyzed the resistome, virulome and mobilome of an MCR-9-producing Enterobacter sp. identified in a muscle sample of seabream (Sparus aurata), collected in a land tank from multitrophic fish farming production. Average Nucleotide Identity analysis identified INSAq77 at the species level as an Enterobacter ludwigii INSAq77 strain that was resistant to chloramphenicol, florfenicol and fosfomycin and was susceptible to all other antibiotics tested. In silico antimicrobial resistance analyses revealed genes conferring in silico resistance to β-lactams (bla_ACT-88), chloramphenicol (catA4-type), fosfomycin (fosA2-type) and colistin (mcr-9.1), as well as several efflux pumps (e.g., oqxAB-type and mar operon). Further bioinformatics analysis revealed five plasmid replicon types, including the IncHI2/HI2A, which are linked to the worldwide dissemination of the mcr-9 gene in different antibiotic resistance reservoirs. The conserved nickel/copper operon rcnR-rcnA-pcoE-ISSgsp1-pcoS-IS903-mcr-9-wbuC was present, which may play a key role in copper tolerance under anaerobic growth and nickel homeostasis. These results highlight that antibiotic resistance in aquaculture are spreading through food, the environment and humans, which places this research in a One Health context. In fact, colistin is used as a last resort for the treatment of serious infections in clinical settings, thus mcr genes may represent a serious threat to human health.

Phenotypic characterization and genome analysis of a novel Salmonella Typhimurium phage having unique tail fiber genes

Salmonella enterica serovar Typhimurium is a foodborne pathogen causing occasional outbreaks of enteric infections in humans. Salmonella has one of the largest pools of temperate phages in its genome that possess evolutionary significance for pathogen. In this study, we characterized a novel temperate phage Salmonella phage BIS20 (BIS20) with unique tail fiber genes. It belongs to the subfamily Peduovirinae genus Eganvirus and infects Salmonella Typhimurium strain (SE-BS17; Acc. NO MZ503545) of poultry origin. Phage BIS20 was viable only at biological pH and temperature ranges (pH7 and 37 °C). Despite being temperate BIS20 significantly slowed down the growth of host strain for 24 h as compared to control (P < 0.009). Phage BIS20 features 29,477-base pair (bp) linear DNA genome with 53% GC content and encodes for 37 putative ORFs. These ORFs have mosaic arrangement as indicated by its ORF similarity to various phages and prophages in NCBI. Genome analysis indicates its similarity to Salmonella enterica serovar Senftenberg prophage (SEStP) sequence (Nucleotide similarity 87.7%) and Escherichia virus 186 (~ 82.4% nucleotide similarity). Capsid genes were conserved however those associated with tail fiber formation and assembly were unique to all members of genus Eganvirus. We found strong evidence of recombination hotspot in tail fiber gene. Our study identifies BIS20 as a new species of genus Eganvirus temperate phages as its maximum nucleotide similarity is 82.4% with any phage in NCBI. Our findings may contribute to understanding of origin of new temperate phages.

Isolation and characterization of phage (ETP-1) specific to multidrug resistant pathogenic Edwardsiella tarda and its in vivo biocontrol efficacy in zebrafish (Danio rerio)

Edwardsiella tarda phage (ETP-1) was isolated from marine fish farm water to characterize its effect against pathogenic multidrug-resistant E. tarda. According to transmission electron microscopy results, ETP-1 is classified as a member of the family Podoviridae. ETP-1 showed MOI dependent E. tarda growth inhibition, a latent period of 60 min, and burst size of 100 PFU per infected cells. In host range tests, five out of eight E. tarda strains were sensitive to ETP-1 which had efficiency of plating index in the range 1-1.28. ETP-1 was stable over a broad range of pH and temperature. The size of the ETP-1 genome was predicted to be approximately 40 kb. Zebrafish exposed to ETP-1 showed no adverse gene responses to the inflammatory mediator cytokines, il1-β, tnf-α, il-6, and il-10, the chemokine, cxcl-8a, and reactive oxygen species, sod-1. When zebrafish were bath exposed to ETP-1 for 12 days and simultaneously challenged with E. tarda (1.08 × 10⁵ CFU fish^-1), the survival rate was higher in phage exposed fish (68%) compared to that of the control (18%) until 4 days post challenge. Our results suggest that ETP-1 can be used as a potential bio-therapeutic candidate to control multi-drug resistant E. tarda infection in aquaculture.

Diversity and Host Specificity Revealed by Biological Characterization and Whole Genome Sequencing of Bacteriophages Infecting Salmonella enterica

Phages infecting members of the opportunistic human pathogen, Salmonella enterica, are widespread in natural environments and offer a potential source of agents that could be used for controlling populations of this bacterium; yet, relatively little is known about these phages. Here we describe the isolation and characterization of 45 phages of Salmonella enterica from disparate geographic locations within British Columbia, Canada. Host-range profiling revealed host-specific patterns of susceptibility and resistance, with several phages identified that have a broad-host range (i.e., able to lyse >40% of bacterial hosts tested). One phage in particular, SE13, is able to lyse 51 out of the 61 Salmonella strains tested. Comparative genomic analyses also revealed an abundance of sequence diversity in the sequenced phages. Alignment of the genomes grouped the phages into 12 clusters with three singletons. Phages within certain clusters exhibited extraordinarily high genome homology (>98% nucleotide identity), yet between clusters, genomes exhibited a span of diversity (<50% nucleotide identity). Alignment of the major capsid protein also supported the clustering pattern observed with alignment of the whole genomes. We further observed associations between genomic relatedness and the site of isolation, as well as genetic elements related to DNA metabolism and host virulence. Our data support the knowledge framework for phage diversity and phage-host interactions that are required for developing phage-based applications for various sectors, including biocontrol, detection and typing.

Genome sequencing and analysis of Salmonella enterica subsp. enterica serovar Stanley UPM 517: Insights on its virulence-associated elements and their potentials as vaccine candidates

Salmonella enterica subsp. enterica serovar Stanley (S. Stanley) is a pathogen that contaminates food, and is related to Salmonella outbreaks in a variety of hosts such as humans and farm animals through products like dairy items and vegetables. Despite the fact that several vaccines of Salmonella strains had been constructed, none of them were developed according to serovar Stanley up to this day. This study presents results of genome sequencing and analysis on our S. Stanley UPM 517 strain taken from fecal swabs of 21-day-old healthy commercial chickens in Perak, Malaysia and used Salmonella enterica subsp. enterica serovar Typhimurium LT2 (S. Typhimurium LT2) as a reference to be compared with. First, sequencing and assembling of the Salmonella Stanley UPM 517 genome into a contiguous form were done. The work was then continued with scaffolding and gap filling. Annotation and alignment of the draft genome was performed with S. Typhimurium LT2. The other elements of virulence estimated in this study included Salmonella pathogenicity islands, resistance genes, prophages, virulence factors, plasmid regions, restriction-modification sites and the CRISPR-Cas system. The S. Stanley UPM 517 draft genome had a length of 4,736,817 bp with 4,730 coding sequence and 58 RNAs. It was discovered via genomic analysis on this strain that there were antimicrobial resistance properties toward a wide variety of antibiotics. Tcf and ste, the two fimbrial virulence clusters related with human and broiler intestinal colonizations which were not found in S. Typhimurium LT2, were atypically discovered in the S. Stanley UPM 517 genome. These clusters are involved in the intestinal colonization of human and broilers, respectively. There were seven Salmonella pathogenicity islands (SPIs) within the draft genome, which contained the virulence factors associated with Salmonella infection (except SPI-14). Five intact prophage regions, mostly comprising of the protein encoding Gifsy-1, Fels-1, RE-2010 and SEN34 prophages, were also encoded in the draft genome. Also identified were Type I–III restriction-modification sites and the CRISPR-Cas system of the Type I–E subtype. As this strain exhibited resistance toward numerous antibiotics, we distinguished several genes that had the potential for removal in the construction of a possible vaccine candidate to restrain and lessen the pervasiveness of salmonellosis and to function as an alternative to antibiotics.

Complete Genome Sequence of Bacteriophage SEN8, a Temperate Phage Isolated from Salmonella enterica subsp. salamae

A temperate phage, SEN8, having a broad activity against pathogenic Salmonella serovars, was isolated from Salmonella enterica subsp. salamae strain Sen8. The complete genome sequence of phage SEN8 was determined (35,203 bp) and showed relatedness to P2-like phages (Salmonella phages Fels-2 and RE-2010).