Comparative genomics reveals close genetic relationships between phages from dairy bacteria and pathogenic Streptococci: evolutionary implications for prophage-host interactions

Complete genome sequencing and comparative genomic analysis of three donkey Streptococcus equi subsp. equi isolates

Introduction

Streptococcus equi subspecies equi (S. equi) is the causative agent of strangles, which is one of the most common and highly contagious respiratory infectious illnesses in horses. Streptococcus equi (S. equi) is a horse-specific pathogen that originated from the closely related zoonotic pathogen Streptococcus equi subspecies zooepidemicus (S. zooepidemicus). Despite decades of research, the movement of genetic material across host-restricted diseases remains a mystery.

Methods

Three S. equi donkey isolates (HTP133, HTP232, and HT1112) were recently isolated from a strangles epidemic on donkey farms in China's Xinjiang Province. In this study, we performed a comprehensive comparative analysis of these isolates using whole genome sequencing and compared them to the published genomic sequences of equine strain S. equi 4047 to uncover evidence of genetic events that shaped the evolution of these donkey S. equi isolates' genomes.

Results

Whole genome sequencing indicated that both strains were closely related, with comparable gene compositions and a high rate of shared core genomes (1788-2004). Our comparative genomic study indicated that the genome structure is substantially conserved across three donkey strains; however, there are several rearrangements and inversions when compared to the horse isolate S. equi 4047. The virulence factors conveyed by genomic islands and prophages, in particular, played a key role in shaping the pathogenic capacity and genetic diversity of these S. equi strains. Furthermore, we discovered that the HT133 isolate had a strong colonization ability and increased motility; the HT1112 isolates had a significantly higher ability for antimicrobial resistance and biofilm formation, and the HT232 isolate gained pathogenic specialization by acquiring a bacteriophage encoding hyaluronate lyase.

Discussion

In summary, our findings show that genetic exchange across S. equi strains influences the development of the donkey S. equi genome, offering important genetic insights for future epidemiological studies of S. equi infection.

Temporal Transcriptional Responses of a Vibrio alginolyticus Strain to Podoviridae Phage HH109 Revealed by RNA-Seq

Vibrio alginolyticus is a common opportunistic pathogen that causes mass mortality in cultured marine animals. Phage HH109 lyses pathogenic V. alginolyticus strain E110 with high efficiency and thus serves as a useful model to understand the dynamic interplay of a phage and its host.

Prophage integration into CRISPR loci enables evasion of antiviral immunity in Streptococcus pyogenes

CRISPR loci are composed of short DNA repeats separated by sequences, known as spacers, that match the genomes of invaders such as phages and plasmids. Spacers are transcribed and processed to generate RNA guides used by CRISPR-associated nucleases to recognize and destroy the complementary nucleic acids of invaders. To counteract this defence, phages can produce small proteins that inhibit these nucleases, termed anti-CRISPRs (Acrs). Here we demonstrate that the ΦAP1.1 temperate phage utilizes an alternative approach to antagonize the type II-A CRISPR response in Streptococcus pyogenes. Immediately after infection, this phage expresses a small anti-CRISPR protein, AcrIIA23, that prevents Cas9 function, allowing ΦAP1.1 to integrate into the direct repeats of the CRISPR locus, neutralizing immunity. However, acrIIA23 is not transcribed during lysogeny and phage integration/excision cycles can result in the deletion and/or transduction of spacers, enabling a complex modulation of the type II-A CRISPR immune response. A bioinformatic search identified prophages integrated not only in the CRISPR repeats, but also the cas genes, of diverse bacterial species, suggesting that prophage disruption of the CRISPR-cas locus is a recurrent mechanism to counteract immunity.

Cryptic prophages within a Streptococcus pyogenes genotype emm4 lineage

The major human pathogen Streptococcus pyogenes shares an intimate evolutionary history with mobile genetic elements, which in many cases carry genes encoding bacterial virulence factors. During recent whole-genome sequencing of a longitudinal sample of S. pyogenes isolates in England, we identified a lineage within emm4 that clustered with the reference genome MEW427. Like MEW427, this lineage was characterized by substantial gene loss within all three prophage regions, compared to MGAS10750 and isolates outside of the MEW427-like lineage. Gene loss primarily affected lysogeny, replicative and regulatory modules, and to a lesser and more variable extent, structural genes. Importantly, prophage-encoded superantigen and DNase genes were retained in all isolates. In isolates where the prophage elements were complete, like MGAS10750, they could be induced experimentally, but not in MEW427-like isolates with degraded prophages. We also found gene loss within the chromosomal island SpyCIM4 of MEW427-like isolates, although surprisingly, the SpyCIM4 element could not be experimentally induced in either MGAS10750-like or MEW427-like isolates. This did not, however, appear to abolish expression of the mismatch repair operon, within which this element resides. The inclusion of further emm4 genomes in our analyses ratified our observations and revealed an international emm4 lineage characterized by prophage degradation. Intriguingly, the USA population of emm4 S. pyogenes appeared to constitute predominantly MEW427-like isolates, whereas the UK population comprised both MEW427-like and MGAS10750-like isolates. The degraded and cryptic nature of these elements may have important phenotypic and fitness ramifications for emm4 S. pyogenes, and the geographical distribution of this lineage raises interesting questions on the population dynamics of the genotype.

Whole genome sequencing and de novo assembly of three virulent Indian isolates of Leptospira

Leptospirosis is a re-emerging bacterial zoonosis caused by pathogenic Leptospira, with a worldwide distribution and becoming a major public health concern. Prophylaxis of this disease is difficult due to several factors such as non-specific variable clinical manifestation, presence of a large number of serovar, species and asymptomatic reservoir hosts, lack of proper diagnostics and vaccines. Despite its global importance and severity of the disease, knowledge about the molecular mechanism of pathogenesis and evolution of pathogenic species of Leptospira remains limited. In this study, we sequenced and analyzed three highly pathogenic species of Indian isolates of Leptospira (interrogans, santarosai, and kirschneri). Additionally, we identified some virulence-related and CRISPR-Cas genes. The virulent analysis showed 232 potential virulence factors encoding proteins in L. interrogans strain Salinem and L. santarosai strain M-4 genome. While the genome of L. kirschneri strain Wumalasena was predicted to encode 198 virulence factor proteins. The variant calling analysis revealed 1151, 19,786, and 22,996 single nucleotide polymorphisms (SNPs) for L. interrogans strain Salinem, L. kirschneri strain Wumalasena and L. santarosai strain M-4, respectively, with a maximum of 5315 missense and 12,221 synonymous mutations for L. santarosai strain M-4. The structural analyses of genomes indicated potential evidence of inversions and structural rearrangment in all three genomes. The availability of these genome sequences and in silico analysis of Leptospira will provide a basis for a deeper understanding of their molecular diversity and pathogenesis mechanism, and further pave a way towards proper management of the disease.

Viral composition and context in metagenomes from biofilm and suspended growth municipal wastewater treatment plants

Wastewater treatment plants (WWTPs) contain high density and diversity of viruses which can significantly impact microbial communities in aquatic systems. While previous studies have investigated viruses in WWTP samples that have been specifically concentrated for viruses and filtered to exclude bacteria, little is known about viral communities associated with bacterial communities throughout wastewater treatment systems. Additionally, differences in viral composition between attached and suspended growth wastewater treatment bioprocesses are not well characterized. Here, shotgun metagenomics was used to analyse wastewater and biomass from transects through two full-scale WWTPs for viral composition and associations with bacterial hosts. One WWTP used a suspended growth activated sludge bioreactor and the other used a biofilm reactor (trickling filter). Myoviridae, Podoviridae and Siphoviridae were the dominant viral families throughout both WWTPs, which are all from the order Caudovirales. Beta diversity analysis of viral sequences showed that samples clustered significantly both by plant and by specific sampling location. For each WWTP, the overall bacterial community structure was significantly different than community structure of bacterial taxa associated with viral sequences. These findings highlight viral community composition in transects through different WWTPs and provide context for dsDNA viral sequences in bacterial communities from these systems.

