Structural and Genetic Analyses Reveal a Key Role in Prophage Excision for the TorI Response Regulator Inhibitor

Comparative genomic analysis of five freshwater cyanophages and reference-guided metagenomic data mining

BACKGROUND

As important producers using photosynthesis on Earth, cyanobacteria contribute to the oxygenation of atmosphere and the primary production of biosphere. However, due to the eutrophication of urban waterbodies and global warming, uncontrollable growth of cyanobacteria usually leads to the seasonal outbreak of cyanobacterial blooms. Cyanophages, a group of viruses that specifically infect and lyse cyanobacteria, are considered as potential environment-friendly agents to control the harmful blooms. Compared to the marine counterparts, only a few freshwater cyanophages have been isolated and genome sequenced to date, largely limiting their characterizations and applications.

RESULTS

Here, we isolated five freshwater cyanophages varying in tail morphology, termed Pam1~Pam5, all of which infect the cyanobacterium Pseudanabaena mucicola Chao 1806 that was isolated from the bloom-suffering Lake Chaohu in Anhui, China. The whole-genome sequencing showed that cyanophages Pam1~Pam5 all contain a dsDNA genome, varying in size from 36 to 142 Kb. Phylogenetic analyses suggested that Pam1~Pam5 possess different DNA packaging mechanisms and are evolutionarily distinct from each other. Notably, Pam1 and Pam5 have lysogeny-associated gene clusters, whereas Pam2 possesses 9 punctuated DNA segments identical to the CRISPR spacers in the host genome. Metagenomic data-based calculation of the relative abundance of Pam1~Pam5 at the Nanfei estuary towards the Lake Chaohu revealed that the short-tailed Pam1 and Pam5 account for the majority of the five cyanophages. Moreover, comparative analyses of the reference genomes of Pam1~Pam5 and previously reported cyanophages enabled us to identify three circular and seven linear contigs of virtual freshwater cyanophages from the metagenomic data of the Lake Chaohu.

CONCLUSIONS

We propose a high-throughput strategy to systematically identify cyanophages based on the currently available metagenomic data and the very limited reference genomes of experimentally isolated cyanophages. This strategy could be applied to mine the complete or partial genomes of unculturable bacteriophages and viruses. Transformation of the synthesized whole genomes of these virtual phages/viruses to proper hosts will enable the rescue of bona fide viral particles and eventually enrich the library of microorganisms that exist on Earth. Video abstract.

Genetic Mining of Newly Isolated Salmophages for Phage Therapy

Salmonella enterica, a Gram-negative zoonotic bacterium, is mainly a food-borne pathogen and the main cause of diarrhea in humans worldwide. The main reservoirs are found in poultry farms, but they are also found in wild birds. The development of antibiotic resistance in S. enterica species raises concerns about the future of efficient therapies against this pathogen and revives the interest in bacteriophages as a useful therapy against bacterial infections. Here, we aimed to decipher and functionally annotate 10 new Salmonella phage genomes isolated in Spain in the light of phage therapy. We designed a bioinformatic pipeline using available building blocks to de novo assemble genomes and perform syntaxic annotation. We then used genome-wide analyses for taxonomic annotation enabled by vContact2 and VICTOR. We were also particularly interested in improving functional annotation using remote homologies detection and comparisons with the recently published phage-specific PHROG protein database. Finally, we searched for useful functions for phage therapy, such as systems encoded by the phage to circumvent cellular defenses with a particular focus on anti-CRISPR proteins. We, thus, were able to genetically characterize nine virulent phages and one temperate phage and identify putative functions relevant to the formulation of phage cocktails for Salmonella biocontrol.

Bioinformatic and experimental characterization of SEN1998: a conserved gene carried by the Enterobacteriaceae-associated ROD21-like family of genomic islands

