Comparative genomic analysis of twelve Streptococcus suis (pro)phages

Isolation of phages infecting the zoonotic pathogen Streptococcus suis reveals novel structural and genomic characteristics

Bacteriophage research has experienced a renaissance in recent years, owing to their therapeutic potential and versatility in biotechnology, particularly in combating antibiotic resistant-bacteria along the farm-to-fork continuum. However, certain pathogens remain underexplored as targets for phage therapy, including the zoonotic pathogen Streptococcus suis which causes infections in pigs and humans. Despite global efforts, the genome of only one infective S. suis phage has been described. Here, we report the isolation of two phages that infect S. suis: Bonnie and Clyde. The phages infect 58 of 100 S. suis strains tested, including representatives of seven different serotypes and thirteen known sequence types from diverse geographical origins. Clyde suppressed bacterial growth in vitro within two multi-strain mixes designed to simulate a polyclonal S. suis infection. Both phages demonstrated stability across various temperatures and pH levels, highlighting their potential to withstand storage conditions and maintain viability in delivery formulations. Genome comparisons revealed that neither phage shares significant nucleotide identity with any cultivated phages in the NCBI database and thereby represent novel species belonging to two distinct novel genera. This study is the first to investigate the adhesion devices of S. suis infecting phages. Structure prediction and analysis of adhesion devices with AlphaFold2 revealed two distinct lineages of S. suis phages: Streptococcus thermophilus-like (Bonnie) and S. suis-like (Clyde). The structural similarities between the adhesion devices of Bonnie and S. thermophilus phages, despite the lack of nucleotide similarity and differing ecological niches, suggest a common ancestor or convergent evolution, highlighting evolutionary links between pathogenic and non-pathogenic streptococcal species. These findings provide valuable insights into the genetic and phenotypic characteristics of phages that can infect S. suis, providing new data for the therapeutic application of phages in a One Health context.

Comparative Genomic Analysis of 66 Bacteriophages Infecting Morganella morganii Strains

Bacteriophages are viruses that specifically target bacteria and play a crucial role in influencing bacterial evolution and the transmission of antibiotic resistance. In this study, we explored the genomic profiles of 66 bacteriophages that infect Morganella morganii, an opportunistic pathogen associated with difficult-to-treat nosocomial and urinary tract infections. Our findings highlight the extraordinary diversity within this phage population, reflected in their genomic features, evolutionary relationships, and potential contributions to bacterial pathogenicity. The 66 phage genomes exhibited diversity in size, spanning from 6 to 115 kilobase pairs, reflecting a heterogeneous genetic material and coding potential. Their guanine-cytosine (G+C) content varied widely, from 43.3% to 64.6%, suggesting diverse evolutionary origins and adaptive strategies. Phylogenetic analysis identified ten distinct evolutionary clusters, some classified as singletons, highlighting unique evolutionary pathways. Several clusters included phages capable of infecting multiple M. morganii strains, indicating a broader host range and the potential for horizontal gene transfer. Genomic analysis also determined a substantial number of hypothetical proteins, underscoring the need for further investigation to clarify their functions. Importantly, we identified a wide array of antibiotic resistance and virulence-associated genes within these phage genomes, illuminating their potential to impact the treatment of M. morganii infections and develop new, more virulent strains. These findings highlight the critical role of phage-mediated gene transfer in shaping bacterial evolution and facilitating the transmission of antibiotic resistance.

Isolation of phages infecting the zoonotic pathogen Streptococcus suis reveals novel structural and genomic characteristics

Bacteriophage research has experienced a renaissance in recent years, owing to their therapeutic potential and versatility in biotechnology, particularly in combating antibiotic resistant-bacteria along the farm-to-fork continuum. However, certain pathogens remain underexplored as targets for phage therapy, including the zoonotic pathogen Streptococcus suis which causes infections in pigs and humans. Despite global efforts, the genome of only one infective S. suis phage has been described. Here, we report the isolation of two phages that infect S. suis: Bonnie and Clyde. The phages infect 58% of 100 S. suis strains tested, including representatives of seven different serotypes and thirteen known sequence types from diverse geographical origins. Clyde suppressed bacterial growth in vitro within two multi-strain mixes designed to simulate a polyclonal S. suis infection. Both phages demonstrated stability across various temperatures and pH levels, highlighting their potential to withstand storage conditions and maintain viability in delivery formulations. Genome comparisons revealed that neither phage shares significant nucleotide identity with any cultivated phages in the NCBI database and thereby represent novel species belonging to two distinct novel genera. This study is the first to investigate the adhesion devices of S. suis infecting phages. Structure prediction and analysis of adhesion devices with AlphaFold2 revealed two distinct lineages of S. suis phages: Streptococcus thermophilus-like (Bonnie) and S. suis-like (Clyde). The structural similarities between the adhesion devices of Bonnie and S. thermophilus phages, despite the lack of nucleotide similarity and differing ecological niches, suggest a common ancestor or convergent evolution, highlighting evolutionary links between pathogenic and non-pathogenic streptococcal species. These findings provide valuable insights into the genetic and phenotypic characteristics of phages that can infect S. suis, providing new data for the therapeutic application of phages in a One Health context.

Comparative Genome Analysis of Two Streptococcus suis Serotype 8 Strains Identifies Two New Virulence-Associated Genes

Streptococcus suis is an important zoonotic pathogen that can cause meningitis and septicemia in swine and humans. Among numerous pathogenic serotypes, S. suis serotype 8 has distinctive characteristics such as a high detection rate and causing multi-host infection. There is no complete genome of serotype 8 strains so far. In this study, the complete genome of two S. suis serotype 8 strains, virulent strain 2018WUSS151 and non-virulent strain WUSS030, were sequenced. Comparative genomic analysis showed that the homology of the two genomes reaches 99.68%, and the main difference is the distinctive prophages. There are 83 genes unique to virulent strain 2018WUSS151, including three putative virulence-associated genes (PVGs). Two PVGs, padR and marR, are passenger genes in ISSsu2 family transposons that are able to form circular DNA intermediates during transposition, indicating the possibility of horizontal transmission among S. suis strains. The deletion mutant of PVGs marR or atpase attenuated the virulence of serotype 2 virulent SC070731 in a mouse infection model, confirming their role in S. suis virulence. These findings contribute to clarifying the genomic characterization of S. suis serotype 8 and S. suis pathogenesis.