Natural Occurrence of Escherichia coli-Infecting Bacteriophages in Clinical Samples

The interaction between bacteriophages, bacteria and the human host as a tripartite system has recently captured attention. The taxonomic diversity of bacteriophages, as a natural parasite of bacteria, still remains obscure in human body biomes, representing a so-called “viral dark matter.” Here, we isolated and characterized coliphages from blood, urine and tracheal aspirates samples collected at a tertiary care hospital in Austria. Phages were more often isolated from blood, followed by urine and tracheal aspirates. Phylogenetic analysis and genome comparisons allowed the identification of phages belonging to the Tunavirinae subfamily, and to the Peduovirus and Tequintavirus genera. Tunavirinae phages cluster together and are found in samples from 14 patients, suggesting their prevalence across a variety of human samples. When compared with other phage genomes, the highest similarity level was at 87.69% average nucleotide identity (ANI), which suggests that these are in fact a newly isolated phage species. Tequintavirus phages share a 95.90% with phage 3_29, challenging the ANI threshold currently accepted to differentiate phage species. The isolated phages appear to be virulent, with the exception of the Peduovirus members, which are integrative and seem to reside as prophages in bacterial genomes.

The Bacteriophages of Streptococcus pyogenes

The bacteriophages of Streptococcus pyogenes (group A streptococcus) play a key role in population shaping, genetic transfer, and virulence of this bacterial pathogen. Lytic phages like A25 can alter population distributions through elimination of susceptible serotypes but also serve as key mediators for genetic transfer of virulence genes and antibiotic resistance via generalized transduction. The sequencing of multiple S. pyogenes genomes has uncovered a large and diverse population of endogenous prophages that are vectors for toxins and other virulence factors and occupy multiple attachment sites in the bacterial genomes. Some of these sites for integration appear to have the potential to alter the bacterial phenotype through gene disruption. Remarkably, the phage-like chromosomal islands (SpyCI), which share many characteristics with endogenous prophages, have evolved to mediate a growth-dependent mutator phenotype while acting as global transcriptional regulators. The diverse population of prophages appears to share a large pool of genetic modules that promotes novel combinations that may help disseminate virulence factors to different subpopulations of S. pyogenes. The study of the bacteriophages of this pathogen, both lytic and lysogenic, will continue to be an important endeavor for our understanding of how S. pyogenes continues to be a significant cause of human disease.

Bacterial 'Grounded' Prophages: Hotspots for Genetic Renovation and Innovation

Bacterial genomes are highly plastic allowing the generation of variants through mutations and acquisition of genetic information. The fittest variants are then selected by the econiche thereby allowing the bacterial adaptation and colonization of the habitat. Larger genomes, however, may impose metabolic burden and hence bacterial genomes are optimized by the loss of frivolous genetic information. The activity of temperate bacteriophages has acute consequences on the bacterial population as well as the bacterial genome through lytic and lysogenic cycles. Lysogeny is a selective advantage as the prophage provides immunity to the lysogen against secondary phage attack. Since the non-lysogens are eliminated by the lytic phages, lysogens multiply and colonize the habitat. Nevertheless, all lysogens have an imminent risk of lytic cycle activation and cell lysis. However, a mutation in the attachment sites or in the genes that encode the specific recombinase responsible for prophage excision could result in 'grounding' of the prophage. Since the lysogens with grounded prophage are immune to respective phage infection as well as dodge the induction of lytic cycle, we hypothesize that the selection of these mutant lysogens is favored relative to their normal lysogenic counterparts. These grounded prophages offer several advantages to the bacterial genome evolution through propensity for genetic variations including inversions, deletions, and insertions via horizontal gene transfer. We propose that the grounded prophages expedite bacterial genome evolution by acting as 'genetic buffer zones' thereby increasing the frequency as well as the diversity of variations on which natural selection favors the beneficial variants. The grounded prophages are also hotspots for horizontal gene transfer wherein several ecologically significant genes such as those involved in stress tolerance, antimicrobial resistance, and novel metabolic pathways, are integrated. Moreover, the high frequency of genetic changes within prophages also allows proportionate probability for the de novo genesis of genetic information. Through sequence analyses of well-characterized E. coli prophages we exemplify various roles of grounded prophages in E. coli ecology and evolution. Therefore, the temperate prophages are one of the most significant drivers of bacterial genome evolution and sites of biogenesis of genetic information.

Role of Temperate Bacteriophage ϕ20617 on Streptococcus thermophilus DSM 20617T Autolysis and Biology

Streptococcus thermophilus DSM 20167^T showed autolytic behavior when cultured in lactose- and sucrose-limited conditions. The amount of cell lysis induced was inversely related to the energetic status of the cells, as demonstrated by exposing cells to membrane-uncoupling and glycolysis inhibitors. Genome sequence analysis of strain DSM 20617^T revealed the presence of a pac-type temperate bacteriophage, designated Φ20617, whose genomic organization and structure resemble those of temperate streptococcal bacteriophages. The prophage integrated at the 3'-end of the gene encoding the glycolytic enzyme enolase (eno), between eno and the lipoteichoic acid synthase-encoding gene ltaS, affecting their transcription. Comparative experiments conducted on the wild-type strain and a phage-cured derivative strain revealed that the cell-wall integrity of the lysogenic strain was compromised even in the absence of detectable cell lysis. More importantly, adhesion to solid surfaces and heat resistance were significantly higher in the lysogenic strain than in the phage-cured derivative. The characterization of the phenotype of a lysogenic S. thermophilus and its phage-cured derivative is relevant to understanding the ecological constraints that drive the stable association between a temperate phage and its bacterial host.

Genomic analysis of Acinetobacter baumannii prophages reveals remarkable diversity and suggests profound impact on bacterial virulence and fitness

The recent nomination by the World Health Organization of Acinetobacter baumannii as the number one priority pathogen for the development of new antibiotics is a direct consequence of its fast evolution of pathogenicity, and in particular of multidrug resistance. While the development of new antibiotics is critical, understanding the mechanisms behind the crescent bacterial antibiotic resistance is equally relevant. Often, resistance and other bacterial virulence elements are contained on highly mobile pieces of DNA that can easily spread to other bacteria. Prophages are one of the mediators of this form of gene transfer, and have been frequently found in bacterial genomes, often offering advantageous features to the host. Here we assess the contribution of prophages for the evolution of A. baumannii pathogenicity. We found prophages to be notably diverse and widely disseminated in A. baumannii genomes. Also remarkably, A. baumannii prophages encode for multiple putative virulence factors that may be implicated in the bacterium’s capacity to colonize host niches, evade the host immune system, subsist in unfavorable environments, and tolerate antibiotics. Overall our results point towards a significant contribution of prophages for the dissemination and evolution of pathogenicity in A. baumannii, and highlight their clinical relevance.

Type II toxin/antitoxin system ParESO /CopASO stabilizes prophage CP4So in Shewanella oneidensis