Genomic islands (GIs) are horizontally transferred elements that shape bacterial genomes and contributes to the adaptation to different environments. Some GIs encode an integrase and a recombination directionality factor (RDF), which are the molecular GI-encoded machinery that promotes the island excision from the chromosome, the first step for the spread of GIs by horizontal transfer. Although less studied, this process can also play a role in the virulence of bacterial pathogens. While the excision of GIs is thought to be similar to that observed in bacteriophages, this mechanism has been only studied in a few families of islands. Here, we aimed to gain a better understanding of the factors involved in the excision of ROD21 a pathogenicity island of the food-borne pathogen Salmonella enterica serovar Enteritidis and the most studied member of the recently described Enterobacteriaceae-associated ROD21-like family of GIs. Using bioinformatic and experimental approaches, we characterized the conserved gene SEN1998, showing that it encodes a protein with the features of an RDF that binds to the regulatory regions involved in the excision of ROD21. While deletion or overexpression of SEN1998 did not alter the expression of the integrase-encoding gene SEN1970, a slight but significant trend was observed in the excision of the island. Surprisingly, we found that the expression of both genes, SEN1998 and SEN1970, were negatively correlated to the excision of ROD21 which showed a growth phase-dependent pattern. Our findings contribute to the growing body of knowledge regarding the excision of GIs, providing insights about ROD21 and the recently described EARL family of genomic islands.

Experimental Evolution of Escherichia coli K-12 at High pH and with RpoS Induction

Experimental evolution of Escherichia coli K-12 W3110 by serial dilutions for 2,200 generations at high pH extended the range of sustained growth from pH 9.0 to pH 9.3. pH 9.3-adapted isolates showed mutations in DNA-binding regulators and envelope proteins. One population showed an IS1 knockout of phoB (encoding the positive regulator of the phosphate regulon). A phoB::kanR knockout increased growth at high pH. phoB mutants are known to increase production of fermentation acids, which could enhance fitness at high pH. Mutations in pcnB [poly(A) polymerase] also increased growth at high pH. Three out of four populations showed deletions of torI, an inhibitor of TorR, which activates expression of torCAD (trimethylamine N-oxide respiration) at high pH. All populations showed point mutations affecting the stationary-phase sigma factor RpoS, either in the coding gene or in genes for regulators of RpoS expression. RpoS is required for survival at extremely high pH. In our microplate assay, rpoS deletion slightly decreased growth at pH 9.1. RpoS protein accumulated faster at pH 9 than at pH 7. The RpoS accumulation at high pH required the presence of one or more antiadaptors that block degradation (IraM, IraD, and IraP). Other genes with mutations after high-pH evolution encode regulators, such as those encoded by yobG (mgrB) (PhoPQ regulator), rpoN (nitrogen starvation sigma factor), malI, and purR, as well as envelope proteins, such as those encoded by ompT and yahO Overall, E. coli evolution at high pH selects for mutations in key transcriptional regulators, including phoB and the stationary-phase sigma factor RpoS.IMPORTANCEEscherichia coli in its native habitat encounters high-pH stress such as that of pancreatic secretions. Experimental evolution over 2,000 generations showed selection for mutations in regulatory factors, such as deletion of the phosphate regulator PhoB and mutations that alter the function of the global stress regulator RpoS. RpoS is induced at high pH via multiple mechanisms.

Regulatory (pan-)genome of an obligate intracellular pathogen in the PVC superphylum

Like other obligate intracellular bacteria, the Chlamydiae feature a compact regulatory genome that remains uncharted owing to poor genetic tractability. Exploiting the reduced number of transcription factors (TFs) encoded in the chlamydial (pan-)genome as a model for TF control supporting the intracellular lifestyle, we determined the conserved landscape of TF specificities by ChIP-Seq (chromatin immunoprecipitation-sequencing) in the chlamydial pathogen Waddlia chondrophila. Among 10 conserved TFs, Euo emerged as a master TF targeting >100 promoters through conserved residues in a DNA excisionase-like winged helix-turn-helix-like (wHTH) fold. Minimal target (Euo) boxes were found in conserved developmentally-regulated genes governing vertical genome transmission (cytokinesis and DNA replication) and genome plasticity (transposases). Our ChIP-Seq analysis with intracellular bacteria not only reveals that global TF regulation is maintained in the reduced regulatory genomes of Chlamydiae, but also predicts that master TFs interpret genomic information in the obligate intracellular α-proteobacteria, including the rickettsiae, from which modern day mitochondria evolved.

