Comparative genomic analysis of 142 bacteriophages infecting Salmonella enterica subsp. enterica

What's on a prophage: analysis of Salmonella spp. prophages identifies a diverse range of cargo with multiple virulence- and metabolism-associated functions

The gain of mobile elements, such as prophages, can introduce cargo to the recipient bacterium that could facilitate its persistence in or expansion to a new environment, such as a host. While previous studies have focused on identifying and characterizing the genetic diversity of prophages, analyses characterizing the cargo that prophages carry have not been extensively explored. We characterized prophage regions from 303 Salmonella spp. genomes (representing 254 unique serovars) to assess the distribution of prophages in diverse Salmonella. On average, prophages accounted for 3.7% (0.1%-8.8%) of the total genomic content of each isolate. Prophage regions annotated as Gifsy 1 and Salmon Fels 1 were the most commonly identified intact prophages, suggesting that they are common throughout the Salmonella genus. Among 21,687 total coding sequences (CDSs) from intact prophage regions in subsp. enterica genomes, 7.5% (median; range: 1.1%-47.6%) were categorized as having a function not related to prophage integration or phage structure, some of which could potentially provide a functional attribute to the host Salmonella cell. These predicted functions could be broadly categorized into CDSs involved in: (i) modification of cell surface structures (i.e., glycosyltransferases); (ii) modulation of host responses (e.g., SodC/SodA, SopE, ArtAB, and typhoid toxin); (iii) conferring resistance to heavy metals and antimicrobials; (iv) metabolism of carbohydrates, amino acids, and nucleotides; and (v) DNA replication, repair, and regulation. Overall, our systematic analysis of prophage cargo highlights a broader role for prophage cargo in influencing the metabolic, virulence, and resistance characteristics of Salmonella.

IMPORTANCE

Lysogenic bacteriophages (phages) can integrate their genome into a bacterial host's genome, potentially introducing genetic elements that can affect the fitness of the host bacterium. The functions of prophage-encoded genes are important to understand as these genes could be mobilized and transferred to a new host. Using a large genomic dataset representing >300 isolates from all known subspecies and species of Salmonella, our study contributes important new findings on the distribution of prophages and the types of cargo that diverse Salmonella prophages carry. We identified a number of coding sequences (CDSs) annotated as having cell surface-modifying attributes, suggesting that prophages may have played an important role in shaping Salmonella's diverse surface antigen repertoire. Furthermore, our characterization of prophages suggests that they play a broader role in facilitating the acquisition and transfer of CDSs associated with metabolism, DNA replication and repair, virulence factors, and to a lesser extent, antimicrobial resistance.

Comparative genomics reveals high genetic similarity among strains of Salmonella enterica serovar Infantis isolated from multiple sources in Brazil

Background

Salmonella enterica serovar Infantis (Salmonella Infantis) is a zoonotic, ubiquitous and foodborne pathogen of worldwide distribution. Despite Brazil's relevance as a major meat exporter, few studies were conducted to characterize strains of this serovar by genomic analyses in this country. Therefore, this study aimed to assess the diversity of 80 Salmonella Infantis strains isolated from veterinary, food and human sources in Brazil between 2013 and 2018 by comparative genomic analyses. Additional genomes of non-Brazilian countries (n = 18) were included for comparison purposes in some analyses.

Methods

Analyses of whole-genome multi-locus sequence typing (wgMLST), using PGAdb-builder, and of fragmented genomes, using Gegenees, were conducted to compare the 80 Brazilian strains to the 18 non-Brazilian genomes. Pangenome analyses and calculations were performed for all Salmonella Infantis genomes analyzed. The presence of prophages was determined using PHASTER for the 80 Brazilian strains. The genome plasticity using BLAST Ring Image Generator (BRIG) and gene synteny using Mauve were evaluated for 20 selected Salmonella Infantis genomes from Brazil and ten from non-Brazilian countries. Unique orthologous protein clusters were searched in ten selected Salmonella Infantis genomes from Brazil and ten from non-Brazilian countries.

Results

wgMLST and Gegenees showed a high genomic similarity among some Brazilian Salmonella Infantis genomes, and also the correlation of some clusters with non-Brazilian genomes. Gegenees also showed an overall similarity >91% among all Salmonella Infantis genomes. Pangenome calculations revealed an open pangenome for all Salmonella Infantis subsets analyzed and a high gene content in the core genomes. Fifteen types of prophages were detected among 97.5% of the Brazilian strains. BRIG and Mauve demonstrated a high structural similarity among the Brazilian and non-Brazilian isolates. Unique orthologous protein clusters related to biological processes, molecular functions, and cellular components were detected among Brazilian and non-Brazilian genomes.