Isolation and characterization of a novel lytic Parabacteroides distasonis bacteriophage φPDS1 from the human gut

The human gut microbiome plays a significant role in health and disease. The viral component (virome) is predominantly composed of bacteriophages (phages) and has received significantly less attention in comparison to the bacteriome. This knowledge gap is largely due to challenges associated with the isolation and characterization of novel gut phages, and bioinformatic hurdles such as the lack of a universal phage marker gene and the absence of sufficient numbers of homologs in viral databases. Here, we describe the isolation from human feces of a novel lytic phage with siphovirus morphology, φPDS1, infecting Parabacteroides distasonis APCS2/PD, and classified within a newly proposed Sagittacolavirus genus. In silico and biological characterization of this phage is presented in this study. Key to the isolation of φPDS1 was the antibiotic-driven selective enrichment of the bacterial host in a fecal fermenter. Despite producing plaques and lacking genes associated with lysogeny, φPDS1 demonstrates the ability to coexist in liquid culture for multiple days without affecting the abundance of its host. Multiple studies have shown that changes in Parabacteroides distasonis abundance can be linked to various disease states, rendering this novel phage-host pair and their interactions of particular interest.

Conjugative transfer of streptococcal prophages harboring antibiotic resistance and virulence genes

Prophages play important roles in the transduction of various functional traits, including virulence factors, but remain debatable in harboring and transmitting antimicrobial resistance genes (ARGs). Herein we characterize a prevalent family of prophages in Streptococcus, designated SMphages, which harbor twenty-five ARGs that collectively confer resistance to ten antimicrobial classes, including vanG-type vancomycin resistance locus and oxazolidinone resistance gene optrA. SMphages integrate into four chromosome attachment sites by utilizing three types of integration modules and undergo excision in response to phage induction. Moreover, we characterize four subtypes of Alp-related surface proteins within SMphages, the lethal effects of which are extensively validated in cell and animal models. SMphages transfer via high-frequency conjugation that is facilitated by integrative and conjugative elements from either donors or recipients. Our findings explain the widespread of SMphages and the rapid dissemination of ARGs observed in members of the Streptococcus genus.

Isolation and characterization of a novel phage belonging to a new genus against Vibrio parahaemolyticus

BACKGROUND

Vibrio parahaemolyticus is a major foodborne pathogen that contaminates aquatic products and causes great economic losses to aquaculture. Because of the emergence of multidrug-resistant V. parahaemolyticus strains, bacteriophages are considered promising agents for their biocontrol as an alternative or supplement to antibiotics. In this study, a lytic vibriophage, vB_VpaM_R16F (R16F), infecting V. parahaemolyticus 1.1997^T was isolated, characterized and evaluated for its biocontrol potential.

METHODS

A vibriophage R16F was isolated from sewage from a seafood market with the double-layer agar method. R16F was studied by transmission electron microscopy, host range, sensitivity of phage particles to chloroform, one-step growth curve and lytic activity. The phage genome was sequenced and in-depth characterized, including phylogenetic and taxonomic analysis.

RESULTS

R16F belongs to the myovirus morphotype and infects V. parahaemolyticus, but not nine other Vibrio spp. As characterized by determining its host range, one-step growth curve, and lytic activity, phage R16F was found to highly effective in lysing host cells with a short latent period (< 10 min) and a small burst size (13 plaque-forming units). R16F has a linear double-stranded DNA with genome size 139,011 bp and a G + C content of 35.21%. Phylogenetic and intergenomic nucleotide sequence similarity analysis revealed that R16F is distinct from currently known vibriophages and belongs to a novel genus. Several genes (e.g., encoding ultraviolet damage endonuclease and endolysin) that may enhance environmental competitiveness were found in the genome of R16F, while no antibiotic resistance- or virulence factor-related gene was detected.

CONCLUSIONS

In consideration of its biological and genetic properties, this newly discovered phage R16F belongs to a novel genus and may be a potential alternate biocontrol agent.

From Farm to Fork: Streptococcus suis as a Model for the Development of Novel Phage-Based Biocontrol Agents

Bacterial infections of livestock threaten the sustainability of agriculture and public health through production losses and contamination of food products. While prophylactic and therapeutic application of antibiotics has been successful in managing such infections, the evolution and spread of antibiotic-resistant strains along the food chain and in the environment necessitates the development of alternative or adjunct preventive and/or therapeutic strategies. Additionally, the growing consumer preference for "greener" antibiotic-free food products has reinforced the need for novel and safer approaches to controlling bacterial infections. The use of bacteriophages (phages), which can target and kill bacteria, are increasingly considered as a suitable measure to reduce bacterial infections and contamination in the food industry. This review primarily elaborates on the recent veterinary applications of phages and discusses their merits and limitations. Furthermore, using Streptococcus suis as a model, we describe the prevalence of prophages and the anti-viral defence arsenal in the genome of the pathogen as a means to define the genetic building blocks that are available for the (synthetic) development of phage-based treatments. The data and approach described herein may provide a framework for the development of therapeutics against an array of bacterial pathogens.

Comprehensive deciphering prophages in genus Acetobacter on the ecology, genomic features, toxin–antitoxin system, and linkage with CRISPR-Cas system

Acetobacter is the predominant microbe in vinegar production, particularly in those natural fermentations that are achieved by complex microbial communities. Co-evolution of prophages with Acetobacter, including integration, release, and dissemination, heavily affects the genome stability and production performance of industrial strains. However, little has been discussed yet about prophages in Acetobacter. Here, prophage prediction analysis using 148 available genomes from 34 Acetobacter species was carried out. In addition, the type II toxin–antitoxin systems (TAs) and CRISPR-Cas systems encoded by prophages or the chromosome were analyzed. Totally, 12,000 prophage fragments were found, of which 350 putatively active prophages were identified in 86.5% of the selected genomes. Most of the active prophages (83.4%) belonged to the order Caudovirales dominated by the families Siphoviridae and Myroviridae prophages (71.4%). Notably, Acetobacter strains survived in complex environments that frequently carried multiple prophages compared with that in restricted habits. Acetobacter prophages showed high genome diversity and horizontal gene transfer across different bacterial species by genomic feature characterization, average nucleotide identity (ANI), and gene structure visualization analyses. About 31.14% of prophages carry type II TAS, suggesting its important role in addiction, bacterial defense, and growth-associated bioprocesses to prophages and hosts. Intriguingly, the genes coding for Cse1, Cse2, Cse3, Cse4, and Cas5e involved in type I-E and Csy4 involved in type I-F CRISPR arrays were firstly found in two prophages. Type II-C CRISPR-Cas system existed only in Acetobacter aceti, while the other Acetobacter species harbored the intact or eroded type I CRISPR-Cas systems. Totally, the results of this study provide fundamental clues for future studies on the role of prophages in the cell physiology and environmental behavior of Acetobacter.