S- and N-Oxide Reductases

Escherichia coli is a versatile facultative anaerobe that can respire on a number of terminal electron acceptors, including oxygen, fumarate, nitrate, and S- and N-oxides. Anaerobic respiration using S- and N-oxides is accomplished by enzymatic reduction of these substrates by dimethyl sulfoxide reductase (DmsABC) and trimethylamine N-oxide reductase (TorCA). Both DmsABC and TorCA are membrane-associated redox enzymes that couple the oxidation of menaquinol to the reduction of S- and N-oxides in the periplasm. DmsABC is membrane bound and is composed of a membrane-extrinsic dimer with a 90.4-kDa catalytic subunit (DmsA) and a 23.1-kDa electron transfer subunit (DmsB). These subunits face the periplasm and are held to the membrane by a 30.8-kDa membrane anchor subunit (DmsC). The enzyme provides the scaffold for an electron transfer relay composed of a quinol binding site, five [4Fe-4S] clusters, and a molybdo-bis(molybdopterin guanine dinucleotide) (present nomenclature: Mo-bis-pyranopterin) (Mo-bisMGD) cofactor. TorCA is composed of a soluble periplasmic subunit (TorA, 92.5 kDa) containing a Mo-bis-MGD. TorA is coupled to the quinone pool via a pentaheme c subunit (TorC, 40.4 kDa) in the membrane. Both DmsABC and TorCA require system-specific chaperones (DmsD or TorD) for assembly, cofactor insertion, and/or targeting to the Tat translocon. In this chapter, we discuss the complex regulation of the dmsABC and torCAD operons, the poorly understood paralogues, and what is known about the assembly and translocation to the periplasmic space by the Tat translocon.

Regulation of CTnDOT conjugative transfer is a complex and highly coordinated series of events

CTnDOT is a 65-kb conjugative transposon that is found in Bacteroides spp., which are one of the more abundant members within the lower human gastrointestinal tract. CTnDOT encodes resistance to the antibiotics erythromycin and tetracycline (Tc). An interesting feature of CTnDOT is that exposure to low levels of Tc induces a cascade of events that ultimately results in CTnDOT conjugative transfer. However, Tc is apparently not a switch that activates transfer but rather a signal that appears to override a series of negative regulators that inhibit premature excision and transfer of CTnDOT. In this minireview, we summarize over 20 years of research that focused on elucidating the highly coordinated regulation of excision, mobilization, and transfer of CTnDOT.

IMPORTANCE

Bacteroides spp. are abundant commensals in the human colon, but they are also considered opportunistic pathogens, as they can cause life-threatening infections if they should escape the colon. Bacteroides spp. are the most common cause of anaerobic infections and are rather difficult to treat due to the prevalence of antibiotic resistance within this genus. Today over 80% of Bacteroides are resistant to tetracycline (Tc), and a study looking at both clinical and community isolates demonstrated that this resistance was specifically due to the conjugative transposon CTnDOT.

The structure of Xis reveals the basis for filament formation and insight into DNA bending within a mycobacteriophage intasome

The recombination directionality factor, Xis, is a DNA bending protein that determines the outcome of integrase-mediated site-specific recombination by redesign of higher-order protein-DNA architectures. Although the attachment site DNA of mycobacteriophage Pukovnik is likely to contain four sites for Xis binding, Xis crystals contain five subunits in the asymmetric unit, four of which align into a Xis filament and a fifth that is generated by an unusual domain swap. Extensive intersubunit contacts stabilize a bent filament-like arrangement with Xis monomers aligned head to tail. The structure implies a DNA bend of ~120°, which is in agreement with DNA bending measured in vitro. Formation of attR-containing intasomes requires only Int and Xis, distinguishing Pukovnik from lambda. Therefore, we conclude that, in Pukovnik, Xis-induced DNA bending is sufficient to promote intramolecular Int-mediated bridges during intasome formation.

