Bhaskar Chandra Mohan Ramisetty
Pavithra Anantharaman Sudhakari- Laboratory of Molecular Biology and Evolution, School of Chemical and Biotechnology, SASTRA Deemed University, Thanjavur, India
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.
Introduction
Genome plasticity is an essential requirement for evolution directed by ‘econiche.’ Bacterial genomes are highly plastic relative to other organisms owing to the diversity and dynamicity of the niche (Dobrindt et al., 2010). The success of bacterial ‘omnipresence’ is predominantly dependent on the propensity and probability of genomic plasticity. The frequent exposure to and the ability to integrate exogenous genetic information enhances their genome plasticity (Hacker and Carniel, 2001; Dobrindt et al., 2010; Darmon and Leach, 2014; D’Souza et al., 2015). The basic principles of bacterial genomic plasticity comprise of (i) acquisition of exogenous genetic information and (ii) deletion of unnecessary genetic information at the population level by the selection imposed by the niche. The niche, including the inter and intraspecies competition, is the prime factor that imposes the selection of beneficial genetic information and optimization of the genomes based on the costs and benefits of each variation (Figure 1). The intraspecies competition achieves deletion of non-beneficial genetic information wherein the individuals with a higher metabolic burden, and slow growth rate are eliminated under natural testing conditions thereby selecting individuals with lesser ‘junk’ DNA (Dobrindt and Hacker, 2001). Furthermore, bacteria with deleterious recombination events or mutations are also eliminated thereby gradually ‘ridding off’ of frivolous genetic information from the population through a process generally referred to as ‘purifying selection.’ In other words, the econiche selects the best ‘composites of genetic information’ of the available ones based on the reproductive potential and stress tolerance imposed by the habitat (Forde et al., 2004).
Figure 1. Bacterial genome evolution. Based on the habitat, bacterial genomes acquire DNA (genetic information) from the environment. Genetic information from the phages, plasmids, and dead organisms maybe integrated into the genome by homologous or non-homologous recombination.
Genome plasticity is predominantly mediated by horizontal gene transfer (HGT) and intragenomic recombinations such as transposition. The three principal HGT mechanisms are (i) transformation (direct uptake of DNA), (ii) transduction (DNA transfer mediated by phages), and (iii) conjugation (DNA transfer through physical contact between two bacteria) (Davison, 1999; von Wintersdorff et al., 2016). Of the three HGT mechanisms, transduction is highly ‘active’ by virtue of viral tropism, and mechanisms of DNA injection into the recipient bacterium. Occasionally, some of the temperate phages reversibly integrate into the host genome by a phenomenon referred to as lysogeny. The integrated phage genome, now referred to as ‘prophage,’ is vertically inherited by the daughter cells; the prophage replicates as a part of the host genome. Eventually, the prophage would switch from lysogenic to lytic cycle often resulting in phage multiplication and host cell death (Gandon, 2016). Bacteriophages, in general, play an important role in the ecology as well as the evolution of bacteria through several types of interactions between the host and phage genomes (Roossinck, 2011; Nasir et al., 2017). In fact, phages and their hosts have coevolved strategies that impact the survival, persistence, and evolution of their respective genomes (Buckling and Brockhurst, 2012; Koskella and Brockhurst, 2014; Nasir et al., 2017). Prophages were shown to influence the microbial community through enhanced recombinations (Nadeem and Wahl, 2017; Braga et al., 2018).
Mutations in one or more genetic elements required for the excision result in the failure of prophage excision from the host genome rendering them ‘grounded,’ also referred to as ‘cryptic’ or defective prophages (Casjens, 2003; Wang et al., 2010). The whole genome sequencing of diverse bacterial species has revealed multiple prophages within each genome. The prevalence of multiple grounded prophages within most bacterial genomes is indicative of eco-evolutionary selection of such genomes. Although the advantages of prophages are being explored (Bobay et al., 2014; Menouni et al., 2015; Gandon, 2016; Howard-Varona et al., 2017), the ecological scenarios and the possible evolutionary impact of these prophages are still unclear. In this manuscript, we discuss the mechanisms of grounding, and propose eco-evolutionary perspectives rationalizing the prevalence and advantages of grounded prophages to the host bacteria. We performed sequence analyses of five grounded prophages (DLP12, e14, Rac, CPZ-55, and Qin) of Escherichia coli to illustrate the distribution and the advantages.
Furthermore, several types of mutations can occur in the genome such as point mutations, duplications, inversions, deletions, etc. Within a population of bacteria, these genetic variations occur at a low frequency. Often, these variations are deleterious: i.e., the benefits caused by the variations are less than the costs incurred. Several parameters, such as nutrition, stress tolerance, competition, and resistance to antibiotics impose stringent selective pressure on the population. Hence, bacteria with deleterious variations usually perish or are outcompeted by the faster growing kin. Bacteria with beneficial genetic variations, such as higher growth rates, novel metabolic pathways, higher tolerance to stress, and higher competitive fitness fare better and are selected.
Each habitat may pose various stresses of varying degrees, allowing the gradual selection of organisms optimal for successful propagation in the prevailing conditions.
Bacterial Genome Plasticity
Genome plasticity drives the successful spatiotemporal propagation of bacteria in terms of ecology as well as evolution. Genetic rearrangements, caused by transformation, transposition, and recombination events (site-specific as well as homologous), enhanced by insertion sequences (IS) elements, integrons, conjugative transposons, plasmids, bacteriophages, pathogenicity islands, and genomic islands result in a gain (or loss) of genes with diverse functions (Dobrindt and Hacker, 2001; Betts et al., 2018). Different bacteria may have different degrees of plasticity (Gaudriault et al., 2006) and, furthermore, regions such as ‘Regions of Genomic Plasticity’ (RGPs) of a bacterial genome are thought to have a higher propensity for genetic rearrangements relative to other core regions. The four different RGPs categories (genomic islands, prophages, regions encoding recombinases, and unclassified) are mosaics consisting of several modules. Each module ranging from 0.5 to 60 kb consists of genes involved in metabolism, intra and intercellular DNA mobility, drug resistance, host/environmental interactions, and antibiotic synthesis (Ogier et al., 2010). Different bacterial species have varying potential to take up exogenous DNA, metabolize it as a nutrient, or recombine the incoming DNA into the genome. The ability to take up exogenous DNA is termed as competence, a complex mechanism of temporal regulation for the uptake as well as the integration of the exogenous DNA (Blackstone and Green, 1999; Palchevskiy and Finkel, 2006; Maughan et al., 2010; Sinha and Redfield, 2012; Mell and Redfield, 2014). Often, the DNA could be catabolized and reused in nucleic acid metabolism or for generation of energy depending on the metabolic state of the cell (Finkel and Kolter, 2001). The DNA that is taken up maybe recombined into the bacterial host genome by homologous recombination or non-homologous mechanisms resulting in transformation (Fischer et al., 2001; Blokesch, 2017; Veening and Blokesch, 2017). Integration of exogenous DNA is the largest contributor of the pan-genome, the cumulative genetic information of all strains of a single bacterial species. The total amount of incoming DNA during the existence of a species in a habitat, theoretically, would imply an enormous genome. However, bacterial populations tend to get rid of larger genomes by intraspecies competition in limiting conditions imposed by the niche. The reduction of genome size is evident from the occurrence of multiple pseudogenes and loss of several genes within the genomes of most bacteria (Lawrence et al., 2001).
A typical bacterial species has expanded the horizons of its habitat through genome plasticity. When non-pathogenic bacteria infect a host, the niche may select the individual bacterium that may have specific mutations/genetic information allowing survival/virulence and eventually evolve into a pathogenic strain (Hogardt et al., 2007). Such mutations confer fitness by the acquisition of virulence genes and loss of antivirulence genes that nullify the action of factors conferring virulence (Dobrindt et al., 2010). For example, the presence of antivirulence genes encoding enzymes involved in arabinose catabolism has reduced the virulence in non-pathogenic Burkholderia thailandensis. This has been experimentally shown wherein B. pseudomallei that harbors arabinose assimilation operon, grown with arabinose as the only carbon source showed Salmonella-like gene cluster down-regulation (TTSS3-Type III secretion system) involved in virulence (Moore et al., 2004). Shigella and enteroinvasive E. coli (EIEC) species have acquired their virulence through structural mutation of the cadA gene. Another mechanism involved in pathoadaptation is the mutation of cis-regulatory elements (CRE). For example, mutational study in Salmonella species, on the role of cis-regulatory elements, revealed that evolutionary mutations in cis-regulatory elements such as SsrB regulatory system enhanced the expression of srfN gene involved in intrahost fitness (Osborne et al., 2009). The enormous genetic variability in Helicobacter pylori, a human pathogen, shows the extent of adaptation brought about by HGT and homologous recombination. However, these sequence variations have occurred without a drastic increase in the genome size of the bacterium (Kraft et al., 2006). Among closely related Salmonella serotypes, extensive genetic diversity was observed in the antigenic determinant genes that encode factors responsible for lipopolysaccharides (LPS), flagella, fimbriae and virulence factors (Fierer and Guiney, 2001). Genome plasticity allows oscillation of acquisition and loss of genetic information from the population as per the change in growth conditions and the habitats.
Prevalence of Prophages in Bacterial Genomes
It is conceivable that the prophages are the largest exogenous contributors of bacterial genetic information, primarily, because of their sheer size. Astonishingly, prophages may constitute 10–20% of a bacterial genome predominantly due to an exposure to a multitude of bacteriophages (Casjens, 2003). Of the completely sequenced E. coli strains, E. coli O157:H7 str. Sakai has the highest number of prophages per genome; it harbors 18 prophage elements constituting 16% of the total genome (Asadulghani et al., 2009). Often, these prophages could be activated by several stresses to undergo lytic cycle. For example, 40 out of 148 Lactobacilli strains, representing 15 species, released phage particles upon mitomycin induction. In particular, 77% of 30 L. Salivarius produced temperate phages, killer particles or defective phages. Temperate prophages were found in 10% of 105 strains of L. bulgaricus and L. lactis (Sechaud et al., 1988). Using VirSorter tool (Roux et al., 2015a), around 5,492 microbial genomes out of the 14,973 publicly available bacterial and archaeal genomes (as on January 2015) had viral-like elements of which 2,445 contained more than one detectable phage elements. 82% of such phage co-infections involved multiple Caudovirales (Roux et al., 2015b). Bacterial genomes may contain intact prophages, LT2 isolate of Salmonella enterica (four prophages), pathogenic Streptococcus pyogenes isolate SF370 (prophage SF370.1) (Canchaya et al., 2002), H. influenza Rd (Flu-Mu) and E. coli sakai (Sp18) or defective prophages like two non-inducible prophages (Lp1 and Lp2) and two prophage remnants (R-Lp3 and R-Lp4) of Lactobacillus plantarum (Ventura et al., 2003). A large number of IS elements that are collectively referred to as minichromosomes comprise highly divergent prophage remnants constituting 22% of Halobacterium sp. NRC-l genome, and code for transposases associated with phage genomes and other mobile genetic elements (Bonneau et al., 2004). A recent web based tool called PhageWeb, for the identification and characterization of prophages within bacterial genomes, may expand the exploration of prophages in a variety of species (Sousa et al., 2018). The prevalence of prophages in bacterial genomes is intriguing, firstly, because of the imminent danger of lytic cycle induction and, secondly, the metabolic burden of replicating and expressing the prophage genetic elements. The eco-evolutionary rationales that operate to sustain prophage rich genomes of bacterial populations in the environment are poorly understood. It is indeed interesting to understand the beneficiaries of lysogeny: the host bacterial genome or the phage genome or both.
Why Are Prophages Prevalent in Bacterial Genomes?
(a) Lysogeny Mechanisms
Lysogeny is a common feature in almost all habitats including the gut microbiome (Howard-Varona et al., 2017; Kim and Bae, 2018). Integration into the host genome is the fundamental step for the establishment of lysogeny. Different phages employ different mechanisms (Figure 2) that broadly fall into (i) site-specific recombination (Landy and Ross, 1977), and (ii) transposition based integration (Allet, 1979; Harshey, 2014). Site-specific recombination is a well-characterized mechanism of phage (specifically lambda bacteriophages) genome integration and is the only one discussed here for brevity. Site-specific recombination requires phage attP and host attB sites that share 15 nucleotides in common, and thus the crossover occurs either within or flanking regions of this core region (Landy and Ross, 1977). Despite defective attB site Int-mediated recombination occurs between attP and secondary attachment sites in the host genome, albeit at a lower frequency (Shimada et al., 1975). The Cre- and Int-dependent recombination between the lox sites are mediated by site-specific recombination system of P1 and λ phages, respectively. The loxP has two 17 bp binding site for Cre, and each of the sites has 13 bp inverted repeat and a 4 bp region which together forms 8 bp spacer region flanked by two 13 bp inverted repeats. The loxP and loxB sites share limited homology, and hence their recombination efficiency is altered with sequence variation. Cre recombinase makes 6 bp long double strand staggered cut in the spacer region, localized between the adjacent loxP sites followed by strand exchange. The process of Cre-mediated cleavage involves phosphodiester bond breakage at 3’ end and covalent interaction of the primed phosphate with the protein (Hoess and Abremski, 1985). A similar rejoining mechanism is mediated by λ Int recombinase (Craig and Nash, 1983). Four recognition sites form a synapse prior to the cleavage, and this process is homology dependent (Ross et al., 1979). In essence, the essentials required for integration or excision are the att sites (on the bacterial as well as the phage genomes) and the specific recombinase that recognizes the concomitant recombinase recognition sites present within the corresponding att sites. Different recombinases, or here integrases, recognize different sequence specificities wherein the recognized sequences may have considerable sequence differences. Cross-activity of a recombinase of one phage with att site of another is proportional to the degree of homology between their integrases and similarity between the core- and arm-type sites in attachment sites of the respective phages (Groth and Calos, 2004).
Figure 2. Mechanisms of bacteriophage genome integration into bacterial chromosome. (A) Site-specific recombination mediated-lysogenization is well characterized for lambdoid phages. Site-specific recombination is catalyzed by the integrase through phage and bacterial attachment sites (attP and attB, respectively) resulting in recombinant attL and attR sites (Landy and Ross, 1977; Menouni et al., 2015). Prophage induction is mediated by excisase and integrase complex. (B) Transposition based integration of Mu phage into the target site in the host genome requires MuA transposase. The integration results in target site duplication (TSD) (Allet, 1979; Harshey, 2014).
(b) Hotspots for Prophage Integration
The fundamental premise is that any genome with deleterious recombination events, concerning the niche conditions, will be lost in the population whereas tolerable recombination events are unaffected. Occasionally, advantageous recombination events would be selected and that would eventually outcompete their kin. Integration of most prophages is by site-specific recombination at att sites or by transposition at a random site. Although there is a significant chance for the occurrence of the att sites at multiple loci of a bacterial genome, many integration events could be deleterious for the host genome because of gene disruptions or dysregulations leading to greater fitness costs. Various analyses of lysogens have shown that prophages are predominantly found at loci closer to non-coding genes (e.g., tRNA genes), a few functional genes, and intergenic regions. Hence, these integration events should be assumed as either tolerable or advantageous in the growth conditions of the context. For example, in Ralstonia solanacearum the attB site for ϕRSA1 prophage integration occurred downstream of the tRNA-Arg gene (Fujiwara et al., 2008) and in E. coli, the relatives of λ phage and P4 phage insert within the tRNA genes (Campbell, 2003). It is rational for a phage with a broader range of host specificity to prefer integration sites that are within protein-coding functional genes as their occurrence is conserved across distantly related species, subject to the tolerable consequences of the recombination (Hill et al., 1989). Based on the composition and length of the integration site, a variable number of sites within a host genome for the phage integration are possible. However, Mu-like prophages show random integration events due to transposition-based integration. A robust analysis of 471 E. coli and Salmonella prophages showed 58 distinctive integration loci for 369 E. coli prophages and 102 loci for S. enterica prophages (Bobay et al., 2013). Based on the empirical evidence of various prophages in numerous bacterial genomes, one may assume that integration of phage genome into that of the bacterial genome at the specified loci are either tolerable or advantageous in the ecological context in which these lysogens would be selected.
(c) Advantages of Lysogeny
The size of the prophages (approximately 45 kb each) and the expression of prophage genes impose a significant metabolic burden on the lysogen growing in stringent natural conditions competing with non-lysogenic peers. It is conceivable that the bacterial cells, which do not harbor any prophages, have a competitive growth advantage relative to lysogens. However, as noted earlier, the prevalence of multiple cryptic prophages indicate that the lysogenic strains had a survival advantage compared to the non-lysogens. Intriguingly, lysogeny, as a trait of the bacterial genome as well as phage genome, must have been selected by virtue of some spatiotemporal advantage to both the host bacterium and the phage (Paterson et al., 2010; Roossinck, 2011). One of the straight forward explanations is that some of the prophages contribute genes of selective advantages, such as antibiotic resistance conferring genes encoding multidrug resistance pump, outer membrane protease, cell division inhibiting factors, small toxic membrane polypeptide, (Wang and Wood, 2016) etc., in the habitat and hence the lysogens would outcompete the non-lysogens. However, harboring large prophages that do not confer a selective advantage is a metabolic burden and hence is a fitness disadvantage to the lysogen. The key to deciphering this paradox—selection despite the metabolic burden—is the vulnerability of non-lysogens to phage-induced lysis while the lysogens are immune to the corresponding phage. Prophages bestow lysogens with immunity against secondary phage infections, referred to as superinfection exclusion (or immunity) (Figure 3). Assuming that the ecological niche of the bacterial population contains the phages, the lysogens have a survival advantage due to the ‘superinfection immunity’ conferred by the corresponding phage. Conferring immunity is ecologically vital for the resident prophage as well as the host bacterium. Lysogeny in conjunction with superinfection immunity confers the selective advantage to the lysogen from phage attack and concomitantly also allows the propagation of the lysogen. Superinfection immunity, as a mechanism to prevent multiple infections of a lysogen by similar phages, is achieved by different mechanisms such as: (i) conformational change in phage receptor sites, (ii) modifications in DNA transfer system, (iii) inhibition of phage lysozyme for penetration, (iv) incoming phage DNA breakdown, and (v) abortive exclusion. The gp15 expression of HK97 had shown resistance to HK97 and its close relative HK75 infections (Cumby et al., 2012). Similar factors include (i) SieA and SieB of Salmonella phage P22 blocks the DNA entry with unknown component of DNA transfer system and leads to abortive infection, respectively (Susskind et al., 1974), and (ii) Sie2009 of Lactococcus lactis phage Tuc2009 (Mahony et al., 2008). Rex system of λ prophage is thought to induce superinfection exclusion via abortive exclusion. RexA (a component of cell cytoplasm) is activated by protein–DNA interaction which then mobilizes RexB, a membrane ion channel protein, thereby altering the membrane potential of the host (Parma et al., 1992). The CI repressor (a positive regulator of its own expression) of the incumbent λ phage inhibits the establishment of superinfecting λ phage upon forming CI repressor-DNA complex and thus preventing the integration and replication of the secondary λ phage (Gottesman and Weisberg, 2004; Fogg et al., 2010). The defective e14 prophage-encoded Lit protease upon its activation by GoI peptide (located within the major head protein) cleaves the elongation factor EF-Tu, thus blocking the late gene expression of T4 phage (McGrath et al., 2002). Upon phages attack on a bacterial population, the only survivors are the lysogens and the spontaneously evolved mutants that are no longer recognized by the phage. Such scenarios would allow the elimination of the non-lysogens resulting in survival and growth of the lysogens due to extensive relaxation of peer competition for resources. Hence, the lysogens and the mutants will have unhindered access to the resources and the ‘spoils’ of the dead kin allowing the repopulation of the habitat, a phenomenon similar to that of bottlenecking.
Figure 3. A hypothetical model to explain the prevalence of cryptic prophages in bacterial genomes. A typical temperate phage injects its genome in the bacterial cells that may undergo either lytic cycle or lysogeny. In the lytic cycle, more copies of phages are produced, the host bacterium is lysed, and a multitude of new phages are released into the surrounding. The lytic cycle repeats with each bacterium, which almost annihilates the bacterial population. In some bacteria, the injected genome undergoes lysogeny by integrating into the bacterial genome. Lysogeny allows vertical propagation of prophage along with that of the bacterial genome. By means of the prophage, the lysogens attain immunity to the infecting secondary phages thereby preventing infection/lysis. Apart from other mutants, the lysogens are the only surviving kin and hence will multiply well due to lack of intraspecies competition. Lysogeny, upon reversal, will ultimately lead to a lytic cycle leading to the death of the host bacterium. In some of the prophages, mutations in the recombinase recognition sites or recombinase genes occur, which prevent the elicitation of the lytic cycle: a phenomenon referred to as ‘grounding of prophages.’ The bacteria with grounded prophages have two advantages; immunity from attack by related secondary phages and irreversible lysogeny that prevents phage-mediated cell death. Eventually, other lysogens undergo lytic cycle upon induction by various stimuli resulting in lysis of the host bacterium. The bacterium with cryptic prophages will survive and have a growth advantage due to reduced competitive pressure from kin.
It is also proposed that the prophage genes could confer a competitive advantage at a population level such as in the case of biofilm formation (Wang et al., 2010). The individual cost of prophage induction is explained as a mechanism of altruism since the induced lysis of the lysogens that benefits their population as a whole. For example, phage-induced lysis of the lysogens in the biofilm further enhances the biofilm formation by providing additional nutrients, bacterial dispersal, and extracellular DNA. Similarly, the bacteria competing for resources produce bacteriocins and require transport protein for their release. In cases where the transport protein is not encoded, like the group B colicins producing operon in E. coli and Shigella, the phage-induced lysis promotes the release of colicin bacteriocin (Obeng et al., 2016). Shiga toxin-producing E. coli are lysogens for lambdoid prophages that also encode Shiga toxin. The prophage induction is followed by the expression of Shiga toxin genes release of the Shiga toxin (Iversen et al., 2015). Shiga toxin targets the eukaryotic protein synthesis of cells that can internalize the toxin. In a small subset of the lysogens the prophage is induced to lytic cycle by stimuli like H2O2 (Shimizu et al., 2009; Los et al., 2010, 2013). The lysis of the bacteria causes the release of Shiga toxin which acts on their predators (protozoan or neutrophils as the case maybe) (Dydecka et al., 2017; Chakraborty et al., 2018). The attenuation or inactivation of the predators is thought to increase the survivability of the rest of the population. It is hypothesized that these Shiga converting prophages is an example of altruism at the population level. However, the concept of altruism is still debatable in the case of bacteria (Ramisetty et al., 2015; Durand et al., 2016; Peeters and de Jonge, 2018). Specifically with free-living bacteria, altruistic behavior would benefit the self and/or kin if only the beneficiaries of the donation are the kin. The ultimate altruism in the form of suicide is detrimental to the unicellular free-living species unless the fraction of donors is very small. Bacteria in nature occur as multispecies colonies/biofilms wherein intraspecies-altruism is highly unlikely because the sacrifice by the donors may benefit the competing species equally. Feeding the competitors is counterproductive to the goals of altruism. In the case of prophage-mediated host lysis, the selective induction of prophage, the mechanisms limiting the induction to a small fraction of the population, and the beneficiaries of the altruistic behavior require more experimental evaluation (Los et al., 2013).
Multiple other ecological advantages of lysogeny were also reported. Phage conversion upon transduction alters the host physiology with respect to metabolism, pathogenicity, and niche adaptation (Nasir et al., 2017). The integrated auxiliary metabolic genes (AMGs) among the marine viromes have been reported to carry essential genes for different nutrient cycles such as dsr and sox genesin sulfur cycle, amok and P-II in nitrogen cycle, photosynthesis (psbA and psbD genes from Prochlorococcus and Synechococcus cyanophages) (Roux et al., 2016), motility proteins (e.g., flaB), and flagellar motor complexes (e.g., motA) (Hurwitz et al., 2015). Prophages increase the survival (Steinberg and Levin, 2007), cell growth (Picozzi et al., 2015), virulence (Waldor and Mekalanos, 1996; Zschöck et al., 2000), promote biofilm formation, phase variation (Kutsukake and Iino, 1980; Silverman and Simon, 1980; Sekulovic and Fortier, 2015), quorum sensing (Hargreaves et al., 2014), confer resistance to antibiotics/environmental stresses (Hyder and Streitfeld, 1978; Lopez et al., 2012; Wipf et al., 2014), and immunotolerance (Farrant et al., 1997). Yet another example for prophage-enhanced environmental stress resistance of the lysogen is phage-regulated acid resistance in E. coli. Comparative analysis of non-lysogen and lysogen of Shiga-toxigenic phage (Φ24B), revealed the upregulation of operons encoding type I fimbriae and the glutamic acid decarboxylase (GAD) acid stress island (Veses-Garcia et al., 2015). It is highly likely that different prophages may have conferred one or more of the above discussed eco-evolutionary advantages. More such rationales could also be expected with the exploration and molecular examination of prophages in many other bacterial genomes. Of all the advantages discussed, immunity to superinfection is the most likely and effective advantage conferred by the prophage.
(d) ‘Grounding’ of Prophages
Similar to transposons, the prophages are also prone to ‘grounding,’ which is to imply irreversible lysogeny due to genetic events that prevent their excision. Multiple trans and cis elements are involved in the site-specific recombination of the phage with the host genome. The site-specific integration of phage (specifically lambdoid bacteriophage) into their specific host genome via host and phage attachment sites (“attB” and “attP” sites, respectively) yields the recombinant “attL” and “attR” sites. attL and attR sites are the essential cis elements for excision of the integrated prophage upon switching from lysogenic to lytic cycle. The prerequisites for the integration and excision of phage genomes are (i) phage-encoded recombinase proteins (integrase and excisase); (ii) conserved prophage flanking “attR” and “attL” excisase recognition sites; and (iii) the expression of accessory host factors such as IHF and factor for inversion stimulation (Fis) proteins. Efficient excision is achieved by the cooperative role of both integrase and excisase. It has been shown that Xis protein promotes excision by recognizing the phage ‘attR’ site in addition to sequence-specific binding of Int or IHF protein (Yin et al., 1985). The assembly of the excisive complex and the rate of prophage curing depends on the amount of Xis expressed in the host. The host Fis enhances the binding specificity of Xis thereby promoting the formation of excisive intasome on attR site (Papagiannis et al., 2007). The phage cII protein, along with IHF, upregulates the expression of Int and late gene repressors and thus promotes lysogeny. The mRNA level of the integrase gene in S. aureus prophage is influenced by the host alternative sigma factor σH (Tao et al., 2010). The Bxb1 prophage of Mycobacterium species encodes an accessory protein, gp47, which is involved in the directionality of Bxb1 integrase for the assembly of excisive intasome (Ghosh et al., 2006). Similarly, gp3 of streptomyces phage ΦBT (Zhang et al., 2013), gp52 of streptococcus bacteriophage ΦJoe, Rv1584c from ΦRv1 phage of Mycobacterium tuberculosis (Bibb et al., 2005), and SprA from SPβ prophage of Bacillus subtilis (Abe et al., 2014) are recombination directionality factors (RDFs). The insertion of the IS element and further rearrangement in lambdoid prophage of Shigella dysenteriae rendered them defective upon deletion of stx region (McDonough and Butterton, 1999). Similarly, the RNA-mediated translational control of int expression by sib results in the formation of a truncated integrase which in turn hampers the excision of defective prophages (Schindler and Echols, 1981).
The intracellular concentration of proteins involved in recombination such as Int, IHF, Xis, and Fis is highly regulated. IHF inhibition of excisive recombination occurs at lower levels of Int protein (Thompson et al., 1986). Hence, the mostly likely events for the grounding of temperate prophages are: (i) defective recombinase proteins (Swalla et al., 2003), (ii) homologous recombination during the super infection of related phages (Canchaya et al., 2003), and (iii) resected att sites. DNA invertase, Xis, and Fis could promote genetic rearrangement like DNA inversion, specifically between the attP site and secondary attachment site within the prophage. Such genetic rearrangements could also affect the expression of phage-encoded proteins and modify the orientation of the attachment sites resulting in the establishment of irreversible lysogeny or ‘grounding’ (Dorgai et al., 1993).
(e) Degeneration of Prophages
The grounded prophages are highly prone to deletions, especially the genes whose products do not confer any advantage or are regressive. Harboring large prophages and expressing the prophage genes is a metabolic burden. Hence, minimizing such genetic elements are beneficial to the host genome. The closely related E. coli K-12 and Shigella flexneri showed variability in their prophage remnants (Brussow et al., 2004). A study on the genome diversity of Salmonella enterica serovar Typhimurium found deletion in gipA gene of L927 and L847 phages. Three complete deletions and four partial deletions in the virulence genes have been found in six of the sequenced strains of Salmonella enterica, emphasizing the diverse ecology of phage types. Study on E. coli O157:H7 EDL933 strain uncovered various types of indels of IS elements in 12 of their prophages (Iguchi et al., 2006). The prophages of E. coli O157 EDL933 have lost their deleterious gene such as λ N gene whose expression would result in host lysis. However, the neutral or beneficial genes for lysogenization like int and other phage structural genes occurred in most of their prophages. It was also found that E. coli K-12 harbored an incomplete Rac prophage, which was intact in EDL933 strain (Lawrence et al., 2001). A simulation study showed that large-scale deletions occurred in the accessory genes such as integrase and cargo genes, rather than the conserved genes involved in repression of lytic genes (Canchaya et al., 2002), replication, expression of capsid proteins, and packaging. For example, tail protein encoding regions are highly conserved over cargo genes in lambdoid prophages; morons 1 and 2 were least conserved in P2 like prophages (Bobay et al., 2014). The low ion-irradiation experiment showed that the IS elements and pseudogenes as a preferential site for deletions (Song and Luo, 2012) while a study on Neisseria meningitides strains showed a IS30 transposase replacement for head and tail morphogenesis genes (Dunning Hotopp et al., 2006). Deletion of inutile parts of the genome confers a growth advantage due to lowered metabolic burden or reduction in the inhibition of growth. Such bacteria will eventually outnumber the counterparts leading to smaller but optimal genomes in the population. However, other views suggest that the degeneration of prophages is driven by the need to eliminate the deleterious genes rather than to preserve chromosome compaction (Lawrence et al., 2001).
What Are the Advantages of Grounded Prophages?
As noted earlier, prophages are dynamic zones within the genome wherein the frequency of genetic variation is very high. There are several disadvantages of bearing the grounded prophages, mostly stemming from the metabolic burden of replicating the large prophage and the expression of prophage genes. However, these disadvantages are subject to contextual economics, i.e., the resources available, growth conditions, and the competition would determine the costs to benefits ratio. The pre-selected lysogens continue to survive unless there are drastic changes in the growth conditions, competition or habitat change. There is also a possibility for the reversal of grounding through homologous recombination of a defective incumbent prophage with a related superinfecting-phage genome, leading to prophage activation and resulting in host lysis. Within a habitat, individual bacteria that may have acquired genetic changes within prophage imposing high costs/benefits are eliminated in the purifying selection. However, this purifying selection based on the prophage would have negligible impact on the species. On the other hand, there is a need to explore various rationales to reason the selection of individuals within a given habitat. There is minimal consistency in the preservation of a prophage. Besides, different prophages are implicated in a plethora of functions to the host genome. Hence, in the following sections we propose various advantages of harboring prophages to explain how prophages may have influenced the bacterial ecological success.
(a) Variations due to Prophage Integration
Every event of integration causes a genetic variation; those that are either tolerable or advantageous are retained, similar to the consequences of transposition. If an integration event were deleterious, the bacterium with such mutations would be eliminated from the population. On the contrary, if a mutation is retained within the population, it must be either tolerable or advantageous. The λ integration at the exuR gene, involved in the regulation of exu regulon, rendered them non-functional in E. coli HfrH 58 strain. It is possible that the expression of exu regulon in the ecological context was not required. However, the functionality of exuR gene is regained upon prophage induction by site-specific recombination (Mata et al., 1978). Insertional inactivation of lipase-encoding gene by L54a and L54b prophages in Staphylococcus aureus PS54 has been reported. However, the curing of L54a prophage has restored the lipase-encoding gene expression but not with L54b curing (Lee and Iandolo, 1985). Similarly, S. aureus MW2 harbors ΦSa3mw prophage whose integration happened in haemolytic toxin gene, β-hemolysin (Hlb). Upon prophage induction, the Hlb gene expression was restored (Katayama et al., 2013). Most of the inversions occur as a result of homologous recombination between prophages (as hot spots for large-scale inversions), IS elements, and RNA-encoding operons (Iguchi et al., 2006; Camilli et al., 2011). Multiple inversions among the Mu-like prophages in various strains are shown to influence the host specificity by altering the expression of various host-cell-receptor recognition genes (Schmucker et al., 1986) and loss of tail genes (Braid et al., 2004). It is indeed difficult to ascertain if any of the phage integrations, solely by virtue of the integration event, has conferred a selective advantage.
(b) Hotspots of Horizontal Gene Transfer
The probability of acquiring exogenous DNA into the genome is dependent on several aspects which include the availability of exogenous DNA, the competence of the host bacterium to take up the DNA, the expression of appropriate recombinases to integrate the DNA, and importantly the costs/benefits of the integration event. There is a lot of evidence that the numerous recombination events in a bacterial genome are associated with their prophages. Multiple prophages encode recombinases, which may aid in integration. However, the prevalence of indels within prophages is not indicative that prophages somehow ‘actively’ acquire genes. Rather, to state accurately, indels may occur randomly in any locus of the bacterial genome, but those recombinant genomes that have resulted in deleterious consequences (high costs) would have been eliminated from the population. Hence, the prevalence of multiple recombination events within the prophages of various strains within a given species implies that these recombinations are either tolerated or advantageous. This aspect could be best explained with the genetic phage mosaicism observed in several prophages. Phages of different groups or broad host range show the shuffling of virulence genes between their genomes (Desiere et al., 2001; Mirold et al., 2001). The number of mosaic phages is higher among gram negative bacteria over gram positive bacteria as the former have a broad host range and are recombinogenic (Brussow and Desiere, 2001).
Lambdoid prophages are capable of supplementing virion proteins to other lambdoid phages (Asadulghani et al., 2009). A study on the genome sequence of HK97, HK022, λ, and P22 prophages has found the mosaic pattern in their functional units such as individual genes, gene clusters, or a portion of protein-encoding genes (Juhala et al., 2000). Recombinases play a vital role in the formation of mosaics among prophages. For example, Rad52-like recombinase promoted the homologous recombination between lambdoid phages and E. coli prophage remnants. The prophage remnants in many cases have been found to complement functional genes involved in the lytic cycle. For example, DLP12 lysis protein essential in E. coli for the biofilm formation is also found active in lambda phages (De Paepe et al., 2014).
Horizontal gene transfer has led to the integration of ‘morons,’ homologs of bacterial genes, into the phage genomes. The moron gp15 of HK97 prophage in Enterobacteriaceae shared no homologs among the other phages but it did with bacterial genes that expressed YebO family of proteins. Other examples include GogB effector protein of Salmonella enteric prophage Gifsy-1 (Coombes et al., 2005) and NelA-like type III effector protein of enterohemorrhagic E. coli O157:H7 prophage BP4795 (Gruenheid et al., 2004). Similarly, oac moron of Pseudomonas phage D3 was homologous to Pseudomonas WbpC proteins involved in LPS biosynthesis (Burrows et al., 1996).
Prophages may confer direct selective advantages to their host. For example, the lysogenization of E. coli by CPS-53 and CP4-57 prophages promote host fitness under oxidative, osmotic, and acidic stresses (Wang et al., 2010). Likewise, RnlA toxin encoding gene in CP4-57 prophage (Koga et al., 2011) and CbtA toxin in CP4-44 (Tan et al., 2011) enhance persister phenotype in E. coli by inhibiting the cell growth. The paralogs of nanS esterase gene, required for the growth of pathogenic bacteria upon exploiting the host sialic acids, are found downstream of stx toxin genes in prophage or prophage remnants of E. coli strain O157:H7 (Rangel et al., 2016).
(c) Hotspots of de novo Genesis
Since the generation of the genetic information and the evolution of the first cell, most of the genetic information must have evolved within cellular entities. Although it is possible for the de novo genesis of genes in bacterial cells, there is minimal tolerance for genetic events within the core gene-rich loci due to high costs and low benefits. The probability and tolerance for recombinations within the prophages serve as a platform for the generation of new genes. Hence, there is a relatively high probability for the de novo formation of new open reading frames (Bartels et al., 2009) and possibly selected if the product of the ORF is beneficial for the cell. With time and generations, the de novo formed gene may accumulate beneficial mutations and eventually be optimized for an appropriate function and regulated expression. For example, 94.2% of the genome of Staphylococcus aureus strain Jevons B contains ϕJB prophage that has 70 predicted ORFs (Bartels et al., 2009) out of which only 21 ORFs were associated with a putative function. The remaining 49 putative proteins have no known function, but structural domains and homologs have been identified for 26 and transmembrane region in 5 ORFs. The remaining proteins were recorded as “hypothetical” which refers to a legitimate ORF but encodes a protein with unknown function or association with any known proteins (Varga et al., 2016). A study on lysogeny module of StB27 prophage among Staphylococcus hominis and Staphylococcus capitis species found that 42% of the conserved proteins have no predicted functions, and 13 proteins shared no homologs. StB27 of coagulase-negative Staphylococci (CoNS) contained ORF3 of unknown function (Deghorain et al., 2012). Several ORFs of unknown function and hypothetical proteins have been found in prophages among E. coli O157 strains (Svab et al., 2013) (Saile et al., 2016). Similarly, numerous ORFs of unknown functions were identified in multiple species such as Helicobacter pylori (Luo et al., 2012), Ralstonia solanacearum strains (Askora et al., 2009), and Pseudomonas aeruginosa (Braid et al., 2004; Wang et al., 2004). P2-like coliphages of Enterobacter have ORFs that have hypothetical protein homologs in plasmids, genomes of other bacterial species (Nilsson et al., 2004). Interestingly, the prophage mEp021 of lysogenic E. coli K-12 harbors ORFs whose expression shows haemolytic activity in the host. The ORF-4 has four initiation codons within the single frame encoding for 83aa (ORF-4.1), 82aa (ORF-4.2), 77aa (ORF-4.3), and 72aa (ORF-4.4). Among these, the expression of ORF-4.3 has an inductive role in haemolytic activity by releasing vesicles containing bacterial protein HlyE and alters the morphology of the host bacterium (Martinez-Penafiel et al., 2012). The above examples reiterate that the great extent of the de novo genesis of novel genes is relative to any other sites on a typical bacterial genome. However, those ORFs whose products are functional and beneficial will be selected based on the costs/benefits and may evolve into beneficial genes with a high degree of propagative potential.
Grounded Prophages of E. coli: Illustrative Analyses
E. coli K12 MG1655 strain has nine cryptic prophages of which DLP12 (20,525 bp) (Lindsey et al., 1989; Toba et al., 2011), Rac (22,890 bp) (Kaiser and Murray, 1979), Qin (20,456 bp) (Kaiser, 1980), e14 (15,204 bp) (Hiom et al., 1991; Mehta et al., 2004), and CPZ-55 (6,802 bp) prophages are well characterized. Several physiological functions such as stress tolerance, biofilm formation, antibiotic resistance, and other advantages to the host genome are attributed to the cryptic prophages in E. coli (Wang et al., 2010). These grounded prophages were suggested as potential drug targets (Wang and Wood, 2016). Analysis of locus-specific prophage sequences in completely sequenced Escherichia and Shigella strains (E. coli MG1655 strain as the reference) illustrates the significance of prophages in bacterial genome evolution and ecology at the level of a species. Examining the alignment of the specific locus in strains with and without prophage led to deduce the probable ancestral type locus and decipher the variations caused due to prophage integration. The tRNA genes are hotspots of phage integration based on sequence homology between the phage genome and tRNA gene (Campbell, 2003) and has been shown in case of integration of DLP12 prophage near the 3′ end of argU gene encoding tRNA-Arg. The probable recombinant attR site is similar to 47 bp of the argU gene, which was located on the other end of the prophage. The site-specific integration of e14 within the icd gene coding for isocitrate dehydrogenase protein (416aa) separated the 163 bp region from 3′ end of icd transcript. The integrated prophage provides the 3’ terminal codon of the icd ORF and thus forming homologous attL and attR sites [Supplementary Table S1 (Datasheet S1)]. Hence, the excision of e14 prophage results in restoration of functional icd gene yet rendering conservative replacement of aspartate by glutamate encoding codon (Hill et al., 1989). The site-specific integration of Rac prophage at the 5′ end of ttcA gene via the attB and attP sites resulted in the recombinant attL and attR sites. Variations in the amino acid sequence of the N-terminus of TtcA protein [a 311aa tRNA-cytidine(32) 2-sulfurtransferase protein involved in post-transcriptional thiolation of tRNA] was observed within the strains harboring Rac prophage. Such variations owing to sequence dissimilarity between the recombinant attR and attL sites reduce the fitness level in carbenicillin positive niche. The excision of Rac enhances biofilm formation, cell lysis, and motility (Hong et al., 2010; Liu et al., 2015). Integration of CPZ-55 prophage has occurred right after the stop codon of the eutA gene coding for ethanolamine utilization protein, causing target site duplication of 8 bp sequence (5′-TCAGGAAG-3′). The intergenic region of 13 bp (5′-TAAGTCGTTCCCT-3′) was conserved. No variation observed in prophage flanking genes, eutB (453 codons) and eutA (467 codons). However, CPZ-55 integration has separated the genes eutA and eutB and thereby disrupting the functionality of eut operon encoding proteins such as EutS, EutQ, EutN, EutD, EutM, EutE, EutH, and others, that are involved in ethanolamine reduction followed by utilization when available as the only nitrogen source. As the initial steps in the utilization of ethanolamine requires the expression of eutA, eutB, and eutC, any alterations as caused by prophage integration would result in the inability of eut- strains to multiply upon utilizing ethanolamine (Soupene et al., 2003) (Figure 4) [Supplementary Table S1 (Datasheet S1)].
Figure 4. Genetic variations caused by prophage integration in E. coli MG1655. The examination of specific integration loci of each of the five prophages (DLP12, e14, Rac, CPZ-55, and Qin) concerning E. coli MG1655 among the strains with and without corresponding prophages uncovered the genetic variations caused post integration event as shown above. (A) Integration of DLP12 via homologous recombination near the 3′ end of the argU gene. The 47 bp region similar to 3′ end of argU gene is located on the other end of the DLP12. (B) Integration of e14 at the 3′ end of icd gene causing the formation of pseudo icd’ (163 bp) gene as a duplicated region. (C) Site-specific integration of Rac at the 5′ end of ttcA gene. The variation in the N-terminus of TtcA protein was observed. (D) The transposon-like integration of CPZ-55 next to the stop codon of the eutA gene causing duplication of 8 bp region. No variations were seen in the bacterial flanking genes (eutB and eutA). (E) The site-specific integration of Qin at the 5′ end of ydfJ gene resulting in promoter loss and truncation (lacking 28 codons at the N terminal) of YdfJ protein.
We analyzed the distribution of the selected cryptic prophages within different Escherichia and Shigella strains varies. Locus-specific sequence analysis of the above mentioned five prophages (loci as per E. coli MG1655) within the completely sequenced Escherichia and Shigella strains (a total of 638 strains) shows that DLP12 is the most prevalent (74%) prophage followed by Rac prophage in 69% of total strains (Figure 5A). Qin prophage is present in 48%, and e14 prophage is encoded by 19%, while only 8% of total strains encoded CPZ-55 prophage. The prevalence in the E. coli genomes indicates that the integration of DLP12 and Rac have happened earlier relative to the other prophages analyzed here. Forty-six strains did not harbor any of the five prophages scored for. The analysis of co-occurrence of prophages within Escherichia and Shigella strains would shed some light on the acquisition of the prophages during the evolution of these bacteria. Hence, we analyzed the co-occurrence of the five prophages (DLP12, e14, Rac, CPZ-55 and Qin) with all possible permutations in 638 completely sequenced and annotated Escherichia and Shigella genomes (available as on March 2018) (Figure 5). Although there are 32 different permutations of the five prophages, only 19 combinations were observed in our analysis [Supplementary Table S3 (Datasheet S2)]. If prophages were acquired sequentially, it would be expected that the various observed permutations would have an incremental pattern of prophages within various genomes. 172 E. coli genomes had DLP12, Rac, and Qin cryptic prophages while 106 genomes had DLP12 prophage and none of the other four prophages. Thirty-nine E. coli genomes including most of the K12 strains had the DLP12, Rac, Qin, e14, and CPZ-55 permutation (Figure 5B). In essence, we did not notice any particular pattern that might imply serial acquisition of cryptic prophages. For example, DLP12 is present with various combinations of other prophages that were analyzed indicating random co-occurrence (Figure 5C). The observed permutations are mostly random which would mean that these cryptic prophages also propagated horizontally among the E. coli strains. Horizontal propagation is common among bacterial strains; however, the mechanism of HGT is most likely through conjugation considering the large size of the cryptic prophages.
Figure 5. Locus-specific distribution of the five prophages (with reference to E. coli MG1655 strain) in 638 completely sequenced E. coli strains. The cryptic prophages (DLP12, e14, Rac, CPZ-55, and Qin) of MG1655 strain were analyzed for their distribution within 638 completely sequenced E. coli strains taking the bacterial and prophage borders sequences as a query in BLASTn. Only cryptic prophages of MG1655 strain with location specificity were analyzed. Other strains may have various other prophages at multiple locations. (A) The prevalence of each cryptic prophage within various Escherichia and Shigella strains. The distribution was determined by curating the BLAST hits obtained in completely sequenced strains showing >85% sequence identity. (B) The prevalence of various permutations of prophage co-occurrence within various Escherichia and Shigella strains. (C) The diversity in the locus-specific prophage distribution across Escherichia and Shigella species. Only representative strains are indicated for illustration purposes.
From our analysis of the 36 genes of Qin prophage from E. coli MG1655 strain, genes like hokD, relBE, cold-shock genes (e.g., cspI, cspB, etc.), and regulators (e.g., dicC, dicB, etc.) were predicted to have been acquired via HGT [Supplementary Table S2 (Datasheet S1)]. In our sequence analyses, we found multiple mutations in intQ’ of most Qin prophage hits, which is indicative that the inactivation of the integrase gene could be a reason for the grounding of the Qin prophage. Although several Qin prophage genes are traditionally associated with general phage genes, numerous genes associated with Qin prophages are obtained by HGT. Several genes from plasmids and bacterial genomes (predominantly from E. coli, E. albertii, and Citrobacter) have integrated into the Qin prophage. The ancestral hokD ORF which is about 70 codons, however, is only 51 codons due to the integration of the relBE Toxin-Antitoxin system which is acquired from plasmids (Ramisetty and Santhosh, 2016). This observation emphasizes cautions to be exercised during characterization of genes that are highly linked to or within prophages because of high probability for genetic alterations. Comparisons of sequences of the same gene across several strains should be made to determine any truncations or modifications induced due to multiple genetic changes.
We then performed a comparative analysis of Qin prophage diversity within various Escherichia strains by estimating the length of prophage in base pairs (distance between the att sites), the number of Qin genes (with reference to MG1655 strain), and the number of de novo ORFs using a semantic script (hypothetical genes and DUFs). The size of the Qin prophage in 189 strains was highly divergent ranging from 11,280 bp in E. coli O6:H6 strain M9682-C1 to 97,450 bp in E. coli O26:H11 str. 11368 DNA (Figure 6). Surprisingly, we noticed a slight variation in Qin prophage size (of difference ∼1,300 bp) even within K-12 strains, which is indicative that prophages are highly prone to mutations. The number of copies of phage genes encoding major and minor capsid proteins, tail measure protein, host specificity protein, transposase, and other phage-related genes increased with the size of the Qin prophage. This observation indicates that cryptic prophages are hotspots for superintegration of phage/plasmid genomes (Lichens-Park et al., 1990). Conceivably, a positive correlation between the size of the prophage and the number of hypothetical proteins can be found. Some of the hypothetical genes could have accumulated as more plasmid/phage genes integrated into the Qin prophage. Nevertheless, there is a high probability for the generation of ORFs de novo due to multiple recombination events (Delaye et al., 2008; Tautz and Domazet-Lošo, 2011). Although the expressivity and functionality of most of the genes that originated through de novo genesis are unclear, they are genes in the making which maybe already functional or are yet to evolve for a beneficial function in the future generations.
Figure 6. The genetic diversity of Qin prophage in various E. coli strains. Correlation between the prophage size and the number of genes carried within the prophage is studied using BLAST tool and curated in excel. For each strain, the length of site- specific Qin prophage and the number of prophage encoding proteins were obtained using BLASTn tool taking gene sequence of each gene as query and hits are curated using excel. The dependence of prophage size and the flux of genes are represented in a scatter plot as shown taking prophage size on primary Y-axis, the number of genes on the prophage as secondary Y-axis and 189 strains on the X-axis. For Qin prophage analyses, the prophage coordinates were used to estimate the prophage size. Having E. coli MG1655 as the reference strain, the number of genes carried on Qin prophage among 189 strains was noted. The number of DUFs (Domain of Unknown Function) and hypothetical proteins in Qin prophage of each strain was fetched using a JavaScript to search the GenBank files of each prophage for “DUF” and “hypothetical proteins.” The graph was plotted to present the diversity of prophage in terms of its size, the number of genes carried, and novel genes in each of the strains.
To predict the possible role of the cryptic prophage in contributing to the ecological advantages to the host bacterium, we manually scouted for genes that are unrelated to the phage/plasmid in the Qin prophages of various E. coli strains [Supplementary Table S3 (Datasheet S2)]. The Qin prophage from strains examined carried diverse genes involved in significant functions such as transportation (e.g., Calcium transporter (chaC) in E. coli strain CRE1540), metabolism (e.g., Acyl-CoA thioesterase in E. coli O111:H-str. 11128 DNA), recombination (e.g., RusA endodeoxyribonuclease in E. coli strain ETEC-2265), virulence (e.g., Shiga toxin Stx2 subunits A and B in E. coli O178:H19 strain 2012C-4431), transcriptional regulation (e.g., LysR transcriptional regulator in E. coli strain CRE1540), and resistance (e.g., Tellurite/Colicin resistance protein in E. coli SE11 DNA) [Supplementary Table S3 (Datasheet S2)]. Using a semantic script, a total of 4,034 hypothetical proteins/DUFs were found in Qin prophages of 189 Escherichia and Shigella genomes indicating an average of 21 hypothetical genes per Qin prophage. Even after discounting 50% of the hypothetical genes for annotation and functionality issues, the high density and diversity of the hypothetical genes per prophage is indicative that grounded prophages are sites of de novo genesis which could expedite genome diversification and evolution (Berglund et al., 2010). In essence, grounded prophages of E. coli are hotspots of multiple recombination events, genetic flux, and de novo genesis of genes which expedite adaptations to various habitats.
Conclusion
The physiology of higher organisms is vastly influenced by their respective microbiomes. Microbiomes, in turn, are influenced by their phages (Scanlan, 2017; De Sordi et al., 2018). Bacteriophages are the major biological drivers of bacterial ecology and evolution through strategies such as symbiosis, dependency, and dormancy (Roossinck, 2011; Nasir et al., 2017). The prevalence of multiple grounded prophages within most bacterial genomes implies that the sequence of events such as lysogenization, grounding of prophages, domestication of phage genes, and the genetic events within the grounded prophages have an evolutionary impact on the selection of bacterial genomes. (1) Lysogeny, as opposed to the lytic cycle, represents an ‘evolutionarily stable strategy’ between bacteria and phages growing in an econiche to allow the existence of both the phage and the host. (2) The selection of lysogens is favored by the elimination of intraspecies competition by the phage and the immunity conferred to lysogens by the prophages. (3) The lysogens with grounded prophages are selected over lysogens (with wild-type prophages) because the latter have an imminent threat of prophage activation followed by lysis. The grounded prophages may provide immunity from specific phage attack, ecologically significant traits, and the evolution of new genes (Figure 7). Many genes that provide stress tolerance, antibiotic resistance, virulence genes, and metabolic genes acquired horizontally are located within prophages. (4) Grounded prophages are zones within genomes that allow multiple recombinations and mutations with minimal deleterious effects on the host physiology. The increased probability of genetic modifications (insertions, deletions, and inversions) within the prophages increases the probability of formation of new ORFs. These de novo genes maybe refined gradually for a function with a selective advantage such as antibiotic resistance and virulence.
Figure 7. A genetic history and advantages of cryptic prophages. The genome of temperate phages integrates with the bacterial chromosome through site-specific recombination between sites similar to those of attP and attB. Occasionally, mutations in the attL, attR, recombinase genes, and genes encoding other essential accessory factors would result in a ‘grounded’ prophage. Eventually but gradually, more mutations, usually deletions, accumulate to ease the metabolic and genetic burden imposed by the cryptic prophage. The cryptic prophage offers several advantages to the host bacterium. The genetic variations caused in the host genome due to integrations, inversions involving prophage sequences, and the action of their recombinases may have a selective advantage. The prophage confers immunity from secondary phage infection and acts as hotspots for gene influx through horizontal gene transfer of important genes such as those involved in metabolism, antibiotic resistance, virulence, and stress tolerance. More importantly, the cryptic prophage is a site for the emergence of new ORFs that may have novel and useful functions.
Since the evolution of the first temperate phages, prophages must have enhanced the rate of evolution in different life forms including higher organisms. Grounded prophages are the most significant drivers of ecology, genetic diversification, and microbial genome evolution. The grounded prophages, and the analogous elements in other cellular life-forms, are the ‘holy grails’ of de novo genesis of genetic information, which must be of immense interest to molecular evolutionists.
Author Contributions
BCMR has conceived the idea, designed the work, performed the analyses, interpreted the results, and wrote the manuscript. PAS has performed the work, interpreted the results, and wrote the manuscript.
Funding
The work in BCMR’s lab is supported by Prof. T. R. Rangarajan grants (internal funding) and the Central Research Facility provided by SASTRA Deemed University, Thanjavur, India.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
The authors thank Sri Uma Aisvarya Balasubramanian for proofreading the manuscript. The authors also thank the reviewers for their valuable suggestions.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2019.00065/full#supplementary-material
References
Abe, K., Kawano, Y., Iwamoto, K., Arai, K., Maruyama, Y., Eichenberger, P., et al. (2014). Developmentally-regulated excision of the SPbeta prophage reconstitutes a gene required for spore envelope maturation in Bacillus subtilis. PLoS Genet. 10:e1004636. doi: 10.1371/journal.pgen.1004636
Allet, B. (1979). Mu insertion duplicates a 5 base pair sequence at the host inserted site. Cell 16, 123–129. doi: 10.1016/0092-8674(79)90193-4
Asadulghani, M., Ogura, Y., Ooka, T., Itoh, T., Sawaguchi, A., Iguchi, A., et al. (2009). The defective prophage pool of Escherichia coli O157: prophage-prophage interactions potentiate horizontal transfer of virulence determinants. PLoS Pathog. 5:e1000408. doi: 10.1371/journal.ppat.1000408
Askora, A., Kawasaki, T., Usami, S., Fujie, M., and Yamada, T. (2009). Host recognition and integration of filamentous phage phiRSM in the phytopathogen, Ralstonia solanacearum. Virology 384, 69–76. doi: 10.1016/j.virol.2008.11.007
Bartels, M. D., Boye, K., Rohde, S. M., Larsen, A. R., Torfs, H., Bouchy, P., et al. (2009). A common variant of staphylococcal cassette chromosome mec type IVa in isolates from Copenhagen, Denmark, is not detected by the BD GeneOhm methicillin-resistant Staphylococcus aureus assay. J. Clin. Microbiol. 47, 1524–1527. doi: 10.1128/JCM.02153-08
Berglund, E. C., Ellegaard, K., Granberg, F., Xie, Z., Maruyama, S., Kosoy, M. Y., et al. (2010). Rapid diversification by recombination in Bartonella grahamii from wild rodents in Asia contrasts with low levels of genomic divergence in Northern Europe and America. Mol. Ecol. 19, 2241–2255. doi: 10.1111/j.1365-294X.2010.04646.x
Betts, A., Gray, C., Zelek, M., Maclean, R. C., and King, K. C. (2018). High parasite diversity accelerates host adaptation and diversification. Science 360, 907–911. doi: 10.1126/science.aam9974
Bibb, L. A., Hancox, M. I., and Hatfull, G. F. (2005). Integration and excision by the large serine recombinase phiRv1 integrase. Mol. Microbiol. 55, 1896–1910. doi: 10.1111/j.1365-2958.2005.04517.x
Blackstone, N. W., and Green, D. R. (1999). The evolution of a mechanism of cell suicide. Bioessays 21, 84–88. doi: 10.1002/(SICI)1521-1878(199901)21:1<84::AID-BIES11>3.0.CO;2-0
Blokesch, M. (2017). In and out-contribution of natural transformation to the shuffling of large genomic regions. Curr. Opin. Microbiol. 38, 22–29. doi: 10.1016/j.mib.2017.04.001
Bobay, L. M., Rocha, E. P., and Touchon, M. (2013). The adaptation of temperate bacteriophages to their host genomes. Mol. Biol. Evol. 30, 737–751. doi: 10.1093/molbev/mss279
Bobay, L. M., Touchon, M., and Rocha, E. P. (2014). Pervasive domestication of defective prophages by bacteria. Proc. Natl. Acad. Sci. U.S.A. 111, 12127–12132. doi: 10.1073/pnas.1405336111
Bonneau, R., Baliga, N. S., Deutsch, E. W., Shannon, P., and Hood, L. (2004). Comprehensive de novo structure prediction in a systems-biology context for the archaea Halobacterium sp. NRC-1. Genome Biol. 5:R52. doi: 10.1186/gb-2004-5-8-r52
Braga, L. P. P., Soucy, S. M., Amgarten, D. E., Da Silva, A. M., and Setubal, J. C. (2018). Bacterial diversification in the light of the interactions with phages: the genetic symbionts and their role in ecological speciation. Front. Ecol. Evol. 6:6. doi: 10.3389/fevo.2018.00006
Braid, M. D., Silhavy, J. L., Kitts, C. L., Cano, R. J., and Howe, M. M. (2004). Complete genomic sequence of bacteriophage B3, a Mu-like phage of Pseudomonas aeruginosa. J. Bacteriol. 186, 6560–6574. doi: 10.1128/JB.186.19.6560-6574.2004
Brody, H., Greener, A., and Hill, C. W. (1985). Excision and reintegration of the Escherichia coli K-12 chromosomal element e14. J. Bacteriol. 161, 1112–1117.
Brody, H., and Hill, C. W. (1988). Attachment site of the genetic element e14. J. Bacteriol. 170, 2040–2044. doi: 10.1128/jb.170.5.2040-2044.1988
Brussow, H., Canchaya, C., and Hardt, W. D. (2004). Phages and the evolution of bacterial pathogens: from genomic rearrangements to lysogenic conversion. Microbiol. Mol. Biol. Rev. 68, 560–602. doi: 10.1128/MMBR.68.3.560-602.2004
Brussow, H., and Desiere, F. (2001). Comparative phage genomics and the evolution of Siphoviridae: insights from dairy phages. Mol. Microbiol. 39, 213–222. doi: 10.1046/j.1365-2958.2001.02228.x
Buckling, A., and Brockhurst, M. (2012). Bacteria-virus coevolution. Adv. Exp. Med. Biol. 751, 347–370. doi: 10.1007/978-1-4614-3567-9_16
Burrows, L. L., Charter, D. F., and Lam, J. S. (1996). Molecular characterization of the Pseudomonas aeruginosa serotype O5 (PAO1) B-band lipopolysaccharide gene cluster. Mol. Microbiol. 22, 481–495. doi: 10.1046/j.1365-2958.1996.1351503.x
Camilli, R., Bonnal, R. J., Del Grosso, M., Iacono, M., Corti, G., Rizzi, E., et al. (2011). Complete genome sequence of a serotype 11A, ST62 Streptococcus pneumoniae invasive isolate. BMC Microbiol. 11:25. doi: 10.1186/1471-2180-11-25
Campbell, A. (2003). Prophage insertion sites. Res. Microbiol. 154, 277–282. doi: 10.1016/S0923-2508(03)00071-8
Canchaya, C., Desiere, F., Mcshan, W. M., Ferretti, J. J., Parkhill, J., and Brussow, H. (2002). Genome analysis of an inducible prophage and prophage remnants integrated in the Streptococcus pyogenes strain SF370. Virology 302, 245–258. doi: 10.1006/viro.2002.1570
Canchaya, C., Proux, C., Fournous, G., Bruttin, A., and Brussow, H. (2003). Prophage genomics. Microbiol. Mol. Biol. Rev. 67, 238–276. doi: 10.1128/MMBR.67.2.238-276.2003
Casjens, S. (2003). Prophages and bacterial genomics: what have we learned so far? Mol. Microbiol. 49, 277–300. doi: 10.1046/j.1365-2958.2003.03580.x
Chakraborty, D., Clark, E., Mauro, S. A., and Koudelka, G. B. (2018). Molecular mechanisms governing “hair-trigger” induction of shiga toxin-encoding prophages. Viruses 10:E228. doi: 10.3390/v10050228
Coombes, B. K., Wickham, M. E., Brown, N. F., Lemire, S., Bossi, L., Hsiao, W. W., et al. (2005). Genetic and molecular analysis of GogB, a phage-encoded type III-secreted substrate in Salmonella enterica serovar typhimurium with autonomous expression from its associated phage. J. Mol. Biol. 348, 817–830. doi: 10.1016/j.jmb.2005.03.024
Craig, N. L., and Nash, H. A. (1983). The mechanism of phage lambda site-specific recombination: site-specific breakage of DNA by Int topoisomerase. Cell 35, 795–803. doi: 10.1016/0092-8674(83)90112-5
Cumby, N., Edwards, A. M., Davidson, A. R., and Maxwell, K. L. (2012). The bacteriophage HK97 gp15 moron element encodes a novel superinfection exclusion protein. J. Bacteriol. 194, 5012–5019. doi: 10.1128/JB.00843-12
Darmon, E., and Leach, D. R. (2014). Bacterial genome instability. Microbiol. Mol. Biol. Rev. 78, 1–39. doi: 10.1128/MMBR.00035-13
Davison, J. (1999). Genetic exchange between bacteria in the environment. Plasmid 42, 73–91. doi: 10.1006/plas.1999.1421
Delaye, L., Deluna, A., Lazcano, A., and Becerra, A. (2008). The origin of a novel gene through overprinting in Escherichia coli. BMC Evol. Biol. 8:31. doi: 10.1186/1471-2148-8-31
De Paepe, M., Hutinet, G., Son, O., Amarir-Bouhram, J., Schbath, S., and Petit, M. A. (2014). Temperate phages acquire DNA from defective prophages by relaxed homologous recombination: the role of Rad52-like recombinases. PLoS Genet. 10:e1004181. doi: 10.1371/journal.pgen.1004181
De Sordi, L., Lourenco, M., and Debarbieux, L. (2018). “I will survive”: a tale of bacteriophage-bacteria coevolution in the gut. Gut Microbes 10, 92–99. doi: 10.1080/19490976.2018.1474322
Deghorain, M., Bobay, L. M., Smeesters, P. R., Bousbata, S., Vermeersch, M., Perez-Morga, D., et al. (2012). Characterization of novel phages isolated in coagulase-negative staphylococci reveals evolutionary relationships with Staphylococcus aureus phages. J. Bacteriol. 194, 5829–5839. doi: 10.1128/JB.01085-12
Desiere, F., Mcshan, W. M., Van Sinderen, D., Ferretti, J. J., and Brussow, H. (2001). Comparative genomics reveals close genetic relationships between phages from dairy bacteria and pathogenic Streptococci: evolutionary implications for prophage-host interactions. Virology 288, 325–341. doi: 10.1006/viro.2001.1085
Dobrindt, U., and Hacker, J. (2001). Whole genome plasticity in pathogenic bacteria. Curr. Opin. Microbiol. 4, 550–557. doi: 10.1016/S1369-5274(00)00250-2
Dobrindt, U., Zdziarski, J., Salvador, E., and Hacker, J. (2010). Bacterial genome plasticity and its impact on adaptation during persistent infection. Int. J. Med. Microbiol. 300, 363–366. doi: 10.1016/j.ijmm.2010.04.010
Dorgai, L., Oberto, J., and Weisberg, R. A. (1993). Xis and Fis proteins prevent site-specific DNA inversion in lysogens of phage HK022. J. Bacteriol. 175, 693–700. doi: 10.1128/jb.175.3.693-700.1993
D’Souza, G., Waschina, S., Kaleta, C., and Kost, C. (2015). Plasticity and epistasis strongly affect bacterial fitness after losing multiple metabolic genes. Evolution 69, 1244–1254. doi: 10.1111/evo.12640
Dunning Hotopp, J. C., Grifantini, R., Kumar, N., Tzeng, Y. L., Fouts, D., Frigimelica, E., et al. (2006). Comparative genomics of Neisseria meningitidis: core genome, islands of horizontal transfer and pathogen-specific genes. Microbiology 152, 3733–3749. doi: 10.1099/mic.0.29261-0
Durand, P. M., Sym, S., and Michod, R. E. (2016). Programmed cell death and complexity in microbial systems. Curr. Biol. 26, R587–R593. doi: 10.1016/j.cub.2016.05.057
Dydecka, A., Bloch, S., Rizvi, A., Perez, S., Nejman-Falenczyk, B., Topka, G., et al. (2017). Bad phages in good bacteria: role of the mysterious orf63 of lambda and shiga toxin-converting Phi24B Bacteriophages. Front. Microbiol. 8:1618. doi: 10.3389/fmicb.2017.01618
Farrant, J. L., Sansone, A., Canvin, J. R., Pallen, M. J., Langford, P. R., Wallis, T. S., et al. (1997). Bacterial copper- and zinc-cofactored superoxide dismutase contributes to the pathogenesis of systemic salmonellosis. Mol. Microbiol. 25, 785–796. doi: 10.1046/j.1365-2958.1997.5151877.x
Fierer, J., and Guiney, D. G. (2001). Diverse virulence traits underlying different clinical outcomes of Salmonella infection. J. Clin. Invest. 107, 775–780. doi: 10.1172/JCI12561
Finkel, S. E., and Kolter, R. (2001). DNA as a nutrient: novel role for bacterial competence gene homologs. J. Bacteriol. 183, 6288–6293. doi: 10.1128/JB.183.21.6288-6293.2001
Fischer, W., Hofreuter, D., and Haas, R. (2001). “Natural transformation, recombination, and repair,” in Helicobacter pylori: Physiology and Genetics, eds H. L. T. Mobley, G. L. Mendz, and S. L. Hazell (Washington, DC: ASM Press).
Fogg, P. C. M., Allison, H. E., Saunders, J. R., and Mccarthy, A. J. (2010). Bacteriophage lambda: a paradigm revisited. J. Virol. 84, 6876–6879. doi: 10.1128/JVI.02177-09
Forde, S. E., Thompson, J. N., and Bohannan, B. J. (2004). Adaptation varies through space and time in a coevolving host-parasitoid interaction. Nature 431, 841–844. doi: 10.1038/nature02906
Fujiwara, A., Kawasaki, T., Usami, S., Fujie, M., and Yamada, T. (2008). Genomic characterization of Ralstonia solanacearum phage phiRSA1 and its related prophage (phiRSX) in strain GMI1000. J. Bacteriol. 190, 143–156. doi: 10.1128/JB.01158-07
Gandon, S. (2016). Why be temperate: lessons from Bacteriophage lambda. Trends Microbiol. 24, 356–365. doi: 10.1016/j.tim.2016.02.008
Gaudriault, S., Duchaud, E., Lanois, A., Canoy, A. S., Bourot, S., Derose, R., et al. (2006). Whole-genome comparison between Photorhabdus strains to identify genomic regions involved in the specificity of nematode interaction. J. Bacteriol. 188, 809–814. doi: 10.1128/JB.188.2.809-814.2006
Ghosh, P., Wasil, L. R., and Hatfull, G. F. (2006). Control of phage Bxb1 excision by a novel recombination directionality factor. PLoS Biol. 4:e186. doi: 10.1371/journal.pbio.0040186
Gottesman, M. E., and Weisberg, R. A. (2004). Little lambda, who made thee? Microbiol. Mol. Biol. Rev. 68, 796–813. doi: 10.1128/MMBR.68.4.796-813.2004
Groth, A. C., and Calos, M. P. (2004). Phage integrases: biology and applications. J. Mol. Biol. 335, 667–678. doi: 10.1016/j.jmb.2003.09.082
Gruenheid, S., Sekirov, I., Thomas, N. A., Deng, W., O’donnell, P., Goode, D., et al. (2004). Identification and characterization of NleA, a non-LEE-encoded type III translocated virulence factor of enterohaemorrhagic Escherichia coli O157:H7. Mol. Microbiol. 51, 1233–1249. doi: 10.1046/j.1365-2958.2003.03911.x
Hacker, J., and Carniel, E. (2001). Ecological fitness, genomic islands and bacterial pathogenicity. A Darwinian view of the evolution of microbes. EMBO Rep. 2, 376–381. doi: 10.1093/embo-reports/kve097
Hargreaves, K. R., Kropinski, A. M., and Clokie, M. R. (2014). What does the talking?: Quorum sensing signalling genes discovered in a bacteriophage genome. PLoS One 9:e85131. doi: 10.1371/journal.pone.0085131
Harshey, R. M. (2014). Transposable phage Mu. Microbiol. Spectr. 2, 669–691. doi: 10.1128/microbiolspec.MDNA3-0007-2014
Hill, C. W., Gray, J. A., and Brody, H. (1989). Use of the isocitrate dehydrogenase structural gene for attachment of e14 in Escherichia coli K-12. J. Bacteriol. 171, 4083–4084. doi: 10.1128/jb.171.7.4083-4084.1989
Hiom, K., Thomas, S. M., and Sedgwick, S. G. (1991). Different mechanisms for SOS induced alleviation of DNA restriction in Escherichia coli. Biochimie 73, 399–405. doi: 10.1016/0300-9084(91)90106-B
Hoess, R. H., and Abremski, K. (1985). Mechanism of strand cleavage and exchange in the Cre-lox site-specific recombination system. J. Mol. Biol. 181, 351–362. doi: 10.1016/0022-2836(85)90224-4
Hogardt, M., Hoboth, C., Schmoldt, S., Henke, C., Bader, L., and Heesemann, J. (2007). Stage-specific adaptation of hypermutable Pseudomonas aeruginosa isolates during chronic pulmonary infection in patients with cystic fibrosis. J. Infect. Dis. 195, 70–80. doi: 10.1086/509821
Hong, S. H., Wang, X., and Wood, T. K. (2010). Controlling biofilm formation, prophage excision and cell death by rewiring global regulator H-NS of Escherichia coli. Microb. Biotechnol. 3, 344–356. doi: 10.1111/j.1751-7915.2010.00164.x
Howard-Varona, C., Hargreaves, K. R., Abedon, S. T., and Sullivan, M. B. (2017). Lysogeny in nature: mechanisms, impact and ecology of temperate phages. ISME J. 11, 1511–1520. doi: 10.1038/ismej.2017.16
Hurwitz, B. L., Brum, J. R., and Sullivan, M. B. (2015). Depth-stratified functional and taxonomic niche specialization in the ’core’ and ’flexible’ Pacific Ocean Virome. ISME J. 9, 472–484. doi: 10.1038/ismej.2014.143
Hyder, S. L., and Streitfeld, M. M. (1978). Transfer of erythromycin resistance from clinically isolated lysogenic strains of Streptococcus pyogenes via their endogenous phage. J. Infect. Dis. 138, 281–286. doi: 10.1093/infdis/138.3.281
Iguchi, A., Iyoda, S., Terajima, J., Watanabe, H., and Osawa, R. (2006). Spontaneous recombination between homologous prophage regions causes large-scale inversions within the Escherichia coli O157:H7 chromosome. Gene 372,199–207. doi: 10.1016/j.gene.2006.01.005
Iversen, H., L’Abée-Lund, T. M., Aspholm, M., Arnesen, L. P., and Lindback, T. (2015). Commensal E. coli Stx2 lysogens produce high levels of phages after spontaneous prophage induction. Front. Cell. Infect. Microbiol. 5:5. doi: 10.3389/fcimb.2015.00005
Juhala, R. J., Ford, M. E., Duda, R. L., Youlton, A., Hatfull, G. F., and Hendrix, R. W. (2000). Genomic sequences of bacteriophages HK97 and HK022: pervasive genetic mosaicism in the lambdoid bacteriophages. J. Mol. Biol. 299, 27–51. doi: 10.1006/jmbi.2000.3729
Kaiser, K. (1980). The origin of Q-independent derivatives of phage lambda. Mol. Gen Genet. 179, 547–554. doi: 10.1007/BF00271744
Kaiser, K., and Murray, N. E. (1979). Physical characterisation of the “Rac prophage” in E. coli K12. Mol. Gen. Genet. 175, 159–174. doi: 10.1007/BF00425532
Katayama, Y., Baba, T., Sekine, M., Fukuda, M., and Hiramatsu, K. (2013). Beta-hemolysin promotes skin colonization by Staphylococcus aureus. J. Bacteriol. 195, 1194–1203. doi: 10.1128/JB.01786-12
Kim, M. S., and Bae, J. W. (2018). Lysogeny is prevalent and widely distributed in the murine gut microbiota. ISME J. 12, 1127–1141. doi: 10.1038/s41396-018-0061-9
Koga, M., Otsuka, Y., Lemire, S., and Yonesaki, T. (2011). Escherichia coli rnlA and rnlB compose a novel toxin-antitoxin system. Genetics 187, 123–130. doi: 10.1534/genetics.110.121798
Koskella, B., and Brockhurst, M. A. (2014). Bacteria-phage coevolution as a driver of ecological and evolutionary processes in microbial communities. FEMS Microbiol. Rev. 38, 916–931. doi: 10.1111/1574-6976.12072
Kraft, C., Stack, A., Josenhans, C., Niehus, E., Dietrich, G., Correa, P., et al. (2006). Genomic changes during chronic Helicobacter pylori infection. J. Bacteriol. 188, 249–254. doi: 10.1128/JB.188.1.249-254.2006
Kutsukake, K., and Iino, T. (1980). A trans-acting factor mediates inversion of a specific DNA segment in flagellar phase variation of Salmonella. Nature 284, 479–481. doi: 10.1038/284479a0
Landy, A., and Ross, W. (1977). Viral integration and excision: structure of the lambda att sites: DNA sequences have been determined for regions involved in lambda site-specific recombination. Science (New York, N.Y.) 197, 1147–1160. doi: 10.1126/science.331474
Lawrence, J. G., Hendrix, R. W., and Casjens, S. (2001). Where are the pseudogenes in bacterial genomes? Trends Microbiol. 9, 535–540. doi: 10.1016/S0966-842X(01)02198-9
Lee, C. Y., and Iandolo, J. J. (1985). Mechanism of bacteriophage conversion of lipase activity in Staphylococcus aureus. J. Bacteriol. 164, 288–293.
Lichens-Park, A., Smith, C. L., and Syvanen, M. (1990). Integration of bacteriophage lambda into the cryptic lambdoid prophages of Escherichia coli. J. Bacteriol. 172, 2201–2208. doi: 10.1128/jb.172.5.2201-2208.1990
Lindsey, D. F., Mullin, D. A., and Walker, J. R. (1989). Characterization of the cryptic lambdoid prophage DLP12 of Escherichia coli and overlap of the DLP12 integrase gene with the tRNA gene argU. J. Bacteriol. 171, 6197–6205. doi: 10.1128/jb.171.11.6197-6205.1989
Liu, X., Li, Y., Guo, Y., Zeng, Z., Li, B., Wood, T. K., et al. (2015). Physiological function of rac prophage during biofilm formation and regulation of rac excision in Escherichia coli K-12. Sci. Rep. 5:16074. doi: 10.1038/srep16074
Lopez, C. A., Winter, S. E., Rivera-Chavez, F., Xavier, M. N., Poon, V., Nuccio, S. P., et al. (2012). Phage-mediated acquisition of a type III secreted effector protein boosts growth of Salmonella by nitrate respiration. MBio 3:e00143-12. doi: 10.1128/mBio.00143-12
Los, J., Los, M., Wegrzyn, A., and Wegrzyn, G. (2013). Altruism of Shiga toxin-producing Escherichia coli: recent hypothesis versus experimental results. Front. Cell. Infect. Microbiol. 2:166. doi: 10.3389/fcimb.2012.00166
Los, J. M., Los, M., Wegrzyn, A., and Wegrzyn, G. (2010). Hydrogen peroxide-mediated induction of the Shiga toxin-converting lambdoid prophage ST2-8624 in Escherichia coli O157:H7. FEMS Immunol. Med. Microbiol. 58, 322–329. doi: 10.1111/j.1574-695X.2009.00644.x
Luo, C. H., Chiou, P. Y., Yang, C. Y., and Lin, N. T. (2012). Genome, integration, and transduction of a novel temperate phage of Helicobacter pylori. J. Virol. 86, 8781–8792. doi: 10.1128/JVI.00446-12
Mahony, J., Mcgrath, S., Fitzgerald, G. F., and Van Sinderen, D. (2008). Identification and characterization of lactococcal-prophage-carried superinfection exclusion genes. Appl. Environ. Microbiol. 74, 6206–6215. doi: 10.1128/AEM.01053-08
Martinez-Penafiel, E., Fernandez-Ramirez, F., Ishida, C., Reyes-Cortes, R., Sepulveda-Robles, O., Guarneros-Pena, G., et al. (2012). Overexpression of Ipe protein from the coliphage mEp021 induces pleiotropic effects involving haemolysis by HlyE-containing vesicles and cell death. Biochimie 94,1262–1273. doi: 10.1016/j.biochi.2012.02.004
Mata, M., Delstanche, M., and Robert-Baudouy, J. (1978). Isolation of specialized transducing bacteriophages carrying the structural genes of the hexuronate system in Escherichia coli K-12: exu region. J. Bacteriol. 133, 549–557.
Maughan, H., Wilson, L. A., and Redfield, R. J. (2010). Bacterial DNA uptake sequences can accumulate by molecular drive alone. Genetics 186, 613–627. doi: 10.1534/genetics.110.119438
McDonough, M. A., and Butterton, J. R. (1999). Spontaneous tandem amplification and deletion of the Shiga toxin operon in Shigella dysenteriae 1. Mol. Microbiol. 34, 1058–1069. doi: 10.1046/j.1365-2958.1999.01669.x
McGrath, S., Fitzgerald, G. F., and Van Sinderen, D. (2002). Identification and characterization of phage-resistance genes in temperate lactococcal bacteriophages. Mol. Microbiol. 43, 509–520. doi: 10.1046/j.1365-2958.2002.02763.x
Mehta, P., Casjens, S., and Krishnaswamy, S. (2004). Analysis of the lambdoid prophage element e14 in the E. coli K-12 genome. BMC Microbiol. 4:4. doi: 10.1186/1471-2180-4-4
Mell, J. C., and Redfield, R. J. (2014). Natural competence and the evolution of DNA uptake specificity. J. Bacteriol. 196, 1471–1483. doi: 10.1128/JB.01293-13
Menouni, R., Hutinet, G., Petit, M. A., and Ansaldi, M. (2015). Bacterial genome remodeling through bacteriophage recombination. FEMS Microbiol. Lett. 362, 1–10. doi: 10.1093/femsle/fnu022
Mirold, S., Rabsch, W., Tschape, H., and Hardt, W. D. (2001). Transfer of the Salmonella type III effector sopE between unrelated phage families. J. Mol. Biol. 312, 7–16. doi: 10.1006/jmbi.2001.4950
Moore, R. A., Reckseidler-Zenteno, S., Kim, H., Nierman, W., Yu, Y., Tuanyok, A., et al. (2004). Contribution of gene loss to the pathogenic evolution of Burkholderia pseudomallei and Burkholderia mallei. Infect. Immun. 72,4172–4187. doi: 10.1128/IAI.72.7.4172-4187.2004
Nadeem, A., and Wahl, L. M. (2017). Prophage as a genetic reservoir: promoting diversity and driving innovation in the host community. Evolution 71,2080–2089. doi: 10.1111/evo.13287
Nasir, A., Kim, K. M., and Caetano-Anolles, G. (2017). Long-term evolution of viruses: a Janus-faced balance. Bioessays 39:1700026. doi: 10.1002/bies.201700026
Nilsson, A. S., Karlsson, J. L., and Haggard-Ljungquist, E. (2004). Site-specific recombination links the evolution of P2-like coliphages and pathogenic enterobacteria. Mol. Biol. Evol. 21, 1–13. doi: 10.1093/molbev/msg223
Obeng, N., Pratama, A. A., and Elsas, J. D. V. (2016). The significance of mutualistic phages for bacterial ecology and evolution. Trends Microbiol. 24, 440–449. doi: 10.1016/j.tim.2015.12.009
Ogier, J. C., Calteau, A., Forst, S., Goodrich-Blair, H., Roche, D., Rouy, Z., et al. (2010). Units of plasticity in bacterial genomes: new insight from the comparative genomics of two bacteria interacting with invertebrates, Photorhabdus and Xenorhabdus. BMC Genomics 11:568. doi: 10.1186/1471-2164-11-568
Osborne, S. E., Walthers, D., Tomljenovic, A. M., Mulder, D. T., Silphaduang, U., Duong, N., et al. (2009). Pathogenic adaptation of intracellular bacteria by rewiring a cis-regulatory input function. Proc. Natl. Acad. Sci. U.S.A. 106, 3982–3987. doi: 10.1073/pnas.0811669106
Palchevskiy, V., and Finkel, S. E. (2006). Escherichia coli competence gene homologs are essential for competitive fitness and the use of DNA as a nutrient. J. Bacteriol. 188, 3902–3910. doi: 10.1128/JB.01974-05
Papagiannis, C. V., Sam, M. D., Abbani, M. A., Yoo, D., Cascio, D., Clubb, R. T., et al. (2007). Fis targets assembly of the Xis nucleoprotein filament to promote excisive recombination by phage lambda. J. Mol. Biol. 367, 328–343. doi: 10.1016/j.jmb.2006.12.071
Parma, D. H., Snyder, M., Sobolevski, S., Nawroz, M., Brody, E., and Gold, L. (1992). The rex system of Bacteriophage lambda: tolerance and altruistic cell death. Genes Dev. 6, 497–510. doi: 10.1101/gad.6.3.497
Paterson, S., Vogwill, T., Buckling, A., Benmayor, R., Spiers, A. J., Thomson, N. R., et al. (2010). Antagonistic coevolution accelerates molecular evolution. Nature 464, 275–278. doi: 10.1038/nature08798
Peeters, S. H., and de Jonge, M. I. (2018). For the greater good: programmed cell death in bacterial communities. Microbiol. Res. 207, 161–169. doi: 10.1016/j.micres.2017.11.016
Picozzi, C., Meissner, D., Chierici, M., Ehrmann, M. A., Vigentini, I., Foschino, R., et al. (2015). Phage-mediated transfer of a dextranase gene in Lactobacillus sanfranciscensis and characterization of the enzyme. Int. J. Food Microbiol. 202, 48–53. doi: 10.1016/j.ijfoodmicro.2015.02.018
Ramisetty, B. C., Natarajan, B., and Santhosh, R. S. (2015). mazEF-mediated programmed cell death in bacteria: “what is this?”. Crit. Rev. Microbiol. 41, 89–100. doi: 10.3109/1040841X.2013.804030
Ramisetty, B. C., and Santhosh, R. S. (2016). Horizontal gene transfer of chromosomal Type II toxin-antitoxin systems of Escherichia coli. FEMS Microbiol. Lett. 363:fnv238. doi: 10.1093/femsle/fnv238
Rangel, A., Steenbergen, S. M., and Vimr, E. R. (2016). Unexpected diversity of Escherichia coli sialate O-acetyl esterase NanS. J. Bacteriol. 198, 2803–2809. doi: 10.1128/JB.00189-16
Roossinck, M. J. (2011). The good viruses: viral mutualistic symbioses. Nat. Rev. Microbiol. 9, 99–108. doi: 10.1038/nrmicro2491
Ross, W., Landy, A., Kikuchi, Y., and Nash, H. (1979). Interaction of int protein with specific sites on lambda att DNA. Cell 18, 297–307. doi: 10.1016/0092-8674(79)90049-7
Roux, S., Brum, J. R., Dutilh, B. E., Sunagawa, S., Duhaime, M. B., Loy, A., et al. (2016). Ecogenomics and potential biogeochemical impacts of globally abundant ocean viruses. Nature 537, 689–693. doi: 10.1038/nature19366
Roux, S., Enault, F., Hurwitz, B. L., and Sullivan, M. B. (2015a). VirSorter: mining viral signal from microbial genomic data. PeerJ 3:e985. doi: 10.7717/peerj.985
Roux, S., Hallam, S. J., Woyke, T., and Sullivan, M. B. (2015b). Viral dark matter and virus-host interactions resolved from publicly available microbial genomes. eLife 4:e08490. doi: 10.7554/eLife.08490
Saile, N., Voigt, A., Kessler, S., Stressler, T., Klumpp, J., Fischer, L., et al. (2016). Escherichia coli O157:H7 strain EDL933 harbors multiple functional prophage-associated genes necessary for the utilization of 5-N-Acetyl-9-O-Acetyl neuraminic acid as a growth substrate. Appl. Environ. Microbiol. 82, 5940–5950. doi: 10.1128/AEM.01671-16
Scanlan, P. D. (2017). Bacteria-bacteriophage coevolution in the human gut: implications for microbial diversity and functionality. Trends Microbiol. 25, 614–623. doi: 10.1016/j.tim.2017.02.012
Schindler, D., and Echols, H. (1981). Retroregulation of the int gene of Bacteriophage lambda: control of translation completion. Proc. Natl. Acad. Sci. U.S.A. 78, 4475–4479. doi: 10.1073/pnas.78.7.4475
Schmucker, R., Ritthaler, W., Stern, B., and Kamp, D. (1986). DNA inversion in bacteriophage Mu: characterization of the inversion site. J. Gen. Virol. 67(Pt 6), 1123–1133. doi: 10.1099/0022-1317-67-6-1123
Sechaud, L., Cluzel, P. J., Rousseau, M., Baumgartner, A., and Accolas, J. P. (1988). Bacteriophages of lactobacilli. Biochimie 70, 401–410. doi: 10.1016/0300-9084(88)90214-3
Sekulovic, O., and Fortier, L. C. (2015). Global transcriptional response of Clostridium difficile carrying the CD38 prophage. Appl. Environ. Microbiol. 81, 1364–1374. doi: 10.1128/AEM.03656-14
Shimada, K., Weisberg, R. A., and Gottesman, M. E. (1975). Prophage lambda at unusual chromosomal locations. III. The components of the secondary attachment sites. J. Mol. Biol. 93, 415–429. doi: 10.1016/0022-2836(75)90237-5
Shimizu, T., Ohta, Y., and Noda, M. (2009). Shiga toxin 2 is specifically released from bacterial cells by two different mechanisms. Infect. Immun. 77, 2813–2823. doi: 10.1128/IAI.00060-09
Silverman, M., and Simon, M. (1980). Phase variation: genetic analysis of switching mutants. Cell 19, 845–854. doi: 10.1016/0092-8674(80)90075-6
Sinha, S., and Redfield, R. J. (2012). Natural DNA uptake by Escherichia coli. PLoS One 7:e35620. doi: 10.1371/journal.pone.0035620
Song, Z., and Luo, L. (2012). Escherichia coli mutants induced by multi-ion irradiation. J. Radiat. Res. 53, 854–859. doi: 10.1093/jrr/rrs061
Soupene, E., Van Heeswijk, W. C., Plumbridge, J., Stewart, V., Bertenthal, D., Lee, H., et al. (2003). Physiological studies of Escherichia coli strain MG1655: growth defects and apparent cross-regulation of gene expression. J. Bacteriol. 185, 5611–5626. doi: 10.1128/JB.185.18.5611-5626.2003
Sousa, A. L. D., Maués, D., Lobato, A., Franco, E. F., Pinheiro, K., Araújo, F., et al. (2018). PhageWeb – Web interface for rapid identification and characterization of prophages in bacterial genomes. Front. Genet. 9:44. doi: 10.3389/fgene.2018.00644
Steinberg, K. M., and Levin, B. R. (2007). Grazing protozoa and the evolution of the Escherichia coli O157:H7 Shiga toxin-encoding prophage. Proc. Biol. Sci. 274, 1921–1929. doi: 10.1098/rspb.2007.0245
Susskind, M. M., Wright, A., and Botstein, D. (1974). Superinfection exclusion by P22 prophage in lysogens of Salmonella typhimurium. IV. Genetics and physiology of sieB exclusion. Virology 62, 367–384. doi: 10.1016/0042-6822(74)90399-7
Svab, D., Horvath, B., Maroti, G., Dobrindt, U., and Toth, I. (2013). Sequence variability of P2-like prophage genomes carrying the cytolethal distending toxin V operon in Escherichia coli O157. Appl. Environ. Microbiol. 79, 4958–4964. doi: 10.1128/AEM.01134-13
Swalla, B. M., Cho, E. H., Gumport, R. I., and Gardner, J. F. (2003). The molecular basis of co-operative DNA binding between lambda integrase and excisionase. Mol. Microbiol. 50, 89–99. doi: 10.1046/j.1365-2958.2003.03687.x
Tan, Q., Awano, N., and Inouye, M. (2011). YeeV is an Escherichia coli toxin that inhibits cell division by targeting the cytoskeleton proteins, FtsZ and MreB. Mol. Microbiol. 79, 109–118. doi: 10.1111/j.1365-2958.2010.07433.x
Tao, L., Wu, X., and Sun, B. (2010). Alternative sigma factor sigmaH modulates prophage integration and excision in Staphylococcus aureus. PLoS Pathog. 6:e1000888. doi: 10.1371/journal.ppat.1000888
Tautz, D., and Domazet-Lošo, T. (2011). The evolutionary origin of orphan genes. Nat. Rev. Genet. 12, 692–702. doi: 10.1038/nrg3053
Thompson, J. F., Waechter-Brulla, D., Gumport, R. I., Gardner, J. F., Moitoso, De Vargas, L., et al. (1986). Mutations in an integration host factor-binding site: effect on lambda site-specific recombination and regulatory implications. J. Bacteriol. 168, 1343–1351. doi: 10.1128/jb.168.3.1343-1351.1986
Toba, F. A., Thompson, M. G., Campbell, B. R., Junker, L. M., Rueggeberg, K. G., and Hay, A. G. (2011). Role of DLP12 lysis genes in Escherichia coli biofilm formation. Microbiology 157, 1640–1650. doi: 10.1099/mic.0.045161-0
Varga, M., Pantucek, R., Ruzickova, V., and Doskar, J. (2016). Molecular characterization of a new efficiently transducing bacteriophage identified in meticillin-resistant Staphylococcus aureus. J. Gen. Virol. 97, 258–268. doi: 10.1099/jgv.0.000329
Veening, J. W., and Blokesch, M. (2017). Interbacterial predation as a strategy for DNA acquisition in naturally competent bacteria. Nat. Rev. Microbiol. 15, 621–629. doi: 10.1038/nrmicro.2017.66
Ventura, M., Canchaya, C., Kleerebezem, M., De Vos, W. M., Siezen, R. J., and Brussow, H. (2003). The prophage sequences of Lactobacillus plantarum strain WCFS1. Virology 316, 245–255. doi: 10.1016/j.virol.2003.08.019
Veses-Garcia, M., Liu, X., Rigden, D. J., Kenny, J. G., Mccarthy, A. J., and Allison, H. E. (2015). Transcriptomic analysis of Shiga-toxigenic bacteriophage carriage reveals a profound regulatory effect on acid resistance in Escherichia coli. Appl. Environ. Microbiol. 81, 8118–8125. doi: 10.1128/AEM.02034-15
von Wintersdorff, C. J., Penders, J., Van Niekerk, J. M., Mills, N. D., Majumder, S., Van Alphen, L. B., et al. (2016). Dissemination of antimicrobial resistance in microbial ecosystems through horizontal gene transfer. Front. Microbiol. 7:173. doi: 10.3389/fmicb.2016.00173
Waldor, M. K., and Mekalanos, J. J. (1996). Lysogenic conversion by a filamentous phage encoding cholera toxin. Science 272, 1910–1914. doi: 10.1126/science.272.5270.1910
Wang, P. W., Chu, L., and Guttman, D. S. (2004). Complete sequence and evolutionary genomic analysis of the Pseudomonas aeruginosa transposable bacteriophage D3112. J. Bacteriol. 186, 400–410. doi: 10.1128/JB.186.2.400-410.2004
Wang, X., Kim, Y., Ma, Q., Hong, S. H., Pokusaeva, K., Sturino, J. M., et al. (2010). Cryptic prophages help bacteria cope with adverse environments. Nat. Commun. 1:147. doi: 10.1038/ncomms1146
Wang, X., and Wood, T. K. (2016). Cryptic prophages as targets for drug development. Drug Resist. Updates 27, 30–38. doi: 10.1016/j.drup.2016.06.001
Wipf, J. R., Schwendener, S., and Perreten, V. (2014). The novel macrolide-lincosamide-streptogramin B resistance gene erm(44) is associated with a prophage in Staphylococcus xylosus. Antimicrob. Agents Chemother. 58,6133–6138. doi: 10.1128/AAC.02949-14
Yin, S., Bushman, W., and Landy, A. (1985). Interaction of the lambda site-specific recombination protein Xis with attachment site DNA. Proc. Natl. Acad. Sci. U.S.A. 82, 1040–1044. doi: 10.1073/pnas.82.4.1040
Zhang, L., Zhu, B., Dai, R., Zhao, G., and Ding, X. (2013). Control of directionality in Streptomyces phage phiBT1 integrase-mediated site-specific recombination. PLoS One 8:e80434. doi: 10.1371/journal.pone.0080434
Zschöck, M., Botzler, D., Blöcher, S., Sommerhäuser, J., and Hamann, H. P. (2000). Detection of genes for enterotoxins (ent) and toxic shock syndrome toxin-1 (tst) in mammary isolates of Staphylococcus aureus by polymerase-chain-reaction. Int. Dairy J. 10, 569–574. doi: 10.1016/S0958-6946(00)00084-4
Keywords: genome plasticity, genome evolution, horizontal gene transfer, bacterial ecology, bacteriophage
Citation: Ramisetty BCM and Sudhakari PA (2019) Bacterial ‘Grounded’ Prophages: Hotspots for Genetic Renovation and Innovation. Front. Genet. 10:65. doi: 10.3389/fgene.2019.00065
Received: 28 August 2018; Accepted: 24 January 2019;
Published: 12 February 2019.
Edited by:
Arshan Nasir, COMSATS Institute of Information Technology, PakistanReviewed by:
Sylwia Bloch, University of Gdańsk, PolandJianhe Sun, Shanghai Jiao Tong University, China
Alicja Wegrzyn, University of Szczecin, Poland
Copyright © 2019 Ramisetty and Sudhakari. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Bhaskar Chandra Mohan Ramisetty, ramisettybcm@biotech.sastra.edu; bcmr25@gmail.com