Prevalence, Diversity and UV-Light Inducibility Potential of Prophages in Bacillus subtilis and Their Possible Roles in Host Properties

Bacillus subtilis is an important bacterial species due to its various industrial, medicinal, and agricultural applications. Prophages are known to play vital roles in host properties. Nevertheless, studies on the prophages and temperate phages of B. subtilis are relatively limited. In the present study, an in silico analysis was carried out in sequenced B. subtilis strains to investigate their prevalence, diversity, insertion sites, and potential roles. In addition, the potential for UV induction and prevalence was investigated. The in silico prophage analysis of 164 genomes of B. subtilis strains revealed that 75.00% of them contained intact prophages that exist as integrated and/or plasmid forms. Comparative genomics revealed the rich diversity of the prophages distributed in 13 main clusters and four distinct singletons. The analysis of the putative prophage proteins indicated the involvement of prophages in encoding the proteins linked to the immunity, bacteriocin production, sporulation, arsenate, and arsenite resistance of the host, enhancing its adaptability to diverse environments. An induction study in 91 B. subtilis collections demonstrated that UV-light treatment was instrumental in producing infective phages in 18.68% of them, showing a wide range of host specificity. The high prevalence and inducibility potential of the prophages observed in this study implies that prophages may play vital roles in the B. subtilis host.

Comparative genomic analysis of dwarf Vibrio myoviruses defines a conserved gene cluster for successful phage infection

Tailed bacteriophages have been at the center of attention, not only for their ability to infect and kill pathogenic bacteria but also due to their peculiar and intriguing complex contractile tail structure. Tailed bacteriophages with contractile tails are known to have a Myoviridae morphotype and are members of the order Caudovirales. Large bacteriophages with a genome larger than 150 kbp have been studied for their ability to use multiple infection and lysis strategies to replicate more efficiently. On the other hand, smaller bacteriophages with fewer genes are represented in the GenBank database in greater numbers, and have several genes with unknown function. Isolation and molecular characterization of a newly reported bacteriophage named Athena1 revealed that it is a strongly lytic bacteriophage with a genome size of 39,826 bp. This prompted us to perform a comparative genomic analysis of Vibrio myoviruses with a genome size of no more than 50 kbp. The results revealed a pattern of genomic organization that includes sets of genes responsible for virion morphogenesis, replication/recombination of DNA, and lysis/lysogeny switching. By studying phylogenetic gene markers, we were able to draw conclusions about evolutionary events that shaped the genomic mosaicism of these phages, pinpointing the importance of a conserved organization of the genomic region encoding the baseplate protein for successful infection of Gram-negative bacteria. In addition, we propose the creation of new genera for dwarf Vibrio myoviruses. Comparative genomics of phages infecting aquatic bacteria could provide information that is useful for combating fish pathogens in aquaculture, using novel strategies.

Genome analysis provides insight into hyper-virulence of Streptococcus suis LSM178, a human strain with a novel sequence type 1005

Streptococcus suis has been well-recognized as a zoonotic pathogen worldwide, and the diversity and unpredictable adaptive potential of sporadic human strains represent a great risk to the public health. In this study, S. suis LSM178, isolated from a patient in contact with pigs and raw pork, was assessed as a hyper-virulent strain and interpreted for the virulence based on its genetic information. The strain was more invasive for Caco-2 cells than two other S. suis strains, SC19 and P1/7. Sequence analysis designated LSM178 with serotype 2 and a novel sequence type 1005. Phylogenetic analysis showed that LSM178 clustered with highly virulent strains including all human strains and epidemic strains. Compared with other strains, these S. suis have the most and the same virulent factors and a type I-89 K pathogenicity island. Further, groups of genes were identified to distinguish these highly virulent strains from other generally virulent strains, emphasizing the key roles of genes modeling transcription, cell barrier, replication, recombination and repair on virulence regulation. Additionally, LSM178 contains a novel prophage conducive potentially to pathogenicity.

Comparative analysis of prophages carried by human and animal-associated Staphylococcus aureus strains spreading across the European regions

Staphylococcus aureus is a major human and animal pathogen although the animal-associated S. aureus can be a potential risk of human zoonoses. Acquisition of phage-related genomic islands determines the S. aureus species diversity. This study characterized and compared the genome architecture, distribution nature, and evolutionary relationship of 65 complete prophages carried by human and animal-associated S. aureus strains spreading across the European regions. The analyzed prophage genomes showed mosaic architecture with extensive variation in genome size. The phylogenetic analyses generated seven clades in which prophages of the animal-associated S. aureus scattered in all the clades. The S. aureus strains with the same SCCmec type, and clonal complex favored the harboring of similar prophage sequences and suggested that the frequency of phage-mediated horizontal gene transfer is higher between them. The presence of various virulence factors in prophages of animal-associated S. aureus suggested that these prophages could have more pathogenic potential than prophages of human-associated S. aureus. This study showed that the S. aureus phages are dispersed among the several S. aureus serotypes and around the European regions. Further, understanding the phage functional genomics is necessary for the phage-host interactions and could be used for tracing the S. aureus strains transmission.

Characterization of Clinical and Carrier Streptococcus agalactiae and Prophage Contribution to the Strain Variability

Streptococcus agalactiae (group B Streptococcus, GBS) represents a leading cause of invasive bacterial infections in newborns and is also responsible for diseases in older and immunocompromised adults. Prophages represent an important factor contributing to the genome plasticity and evolution of new strains. In the present study, prophage content was analyzed in human GBS isolates. Thirty-seven prophages were identified in genomes of 20 representative sequenced strains. On the basis of the sequence comparison, we divided the prophages into eight groups named A–H. This division also corresponded to the clustering of phage integrase, even though several different integration sites were observed in some relative prophages. Next, PCR method was used for detection of the prophages in 123 GBS strains from adult hospitalized patients and from pregnancy screening. At least one prophage was present in 105 isolates (85%). The highest prevalence was observed for prophage group A (71%) and satellite prophage group B (62%). Other groups were detected infrequently (1–6%). Prophage distribution did not differ between clinical and screening strains, but it was unevenly distributed in MLST (multi locus sequence typing) sequence types. High content of full-length and satellite prophages detected in present study implies that prophages could be beneficial for the host bacterium and could contribute to evolution of more adapted strains.