Toxin/antitoxin (TA) loci are commonly found in mobile genetic elements such as plasmids and prophages. However, the physiological functions of these TA loci in prophages and cross-regulation among these TA loci remain largely unexplored. Here, we characterized a newly discovered type II TA pair, ParE_SO /CopA_SO , in the CP4So prophage in Shewanella oneidensis. We demonstrated that ParE_SO /CopA_SO plays a critical role in the maintenance of CP4So in host cells after its excision. The toxin ParE_SO inhibited cell growth, resulting in filamentous growth and eventually cell death. The antitoxin CopA_SO neutralized the toxicity of ParE_SO through direct protein-protein interactions and repressed transcription of the TA operon by binding to a DNA motif in the promoter region containing two inverted repeats [5'-GTANTAC (N)₃ GTANTAC-3']. CopA_SO also repressed transcription of another TA system PemK_SO /PemI_SO in megaplasmid pMR-1 of S. oneidensis through binding to a highly similar DNA motif in its promoter region. CopA_SO homologs are widely spread in Shewanella and other Proteobacteria, either as a component of a TA pair or as orphan antitoxins. Our study thus illustrated the cross-regulation of the TA systems in different mobile genetic elements and expanded our understanding of the physiological function of TA systems.

Why Bacteriophage Encode Exotoxins and other Virulence Factors

This study considers gene location within bacteria as a function of genetic element mobility. Our emphasis is on prophage encoding of bacterial virulence factors (VFs). At least four mechanisms potentially contribute to phage encoding of bacterial VFs: (i) Enhanced gene mobility could result in greater VF gene representation within bacterial populations. We question, though, why certain genes but not others might benefit from this mobility. (ii) Epistatic interactions—between VF genes and phage genes that enhance VF utility to bacteria—could maintain phage genes via selection acting on individual, VF-expressing bacteria. However, is this mechanism sufficient to maintain the rest of phage genomes or, without gene co-regulation, even genetic linkage between phage and VF genes? (iii) Phage could amplify VFs during disease progression by carrying them to otherwise commensal bacteria colocated within the same environment. However, lytic phage kill bacteria, thus requiring assumptions of inclusive fitness within bacterial populations to explain retention of phage-mediated VF amplification for the sake of bacterial utility. Finally, (iv) phage-encoded VFs could enhance phage Darwinian fitness, particularly by acting as ecosystem-modifying agents. That is, VF-supplied nutrients could enhance phage growth by increasing the density or by improving the physiology of phage-susceptible bacteria. Alternatively, VF-mediated break down of diffusion-inhibiting spatial structure found within the multicellular bodies of host organisms could augment phage dissemination to new bacteria or to environments. Such phage-fitness enhancing mechanisms could apply particularly given VF expression within microbiologically heterogeneous environments, ie, ones where phage have some reasonable potential to acquire phage-susceptible bacteria.

Global Survey and Genome Exploration of Bacteriophages Infecting the Lactic Acid Bacterium Streptococcus thermophilus

Despite the persistent and costly problem caused by (bacterio)phage predation of Streptococcus thermophilus in dairy plants, DNA sequence information relating to these phages remains limited. Genome sequencing is necessary to better understand the diversity and proliferative strategies of virulent phages. In this report, whole genome sequences of 40 distinct bacteriophages infecting S. thermophilus were analyzed for general characteristics, genomic structure and novel features. The bacteriophage genomes display a high degree of conservation within defined groupings, particularly across the structural modules. Supporting this observation, four novel members of a recently discovered third group of S. thermophilus phages (termed the 5093 group) were found to be conserved relative to both phage 5093 and to each other. Replication modules of S. thermophilus phages generally fall within two main groups, while such phage genomes typically encode one putative transcriptional regulator. Such features are indicative of widespread functional synteny across genetically distinct phage groups. Phage genomes also display nucleotide divergence between groups, and between individual phages of the same group (within replication modules and at the 3′ end of the lysis module)—through various insertions and/or deletions. A previously described multiplex PCR phage detection system was updated to reflect current knowledge on S. thermophilus phages. Furthermore, the structural protein complement as well as the antireceptor (responsible for the initial attachment of the phage to the host cell) of a representative of the 5093 group was defined. Our data more than triples the currently available genomic information on S. thermophilus phages, being of significant value to the dairy industry, where genetic knowledge of lytic phages is crucial for phage detection and monitoring purposes. In particular, the updated PCR detection methodology for S. thermophilus phages is highly useful in monitoring particular phage group(s) present in a given whey sample. Studies of this nature therefore not only provide information on the prevalence and associated threat of known S. thermophilus phages, but may also uncover newly emerging and genomically distinct phages infecting this dairy starter bacterium.

Identification and Analysis of a Novel Group of Bacteriophages Infecting the Lactic Acid Bacterium Streptococcus thermophilus

We present the complete genome sequences of four members of a novel group of phages infecting Streptococcus thermophilus, designated here as the 987 group. Members of this phage group appear to have resulted from genetic exchange events, as evidenced by their "hybrid" genomic architecture, exhibiting DNA sequence relatedness to the morphogenesis modules of certain P335 group Lactococcus lactis phages and to the replication modules of S. thermophilus phages. All four identified members of the 987 phage group were shown to elicit adsorption affinity to both their cognate S. thermophilus hosts and a particular L. lactis starter strain. The receptor binding protein of one of these phages (as a representative of this novel group) was defined using an adsorption inhibition assay. The emergence of a novel phage group infecting S. thermophilus highlights the continuous need for phage monitoring and development of new phage control measures.

IMPORTANCE

Phage predation of S. thermophilus is an important issue for the dairy industry, where viral contamination can lead to fermentation inefficiency or complete fermentation failure. Genome information and phage-host interaction studies of S. thermophilus phages, particularly those emerging in the marketplace, are an important part of limiting the detrimental impact of these viruses in the dairy environment.

Bacteriophage Therapy: Advances in Formulation Strategies and Human Clinical Trials

Recently, a number of phage therapy phase I and II safety trials have been concluded, showing no notable safety concerns associated with the use of phage. Though hurdles for efficient treatment remain, these trials hold promise for future phase III clinical trials. Interestingly, most phage formulations used in these clinical trials are straightforward phage suspensions, and not much research has focused on the processing of phage cocktails in specific pharmaceutical dosage forms. Additional research on formulation strategies and the stability of phage-based drugs will be of key importance, especially with phage therapy advancing toward phase III clinical trials.

Isolation, genome sequencing and functional analysis of two T7-like coliphages of avian pathogenic Escherichia coli

Avian pathogenic Escherichia coli (APEC) causes colibacillosis, which results in significant economic losses to the poultry industry worldwide. Due to the drug residues and increased antibiotic resistance caused by antibiotic use, bacteriophages and other alternative therapeutic agents are expected to control APEC infection in poultry. Two APEC phages, named P483 and P694, were isolated from the feces from the farmers market in China. We then studied their biological properties, and carried out high-throughput genome sequencing and homology analyses of these phages. Assembly results of high-throughput sequencing showed that the structures of both P483 and P694 genomes consist of linear and double-stranded DNA. Results of the electron microscopy and homology analysis revealed that both P483 and P694 belong to T7-like virus which is a member of the Podoviridae family of the Caudovirales order. Comparative genomic analysis showed that most of the predicted proteins of these two phages showed strongest sequence similarity to the Enterobacteria phages BA14 and 285P, Erwinia phage FE44, and Kluyvera phage Kvp1; however, some proteins such as gp0.6a, gp1.7 and gp17 showed lower similarity (<85%) with the homologs of other phages in the T7 subgroup. We also found some unique characteristics of P483 and P694, such as the two types of the genes of P694 and no lytic activity of P694 against its host bacteria in liquid medium. Our results serve to further our understanding of phage evolution of T7-like coliphages and provide the potential application of the phages as therapeutic agents for the treatment of diseases.

Bacteriophages affect evolution of bacterial communities in spatially distributed habitats: a simulation study

Background

Bacteriophages are known to be one of the driving forces of bacterial evolution. Besides promoting horizontal transfer of genes between cells, they may induce directional selection of cells (for instance, according to more or less resistance to phage infection). Switching between lysogenic and lytic pathways results in various types of (co)evolution in host-phage systems. Spatial (more generally, ecological) organization of the living environment is another factor affecting evolution. In this study, we have simulated and analyzed a series of computer models of microbial communities evolving in spatially distributed environments under the pressure of phage infection.

