Genome landscapes and bacteriophage codon usage

Across all kingdoms of biological life, protein-coding genes exhibit unequal usage of synonmous codons. Although alternative theories abound, translational selection has been accepted as an important mechanism that shapes the patterns of codon usage in prokaryotes and simple eukaryotes. Here we analyze patterns of codon usage across 74 diverse bacteriophages that infect E. coli, P. aeruginosa and L. lactis as their primary host. We introduce the concept of a `genome landscape,' which helps reveal non-trivial, long-range patterns in codon usage across a genome. We develop a series of randomization tests that allow us to interrogate the significance of one aspect of codon usage, such a GC content, while controlling for another aspect, such as adaptation to host-preferred codons. We find that 33 phage genomes exhibit highly non-random patterns in their GC3-content, use of host-preferred codons, or both. We show that the head and tail proteins of these phages exhibit significant bias towards host-preferred codons, relative to the non-structural phage proteins. Our results support the hypothesis of translational selection on viral genes for host-preferred codons, over a broad range of bacteriophages.Comment: 9 Color Figures, 5 Tables, 53 Reference

Evolutionary Genomics of a Temperate Bacteriophage in an Obligate Intracellular Bacteria (Wolbachia)

Genome evolution of bacteria is usually influenced by ecology, such that bacteria with a free-living stage have large genomes and high rates of horizontal gene transfer, while obligate intracellular bacteria have small genomes with typically low amounts of gene exchange. However, recent studies indicate that obligate intracellular species that host-switch frequently harbor agents of horizontal transfer such as mobile elements. For example, the temperate double-stranded DNA bacteriophage WO in Wolbachia persistently transfers between bacterial coinfections in the same host. Here we show that despite the phage's rampant mobility between coinfections, the prophage's genome displays features of constraint related to its intracellular niche. First, there is always at least one intact prophage WO and usually several degenerate, independently-acquired WO prophages in each Wolbachia genome. Second, while the prophage genomes are modular in composition with genes of similar function grouping together, the modules are generally not interchangeable with other unrelated phages and thus do not evolve by the Modular Theory. Third, there is an unusual core genome that strictly consists of head and baseplate genes; other gene modules are frequently deleted. Fourth, the prophage recombinases are diverse and there is no conserved integration sequence. Finally, the molecular evolutionary forces acting on prophage WO are point mutation, intragenic recombination, deletion, and purifying selection. Taken together, these analyses indicate that while lateral transfer of phage WO is pervasive between Wolbachia with occasional new gene uptake, constraints of the intracellular niche obstruct extensive mixture between WO and the global phage population. Although the Modular Theory has long been considered the paradigm of temperate bacteriophage evolution in free-living bacteria, it appears irrelevant in phages of obligate intracellular bacteria

Population Genomics and Phylogeography of an Australian Dairy Factory Derived Lytic Bacteriophage

In this study, we present the full genomic sequences and evolutionary analyses of a serially sampled population of 28 Lactococcus lactis–infecting phage belonging to the 936-like group in Australia. Genome sizes were consistent with previously available genomes ranging in length from 30.9 to 32.1 Kbp and consisted of 55–65 open reading frames. We analyzed their genetic diversity and found that regions of high diversity are correlated with high recombination rate regions (P value = 0.01). Phylogenetic inference showed two major clades that correlate well with known host range. Using the extended Bayesian Skyline model, we found that population size has remained mostly constant through time. Moreover, the dispersion pattern of these genomes is in agreement with human-driven dispersion as suggested by phylogeographic analysis. In addition, selection analysis found evidence of positive selection on codon positions of the Receptor Binding Protein (RBP). Likewise, positively selected sites in the RBP were located within the neck and head region in the crystal structure, both known determinants of host range. Our study demonstrates the utility of phylogenetic methods applied to whole genome data collected from populations of phage for providing insights into applied microbiology

Virulence related sequences: insights provided by comparative genomics of Streptococcus uberis of differing virulence

Background: Streptococcus uberis, a Gram-positive, catalase-negative member of the family Streptococcaceae is an important environmental pathogen responsible for a significant proportion of subclinical and clinical bovine intramammary infections. Currently, the genome of only a single reference strain (0140J) has been described. Here we present a comparative analysis of complete draft genome sequences of an additional twelve S. uberis strains. Results: Pan and core genome analysis revealed the core genome common to all strains to be 1,550 genes in 1,509 orthologous clusters, complemented by 115-246 accessory genes present in one or more S. uberis strains but absent in the reference strain 0140J. Most of the previously predicted virulent genes were present in the core genome of all 13 strains but gene gain/loss was observed between the isolates in CDS associated with clustered regularly interspaced short palindromic repeats (CRISPRs), prophage and bacteriocin production. Experimental challenge experiments confirmed strain EF20 as non-virulent; only able to infect in a transient manner that did not result in clinical mastitis. Comparison of the genome sequence of EF20 with the validated virulent strain 0140J identified genes associated with virulence, however these did not relate clearly with clinical/non-clinical status of infection. Conclusion: The gain/loss of mobile genetic elements such as CRISPRs and prophage are a potential driving force for evolutionary change. This first “whole-genome” comparison of strains isolated from clinical vs non-clinical intramammary infections including the type virulent vs non-virulent strains did not identify simple gene gain/loss rules that readily explain, or be confidently associated with, differences in virulence. This suggests that a more complex dynamic determines infection potential and clinical outcome not simply gene content

Association of a Bacteriophage with Meningococcal Disease in Young Adults

Despite being the agent of life-threatening meningitis, Neisseria meningitidis is usually carried asymptomatically in the nasopharynx of humans and only occasionally causes disease. The genetic bases for virulence have not been entirely elucidated and the search for new virulence factors in this species is hampered by the lack of an animal model representative of the human disease. As an alternative strategy we employ a molecular epidemiological approach to establish a statistical association of a candidate virulence gene with disease in the human population. We examine the distribution of a previously-identified genetic element, a temperate bacteriophage, in 1288 meningococci isolated from cases of disease and asymptomatic carriage. The phage was over-represented in disease isolates from young adults indicating that it may contribute to invasive disease in this age group. Further statistical analysis indicated that between 20% and 45% of the pathogenic potential of the five most common disease-causing meningococcal groups was linked to the presence of the phage. In the absence of an animal model of human disease, this molecular epidemiological approach permitted the estimation of the influence of the candidate virulence factor. Such an approach is particularly valuable in the investigation of exclusively human diseases

Diverse Temperate Bacteriophage Carriage in Clostridium difficile 027 Strains

The hypervirulent Clostridium difficile ribotype 027 can be classified into subtypes, but it unknown if these differ in terms of severity of C. difficile infection (CDI). Genomic studies of C. difficile 027 strains have established that they are rich in mobile genetic elements including prophages. This study combined physiological studies, electron microscopy analysis and molecular biology to determine the potential role of temperate bacteriophages in disease and diversity of C. difficile 027.We induced prophages from 91 clinical C. difficile 027 isolates and used transmission electron microscopy and pulsed-field gel electrophoresis to characterise the bacteriophages present. We established a correlation between phage morphology and subtype. Morphologically distinct tailed bacteriophages belonging to Myoviridae and Siphoviridae were identified in 63 and three isolates, respectively. Dual phage carriage was observed in four isolates. In addition, there were inducible phage tail-like particles (PT-LPs) in all isolates. The capacity of two antibiotics mitomycin C and norfloxacin to induce prophages was compared and it was shown that they induced specific prophages from C. difficile isolates. A PCR assay targeting the capsid gene of the myoviruses was designed to examine molecular diversity of C. difficile myoviruses. Phylogenetic analysis of the capsid gene sequences from eight ribotypes showed that all sequences found in the ribotype 027 isolates were identical and distinct from other C. difficile ribotypes and other bacteria species.A diverse set of temperate bacteriophages are associated with C. difficile 027. The observed correlation between phage carriage and the subtypes suggests that temperate bacteriophages contribute to the diversity of C. difficile 027 and may play a role in severity of disease associated with this ribotype. The capsid gene can be used as a tool to identify C. difficile myoviruses present within bacterial genomes

Diversity of Prophage DNA Regions of Streptococcus agalactiae Clonal Lineages from Adults and Neonates with Invasive Infectious Disease

The phylogenetic position and prophage DNA content of the genomes of 142 S. agalactiae (group-B streptococcus, GBS) isolates responsible for bacteremia and meningitis in adults and neonates were studied and compared. The distribution of the invasive isolates between the various serotypes, sequence types (STs) and clonal complexes (CCs) differed significantly between adult and neonatal isolates. Use of the neighbor-net algorithm with the PHI test revealed evidence for recombination in the population studied (PHI, P = 2.01×10−6), and the recombination-mutation ratio (R/M) was 6∶7. Nevertheless, the estimated R/M ratio differed between CCs. Analysis of the prophage DNA regions of the genomes of the isolates assigned 90% of the isolates to five major prophage DNA groups: A to E. The mean number of prophage DNA fragments amplified per isolate varied from 2.6 for the isolates of prophage DNA group E to 4.0 for the isolates of prophage DNA group C. The isolates from adults and neonates with invasive diseases were distributed differently between the various prophage DNA groups (P<0.00001). Group C prophage DNA fragments were found in 52% of adult invasive isolates, whereas 74% of neonatal invasive isolates had prophage DNA fragments of groups A and B. Differences in prophage DNA content were also found between serotypes, STs and CCs (P<0.00001). All the ST-1 and CC1 isolates, mostly of serotype V, belonged to the prophage DNA group C, whereas 84% of the ST-17 and CC17 isolates, all of serotype III, belonged to prophage DNA groups A and B. These data indicate that the transduction mechanisms, i.e., gene transfer from one bacterium to another by a bacteriophage, underlying genetic recombination in S. agalactiae species, are specific to each intraspecies lineage and population of strains responsible for invasive diseases in adults and neonates

Influence of the Escherichia coli oxyR gene function on λ prophage maintenance

In Escherichia coli hosts, hydrogen peroxide is one of the factors that may cause induction of λ prophage. Here, we demonstrate that H2O2-mediated λ prophage induction is significantly enhanced in the oxyR mutant host. The mRNA levels for cI gene expression were increased in a λ lysogen in the presence of H2O2. On the other hand, stimulation of the pM promoter by cI857 overproduced from a multicopy plasmid was decreased in the ΔoxyR mutant in the presence of H2O2 but not under normal growth conditions. The purified OxyR protein did bind specifically to the pM promoter region. This binding impaired efficiency of interaction of the cI protein with the OR3 site, while stimulating such a binding to OR2 and OR1 sites, in the regulatory region of the pM promoter. We propose that changes in cI gene expression, perhaps in combination with moderately induced SOS response, may be responsible for enhanced λ prophage induction by hydrogen peroxide in the oxyR mutant. Therefore, OxyR seems to be a factor stimulating λ prophage maintenance under conditions of oxidative stress. This proposal is discussed in the light of efficiency of induction of lambdoid prophages bearing genes coding for Shiga toxins

Role of Intraspecies Recombination in the Spread of Pathogenicity Islands within the Escherichia coli Species

Horizontal gene transfer is a key step in the evolution of bacterial pathogens. Besides phages and plasmids, pathogenicity islands (PAIs) are subjected to horizontal transfer. The transfer mechanisms of PAIs within a certain bacterial species or between different species are still not well understood. This study is focused on the High-Pathogenicity Island (HPI), which is a PAI widely spread among extraintestinal pathogenic Escherichia coli and serves as a model for horizontal transfer of PAIs in general. We applied a phylogenetic approach using multilocus sequence typing on HPI-positive and -negative natural E. coli isolates representative of the species diversity to infer the mechanism of horizontal HPI transfer within the E. coli species. In each strain, the partial nucleotide sequences of 6 HPI–encoded genes and 6 housekeeping genes of the genomic backbone, as well as DNA fragments immediately upstream and downstream of the HPI were compared. This revealed that the HPI is not solely vertically transmitted, but that recombination of large DNA fragments beyond the HPI plays a major role in the spread of the HPI within E. coli species. In support of the results of the phylogenetic analyses, we experimentally demonstrated that HPI can be transferred between different E. coli strains by F-plasmid mediated mobilization. Sequencing of the chromosomal DNA regions immediately upstream and downstream of the HPI in the recipient strain indicated that the HPI was transferred and integrated together with HPI–flanking DNA regions of the donor strain. The results of this study demonstrate for the first time that conjugative transfer and homologous DNA recombination play a major role in horizontal transfer of a pathogenicity island within the species E. coli

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