Characterization of the Bacteriophage-Derived Endolysins PlySs2 and PlySs9 with In Vitro Lytic Activity against Bovine Mastitis Streptococcus uberis

Bovine mastitis, an infection of the cow’s mammary gland, is frequently caused by Streptococcus uberis and causes major economic losses in the dairy industry. The intramammary administration of antibiotics currently remains the predominant preventive and therapeutic measure. These antimicrobial compounds, of which some are considered critical in human health care, are frequently applied as dry therapy resulting in their consistent overuse. Therefore, the use of antibiotics in the dairy sector is being questioned. We here identified two endolysins, i.e., PlySs2 and PlySs9, respectively derived from Streptococcus suis serotype-2 and -9 prophages, with lytic activity against S. uberis in an in vitro setting. Both endolysins gave clear lysis zones in spot-on-plate assays and caused a reduction of the optical density in a turbidity reduction assay. In depth characterization identified PlySs9 as the more potent endolysin over PlySs2 with a lower MIC value and about one additional log of killing. PlySs2 and PlySs9 were challenged to a panel of subclinical and clinical S. uberis milk isolates and were both able to lyse all strains tested. Molecular dissection of these endolysins in catalytic and cell wall binding subdomains resulted in major loss of killing and binding activity, respectively. Taken together, we here propose PlySs2 and PlySs9 as candidate compounds to the current antimicrobial arsenal known against bovine mastitis-causing S. uberis as future add-on or replacement strategy to the currently used intramammary antibiotics.

The Isolation and Genome Sequencing of Five Novel Bacteriophages From the Rumen Active Against Butyrivibrio fibrisolvens

Although the prokaryotic communities of the rumen microbiome are being uncovered through genome sequencing, little is known about the resident viral populations. Whilst temperate phages can be predicted as integrated prophages when analyzing bacterial and archaeal genomes, the genetics underpinning lytic phages remain poorly characterized. To the five genomes of bacteriophages isolated from rumen-associated samples sequenced and analyzed previously, this study adds a further five novel genomes and predictions gleaned from them to further the understanding of the rumen phage population. Lytic bacteriophages isolated from fresh ovine and bovine fecal and rumen fluid samples were active against the predominant fibrolytic ruminal bacterium Butyrivibrio fibrisolvens. The double stranded DNA genomes were sequenced and reconstructed into single circular complete contigs. Based on sequence similarity and genome distances, the five phages represent four species from three separate genera, consisting of: (1) Butyrivibrio phages Arian and Bo-Finn; (2) Butyrivibrio phages Idris and Arawn; and (3) Butyrivibrio phage Ceridwen. They were predicted to all belong to the Siphoviridae family, based on evidence in the genomes such as size, the presence of the tail morphogenesis module, genes that share similarity to those in other siphovirus isolates and phylogenetic analysis using phage proteomes. Yet, phylogenomic analysis and sequence similarity of the entire phage genomes revealed that these five phages are unique and novel. These phages have only been observed undergoing the lytic lifecycle, but there is evidence in the genomes of phages Arawn and Idris for the potential to be temperate. However, there is no evidence in the genome of the bacterial host Butyrivibrio fibrisolvens of prophage genes or genes that share similarity with the phage genomes.

Characterization of Bacteriophage Peptides of Pathogenic Streptococcus by LC-ESI-MS/MS: Bacteriophage Phylogenomics and Their Relationship to Their Host

The present work focuses on LC-ESI-MS/MS (liquid chromatography-electrospray ionization-tandem mass spectrometry) analysis of phage-origin tryptic digestion peptides from mastitis-causing Streptococcus spp. isolated from milk. A total of 2,546 non-redundant peptides belonging to 1,890 proteins were identified and analyzed. Among them, 65 phage-origin peptides were determined as specific Streptococcus spp. peptides. These peptides belong to proteins such as phage repressors, phage endopeptidases, structural phage proteins, and uncharacterized phage proteins. Studies involving bacteriophage phylogeny and the relationship between phages encoding the peptides determined and the bacteria they infect were also performed. The results show how specific peptides are present in closely related phages, and a link exists between bacteriophage phylogeny and the Streptococcus spp. they infect. Moreover, the phage peptide M∗ATNLGQAYVQIM∗PSAK is unique and specific for Streptococcus agalactiae. These results revealed that diagnostic peptides, among others, could be useful for the identification and characterization of mastitis-causing Streptococcus spp., particularly peptides that belong to specific functional proteins, such as phage-origin proteins, because of their specificity to bacterial hosts.

Comparative Genome Analysis of a Pathogenic Erysipelothrix rhusiopathiae Isolate WH13013 from Pig Reveals Potential Genes Involve in Bacterial Adaptions and Pathogenesis

Erysipelothrix rhusiopathiae is a common pathogen responsible for pig erysipelas. However, the molecular basis for the pathogenesis of E. rhusiopathiae remains to be elucidated. In this study, the complete genome sequence of the E. rhusiopathiae strain WH13013, a pathogenic isolate from a diseased pig, was generated using a combined strategy of PacBio RSII and Illumina sequencing technologies. The strategy finally generated a single circular chromosome of approximately 1.78 Mb in size for the complete genome of WH13013, with an average GC content of 36.49%. The genome of WH13013 encoded 1633 predicted proteins, 55 tRNAs, as well as 15 rRNAs. It contained four genomic islands and several resistance-associated genes were identified within these islands. Phylogenetic analysis revealed that WH13013 was close to many other sequenced E. rhusiopathiae virulent strains. The comprehensive comparative analysis of eight E. rhusiopathiae virulent strains, including WH13013, identified a total of 1184 core genes. A large proportion (approximately 75.31%) of these core genes participated in nutrition and energy uptake and metabolism as well as the other bioactivities that are necessary for bacterial survival and adaption. The core genes also contained those encoding proteins participating in the biosynthesis and/or the components of the proposed virulence factors of E. rhusiopathiae, including the capsule (cpsA, cpsB, cpsC), neuraminidase (nanH), hyaluronidase (hylA, hylB, hylC), and surface proteins (spaA, rspA, rspB). The obtaining of the complete genome sequence of this virulent strain, WH13013, and this comprehensive comparative genome analysis will help in further studies of the genetic basis of the pathogenesis of E. rhusiopathiae.