Transcription termination controls prophage maintenance in Escherichia coli genomes

Prophages represent a large fraction of prokaryotic genomes and often provide new functions to their hosts, in particular virulence and fitness. How prokaryotic cells maintain such gene providers is central for understanding bacterial genome evolution by horizontal transfer. Prophage excision occurs through site-specific recombination mediated by a prophage-encoded integrase. In addition, a recombination directionality factor (or excisionase) directs the reaction toward excision and prevents the phage genome from being reintegrated. In this work, we describe the role of the transcription termination factor Rho in prophage maintenance through control of the synthesis of transcripts that mediate recombination directionality factor expression and, thus, excisive recombination. We show that Rho inhibition by bicyclomycin allows for the expression of prophage genes that lead to excisive recombination. Thus, besides its role in the silencing of horizontally acquired genes, Rho also maintains lysogeny of defective and functional prophages.

Insights into the functions of a prophage recombination directionality factor

Recombination directionality factors (RDFs), or excisionases, are essential players of prophage excisive recombination. Despite the essentially catalytic role of the integrase in both integrative and excisive recombination, RDFs are required to direct the reaction towards excision and to prevent re-integration of the prophage genome when entering a lytic cycle. KplE1, HK620 and numerous (pro)phages that integrate at the same site in enterobacteria genomes (such as the argW tRNA gene) all share a highly conserved recombination module. This module comprises the attL and attR recombination sites and the RDF and integrase genes. The KplE1 RDF was named TorI after its initial identification as a negative regulator of the tor operon. However, it was characterized as an essential factor of excisive recombination. In this study, we designed an extensive random mutagenesis protocol of the torI gene and identified key residues involved in both functions of the TorI protein. We show that, in addition to TorI-TorR protein-protein interaction, TorI interacts in solution with the IntS integrase. Moreover, in vitro, TorR and IntS appear to compete for TorI binding. Finally, our mutagenesis results suggest that the C-terminal part of the TorI protein is dedicated to protein-protein interactions with both proteins TorR and IntS.

DnaJ (Hsp40 protein) binding to folded substrate impacts KplE1 prophage excision efficiency

Temperate phages mediate gene transfer and can modify the properties of their host organisms through the acquisition of novel genes, a process called lysogeny. The KplE1 prophage is one of the 10 prophage regions in Escherichia coli K12 MG1655. KplE1 is defective for lysis but fully competent for site-specific recombination. The TorI recombination directionality factor is strictly required for prophage excision from the host genome. We have previously shown that DnaJ promotes KplE1 excision by increasing the affinity of TorI for its site-specific recombination DNA target. Here, we provide evidence of a direct association between TorI and DnaJ using in vitro cross-linking assays and limited proteolysis experiments that show that this interaction allows both proteins to be transiently protected from trypsin digestion. Interestingly, NMR titration experiments showed that binding of DnaJ involves specific regions of the TorI structure. These regions, mainly composed of α-helices, are located on a surface opposite the DNA-binding site. Taken together, we propose that DnaJ, without the aid of DnaK/GrpE, is capable of increasing the efficiency of KplE1 excision by causing a conformational stabilization that allows TorI to adopt a more favorable conformation for binding to its specific DNA target.

Chaperone-assisted excisive recombination, a solitary role for DnaJ (Hsp40) chaperone in lysogeny escape