Conclusion

The results presented using different genomic approaches emphasized the significant genomic similarity among Brazilian Salmonella Infantis genomes analyzed, suggesting wide distribution of closely related genotypes among diverse sources in Brazil. The data generated contributed to novel information regarding the genomic diversity of Brazilian and non-Brazilian Salmonella Infantis in comparison. The different genetically related subtypes of Salmonella Infantis from Brazil can either occur exclusively within the country, or also in other countries, suggesting that some exportation of the Brazilian genotypes may have already occurred.

Bacteriophages in nature: recent advances in research tools and diverse environmental and biotechnological applications

Bacteriophages infect and replicate within bacteria and play a key role in the environment, particularly in microbial ecosystems and bacterial population dynamics. The increasing recognition of their significance stems from their wide array of environmental and biotechnological uses, which encompass the mounting issue of antimicrobial resistance (AMR). Beyond their therapeutic potential in combating antibiotic-resistant infections, bacteriophages also find vast applications such as water quality monitoring, bioremediation, and nutrient cycling within environmental sciences. Researchers are actively involved in isolating and characterizing bacteriophages from different natural sources to explore their applications. Gaining insights into key aspects such as the life cycle of bacteriophages, their host range, immune interactions, and physical stability is vital to enhance their application potential. The establishment of diverse phage libraries has become indispensable to facilitate their wide-ranging uses. Consequently, numerous protocols, ranging from traditional to cutting-edge techniques, have been developed for the isolation, detection, purification, and characterization of bacteriophages from diverse environmental sources. This review offers an exploration of tools, delves into the methods of isolation, characterization, and the extensive environmental applications of bacteriophages, particularly in areas like water quality assessment, the food sector, therapeutic interventions, and the phage therapy in various infections and diseases.

Comparative Analyses of Bacteriophage Genomes

Bacterial viruses (bacteriophages or phages) are the most abundant and diverse biological entities on Earth. There is a renewed worldwide interest in phage-centered research motivated by their enormous potential as antimicrobials to cope with multidrug-resistant pathogens. An ever-growing number of complete phage genomes are becoming available, derived either from newly isolated phages (cultivated phages) or recovered from metagenomic sequencing data (uncultivated phages). Robust comparative analysis is crucial for a comprehensive understanding of genotypic variations of phages and their related evolutionary processes, and to investigate the interaction mechanisms between phages and their hosts. In this chapter, we present a protocol for phage comparative genomics employing tools selected out of the many currently available, focusing on complete genomes of phages classified in the class Caudoviricetes. This protocol provides accurate identification of similarities, differences, and patterns among new and previously known complete phage genomes as well as phage clustering and taxonomic classification.

Viromes vs. mixed community metagenomes: choice of method dictates interpretation of viral community ecology

Background

Viruses, the majority of which are uncultivated, are among the most abundant biological entities on Earth. From altering microbial physiology to driving community dynamics, viruses are fundamental members of microbiomes. While the number of studies leveraging viral metagenomics (viromics) for studying uncultivated viruses is growing, standards for viromics research are lacking. Viromics can utilize computational discovery of viruses from total metagenomes of all community members (hereafter metagenomes) or use physical separation of virus-specific fractions (hereafter viromes). However, differences in the recovery and interpretation of viruses from metagenomes and viromes obtained from the same samples remain understudied.

Results

Here, we compare viral communities from paired viromes and metagenomes obtained from 60 diverse samples across human gut, soil, freshwater, and marine ecosystems. Overall, viral communities obtained from viromes were more abundant and species rich than those obtained from metagenomes, although there were some exceptions. Despite this, metagenomes still contained many viral genomes not detected in viromes. We also found notable differences in the predicted lytic state of viruses detected in viromes vs metagenomes at the time of sequencing. Other forms of variation observed include genome presence/absence, genome quality, and encoded protein content between viromes and metagenomes, but the magnitude of these differences varied by environment.

Conclusions

Overall, our results show that the choice of method can lead to differing interpretations of viral community ecology. We suggest that the choice of whether to target a metagenome or virome to study viral communities should be dependent on the environmental context and ecological questions being asked. However, our overall recommendation to researchers investigating viral ecology and evolution is to pair both approaches to maximize their respective benefits.

Evaluation of metagenomic assembly methods for the detection and characterization of antimicrobial resistance determinants and associated mobilizable elements