phiD12-Like Livestock-Associated Prophages Are Associated With Novel Subpopulations of Streptococcus agalactiae Infecting Neonates

Group B Streptococcus (GBS) is a major cause of invasive disease in neonates worldwide. Monitoring data have revealed a continuing trend toward an increase in neonatal GBS infections, despite the introduction of preventive measures. We investigated this trend, by performing the first ever characterization of the prophage content for 106 GBS strains causing neonatal infections between 2002 and 2018. We determined whether the genome of each strain harbored prophages, and identified the insertion site of each of the prophages identified. We found that 71.7% of the strains carried at least one prophage, and that prophages genetically similar to livestock-associated phiD12, carrying genes potentially involved in GBS pathogenesis (e.g., genes encoding putative virulence factors and factors involved in biofilm formation, bacterial persistence, or adaptation to challenging environments) predominated. The phiD12-like prophages were (1) associated with CC17 and 1 strains (p = 0.002), (2) more frequent among strains recovered during the 2011–2018 period than among those from 2002–2010 (p < 0.001), and (3) located at two major insertion sites close to bacterial genes involved in host adaptation and colonization. Our data provide evidence for a recent increase in lysogeny in GBS, characterized by the acquisition, within the genome, of genetic features typical of animal-associated mobile genetic elements by GBS strains causing neonatal infection. We suggest that lysogeny and phiD12-like prophage genetic elements may have conferred an advantage on GBS strains for adaptation to or colonization of the maternal vaginal tract, or for pathogenicity, and that these factors are currently playing a key role in the increasing ability of GBS strains to infect neonates.

Bacterial and Archaeal Viruses of Himalayan Hot Springs at Manikaran Modulate Host Genomes

Hot spring-associated viruses, particularly the archaeal viruses, remain under-examined compared to bacteriophages. Previous metagenomic studies of the Manikaran hot springs in India suggested an abundance of viral DNA, which prompted us to examine the virus–host (bacterial and archaeal) interactions in sediment and microbial mat samples collected from the thermal discharges. Here, we characterize the viruses (both bacterial and archaeal) from this Himalayan hot spring using both metagenomics assembly and electron microscopy. We utilized four shotgun samples from sediment (78–98°C) and two from microbial mats (50°C) to reconstruct 65 bacteriophage genomes (24–200 kb). We also identified 59 archaeal viruses that were notably abundant across the sediment samples. Whole-genome analyses of the reconstructed bacteriophage genomes revealed greater genomic conservation in sediments (65%) compared to microbial mats (49%). However, a minimal phage genome was still maintained across both sediment and microbial mats suggesting a common origin. To complement the metagenomic data, scanning-electron and helium-ion microscopy were used to reveal diverse morphotypes of Caudovirales and archaeal viruses. The genome level annotations provide further evidence for gene-level exchange between virus and host in these hot springs, and augments our knowledgebase for bacteriophages, archaeal viruses and Clustered Regularly Interspaced Short Palindromic Repeat cassettes, which provide a critical resource for studying viromes in extreme natural environments.

Two novel temperate bacteriophages infecting Streptococcus pyogenes: Their genomes, morphology and stability

Only 3% of phage genomes in NCBI nucleotide database represent phages that are active against Streptococcus sp. With the aim to increase general awareness of phage diversity, we isolated two bacteriophages, Str01 and Str03, active against health-threatening Group A Streptococcus (GAS). Both phages are members of the Siphoviridae, but their analysis revealed that Str01 and Str03 do not belong to any known genus. We identified their structural proteins based on LC-ESI29 MS/MS and list their basic thermal stability and physico-chemical features including optimum pH. Annotated genomic sequences of the phages are deposited in GenBank (NCBI accession numbers KY349816 and KY363359, respectively).

Csl2, a novel chimeric bacteriophage lysin to fight infections caused by Streptococcus suis, an emerging zoonotic pathogen

Streptococcus suis is a Gram-positive bacterium that infects humans and various animals, causing human mortality rates ranging from 5 to 20%, as well as important losses for the swine industry. In addition, there is no effective vaccine for S. suis and isolates with increasing antibiotic multiresistance are emerging worldwide. Facing this situation, wild type or engineered bacteriophage lysins constitute a promising alternative to conventional antibiotics. In this study, we have constructed a new chimeric lysin, Csl2, by fusing the catalytic domain of Cpl-7 lysozyme to the CW_7 repeats of LySMP lysin from an S. suis phage. Csl2 efficiently kills different S. suis strains and shows noticeable activity against a few streptococci of the mitis group. Specifically, 15 µg/ml Csl2 killed 4.3 logs of S. suis serotype 2 S735 strain in 60 min, in a buffer containing 150 mM NaCl and 10 mM CaCl₂, at pH 6.0. We have set up a protocol to form a good biofilm with the non-encapsulated S. suis mutant strain BD101, and the use of 30 µg/ml Csl2 was enough for dispersing such biofilms and reducing 1–2 logs the number of planktonic bacteria. In vitro results have been validated in an adult zebrafish model of infection.

Comparative analysis of prophages in Streptococcus mutans genomes

Prophages have been considered genetic units that have an intimate association with novel phenotypic properties of bacterial hosts, such as pathogenicity and genomic variation. Little is known about the genetic information of prophages in the genome of Streptococcus mutans, a major pathogen of human dental caries. In this study, we identified 35 prophage-like elements in S. mutans genomes and performed a comparative genomic analysis. Comparative genomic and phylogenetic analyses of prophage sequences revealed that the prophages could be classified into three main large clusters: Cluster A, Cluster B, and Cluster C. The S. mutans prophages in each cluster were compared. The genomic sequences of phismuN66-1, phismuNLML9-1, and phismu24-1 all shared similarities with the previously reported S. mutans phages M102, M102AD, and ϕAPCM01. The genomes were organized into seven major gene clusters according to the putative functions of the predicted open reading frames: packaging and structural modules, integrase, host lysis modules, DNA replication/recombination modules, transcriptional regulatory modules, other protein modules, and hypothetical protein modules. Moreover, an integrase gene was only identified in phismuNLML9-1 prophages.