Temperate bacteriophage lytic development is intrinsically related to the stress response in particular at the DNA replication and virion maturation steps. Alternatively, temperate phages become lysogenic and integrate their genome into the host chromosome. Under stressful conditions, the prophage resumes a lytic development program, and the phage DNA is excised before being replicated. The KplE1 defective prophage of Escherichia coli K12 constitutes a model system because it is fully competent for integrative as well as excisive recombination and presents an atypical recombination module, which is conserved in various phage genomes. In this work, we identified the host-encoded stress-responsive molecular chaperone DnaJ (Hsp40) as an active participant in KplE1 prophage excision. We first show that the recombination directionality factor TorI of KplE1 specifically interacts with DnaJ. In addition, we found that DnaJ dramatically enhances both TorI binding to its DNA target and excisive recombination in vitro. Remarkably, such stimulatory effect by DnaJ was performed independently of its DnaK chaperone partner and did not require a functional DnaJ J-domain. Taken together, our results underline a novel and unsuspected functional interaction between the generic host stress-regulated chaperone and temperate bacteriophage lysogenic development.

Characterization of the relationship between integrase, excisionase and antirepressor activities associated with a superinfecting Shiga toxin encoding bacteriophage

Shigatoxigenic Escherichia coli emerged as new food borne pathogens in the early 1980s, primarily driven by the dispersal of Shiga toxin-encoding lambdoid bacteriophages. At least some of these Stx phages display superinfection phenotypes, which differ significantly from lambda phage itself, driving through in situ recombination further phage evolution, increasing host range and potentially increasing the host's pathogenic profile. Here, increasing levels of Stx phage Φ24(B) integrase expression in multiple lysogen cultures are demonstrated along with apparently negligible repression of integrase expression by the cognate CI repressor. The Φ24(B) int transcription start site and promoter region were identified and found to differ from in silico predictions. The unidirectional activity of this integrase was determined in an in situ, inducible tri-partite reaction. This indicated that Φ24(B) must encode a novel directionality factor that is controlling excision events during prophage induction. This excisionase was subsequently identified and characterized through complementation experiments. In addition, the previous proposal that a putative antirepressor was responsible for the lack of immunity to superinfection through inactivation of CI has been revisited and a new hypothesis involving the role of this protein in promoting efficient induction of the Φ24(B) prophage is proposed.

Tight regulation of the intS gene of the KplE1 prophage: a new paradigm for integrase gene regulation

Temperate phages have the ability to maintain their genome in their host, a process called lysogeny. For most, passive replication of the phage genome relies on integration into the host's chromosome and becoming a prophage. Prophages remain silent in the absence of stress and replicate passively within their host genome. However, when stressful conditions occur, a prophage excises itself and resumes the viral cycle. Integration and excision of phage genomes are mediated by regulated site-specific recombination catalyzed by tyrosine and serine recombinases. In the KplE1 prophage, site-specific recombination is mediated by the IntS integrase and the TorI recombination directionality factor (RDF). We previously described a sub-family of temperate phages that is characterized by an unusual organization of the recombination module. Consequently, the attL recombination region overlaps with the integrase promoter, and the integrase and RDF genes do not share a common activated promoter upon lytic induction as in the lambda prophage. In this study, we show that the intS gene is tightly regulated by its own product as well as by the TorI RDF protein. In silico analysis revealed that overlap of the attL region with the integrase promoter is widely encountered in prophages present in prokaryotic genomes, suggesting a general occurrence of negatively autoregulated integrase genes. The prediction that these integrase genes are negatively autoregulated was biologically assessed by studying the regulation of several integrase genes from two different Escherichia coli strains. Our results suggest that the majority of tRNA-associated integrase genes in prokaryotic genomes could be autoregulated and that this might be correlated with the recombination efficiency as in KplE1. The consequences of this unprecedented regulation for excessive recombination are discussed.

Protein binding sites involved in the assembly of the KplE1 prophage intasome