Antimicrobial resistance genes (ARGs) can be transferred between members of a bacterial population by mobile genetic elements (MGE). Understanding the risk of these transfer events is important in monitoring and predicting antimicrobial resistance (AMR), especially in the context of a One Health Continuum. However, there is no universally accepted method for detection of ARGs and MGEs, and especially for determining their linkages. This study used publicly available shotgun metagenomic DNA short-read (Illumina, 100 bp paired-end) sequence data from samples across the One Health Continuum (including beef cattle composite feces from feedlots, catch basin water at feedlots, agricultural soil from feedlot manured surrounding fields, and urban/municipal sewage influent from two municipal wastewater treatment plants) to develop a workflow to identify and associate ARGs and MGEs. ARG- and MGE-based targeted-assemblies with available short-read data were unable to meet this analysis goal. In contrast, de novo assembly of contigs provided enough sequence context to associate ARGs and MGEs, without compromising discovery rate. However, to estimate the relative abundance of these elements, unassembled sequence data must still be used.

The rare Salmonella enterica serovar Isangi: genomic characterization of the antimicrobial resistance, virulence potential and epidemiology of Brazilian strains in comparison to global isolates

Introduction. Salmonella enterica serovar Isangi (S. Isangi) is a rare non-typhoidal serovar, related to invasive nosocomial infections in various countries and to increasing antimicrobial resistance rates.Gap statement. Despite existing reports on S. Isangi, there is a lack of information of specific traits regarding this serovar, which could be improved through genomic analyses.Aim. Our goals were to characterize the antimicrobial resistance, virulence potential and genomic relatedness of 11 S. Isangi strains from Brazil in comparison to 185 genomes of global isolates using whole-genome sequencing (WGS) data.Methodology. Phenotypic resistance was determined by disc-diffusion. The search for resistance genes, plasmids, prophages, Salmonella pathogenicity islands (SPIs) and virulence genes, plus multi-locus sequence typing (MLST) and core-genome MLST (cgMLST) were performed using WGS.Results. Brazilian S. Isangi strains showed phenotypic resistance to nalidixic acid, ciprofloxacin and streptomycin, and harboured antimicrobial resistance [qnrB19, aac(6')-Iaa, mdsAB] and heavy metal tolerance (arsD, golST) genes. Col(pHAD28) and IncFII(S) plasmids, virulence genes related to adherence, macrophage induction, magnesium uptake, regulation and type III secretion systems, 12 SPIs and eight prophages were detected. The 185 additional global genomes analysed harboured resistance genes against 11 classes of antimicrobial compounds, 22 types of plasmids, 32 prophages, 14 SPIs, and additional virulence genes related to serum resistance, stress adaptation and toxins. Sequence type (ST)216 was assigned to genomes from Brazil and other countries, while ST335 was the most frequent ST, especially among South African genomes. cgMLST showed that Brazilian genomes were more closely related to genomes from European and African countries, the USA and Taiwan, while the majority of South African genomes were more closely related among each other.Conclusion. The presence of S. Isangi strains from Brazil and different countries showing a close genomic correlation, antimicrobial resistance profiles to drugs used in human therapy and a large number of virulence determinants reinforced the need for stronger initiatives to monitor rare non-typhoidal Salmonella serovars such as S. Isangi in order to prevent its dissemination among human and non-human sources.

Application of prophage sequence analysis to investigate a disease outbreak involving Salmonella Adjame, a rare serovar and implications for the population structure

Introduction

Outbreak investigation of foodborne salmonellosis is hindered when the food source is contaminated by multiple strains of Salmonella, creating difficulties matching an incriminated organism recovered from patients with the specific strain in the suspect food. An outbreak of the rare Salmonella Adjame was caused by multiple strains of the organism as revealed by single-nucleotide polymorphism (SNP) variation. The use of highly discriminatory prophage analysis to characterize strains of Salmonella should enable a more precise strain characterization and aid the investigation of foodborne salmonellosis.

Methods

We have carried out genomic analysis of S. Adjame strains recovered during the course of a recent outbreak and compared them with other strains of the organism (n = 38 strains), using SNPs to evaluate strain differences present in the core genome, and prophage sequence typing (PST) to evaluate the accessory genome. Phylogenetic analyses were performed using both total prophage content and conserved prophages.

Results

The PST analysis of the S. Adjame isolates showed a high degree of strain heterogeneity. We observed small clusters made up of 2-6 isolates (n = 27) and singletons (n = 11) in stark contrast with the three clusters observed by SNP analysis. In total, we detected 24 prophages of which only four were highly prevalent, namely: Entero_p88 (36/38 strains), Salmon_SEN34 (35/38 strains), Burkho_phiE255 (33/38 strains) and Edward_GF (28/38 strains). Despite the marked strain diversity seen with prophage analysis, the distribution of the four most common prophages matched the clustering observed using core genome.

Discussion

Mutations in the core and accessory genomes of S. Adjame have shed light on the evolutionary relationships among the Adjame strains and demonstrated a convergence of the variations observed in both fractions of the genome. We conclude that core and accessory genomes analyses should be adopted in foodborne bacteria outbreak investigations to provide a more accurate strain description and facilitate reliable matching of isolates from patients and incriminated food sources. The outcomes should translate to a better understanding of the microbial population structure and an 46 improved source attribution in foodborne illnesses.

Genomic diversity and comprehensive taxonomical classification of 61 Bacillus subtilis group member infecting bacteriophages, and the identification of ortholog taxonomic signature genes

BACKGROUND

Despite the applications of Bacillus subtilis group species in various sectors, limited information is available regarding their phages. Here, 61 B. subtilis group species-infecting phages (BSPs) were studied for their taxonomic classification considering the genome-size, genomic diversity, and the host, followed by the identification of orthologs taxonomic signature genes.

RESULTS

BSPs have widely ranging genome sizes that can be bunched into groups to demonstrate correlations to family and subfamily classifications. Comparative analysis re-confirmed the existing, BSPs-containing 14 genera and 21 species and displayed inter-genera similarities within existing subfamilies. Importantly, it also revealed the need for the creation of new taxonomic classifications, including 28 species, nine genera, and two subfamilies (New subfamily1 and New subfamily2) to accommodate inter-genera relatedness. Following pangenome analysis, no ortholog shared by all BSPs was identified, while orthologs, namely, the tail fibers/spike proteins and poly-gamma-glutamate hydrolase, that are shared by more than two-thirds of the BSPs were identified. More importantly, major capsid protein (MCP) type I, MCP type II, MCP type III and peptidoglycan binding proteins that are distinctive orthologs for Herelleviridae, Salasmaviridae, New subfamily1, and New subfamily2, respectively, were identified and analyzed which could serve as signatures to distinguish BSP members of the respective taxon.

CONCLUSIONS

In this study, we show the genomic diversity and propose a comprehensive classification of 61 BSPs, including the proposition for the creation of two new subfamilies, followed by the identification of orthologs taxonomic signature genes, potentially contributing to phage taxonomy.

Prophage-mediated genome differentiation of the Salmonella Derby ST71 population

Although Salmonella Derby ST71 strains have been recognized as poultry-specific by previous studies, multiple swine-associated S. Derby ST71 strains were identified in this long-term, multi-site epidemic study. Here, 15 representative swine-associated S. Derby ST71 strains were sequenced and compared with 65 (one swine-associated and 64 poultry-associated) S. Derby ST71 strains available in the NCBI database at a pangenomic level through comparative genomics analysis to identify genomic features related to the differentiation of swine-associated strains and previously reported poultry-associated strains. The distribution patterns of known Salmonella pathogenicity islands (SPIs) and virulence factor (VF) encoding genes were not capable of differentiating between the two strain groups. The results demonstrated that the S. Derby ST71 population harbours an open pan-genome, and swine-associated ST71 strains contain many more genes than the poultry-associated strains, mainly attributed to the prophage sequence contents in the genomes. The numbers of prophage sequences identified in the swine-associated strains were higher than those in the poultry-associated strains. Prophages specifically harboured by the swine-associated strains were found to contain genes that facilitate niche adaptation for the bacterial hosts. Gene deletion experiments revealed that the dam gene specifically present in the prophage of the swine-associated strains is important for S. Derby to adhere onto the host cells. This study provides novel insights into the roles of prophages during the genome differentiation of Salmonella.

Screening of Potential Vibrio cholerae Bacteriophages for Cholera Therapy: A Comparative Genomic Approach

Cholera continues to be a major burden for developing nations, especially where sanitation, quality of water supply, and hospitalization have remained an issue. Recently, growing antimicrobial-resistant strains of Vibrio cholerae underscores alternative therapeutic strategies for cholera. Bacteriophage therapy is considered one of the best alternatives for antibiotic treatment. For the identification of potential therapeutic phages for cholera, we have introduced a comprehensive comparative analysis of whole-genome sequences of 86 Vibrio cholerae phages. We have witnessed extensive variation in genome size (ranging from 33 to 148 kbp), GC (G + C) content (varies from 34.5 to 50.8%), and the number of proteins (ranging from 15 to 232). We have identified nine clusters and three singletons using BLASTn, confirmed by nucleotide dot plot and sequence identity. A high degree of sequence and functional similarities in both the genomic and proteomic levels have been observed within the clusters. Evolutionary analysis confirms that phages are conserved within the clusters but diverse between the clusters. For each therapeutic phage, the top 2 closest phages have been identified using a system biology approach and proposed as potential therapeutic phages for cholera. This method can be applied for the classification of the newly isolated Vibrio cholerae phage. Furthermore, this systematic approach might be useful as a model for screening potential therapeutic phages for other bacterial diseases.

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.

Phenotypic and Genomic Comparison of Klebsiella pneumoniae Lytic Phages: vB_KpnM-VAC66 and vB_KpnM-VAC13

Klebsiella pneumoniae is a human pathogen that worsens the prognosis of many immunocompromised patients. Here, we annotated and compared the genomes of two lytic phages that infect clinical strains of K. pneumoniae (vB_KpnM-VAC13 and vB_KpnM-VAC66) and phenotypically characterized vB_KpnM-VAC66 (time of adsorption of 12 min, burst size of 31.49 ± 0.61 PFU/infected cell, and a host range of 20.8% of the tested strains). Transmission electronic microscopy showed that vB_KpnM-VAC66 belongs to the Myoviridae family. The genomic analysis of the phage vB_KpnM-VAC66 revealed that its genome encoded 289 proteins. When compared to the genome of vB_KpnM-VAC13, they showed a nucleotide similarity of 97.56%, with a 93% of query cover, and the phylogenetic study performed with other Tevenvirinae phages showed a close common ancestor. However, there were 21 coding sequences which differed. Interestingly, the main differences were that vB_KpnM-VAC66 encoded 10 more homing endonucleases than vB_KpnM-VAC13, and that the nucleotidic and amino-acid sequences of the L-shaped tail fiber protein were highly dissimilar, leading to different three-dimensional protein predictions. Both phages differed significantly in their host range. These viruses may be useful in the development of alternative therapies to antibiotics or as a co-therapy increasing its antimicrobial potential, especially when addressing multidrug resistant (MDR) pathogens.

The genetic basis of phage susceptibility, cross-resistance and host-range in Salmonella

Though bacteriophages (phages) are known to play a crucial role in bacterial fitness and virulence, our knowledge about the genetic basis of their interaction, cross-resistance and host-range is sparse. Here, we employed genome-wide screens in Salmonella enterica serovar Typhimurium to discover host determinants involved in resistance to eleven diverse lytic phages including four new phages isolated from a therapeutic phage cocktail. We uncovered 301 diverse host factors essential in phage infection, many of which are shared between multiple phages demonstrating potential cross-resistance mechanisms. We validate many of these novel findings and uncover the intricate interplay between RpoS, the virulence-associated general stress response sigma factor and RpoN, the nitrogen starvation sigma factor in phage cross-resistance. Finally, the infectivity pattern of eleven phages across a panel of 23 genome sequenced Salmonella strains indicates that additional constraints and interactions beyond the host factors uncovered here define the phage host range.

Isolation and Characterization of Bacteriophage ZCSE6 against Salmonella spp.: Phage Application in Milk

Food safety is very important in the food industry as most pathogenic bacteria can cause food-borne diseases and negatively affect public health. In the milk industry, contamination with Salmonella has always been a challenge, but the risks have dramatically increased as almost all bacteria now show resistance to a wide range of commercial antibiotics. This study aimed to isolate a bacteriophage to be used as a bactericidal agent against Salmonella in milk and dairy products. Here, phage ZCSE6 has been isolated from raw milk sample sand molecularly and chemically characterized. At different multiplicities of infection (MOIs) of 0.1, 0.01, and 0.001, the phage–Salmonella interaction was studied for 6 h at 37 °C and 24 h at 8 °C. In addition, ZCSE6 was tested against Salmonella contamination in milk to examine its lytic activity for 3 h at 37 °C. The results showed that ZCSE6 has a small genome size (<48.5 kbp) and belongs to the Siphovirus family. Phage ZCSE6 revealed a high thermal and pH stability at various conditions that mimic milk manufacturing and supply chain conditions. It also demonstrated a significant reduction in Salmonella concentration in media at various MOIs, with higher bacterial eradication at higher MOI. Moreover, it significantly reduced Salmonella growth (MOI 1) in milk, manifesting a 1000-fold decrease in bacteria concentration following 3 h incubation at 37 °C. The results highlighted the strong ability of ZCSE6 to kill Salmonella and control its growth in milk. Thus, ZCSE6 is recommended as a biocontrol agent in milk to limit bacterial growth and increase the milk shelf-life.