Phenetic Comparison of Prokaryotic Genomes Using k-mers

Bacterial genomics studies are getting more extensive and complex, requiring new ways to envision analyses. Using the Ray Surveyor software, we demonstrate that comparison of genomes based on their k-mer content allows reconstruction of phenetic trees without the need of prior data curation, such as core genome alignment of a species. We validated the methodology using simulated genomes and previously published phylogenomic studies of Streptococcus pneumoniae and Pseudomonas aeruginosa. We also investigated the relationship of specific genetic determinants with bacterial population structures. By comparing clusters from the complete genomic content of a genome population with clusters from specific functional categories of genes, we can determine how the population structures are correlated. Indeed, the strain clustering based on a subset of k-mers allows determination of its similarity with the whole genome clusters. We also applied this methodology on 42 species of bacteria to determine the correlational significance of five important bacterial genomic characteristics. For example, intrinsic resistance is more important in P. aeruginosa than in S. pneumoniae, and the former has increased correlation of its population structure with antibiotic resistance genes. The global view of the pangenome of bacteria also demonstrated the taxa-dependent interaction of population structure with antibiotic resistance, bacteriophage, plasmid, and mobile element k-mer data sets.

Analysis of the prophages carried by human infecting isolates provides new insight into the evolution of Group B Streptococcus species

OBJECTIVES

Group B Streptococcus (GBS) emerged in the 1970s as a major cause of neonatal infections, and has been increasingly associated with infections in adults since the 1990s. Prophages have been suspected to have driven these epidemiological trends. We have characterized the prophages harboured by 275 human GBS isolates belonging to the major lineages.

METHODS

We applied whole genome sequencing (WGS) to 14 isolates representative of the diversity within GBS species, located and identified their prophages. Using prediction tools, we searched for prophage elements potentially involved with the ability of GBS to infect humans. Using the data obtained by WGS, we designed a PCR-based tool and studied the prophage content of 275 isolates.

RESULTS

WGS of the 14 isolates revealed 22 prophages (i) distributed into six groups (A-F), (ii) similar to phages and prophages from GBS and non-GBS streptococci recovered from livestock, and (iii) carrying genes encoding factors previously associated with host adaptation and virulence. PCR-based detection of prophages revealed the presence of at least one prophage in 72.4% of the 275 isolates and a significant association between neonatal infecting isolates and prophages C, and between adult infecting isolates and prophages A.

CONCLUSIONS

Our results suggest that prophages (possibly animal-associated) have conditioned bacterial adaptation and ability to cause infections in neonates and adults, and support a role of lysogeny with the emergence of GBS as a pathogen in human.

Genome Plasticity and Polymorphisms in Critical Genes Correlate with Increased Virulence of Dutch Outbreak-Related Coxiella burnetii Strains

Coxiella burnetii is an obligate intracellular bacterium and the etiological agent of Q fever. During 2007-2010 the largest Q fever outbreak ever reported occurred in The Netherlands. It is anticipated that strains from this outbreak demonstrated an increased zoonotic potential as more than 40,000 individuals were assumed to be infected. The acquisition of novel genetic factors by these C. burnetii outbreak strains, such as virulence-related genes, has frequently been proposed and discussed, but is not proved yet. In the present study, the whole genome sequence of several Dutch strains (CbNL01 and CbNL12 genotypes), a few additionally selected strains from different geographical locations and publicly available genome sequences were used for a comparative bioinformatics approach. The study focuses on the identification of specific genetic differences in the outbreak related CbNL01 strains compared to other C. burnetii strains. In this approach we investigated the phylogenetic relationship and genomic aspects of virulence and host-specificity. Phylogenetic clustering of whole genome sequences showed a genotype-specific clustering that correlated with the clustering observed using Multiple Locus Variable-number Tandem Repeat Analysis (MLVA). Ortholog analysis on predicted genes and single nucleotide polymorphism (SNP) analysis of complete genome sequences demonstrated the presence of genotype-specific gene contents and SNP variations in C. burnetii strains. It also demonstrated that the currently used MLVA genotyping methods are highly discriminatory for the investigated outbreak strains. In the fully reconstructed genome sequence of the Dutch outbreak NL3262 strain of the CbNL01 genotype, a relatively large number of transposon-linked genes were identified as compared to the other published complete genome sequences of C. burnetii. Additionally, large numbers of SNPs in its membrane proteins and predicted virulence-associated genes were identified in all Dutch outbreak strains compared to the NM reference strain and other strains of the CbNL12 genotype. The presence of large numbers of transposable elements and mutated genes, thereof most likely resulted in high level of genome rearrangements and genotype-specific pathogenicity of outbreak strains. Thus, the epidemic potential of Dutch outbreak strains could be linked to increased genome plasticity and mutations in critical genes involved in virulence and the evasion of the host immune system.

A Glimpse into the World of Integrative and Mobilizable Elements in Streptococci Reveals an Unexpected Diversity and Novel Families of Mobilization Proteins

Recent analyses of bacterial genomes have shown that integrated elements that transfer by conjugation play an essential role in horizontal gene transfer. Among these elements, the integrative and mobilizable elements (IMEs) are known to encode their own excision and integration machinery, and to carry all the sequences or genes necessary to hijack the mating pore of a conjugative element for their own transfer. However, knowledge of their prevalence and diversity is still severely lacking. In this work, an extensive analysis of 124 genomes from 27 species of Streptococcus reveals 144 IMEs. These IMEs encode either tyrosine or serine integrases. The identification of IME boundaries shows that 141 are specifically integrated in 17 target sites. The IME-encoded relaxases belong to nine superfamilies, among which four are previously unknown in any mobilizable or conjugative element. A total of 118 IMEs are found to encode a non-canonical relaxase related to rolling circle replication initiators (belonging to the four novel families or to MobT). Surprisingly, among these, 83 encode a TcpA protein (i.e., a non-canonical coupling protein (CP) that is more closely related to FtsK than VirD4) that was not previously known to be encoded by mobilizable elements. Phylogenetic analyses reveal not only many integration/excision module replacements but also losses, acquisitions or replacements of TcpA genes between IMEs. This glimpse into the still poorly known world of IMEs reveals that mobilizable elements have a very high prevalence. Their diversity is even greater than expected, with most encoding a CP and/or a non-canonical relaxase.