The organization of the recombination regions of the KplE1 prophage in Escherichia coli K12 differs from that observed in the lambda prophage. Indeed, the binding sites characterized for the IntS integrase, the TorI recombination directionality factor (RDF) and the integration host factor (IHF) vary in number, spacing and orientation on the attL and attR regions. In this paper, we performed site-directed mutagenesis of the recombination sites to decipher if all sites are essential for the site-specific recombination reaction and how the KplE1 intasome is assembled. We also show that TorI and IntS form oligomers that are stabilized in the presence of their target DNA. Moreover, we found that IHF is the only nucleoid associated protein (NAP) involved in KplE1 recombination, although it is dispensable. This is consistent with the presence of only one functional IHF site on attR and none on attL.

A comparative analysis of the bifunctional Cox proteins of two heteroimmune P2-like phages with different host integration sites

The Cox protein of the coliphage P2 is multifunctional; it acts as a transcriptional repressor of the Pc promoter, as a transcriptional activator of the P(LL) promoter of satellite phage P4, and as a directionality factor for site-specific recombination. The Cox proteins constitute a unique group of directionality factors since they couple the developmental switch with the integration or excision of the phage genome. In this work, the DNA binding characteristics of the Cox protein of WPhi, a P2-related phage, are compared with those of P2 Cox. P2 Cox has been shown to recognize a 9 bp sequence, repeated at least 6 times in different targets. In contrast to P2 Cox, WPhi Cox binds with a strong affinity to the early control region that contains an imperfect direct repeat of 12 nucleotides. The removal of one of the repeats has drastic effects on the capacity of WPhi to bind to the Pe-Pc region. Again in contrast to P2 Cox, WPhi Cox has a lower affinity to attP compared to the Pe-Pc region, and a repeat of 9 bp can be found that has 5 bp in common with the repeat in the Pe-Pc region. WPhi Cox, however, is essential for excisive recombination in vitro. WPhi Cox, like P2 Cox, binds cooperatively with integrase to attP. Both Cox proteins induce a strong bend in their DNA targets upon binding.

Control and regulation of KplE1 prophage site-specific recombination: a new recombination module analyzed

KplE1 is one of the 10 prophage regions of Escherichia coli K12, located at 2464 kb on the chromosome. KplE1 is defective for lysis, but it is fully competent for excisive recombination. In this study, we have mapped the binding sites of the recombination proteins, namely IntS, TorI, and IHF on attL and attR, and the organization of these sites suggests that the intasome is architecturally different from the lambda canonical form. We also measured the relative contribution of these proteins to both excisive and integrative recombination by using a quantitative in vitro assay. These experiments show a requirement of the TorI excisionase for excisive recombination and of the IntS integrase for both integration and excision. Moreover, we observed a strong influence of the supercoiled state of the substrates. The KplE1 recombination module, composed of the integrase and excisionase genes together with the attL and attR DNA regions, is highly similar to that of several phages infecting various E. coli strains as well as Shigella flexneri and Shigella sonnei. The in vitro recombination data reveal that HK620 and KplE1 att sequences are exchangeable. This study thus defines a new site-specific recombination module, and implications for the mechanism and regulation of recombination are discussed.

Evolution of immunity and host chromosome integration site of P2-like coliphages

The amount and distribution of variation in the genomic region containing the genes in the lytic-lysogenic genetic switch and the sequence that determines the integration site into the host chromosome were analyzed for 38 P2-like phages from Escherichia coli. The genetic switch consists of two convergent mutually exclusive promoters, Pe and Pc, and two repressors, C and Cox. The immunity repressor C blocks the early Pe promoter, leading to the establishment of lysogeny. The Cox repressor blocks expression of Pc, allowing lytic growth. Phylogenetic analyses showed that the C and Cox proteins were distributed into seven distinct classes. The phylogenetic relationship differed between the two proteins, and we showed that homologous recombination plays a major role in creating alterations in the genetic switch, leading to new immunity classes. Analyses of the host integration site for these phages resulted in the discovery of a previously unknown site, and there were at least four regular integration sites. Interestingly, we found no case where phages of the same immunity class had different host attachment sites. The evolution of immunity and integration sites is complex, since it involves interactions both between the phages themselves and between phages and hosts, and often, both regulatory proteins and target DNA must change.