Results

We modeled evolving microbial communities living in spatially distributed flowing environments. Non-specific nutrient supplied in the only spatial direction, resulting in its non-uniform distribution in environment. We varied the time and the location of initial phage infestation of cells as well as switched chemotaxis on and off. Simulations were performed with the Haploid evolutionary constructor software (http://evol-constructor.bionet.nsc.ru/).

Conclusion

Simulations have shown that the spatial location of initial phage invasion may lead to different evolutionary scenarios. Phage infection decreases the speciation rate by more than one order as far as intensified selection blocks the origin of novel viable populations/species, which could carve out potential ecological niches. The dependence of speciation rate on the invasion node location varied on the time of invasion. Speciation rate was found to be lower when the phage invaded fully formed community of sedentary cells (at middle and late times) at the species-rich regions. This is especially noticeable in the case of late-time invasion.

Our simulation study has shown that phage infection affects evolution of microbial community slowing down speciation and stabilizing the system as a whole. This influencing varied in its efficiency depending on spatially-ecological factors as well as community state at the moment of phage invasion.

Electronic supplementary material

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

Targeted Curing of All Lysogenic Bacteriophage from Streptococcus pyogenes Using a Novel Counter-selection Technique

Streptococcus pyogenes is a human commensal and a bacterial pathogen responsible for a wide variety of human diseases differing in symptoms, severity, and tissue tropism. The completed genome sequences of >37 strains of S. pyogenes, representing diverse disease-causing serotypes, have been published. The greatest genetic variation among these strains is attributed to numerous integrated prophage and prophage-like elements, encoding several virulence factors. A comparison of isogenic strains, differing in prophage content, would reveal the effects of these elements on streptococcal pathogenesis. However, curing strains of prophage is often difficult and sometimes unattainable. We have applied a novel counter-selection approach to identify rare S. pyogenes mutants spontaneously cured of select prophage. To accomplish this, we first inserted a two-gene cassette containing a gene for kanamycin resistance (Kan^R) and the rpsL wild-type gene, responsible for dominant streptomycin sensitivity (Sm^S), into a targeted prophage on the chromosome of a streptomycin resistant (Sm^R) mutant of S. pyogenes strain SF370. We then applied antibiotic counter-selection for the re-establishment of the Kan^S/Sm^R phenotype to select for isolates cured of targeted prophage. This methodology allowed for the precise selection of spontaneous phage loss and restoration of the natural phage attB attachment sites for all four prophage-like elements in this S. pyogenes chromosome. Overall, 15 mutants were constructed that encompassed every permutation of phage knockout as well as a mutant strain, named CEM1ΔΦ, completely cured of all bacteriophage elements (a ~10% loss of the genome); the only reported S. pyogenes strain free of prophage-like elements. We compared CEM1ΔΦ to the WT strain by analyzing differences in secreted DNase activity, as well as lytic and lysogenic potential. These mutant strains should allow for the direct examination of bacteriophage relationships within S. pyogenes and further elucidate how the presence of prophage may affect overall streptococcal survival, pathogenicity, and evolution.

Elimination of Chromosomal Island SpyCIM1 from Streptococcus pyogenes Strain SF370 Reverses the Mutator Phenotype and Alters Global Transcription

Streptococcus pyogenes chromosomal island M1 (SpyCIM1) integrates by site-specific recombination into the 5’ end of DNA mismatch repair (MMR) gene mutL in strain SF370SmR, blocking transcription of it and the downstream operon genes. During exponential growth, SpyCIM1 excises from the chromosome and replicates as an episome, restoring mutL transcription. This process is reversed in stationary phase with SpyCIM1 re-integrating into mutL, returning the cells to a mutator phenotype. Here we show that elimination of SpyCIM1 relieves this mutator phenotype. The downstream MMR operon genes, multidrug efflux pump lmrP, Holliday junction resolution helicase ruvA, and DNA base excision repair glycosylase tag, are also restored to constitutive expression by elimination of SpyCIM1. The presence of SpyCIM1 alters global transcription patterns in SF370SmR. RNA sequencing (RNA-Seq) demonstrated that loss of SpyCIM1 in the SpyCIM1 deletion mutant, CEM1Δ4, impacted the expression of over 100 genes involved in virulence and metabolism both in early exponential phase, when the SpyCIM1 is episomal, as well as at the onset of stationary phase, when SpyCIM1 has reintegrated into mutL. Among these changes, the up-regulation of the genes for the antiphagocytic M protein (emm1), streptolysin O (slo), capsule operon (hasABC), and streptococcal pyrogenic exotoxin (speB), are particularly notable. The expression pattern of the MMR operon confirmed our earlier observations that these genes are transcribed in early exponential phase but silenced as stationary phase is approached. Thus, the direct role of SpyCIM1 in causing the mutator phenotype is confirmed, and further, its influence upon the biology of S. pyogenes was found to impact multiple genes in addition to the MMR operon, which is a novel function for a mobile genetic element. We suggest that such chromosomal islands are a remarkable evolutionary adaptation to promote the survival of its S. pyogenes host cell in changing environments.

Identification and characterization of cis- and trans-acting elements involved in prophage induction in Streptococcus thermophilus J34

The genetic switch region of temperate Streptococcus thermophilus phage TP-J34 contains two divergently oriented promoters and several predicted operator sites. It separates lytic cycle-promoting genes from those promoting lysogeny. A polycistronic transcript comprises the genes coding for repressor Crh, metalloproteinase-motif protein Rir and superinfection exclusion lipoprotein Ltp. Weak promoters effecting monocistronic transcripts were localized for ltp and int (encoding integrase) by Northern blot and 5'-RACE-PCR. These transcripts appeared in lysogenic as well as lytic state. A polycistronic transcript comprising genes coh (encoding Cro homolog), ant (encoding putative antirepressor), orf7, orf8 and orf9 was only detected in the lytic state. Four operator sites, of which three were located in the intergenic regions between crh and coh, and one between coh and ant, were identified by competition electromobility shift assays. Cooperative binding of Crh to two operator sites immediately upstream of coh could be demonstrated. Coh was shown to bind to the operator closest to crh only. Oligomerization was proven by cross-linking Crh by glutaraldehyde. Knock-out of rir revealed a key role in prophage induction. Rir and Crh were shown to form a complex in solution and Rir prevented binding of Crh to its operator sites.

Molecular epidemiology and genomics of group A Streptococcus

Streptococcus pyogenes (group A Streptococcus; GAS) is a strict human pathogen with a very high prevalence worldwide. This review highlights the genetic organization of the species and the important ecological considerations that impact its evolution. Recent advances are presented on the topics of molecular epidemiology, population biology, molecular basis for genetic change, genome structure and genetic flux, phylogenomics and closely related streptococcal species, and the long- and short-term evolution of GAS. The application of whole genome sequence data to addressing key biological questions is discussed.

Life-style and genome structure of marine Pseudoalteromonas siphovirus B8b isolated from the northwestern Mediterranean Sea

Marine viruses (phages) alter bacterial diversity and evolution with impacts on marine biogeochemical cycles, and yet few well-developed model systems limit opportunities for hypothesis testing. Here we isolate phage B8b from the Mediterranean Sea using Pseudoalteromonas sp. QC-44 as a host and characterize it using myriad techniques. Morphologically, phage B8b was classified as a member of the Siphoviridae family. One-step growth analyses showed that this siphovirus had a latent period of 70 min and released 172 new viral particles per cell. Host range analysis against 89 bacterial host strains revealed that phage B8b infected 3 Pseudoalteromonas strains (52 tested, >99.9% 16S rRNA gene nucleotide identity) and 1 non-Pseudoaltermonas strain belonging to Alteromonas sp. (37 strains from 6 genera tested), which helps bound the phylogenetic distance possible in a phage-mediated horizontal gene transfer event. The Pseudoalteromonas phage B8b genome size was 42.7 kb, with clear structural and replication modules where the former were delineated leveraging identification of 16 structural genes by virion structural proteomics, only 4 of which had any similarity to known structural proteins. In nature, this phage was common in coastal marine environments in both photic and aphotic layers (found in 26.5% of available viral metagenomes), but not abundant in any sample (average per sample abundance was 0.65% of the reads). Together these data improve our understanding of siphoviruses in nature, and provide foundational information for a new ‘rare virosphere’ phage–host model system.

Lytic activity by temperate phages of Pseudomonas aeruginosa in long-term cystic fibrosis chronic lung infections

Pseudomonas aeruginosa is the most common bacterial pathogen infecting the lungs of cystic fibrosis (CF) patients. The transmissible Liverpool epidemic strain (LES) harbours multiple inducible prophages (LESϕ2; LESϕ3; LESϕ4; LESϕ5; and LESϕ6), some of which are known to confer a competitive advantage in an in vivo rat model of chronic lung infection. We used quantitative PCR (Q-PCR) to measure the density and dynamics of all five LES phages in the sputa of 10 LES-infected CF patients over a period of 2 years. In all patients, the densities of free-LES phages were positively correlated with the densities of P. aeruginosa, and total free-phage densities consistently exceeded bacterial host densities 10–100-fold. Further, we observed a negative correlation between the phage-to-bacterium ratio and bacterial density, suggesting a role for lysis by temperate phages in regulation of the bacterial population densities. In 9/10 patients, LESϕ2 and LESϕ4 were the most abundant free phages, which reflects the differential in vitro induction properties of the phages. These data indicate that temperate phages of P. aeruginosa retain lytic activity after prolonged periods of chronic infection in the CF lung, and suggest that temperate phage lysis may contribute to regulation of P. aeruginosa density in vivo.

Chromosomal islands of Streptococcus pyogenes and related streptococci: molecular switches for survival and virulence

Streptococcus pyogenes is a significant pathogen of humans, annually causing over 700,000,000 infections and 500,000 deaths. Virulence in S. pyogenes is closely linked to mobile genetic elements like phages and chromosomal islands (CI). S. pyogenes phage-like chromosomal islands (SpyCI) confer a complex mutator phenotype on their host. SpyCI integrate into the 5′ end of DNA mismatch repair (MMR) gene mutL, which also disrupts downstream operon genes lmrP, ruvA, and tag. During early logarithmic growth, SpyCI excise from the bacterial chromosome and replicate as episomes, relieving the mutator phenotype. As growth slows and the cells enter stationary phase, SpyCI reintegrate into the chromosome, again silencing the MMR operon. This system creates a unique growth-dependent and reversible mutator phenotype. Additional CI using the identical attachment site in mutL have been identified in related species, including Streptococcus dysgalactiae subsp. equisimilis, Streptococcus anginosus, Streptococcus intermedius, Streptococcus parauberis, and Streptococcus canis. These CI have small genomes, which range from 13 to 20 kB, conserved integrase and DNA replication genes, and no identifiable genes encoding capsid proteins. SpyCI may employ a helper phage for packaging and dissemination in a fashion similar to the Staphylococcus aureus pathogenicity islands (SaPI). Outside of the core replication and integration genes, SpyCI and related CI show considerable diversity with the presence of many indels that may contribute to the host cell phenotype or fitness. SpyCI are a subset of a larger family of streptococcal CI who potentially regulate the expression of other host genes. The biological and phylogenetic analysis of streptococcal chromosomal islands provides important clues as to how these chromosomal islands help S. pyogenes and other streptococcal species persist in human populations in spite of antibiotic therapy and immune challenges.

Molecular bases and role of viruses in the human microbiome

Viruses are dependent biological entities that interact with the genetic material of most cells on the planet, including the trillions within the human microbiome. Their tremendous diversity renders analysis of human viral communities ("viromes") to be highly complex. Because many of the viruses in humans are bacteriophage, their dynamic interactions with their cellular hosts add greatly to the complexities observed in examining human microbial ecosystems. We are only beginning to be able to study human viral communities on a large scale, mostly as a result of recent and continued advancements in sequencing and bioinformatic technologies. Bacteriophage community diversity in humans not only is inexorably linked to the diversity of their cellular hosts but also is due to their rapid evolution, horizontal gene transfers, and intimate interactions with host nucleic acids. There are vast numbers of observed viral genotypes on many body surfaces studied, including the oral, gastrointestinal, and respiratory tracts, and even in the human bloodstream, which previously was considered a purely sterile environment. The presence of viruses in blood suggests that virome members can traverse mucosal barriers, as indeed these communities are substantially altered when mucosal defenses are weakened. Perhaps the most interesting aspect of human viral communities is the extent to which they can carry gene functions involved in the pathogenesis of their hosts, particularly antibiotic resistance. Persons in close contact with each other have been shown to share a fraction of oral virobiota, which could potentially have important implications for the spread of antibiotic resistance to healthy individuals. Because viruses can have a large impact on ecosystem dynamics through mechanisms such as the transfers of beneficial gene functions or the lysis of certain populations of cellular hosts, they may have both beneficial and detrimental roles that affect human health, including improvements in microbial resilience to disturbances, immune evasion, maintenance of physiologic processes, and altering the microbial community in ways that promote or prevent pathogen colonization.

Characterization and whole genome sequences of the Rhodococcus bacteriophages RGL3 and RER2

We report here the isolation and genome sequences of two novel phages, lytic for Rhodococcus and Nocardia species. Named RER2 and RGL3, both are members of the family Siphoviridae, and each possesses a novel genome of 46,586 bp and 48,072 bp, respectively. RER2 and RGL3 phages share a modular genome organization, as seen in other sequenced Siphoviridae phage genomes, and appear to share a common evolutionary origin. The genomes of these phages share no similarity with other Rhodococcus or Nocardia phages but are related to Mycobacterium phages. The data presented here extend our understanding of Rhodococcus phage genomics.

Whole genome analysis of Leptospira licerasiae provides insight into leptospiral evolution and pathogenicity

The whole genome analysis of two strains of the first intermediately pathogenic leptospiral species to be sequenced (Leptospira licerasiae strains VAR010 and MMD0835) provides insight into their pathogenic potential and deepens our understanding of leptospiral evolution. Comparative analysis of eight leptospiral genomes shows the existence of a core leptospiral genome comprising 1547 genes and 452 conserved genes restricted to infectious species (including L. licerasiae) that are likely to be pathogenicity-related. Comparisons of the functional content of the genomes suggests that L. licerasiae retains several proteins related to nitrogen, amino acid and carbohydrate metabolism which might help to explain why these Leptospira grow well in artificial media compared with pathogenic species. L. licerasiae strains VAR010^T and MMD0835 possess two prophage elements. While one element is circular and shares homology with LE1 of L. biflexa, the second is cryptic and homologous to a previously identified but unnamed region in L. interrogans serovars Copenhageni and Lai. We also report a unique O-antigen locus in L. licerasiae comprised of a 6-gene cluster that is unexpectedly short compared with L. interrogans in which analogous regions may include >90 such genes. Sequence homology searches suggest that these genes were acquired by lateral gene transfer (LGT). Furthermore, seven putative genomic islands ranging in size from 5 to 36 kb are present also suggestive of antecedent LGT. How Leptospira become naturally competent remains to be determined, but considering the phylogenetic origins of the genes comprising the O-antigen cluster and other putative laterally transferred genes, L. licerasiae must be able to exchange genetic material with non-invasive environmental bacteria. The data presented here demonstrate that L. licerasiae is genetically more closely related to pathogenic than to saprophytic Leptospira and provide insight into the genomic bases for its infectiousness and its unique antigenic characteristics.

Leptospirosis is one of the most common diseases transmitted by animals worldwide and is important because it is a major cause of febrile illness in tropical areas and also occurs in epidemic form associated with natural disasters and flooding. The mechanisms through which Leptospira cause disease are not well understood. In this study we have sequenced the genomes of two strains of Leptospira licerasiae isolated from a person and a marsupial in the Peruvian Amazon. These strains were thought to be able to cause only mild disease in humans. We have compared these genomes with other leptospires that can cause severe illness and death and another leptospire that does not infect humans or animals. These comparisons have allowed us to demonstrate similarities among the disease-causing Leptospira. Studying genes that are common among infectious strains will allow us to identify genetic factors necessary for infecting, causing disease and determining the severity of disease. We have also found that L. licerasiae seems to be able to uptake and incorporate genetic information from other bacteria found in the environment. This information will allow us to begin to understand how Leptospira species have evolved.

Novel Podoviridae family bacteriophage infecting Weissella cibaria isolated from Kimchi

The first complete genome sequence of a phage infecting Weissella cibaria (Weissella kimchii) is presented. The bacteriophage YS61 was isolated from kimchi, a Korean fermented vegetable dish. Bacteriophages are recognized as a serious problem in industrial fermentations; however, YS61 differed from many virulent phages associated with food fermentations since it was difficult to propagate and was very susceptible to resistance development. Sequence analysis revealed that YS61 resembles Podoviridae of the subfamily Picovirinae. Within the subfamily Picovirinae, the 29-like phages have been extensively studied, and their terminal protein-primed DNA replication is well characterized. Our data strongly suggest that YS61 also replicates by a protein-primed mechanism. Weissella phage YS61 is, however, markedly different from members of the Picovirinae with respect to genome size and morphology. Picovirinae are characterized by small (approximately 20-kb) genomes which contrasts with the 33,594-bp genome of YS61. Based on electron microscopy analysis, YS61 was classified as a member of the Podoviridae of morphotype C2, similar to the 29-like phages, but its capsid dimensions are significantly larger than those reported for these phages. The novelty of YS61 was also emphasized by the low number of open reading frames (ORFs) showing significant similarity to database sequences. We propose that the bacteriophage YS61 should represent a new subfamily within the family Podoviridae.

Tracing lifestyle adaptation in prokaryotic genomes

Lifestyle adaptation of microbes due to changes in their ecological niches or acquisition of new environments is a major driving force for genetic changes in their respective genomes. Moving into more specialized niches often results in the acquisition of new gene sets via horizontal gene transfer to utilize previously unavailable metabolites, while genetic ballast is shed by gene loss and/or gene inactivation. In some cases, larger genome rearrangements can be observed, such as the incorporation of whole genetic islands, providing a range of new phenotypic capabilities. Until recently these changes could not be comprehensively followed and identified due to the lack of complete microbial genome sequences. The advent of high-throughput DNA sequencing has dramatically changed the scientific landscape and today microbial genomes have become increasingly abundant. Currently, more than 2,900 genomes are published and more than 11,000 genome projects are listed in the Genomes Online Database. Although this wealth of information provides many new opportunities to assess microbial functionality, it also creates a new array of challenges when a comparison between multiple microbial genomes is required. Here, functional genome distribution (FGD) is introduced, analyzing the diversity between microbes based on their predicted ORFeome. FGD is therefore a comparative genomics approach, emphasizing the assessments of gene complements. To further facilitate the comparison between two or more genomes, degrees of amino-acid similarities between ORFeomes can be visualized in the Artemis comparison tool, graphically depicting small and large scale genome rearrangements, insertion and deletion events, and levels of similarity between individual open reading frames. FGD provides a new tool for comparative microbial genomics and the interpretation of differences in the genetic makeup of bacteria.

Identification of Rgg binding sites in the Streptococcus pyogenes chromosome

Streptococcus pyogenes Rgg is a regulatory protein that controls the transcription of 588 genes in strain NZ131 during the post-exponential phase of growth, including the virulence-associated genes encoding the extracellular SpeB protease, pullulanase A (PulA), and two extracellular nucleases (SdaB and Spd-3). Rgg binds to DNA proximally to the speB promoter (PspeB) to activate transcription; however, it is not known if Rgg binds to the promoters of other genes to influence expression, or if the perturbation of other global regulons accounts for the genome-wide changes in expression associated with the mutant. To address this issue, chromatin immunoprecipitation followed by DNA microarray analysis (ChIP-chip) was used to identify the DNA binding sites of Rgg. Rgg bound to 65 sites in the chromosome. Thirty-five were within noncoding DNA, and 43% of these were adjacent to genes previously identified as regulated by Rgg. Electrophoretic mobility shift assays were used to assess the binding of Rgg to a subset of sites bound in vivo, including the noncoding DNA upstream of speB, the genes encoding PulA, Spd-3, and a transcriptional regulator (SPY49_1113), and prophage-associated genes encoding a putative integrase (SPY49_0746) and a surface antigen (SPY49_0396). Rgg bound to all target DNAs in vitro, consistent with the in vivo results. Finally, analyses with a transcriptional reporter system showed that the DNA bound by Rgg contained an active promoter that was regulated by Rgg. Overall, the results indicate that Rgg binds specifically to multiple sites in the chromosome, including prophage DNA, to influence gene expression.

Prophage spontaneous activation promotes DNA release enhancing biofilm formation in Streptococcus pneumoniae

Streptococcus pneumoniae (pneumococcus) is able to form biofilms in vivo and previous studies propose that pneumococcal biofilms play a relevant role both in colonization and infection. Additionally, pneumococci recovered from human infections are characterized by a high prevalence of lysogenic bacteriophages (phages) residing quiescently in their host chromosome. We investigated a possible link between lysogeny and biofilm formation. Considering that extracellular DNA (eDNA) is a key factor in the biofilm matrix, we reasoned that prophage spontaneous activation with the consequent bacterial host lysis could provide a source of eDNA, enhancing pneumococcal biofilm development. Monitoring biofilm growth of lysogenic and non-lysogenic pneumococcal strains indicated that phage-infected bacteria are more proficient at forming biofilms, that is their biofilms are characterized by a higher biomass and cell viability. The presence of phage particles throughout the lysogenic strains biofilm development implicated prophage spontaneous induction in this effect. Analysis of lysogens deficient for phage lysin and the bacterial major autolysin revealed that the absence of either lytic activity impaired biofilm development and the addition of DNA restored the ability of mutant strains to form robust biofilms. These findings establish that limited phage-mediated host lysis of a fraction of the bacterial population, due to spontaneous phage induction, constitutes an important source of eDNA for the S. pneumoniae biofilm matrix and that this localized release of eDNA favors biofilm formation by the remaining bacterial population.

Genome analysis of Lactobacillus fermentum temperate bacteriophage ФPYB5

The complete genomic sequence of Lactobacillus fermentum temperate bacteriophage ФPYB5 was determined. The phage possesses a linear, double-stranded DNA genome of 32,847 bp with a G+C content of 45.21%. A total of 46 putative open reading frames (ORFs) were identified. On the basis of homology comparisons, 25 ORFs could be assigned putative functions. The genome of bacteriophage ФPYB5 is highly modular with functionally related genes clustered together. Genome DNA of temperate bacteriophage ФPYB5, induced from heterofermentative lactic acid bacteria, showed to be closely related to that of the prophage of heterofermentative Lactobacillus reuteri 100-23 and heterofermetative Leuconostoc oenos bacteriophage 10MC in an evolutionary aspect.

Molecular characterization and prophage DNA contents of Streptococcus agalactiae strains isolated from adult skin and osteoarticular infections

Skin and osteoarticular infections (SKI and OAI, respectively) account for almost one-third of Streptococcus agalactiae infections in nonpregnant adults. We evaluated the genetic diversity and phylogeny of 58 S. agalactiae strains responsible for adult SKI or OAI and of 61 S. agalactiae strains from cases of adult human colonization (HCol) by serotyping and multilocus sequence typing (MLST). We also assessed the prophage DNA content of the genomes of these strains by a PCR-based method. We found that 63% of SKI and 56% of OAI occurred in people aged 55 years and over. Overall, 71% of SKI strains were of serotype Ia or V, and 91% of OAI strains were of serotype Ia, III, or V. Strains of clonal complexes 1 and 23 (CC1 and CC23) were associated with 79% of SKI cases and 62% of OAI cases. Seven groups of strains, groups A, B, C, D, E, F, and G, were obtained by performing a hierarchical analysis on the basis of prophage DNA-PCR data. We found that 85% of CC1 strains clustered in DNA prophage group D, the group with the highest prophage DNA content (average, 4.4; average of absolute deviations [AVEDEV], 0.9). The CC23 strains displayed the greatest diversity in prophage DNA fragment content, but 47% of CC23 strains clustered in group B, which also had a high average prophage DNA content per strain (average, 2.3; AVEDEV, 0.6). Many (65%) of the OAI strains were in prophage DNA group D, whereas 83% of the SKI strains were in prophage DNA groups B and D. These data suggest that S. agalactiae strains from CC1 and CC23 may be subject to particular transduction mechanisms in gene recombination, rendering them particularly capable of invading the skin, bone, or joints in adults.

The source of self: genetic parasites and the origin of adaptive immunity

Stable colonization of the host by viruses (genetic parasites) can alter the systems of host identity and provide immunity against related viruses. To attain the needed stability, some viruses of prokaryotes (P1 phage) use a strategy called an addiction module. The linked protective and destructive gene functions of an addiction module insures both virus persistence but will also destroy cells that interrupt this module and thereby prevent infection by competitors. Previously, I have generalized this concept to also include persistent and lytic states of virus infection, which can be considered as a virus addiction module. Such states often involve defective viruses. In this report, I examine the origin of the adaptive immune system from the perspective of a virus addiction module. The likely role of both endogenous and exogenous retroviruses, DNA viruses, and their defective elements is considered in the origin of all the basal components of adaptive immunity (T-cell receptor, RAG-mediated gene rearrangement, clonal lymphocyte proliferation, antigen surface presentation, apoptosis, and education of immune cells). It is concluded that colonization by viruses and their defectives provides a more coherent explanation for the origin of adaptive immunity.

Viruses are essential agents within the roots and stem of the tree of life

In contrast with former definitions of life limited to membrane-bound cellular life forms which feed, grow, metabolise and replicate (i) a role of viruses as genetic symbionts, (ii) along with peripheral phenomena such as cryptobiosis and (iii) the horizontal nature of genetic information acquisition and processing broaden our view of the tree of life. Some researchers insist on the traditional textbook conviction of what is part of the community of life. In a recent review [Moreira, D., Lopez-Garcia, P., 2009. Ten reasons to exclude viruses from the tree of life. Nat. Rev. Microbiol. 7, 306-311.] they assemble four main arguments which should exclude viruses from the tree of life because of their inability to self-sustain and self-replicate, their polyphyly, the cellular origin of their cell-like genes and the volatility of their genomes. In this article we will show that these features are not coherent with current knowledge about viruses but that viral agents play key roles within the roots and stem of the tree of life.

Comparative genomic analysis of ten Streptococcus pneumoniae temperate bacteriophages

Streptococcus pneumoniae is an important human pathogen that often carries temperate bacteriophages. As part of a program to characterize the genetic makeup of prophages associated with clinical strains and to assess the potential roles that they play in the biology and pathogenesis in their host, we performed comparative genomic analysis of 10 temperate pneumococcal phages. All of the genomes are organized into five major gene clusters: lysogeny, replication, packaging, morphogenesis, and lysis clusters. All of the phage particles observed showed a Siphoviridae morphology. The only genes that are well conserved in all the genomes studied are those involved in the integration and the lysis of the host in addition to two genes, of unknown function, within the replication module. We observed that a high percentage of the open reading frames contained no similarities to any sequences catalogued in public databases; however, genes that were homologous to known phage virulence genes, including the pblB gene of Streptococcus mitis and the vapE gene of Dichelobacter nodosus, were also identified. Interestingly, bioinformatic tools showed the presence of a toxin-antitoxin system in the phage phiSpn_6, and this represents the first time that an addition system in a pneumophage has been identified. Collectively, the temperate pneumophages contain a diverse set of genes with various levels of similarity among them.

Identification and characterization of lactococcal-prophage-carried superinfection exclusion genes

Superinfection exclusion (Sie) proteins are prophage-encoded phage resistance systems. In this study, genes encoding Sie systems were identified on the genomes of Lactococcus lactis subsp. cremoris MG1363 and L. lactis subsp. lactis IL1403. These Sie systems are genetically distinct and yet were shown to act specifically against a particular subset of the 936 phage group. Each of the systems allows normal phage adsorption while affecting plasmid transduction and intracellular phage DNA replication, which points to the blocking of phage DNA injection as their common mode of action. Sie-specifying genes found on the MG1363 prophages are also present in various lactococcal strains, whereas the prophage-encoded Sie systems of IL1403 do not appear to be as widely disseminated.

Phage-associated mutator phenotype in group A streptococcus

Defects in DNA mismatch repair (MMR) occur frequently in natural populations of pathogenic and commensal bacteria, resulting in a mutator phenotype. We identified a unique genetic element in Streptococcus pyogenes strain SF370 that controls MMR via a dynamic process of prophage excision and reintegration in response to growth. In S. pyogenes, mutS and mutL are organized on a polycistronic mRNA under control of a common promoter. Prophage SF370.4 is integrated between the two genes, blocking expression of the downstream gene (mutL) and resulting in a mutator phenotype. However, in rapidly growing cells the prophage excises and replicates as an episome, allowing mutL to be expressed. Excision of prophage SF370.4 and expression of MutL mRNA occur simultaneously during early logarithmic growth when cell densities are low; this brief window of MutL gene expression ends as the cell density increases. However, detectable amounts of MutL protein remain in the cell until the onset of stationary phase. Thus, MMR in S. pyogenes SF370 is functional in exponentially growing cells but defective when resources are limiting. The presence of a prophage integrated into the 5' end of mutL correlates with a mutator phenotype (10(-7) to 10(-8) mutation/generation, an approximately a 100-fold increase in the rate of spontaneous mutation compared with prophage-free strains [10(-9) to 10(-10) mutation/generation]). Such genetic elements may be common in S. pyogenes since 6 of 13 completed genomes have related prophages, and a survey of 100 strains found that about 20% of them are positive for phages occupying the SF370.4 attP site. The dynamic control of a major DNA repair system by a bacteriophage is a novel method for achieving the mutator phenotype and may allow the organism to respond rapidly to a changing environment while minimizing the risks associated with long-term hypermutability.

Plasmid transfer via transduction from Streptococcus thermophilus to Lactococcus lactis

Using Streptococcus thermophilus phages, plasmid transduction in Lactococcus lactis was demonstrated. The transduction frequencies were 4 orders of magnitude lower in L. lactis than in S. thermophilus. These results are the first evidence that there is phage-mediated direct transfer of DNA from S. thermophilus to L. lactis. The implications of these results for phage evolution are discussed.

Phage_Finder: automated identification and classification of prophage regions in complete bacterial genome sequences

Phage_Finder, a heuristic computer program, was created to identify prophage regions in completed bacterial genomes. Using a test dataset of 42 bacterial genomes whose prophages have been manually identified, Phage_Finder found 91% of the regions, resulting in 7% false positive and 9% false negative prophages. A search of 302 complete bacterial genomes predicted 403 putative prophage regions, accounting for 2.7% of the total bacterial DNA. Analysis of the 285 putative attachment sites revealed tRNAs are targets for integration slightly more frequently (33%) than intergenic (31%) or intragenic (28%) regions, while tmRNAs were targeted in 8% of the regions. The most popular tRNA targets were Arg, Leu, Ser and Thr. Mapping of the insertion point on a consensus tRNA molecule revealed novel insertion points on the 5′ side of the D loop, the 3′ side of the anticodon loop and the anticodon. A novel method of constructing phylogenetic trees of phages and prophages was developed based on the mean of the BLAST score ratio (BSR) of the phage/prophage proteomes. This method verified many known bacteriophage groups, making this a useful tool for predicting the relationships of prophages from bacterial genomes.

Analysis of the collar-whisker structure of temperate lactococcal bacteriophage TP901-1

Proteins homologous to the protein NPS (neck passage structure) are widespread among lactococcal phages. We investigated the hypothesis that NPS is involved in the infection of phage TP901-1 by analysis of an NPS- mutant. NPS was determined to form a collar-whisker complex but was shown to be nonessential for infection, phage assembly, and stability.

Comparative genomics and transcriptional analysis of prophages identified in the genomes of Lactobacillus gasseri, Lactobacillus salivarius, and Lactobacillus casei

Lactobacillus gasseri ATCC 33323, Lactobacillus salivarius subsp. salivarius UCC 118, and Lactobacillus casei ATCC 334 contain one (LgaI), four (Sal1, Sal2, Sal3, Sal4), and one (Lca1) distinguishable prophage sequences, respectively. Sequence analysis revealed that LgaI, Lca1, Sal1, and Sal2 prophages belong to the group of Sfi11-like pac site and cos site Siphoviridae, respectively. Phylogenetic investigation of these newly described prophage sequences revealed that they have not followed an evolutionary development similar to that of their bacterial hosts and that they show a high degree of diversity, even within a species. The attachment sites were determined for all these prophage elements; LgaI as well as Sal1 integrates in tRNA genes, while prophage Sal2 integrates in a predicted arginino-succinate lyase-encoding gene. In contrast, Lca1 and the Sal3 and Sal4 prophage remnants are integrated in noncoding regions in the L. casei ATCC 334 and L. salivarius UCC 118 genomes. Northern analysis showed that large parts of the prophage genomes are transcriptionally silent and that transcription is limited to genome segments located near the attachment site. Finally, pulsed-field gel electrophoresis followed by Southern blot hybridization with specific prophage probes indicates that these prophage sequences are narrowly distributed within lactobacilli.

Detailed genomic analysis of the Wbeta and gamma phages infecting Bacillus anthracis: implications for evolution of environmental fitness and antibiotic resistance

Phage-mediated lysis has been an essential laboratory tool for rapidly identifying Bacillus anthracis for more than 40 years, relying on the gamma phage derivative of a Bacillus cereus prophage called W. The complete genomic sequences of the temperate W phage, referred to as Wbeta, and its lytic variant gamma were determined and found to encode 53 open reading frames each, spanning 40,864 bp and 37,373 bp, respectively. Direct comparison of the genomes showed that gamma evolved through mutations at key loci controlling host recognition, lysogenic growth, and possibly host phenotypic modification. Included are a cluster of point mutations at the gp14 tail fiber locus of gamma, encoding a protein that, when fused to green fluorescent protein, binds specifically to B. anthracis. A large 2,003-bp deletion was also identified at the gamma lysogeny module, explaining its shift from a temperate to a lytic lifestyle. Finally, evidence of recombination was observed at a dicistronic Wbeta locus, encoding putative bacterial cell surface-modifying proteins, replaced in gamma with a locus, likely obtained from a B. anthracis prophage, encoding demonstrable fosfomycin resistance. Reverse transcriptase PCR analysis confirmed strong induction at the dicistronic Wbeta locus and at four other phage loci in B. anthracis and/or B. cereus lysogens. In all, this study represents the first genomic and functional description of two historically important phages and is part of a broader investigation into contributions of phage to the B. anthracis life cycle. Initial findings suggest that lysogeny of B. anthracis promotes ecological adaptation, rather than virulence, as with other gram-positive pathogens.

Identification of the lower baseplate protein as the antireceptor of the temperate lactococcal bacteriophages TP901-1 and Tuc2009

The first step in the infection process of tailed phages is recognition and binding to the host receptor. This interaction is mediated by the phage antireceptor located in the distal tail structure. The temperate Lactococcus lactis phage TP901-1 belongs to the P335 species of the Siphoviridae family, which also includes the related phage Tuc2009. The distal tail structure of TP901-1 is well characterized and contains a double-disk baseplate and a central tail fiber. The structural tail proteins of TP901-1 and Tuc2009 are highly similar, but the phages have different host ranges and must therefore encode different antireceptors. In order to identify the antireceptors of TP901-1 and Tuc2009, a chimeric phage was generated in which the gene encoding the TP901-1 lower baseplate protein (bppL(TP901-1)) was exchanged with the analogous gene (orf53(2009)) of phage Tuc2009. The chimeric phage (TP901-1C) infected the Tuc2009 host strain efficiently and thus displayed an altered host range compared to TP901-1. Genomic analysis and sequencing verified that TP901-1C is a TP901-1 derivative containing the orf53(2009) gene in exchange for bppL(TP901-1); however, a new sequence in the late promoter region was also discovered. Protein analysis confirmed that TP901-1C contains ORF53(2009) and not the lower baseplate protein BppL(TP901-1), and it was concluded that BppL(TP901-1) and ORF53(2009) constitute antireceptor proteins of TP901-1 and Tuc2009, respectively. Electron micrographs revealed altered baseplate morphology of TP901-1C compared to that of the parental phage.

Mosaic prophages with horizontally acquired genes account for the emergence and diversification of the globally disseminated M1T1 clone of Streptococcus pyogenes

The recrudescence of severe invasive group A streptococcal (GAS) diseases has been associated with relatively few strains, including the M1T1 subclone that has shown an unprecedented global spread and prevalence and high virulence in susceptible hosts. To understand its unusual epidemiology, we aimed to identify unique genomic features that differentiate it from the fully sequenced M1 SF370 strain. We constructed DNA microarrays from an M1T1 shotgun library and, using differential hybridization, we found that both M1 strains are 95% identical and that the 5% unique M1T1 clone sequences more closely resemble sequences found in the M3 strain, which is also associated with severe disease. Careful analysis of these unique sequences revealed three unique prophages that we named M1T1.X, M1T1.Y, and M1T1.Z. While M1T1.Y is similar to phage 370.3 of the M1-SF370 strain, M1T1.X and M1T1.Z are novel and encode the toxins SpeA2 and Sda1, respectively. The genomes of these prophages are highly mosaic, with different segments being related to distinct streptococcal phages, suggesting that GAS phages continue to exchange genetic material. Bioinformatic and phylogenetic analyses revealed a highly conserved open reading frame (ORF) adjacent to the toxins in 18 of the 21 toxin-carrying GAS prophages. We named this ORF paratox, determined its allelic distribution among different phages, and found linkage disequilibrium between particular paratox alleles and specific toxin genes, suggesting that they may move as a single cassette. Based on the conservation of paratox and other genes flanking the toxins, we propose a recombination-based model for toxin dissemination among prophages. We also provide evidence that a minor population of the M1T1 clonal isolates have exchanged their virulence module on phage M1T1.Y, replacing it with a different module identical to that found on a related M3 phage. Taken together, the data demonstrate that mosaicism of the GAS prophages has contributed to the emergence and diversification of the M1T1 subclone.