Experimental pathogenicity and complete genome characterization of a pig origin Pasteurella multocida serogroup F isolate HN07

Pasteurella multocida serotype F isolates are predominately prevalent in avian hosts, but rarely seen in pigs. However, we isolated several strains of P. multocida serotype F from clinical samples of pigs in China. To understand the pathogenicity of these strains, one of the serotype F isolates designated HN07, was used to challenge experimental chickens, as P. multocida of this serotype is predominately prevalent in avian hosts. However, strain HN07 could not resulted in significant clinical signs in experimental chickens even at an infective dose of ∼10⁹ CFU, suggesting the isolate was avirulent to chickens and therefore raising the possibility that the porcine serotype F isolate is not transmitted by chickens. We then used HN07 to challenge experimental pigs, as this strain was isolated from pigs. As expected, the strain led to the clinical signs and the pathological lesions in experimental pigs that are similar to the pasteurellosis disease. We then determined the complete genome sequence of the pig origin serogroup F isolate HN07 for the first time. Genome comparison between HN07 and the avian serotype F P. multocida Pm70 identified a novel integrative conjugative element (ICE) ICEpmcn07 which was likely to harbor a series of genes responsible for a putative type IV secretion system (T4SS) in HN07. This is the first time that we determined an ICE carrying a T4SS in P. multocida. Besides, comparative analysis also defined a number of virulence-associated genes in HN07 but absent in Pm70 which may have a contribution to the pathogenicity of the strain. This is the first report of the pathogenicity and genome characterization of a pig origin Pasteurella multocida serogroup F isolate. The pathogenic and genomic definition of the pig origin P. multocida serogroup F in our study would have significance on the pathogenesis and genetic diversity and virulence variability of P. multocida.

A novel prophage lysin Ply5218 with extended lytic activity and stability against Streptococcus suis infection

Streptococcus suis (S. suis) is an emerging zoonotic agent that exhibits high level resistance to classic antibiotics and a heavy burden in the swine industry. Therefore alternative antibacterial agents need to be developed. A novel endolysin derived from the S. suis temperate phage phi5218, termed Ply5218, was identified. The minimum inhibitory concentration (MIC) of Ply5218 was 2.5 μg ml(-1) against S. suis strain HA9801, an activity many times greater than the lysins reported previously (MIC of LY7917 and Ply30 against HA9801 were 80 and 64 μg ml(-1), respectively). Ply5218 at 10 μg ml(-1) in vitro exerted broad antibacterial activities against S. suis strains with OD600 ratios decreased from 1 to <0.2 within 1 h. Moreover, Ply5218 showed favorable thermal stability. It was stable at 50°C >30 min, 4°C >30 days, -80°C >7 months, and >60% of the enzyme activity remained after 5 min pre-incubation at 70°C. In vivo, a 0.2 mg dose of Ply5218 protected 90% (9/10) of mice after infection with S. suis HA9801. Finally, Ply5218 maintained high antibacterial activity in some bio-matrices, such as culture media and milk. The data indicate that Ply5218 has all the characteristics to be an effective therapeutic agent against multiple S. suis infections.

Comparative analysis of prophage-like elements in Helicobacter sp. genomes

Prophages are regarded as one of the factors underlying bacterial virulence, genomic diversification, and fitness, and are ubiquitous in bacterial genomes. Information on Helicobacter sp. prophages remains scarce. In this study, sixteen prophages were identified and analyzed in detail. Eight of them are described for the first time. Based on a comparative genomic analysis, these sixteen prophages can be classified into four different clusters. Phylogenetic relationships of Cluster A Helicobacter prophages were investigated. Furthermore, genomes of Helicobacter prophages from Clusters B, C, and D were analyzed. Interestingly, some putative antibiotic resistance proteins and virulence factors were associated with Helicobacter prophages.

A Zebrafish Larval Model to Assess Virulence of Porcine Streptococcus suis Strains

Streptococcus suis is an encapsulated Gram-positive bacterium, and the leading cause of sepsis and meningitis in young pigs resulting in considerable economic losses in the porcine industry. It is also considered an emerging zoonotic agent. In the environment, both avirulent and virulent strains occur in pigs, and virulent strains appear to cause disease in both humans and pigs. There is a need for a convenient, reliable and standardized animal model to assess S. suis virulence. A zebrafish (Danio rerio) larvae infection model has several advantages, including transparency of larvae, low cost, ease of use and exemption from ethical legislation up to 6 days post fertilization, but has not been previously established as a model for S. suis. Microinjection of different porcine strains of S. suis in zebrafish larvae resulted in highly reproducible dose- and strain-dependent larval death, strongly correlating with presence of the S. suis capsule and to the original virulence of the strain in pigs. Additionally we compared the virulence of the two-component system mutant of ciaRH, which is attenuated for virulence in both mice and pigs in vivo. Infection of larvae with the ΔciaRH strain resulted in significantly higher survival rate compared to infection with the S10 wild-type strain. Our data demonstrate that zebrafish larvae are a rapid and reliable model to assess the virulence of clinical porcine S. suis isolates.

Genomic characterization of Pasteurella multocida HB01, a serotype A bovine isolate from China

Pasteurella multocida infects various domestic and feral animals, generally causing clinical disease. To investigate P. multocida disease in cattle, we sequenced the complete genome of P. multocida HB01 (GenBank accession CP006976), a serotype A organism isolated from a cow in China. The genome is composed of a single circular chromosome of 2,416,068 base pairs containing 2212 protein-coding sequences, 6 ribosomal rRNA operons, and 56 tRNA genes. The present study confirms that P. multocida HB01 possesses a more complete metabolic pathway with an intact trichloroacetic acid cycle for anabolism compared with A. pleuropneumoniae and Haemophilus parasuis. This is the first time that this metabolic mechanism of P. multocida has been described. We also identified a full spectrum of genes related to known virulence factors of P. multocida. The differences in virulence factors between strains of different serotypes and origins were also compared. This comprehensive comparative genome analysis will help in further studies of the metabolic pathways, genetic basis of serotype, and virulence of P. multocida.

Novel Moraxella catarrhalis prophages display hyperconserved non-structural genes despite their genomic diversity

Background

Moraxella catarrhalis is an important pathogen that often causes otitis media in children, a disease that is not currently vaccine preventable. Asymptomatic colonisation of the human upper respiratory tract is common and lack of clearance by the immune system is likely due to the emergence of seroresistant genetic lineages. No active bacteriophages or prophages have been described in this species. This study was undertaken to identify and categorise prophages in M. catarrhalis, their genetic diversity and the relationship of such diversity with the host-species phylogeny.

Results

This study presents a comparative analysis of 32 putative prophages identified in 95 phylogenetically variable, newly sequenced M. catarrhalis genomes. The prophages were genotypically classified into four diverse clades. The genetic synteny of each clade is similar to the group 1 phage family Siphoviridae, however, they form genotypic clusters that are distinct from other members of this family. No core genetic sequences exist across the 32 prophages despite clades 2, 3, and 4 sharing the most sequence identity. The analysis of non-structural prophage genes (coding the integrase, and terminase), and portal gene showed that the respective genes were identical for clades 2, 3, and 4, but unique for clade 1. Empirical analysis calculated that these genes are unexpectedly hyperconserved, under purifying selection, suggesting a tightly regulated functional role. As such, it is improbable that the prophages are decaying remnants but stable components of a fluctuating, flexible and unpredictable system ultimately maintained by functional constraints on non-structural and packaging genes. Additionally, the plate encoding genes were well conserved across all four prophage clades, and the tail fibre genes, commonly responsible for receptor recognition, were clustered into three major groups distributed across the prophage clades. A pan-genome of 283,622 bp was identified, and the prophages were mapped onto the diverse M. catarrhalis multi-locus sequence type (MLST) backbone.

Conclusion

This study has provided the first evidence of putatively mobile prophages in M. catarrhalis, identifying a diverse and fluctuating system dependent on the hyperconservation of a few key, non-structural genes. Some prophages harbour virulence-related genes, and potentially influence the physiology and virulence of M. catarrhalis. Importantly our data will provide supporting information on the identification of novel prophages in other species by adding greater weight to the identification of non-structural genes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2104-1) contains supplementary material, which is available to authorized users.

Prophage lysin Ply30 protects mice from Streptococcus suis and Streptococcus equi subsp. zooepidemicus infections

Streptococcus suis and Streptococcus equi subsp. zooepidemicus are capable of infecting humans and various animals, causing significant problems for the worldwide swine industry. As antibiotic resistance has increased, lysosomal enzymes encoded by phages have shown potential for use against pathogenic bacteria. In this study, a novel bacteriophage lysin, Ply30, encoded by the S. suis prophage phi30c, was recombinantly expressed and purified. Ply30 showed high bacteriolysis activity on S. suis and S. equi subsp. zooepidemicus in vitro. The ratio of the optical density at 600 nm (OD600) with treatment versus the OD600 with no treatment for most tested S. suis and S. equi subsp. zooepidemicus strains decreased from 1 to <0.3 and <0.5, respectively, within 1 h. The results of plate viability assays showed that treated bacteria suffered a 1- to 2-log decrease in CFU within 1 h. The optimal concentration of Ply30 was 50 μg/ml, and the optimal pH was 7. Moreover, Ply30 maintained high activity over a wide pH range (pH 6 to 10). The MICs of Ply30 against Streptococcus strains ranged from 16 to 512 μg/ml. In vivo, a 2-mg dose of Ply30 protected 90% (9/10 mice) of mice from infection with S. equi subsp. zooepidemicus and 80% (8/10 mice) of mice from infection with S. suis. Seven days after lysin Ply30 treatment, bacterial loads were significantly decreased in all tested organs and blood compared with those at 1 h postinfection without Ply30 treatment. Ply30 showed in vitro and in vivo antimicrobial efficiency and protected mice against two kinds of bacterial infections, indicating that Ply30 may be an effective therapeutic against streptococci.

Comparative genomic and morphological analyses of Listeria phages isolated from farm environments

The genus Listeria is ubiquitous in the environment and includes the globally important food-borne pathogen Listeria monocytogenes. While the genomic diversity of Listeria has been well studied, considerably less is known about the genomic and morphological diversity of Listeria bacteriophages. In this study, we sequenced and analyzed the genomes of 14 Listeria phages isolated mostly from New York dairy farm environments as well as one related Enterococcus faecalis phage to obtain information on genome characteristics and diversity. We also examined 12 of the phages by electron microscopy to characterize their morphology. These Listeria phages, based on gene orthology and morphology, together with previously sequenced Listeria phages could be classified into five orthoclusters, including one novel orthocluster. One orthocluster (orthocluster I) consists of large genome (~135-kb) myoviruses belonging to the genus “Twort-like viruses,” three orthoclusters (orthoclusters II to IV) contain small-genome (36- to 43-kb) siphoviruses with icosahedral heads, and the novel orthocluster V contains medium-sized-genome (~66-kb) siphoviruses with elongated heads. A novel orthocluster (orthocluster VI) of E. faecalis phages, with medium-sized genomes (~56 kb), was identified, which grouped together and shares morphological features with the novel Listeria phage orthocluster V. This new group of phages (i.e., orthoclusters V and VI) is composed of putative lytic phages that may prove to be useful in phage-based applications for biocontrol, detection, and therapeutic purposes.

Cell-free propagation of Coxiella burnetii does not affect its relative virulence

Q fever is caused by the obligate intracellular bacterium Coxiella burnetii. In vitro growth of the bacterium is usually limited to viable eukaryotic host cells imposing experimental constraints for molecular studies, such as the identification and characterisation of major virulence factors. Studies of pathogenicity may benefit from the recent development of an extracellular growth medium for C. burnetii. However, it is crucial to investigate the consistency of the virulence phenotype of strains propagated by the two fundamentally different culturing systems. In the present study, we assessed the viability of C. burnetii and the lipopolysaccaride (LPS) encoding region of the bacteria in both culture systems as indirect but key parameters to the infection potential of C. burnetii. Propidium monoazide (PMA) treatment-based real-time PCR was used for enumeration of viable C. burnetii which were validated by fluorescent infectious focus forming unit counting assays. Furthermore, RNA isolated from C. burnetiipropagated in both the culture systems was examined for LPS-related gene expression. All thus far known LPS-related genes were found to be expressed in early passages in both culturing systems indicating the presence of predominantly the phase I form of C. burnetii. Finally, we used immune-competent mice to provide direct evidence, that the relative virulence of different C. burnetii strains is essentially the same for both axenic and cell-based methods of propagation.