Characterization of a new plasmid-like prophage in a pandemic Vibrio parahaemolyticus O3:K6 strain

Genome-based characterization of a novel prophage of Vibrio parahaemolyticus, VPS05ph1, a novel member of Peduoviridae

Vibrio parahaemolyticus is a globally important bacterium related to climate warming and health threat to human and marine animals. Yet, there is limited knowledge about its polylysogeny harboring multiple prophages and the genetic information. In this study, two prophages (VPS05ph1 and VPS05ph2) were identified in a V. parahaemolyticus isolate through genomic and transcriptional analyses. Both prophages were determined as HP1-like phages, located in a novel phylogenetic lineage of Peduoviridae. They shared a moderate genome-wide sequence similarity with each other and high synteny with the closest relatives, but showed low identities to the repressor counterparts of the representative phages within the family. In addition, no bacterial virulence genes, antibiotic resistance genes and known phage-encoded lytic proteins were identified on both prophage genomes. Moreover, the V. parahaemolyticus isolate was induced with mitomycin, which caused aberrant cellular morphology and nonviability of bacterial cells and excision of prophage VPS05ph1, accompanied by the respective inhibition and promotion of transcriptions of the cI-like and cox-like regulator genes for phage decision making. Results in this study provide the genetic context of polylysogeny in the V. parahaemolyticus isolate, support the diversity and prevalence of HP1-like phages in vibrios, and promote to explore interactions between the HP1-like prophage and its vibrio host.

Small RNAs direct attack and defense mechanisms in a quorum sensing phage and its host

Many, if not all, bacteria use quorum sensing (QS) to control collective behaviors, and more recently, QS has also been discovered in bacteriophages (phages). Phages can produce communication molecules of their own, or "listen in" on the host's communication processes, to switch between lytic and lysogenic modes of infection. Here, we study the interaction of Vibrio cholerae with the lysogenic phage VP882, which is activated by the QS molecule DPO. We discover that induction of VP882 results in the binding of phage transcripts to the major RNA chaperone Hfq, which in turn outcompetes and downregulates host-encoded small RNAs (sRNAs). VP882 itself also encodes Hfq-binding sRNAs, and we demonstrate that one of these sRNAs, named VpdS, promotes phage replication by regulating host and phage mRNA levels. We further show that host-encoded sRNAs can antagonize phage replication by downregulating phage mRNA expression and thus might be part of the host's phage defense arsenal.

The Novel Yersinia enterocolitica Telomere Phage vB_YenS_P840 Is Closely Related to PY54, but Reveals Some Striking Differences

Telomere phages are a small group of temperate phages, whose prophages replicate as a linear plasmid with covalently closed ends. They have been isolated from some Enterobacteriaceae and from bacterial species living in aquatic environments. Phage PY54 was the first Yersinia (Y.) enterocolitica telomere phage isolated from a nonpathogenic O:5 strain, but recently a second telomeric Yersinia phage (vB_YenS_P840) was isolated from a tonsil of a wild boar in Germany. Both PY54 and vB_YenS_P840 (P840) have a siphoviridal morphology and a similar genome organization including the primary immunity region immB and telomere resolution site telRL. However, whereas PY54 only possesses one prophage repressor for the lysogenic cycle, vB_YenS_P840 encodes two. The telRL region of this phage was shown to be processed by the PY54 protelomerase under in vivo conditions, but unlike with PY54, a flanking inverted repeat was not required for processing. A further substantial difference between the phages is their host specificity. While PY54 infects Y. enterocolitica strains belonging to the serotypes O:5 and O:5,27, vB_YenS_P840 exclusively lyses O:3 strains. As the tail fiber and tail fiber assembly proteins of the phages differ significantly, we introduced the corresponding genes of vB_YenS_P840 by transposon mutagenesis into the PY54 genome and isolated several mutants that were able to infect both serotypes, O:5,27 and O:3.

Natural silencing of quorum-sensing activity protects Vibrio parahaemolyticus from lysis by an autoinducer-detecting phage

Quorum sensing (QS) is a chemical communication process that bacteria use to track population density and orchestrate collective behaviors. QS relies on the production, accumulation, and group-wide detection of extracellular signal molecules called autoinducers. Vibriophage 882 (phage VP882), a bacterial virus, encodes a homolog of the Vibrio QS receptor-transcription factor, called VqmA, that monitors the Vibrio QS autoinducer DPO. Phage VqmA binds DPO at high host-cell density and activates transcription of the phage gene qtip. Qtip, an antirepressor, launches the phage lysis program. Phage-encoded VqmA when bound to DPO also manipulates host QS by activating transcription of the host gene vqmR. VqmR is a small RNA that controls downstream QS target genes. Here, we sequence Vibrio parahaemolyticus strain O3:K6 882, the strain from which phage VP882 was initially isolated. The chromosomal region normally encoding vqmR and vqmA harbors a deletion encompassing vqmR and a portion of the vqmA promoter, inactivating that QS system. We discover that V. parahaemolyticus strain O3:K6 882 is also defective in its other QS systems, due to a mutation in luxO, encoding the central QS transcriptional regulator LuxO. Both the vqmR-vqmA and luxO mutations lock V. parahaemolyticus strain O3:K6 882 into the low-cell density QS state. Reparation of the QS defects in V. parahaemolyticus strain O3:K6 882 promotes activation of phage VP882 lytic gene expression and LuxO is primarily responsible for this effect. Phage VP882-infected QS-competent V. parahaemolyticus strain O3:K6 882 cells lyse more rapidly and produce more viral particles than the QS-deficient parent strain. We propose that, in V. parahaemolyticus strain O3:K6 882, constitutive maintenance of the low-cell density QS state suppresses the launch of the phage VP882 lytic cascade, thereby protecting the bacterial host from phage-mediated lysis.

Natural Silencing of Quorum-Sensing Activity Protects Vibrio parahaemolyticus from Lysis by an Autoinducer-Detecting Phage

Quorum sensing (QS) is a chemical communication process that bacteria use to track population density and orchestrate collective behaviors. QS relies on the production, accumulation, and group-wide detection of extracellular signal molecules called autoinducers. Vibriophage 882 (phage VP882), a bacterial virus, encodes a homolog of the Vibrio QS receptor-transcription factor, called VqmA, that monitors the Vibrio QS autoinducer DPO. Phage VqmA binds DPO at high host-cell density and activates transcription of the phage gene qtip . Qtip, an antirepressor, launches the phage lysis program. Phage-encoded VqmA when bound to DPO also manipulates host QS by activating transcription of the host gene vqmR . VqmR is a small RNA that controls downstream QS target genes. Here, we sequence Vibrio parahaemolyticus strain O3:K6 882, the strain from which phage VP882 was initially isolated. The chromosomal region normally encoding vqmR and vqmA harbors a deletion encompassing vqmR and a portion of the vqmA promoter, inactivating that QS system. We discover that V. parahaemolyticus strain O3:K6 882 is also defective in its other QS systems, due to a mutation in luxO , encoding the central QS transcriptional regulator LuxO. Both the vqmR-vqmA and luxO mutations lock V. parahaemolyticus strain O3:K6 882 into the low-cell density QS state. Reparation of the QS defects in V. parahaemolyticus strain O3:K6 882 promotes activation of phage VP882 lytic gene expression and LuxO is primarily responsible for this effect. Phage VP882-infected QS-competent V. parahaemolyticus strain O3:K6 882 cells lyse more rapidly and produce more viral particles than the QS-deficient parent strain. We propose that, in V. parahaemolyticus strain O3:K6 882, constitutive maintenance of the low-cell density QS state suppresses the launch of the phage VP882 lytic cascade, thereby protecting the bacterial host from phage-mediated lysis.

Homologs of the Escherichia coli F Element Protein TraR, Including Phage Lambda Orf73, Directly Reprogram Host Transcription

Bacterial cells and their associated plasmids and bacteriophages encode numerous small proteins of unknown function. One example, the 73-amino-acid protein TraR, is encoded by the transfer operon of the conjugative F plasmid of Escherichia coli. TraR is a distant homolog of DksA, a protein found in almost all proteobacterial species that is required for ppGpp to regulate transcription during the stringent response. TraR and DksA increase or decrease transcription initiation depending on the kinetic features of the promoter by binding directly to RNA polymerase without binding to DNA. Unlike DksA, whose full activity requires ppGpp as a cofactor, TraR is fully active by itself and unaffected by ppGpp. TraR belongs to a family of divergent proteins encoded by proteobacterial bacteriophages and other mobile elements. Here, we experimentally addressed whether other members of the TraR family function like the F element-encoded TraR. Purified TraR and all 5 homologs that were examined bound to RNA polymerase, functioned at lower concentrations than DksA, and complemented a dksA-null strain for growth on minimal medium. One of the homologs, λ Orf73, encoded by bacteriophage lambda, was examined in greater detail. λ Orf73 slowed host growth and increased phage burst size. Mutational analysis suggested that λ Orf73 and TraR have a similar mechanism for inhibiting rRNA and r-protein promoters. We suggest that TraR and its homologs regulate host transcription to divert cellular resources to phage propagation or conjugation without induction of ppGpp and a stringent response.

Pyocin-mediated antagonistic interactions in Pseudomonas spp. isolated in James Ross Island, Antarctica

Interactions within bacterial communities are frequently mediated by the production of antimicrobial agents. Despite the increasing interest in research of new antimicrobials, studies describing antagonistic interactions among cold-adapted microorganisms are still rare. Our study assessed the antimicrobial interactions of 36 Antarctic Pseudomonas spp. and described the genetic background of these interactions in selected strains. The overall bacteriocinogeny was greater compared to mesophilic Pseudomonas non-aeruginosa species. R-type tailocins were detected on transmission electron micrographs in 16 strains (44.4%); phylogenetic analysis of the corresponding gene clusters revealed that the P. prosekii CCM 8878 tailocin was related to the Rp3 group, whereas the tailocin in Pseudomonas sp. CCM 8880 to the Rp4 group. Soluble antimicrobials were produced by eight strains (22.-2%); gene mining found pyocin L homologues in the genomes of P. prosekii CCM 8881 and CCM 8879 and pyocin S9-like homologues in P. prosekii CCM 8881 and Pseudomonas sp. CCM 8880. Analysis of secretomes confirmed the production of all S- and L-type pyocin genes. Our results suggest that bacteriocin-based inhibition plays an important role in interactions among Antarctic soil bacteria, and these native, cold-adapted microorganisms could be a promising source of new antimicrobials.

Adaptations of Vibrio parahaemolyticus to Stress During Environmental Survival, Host Colonization, and Infection

Vibrio parahaemolyticus (Vp) is an aquatic Gram-negative bacterium that may infect humans and cause gastroenteritis and wound infections. The first pandemic of Vp associated infection was caused by the serovar O3:K6 and epidemics caused by the other serovars are increasingly reported. The two major virulence factors, thermostable direct hemolysin (TDH) and/or TDH-related hemolysin (TRH), are associated with hemolysis and cytotoxicity. Vp strains lacking tdh and/or trh are avirulent and able to colonize in the human gut and cause infection using other unknown factors. This pathogen is well adapted to survive in the environment and human host using several genetic mechanisms. The presence of prophages in Vp contributes to the emergence of pathogenic strains from the marine environment. Vp has two putative type-III and type-VI secretion systems (T3SS and T6SS, respectively) located on both the chromosomes. T3SS play a crucial role during the infection process by causing cytotoxicity and enterotoxicity. T6SS contribute to adhesion, virulence associated with interbacterial competition in the gut milieu. Due to differential expression, type III secretion system 2 (encoded on chromosome-2, T3SS2) and other genes are activated and transcribed by interaction with bile salts within the host. Chromosome-1 encoded T6SS1 has been predominantly identified in clinical isolates. Acquisition of genomic islands by horizontal gene transfer provides enhanced tolerance of Vp toward several antibiotics and heavy metals. Vp consists of evolutionarily conserved targets of GTPases and kinases. Expression of these genes is responsible for the survival of Vp in the host and biochemical changes during its survival. Advanced genomic analysis has revealed that various genes are encoded in Vp pathogenicity island that control and expression of virulence in the host. In the environment, the biofilm gene expression has been positively correlated to tolerance toward aerobic, anaerobic, and micro-aerobic conditions. The genetic similarity analysis of toxin/antitoxin systems of Escherichia coli with VP genome has shown a function that could induce a viable non-culturable state by preventing cell division. A better interpretation of the Vp virulence and other mechanisms that support its environmental fitness are important for diagnosis, treatment, prevention and spread of infections. This review identifies some of the common regulatory pathways of Vp in response to different stresses that influence its survival, gut colonization and virulence.

Pathogenic Vibrio Species Are Associated with Distinct Environmental Niches and Planktonic Taxa in Southern California (USA) Aquatic Microbiomes

Interactions between vibrio bacteria and the planktonic community impact marine ecology and human health. Many coastal Vibrio spp. can infect humans, representing a growing threat linked to increasing seawater temperatures. Interactions with eukaryotic organisms may provide attachment substrate and critical nutrients that facilitate the persistence, diversification, and spread of pathogenic Vibrio spp. However, vibrio interactions with planktonic organisms in an environmental context are poorly understood. We quantified the pathogenic Vibrio species V. cholerae, V. parahaemolyticus, and V. vulnificus monthly for 1 year at five sites and observed high abundances, particularly during summer months, with species-specific temperature and salinity distributions. Using metabarcoding, we established a detailed profile of both prokaryotic and eukaryotic coastal microbial communities. We found that pathogenic Vibrio species were frequently associated with distinct eukaryotic amplicon sequence variants (ASVs), including diatoms and copepods. Shared environmental conditions, such as high temperatures and low salinities, were associated with both high concentrations of pathogenic vibrios and potential environmental reservoirs, which may influence vibrio infection risks linked to climate change and should be incorporated into predictive ecological models and experimental laboratory systems. IMPORTANCE Many species of coastal vibrio bacteria can infect humans, representing a growing health threat linked to increasing seawater temperatures. However, their interactions with surrounding microbes in the environment, especially eukaryotic organisms that may provide nutrients and attachment substrate, are poorly understood. We quantified three pathogenic Vibrio species monthly for a duration of 1 year, finding that all three species were abundant and exhibited species-specific temperature and salinity distributions. Using metabarcoding, we investigated associations between these pathogenic species and prokaryotic and eukaryotic microbes, revealing genus and amplicon sequence variant (ASV)-specific relationships with potential functional implications. For example, pathogenic species were frequently associated with chitin-producing eukaryotes, such as diatoms in the genus Thalassiosira and copepods. These associations between high concentrations of pathogenic vibrios and potential environmental reservoirs should be considered when predicting infection risk and developing ecologically relevant model systems.

Bacteria have numerous distinctive groups of phage-plasmids with conserved phage and variable plasmid gene repertoires

Plasmids and temperate phages are key contributors to bacterial evolution. They are usually regarded as very distinct. However, some elements, termed phage–plasmids, are known to be both plasmids and phages, e.g. P1, N15 or SSU5. The number, distribution, relatedness and characteristics of these phage–plasmids are poorly known. Here, we screened for these elements among ca. 2500 phages and 12000 plasmids and identified 780 phage–plasmids across very diverse bacterial phyla. We grouped 92% of them by similarity of gene repertoires to eight defined groups and 18 other broader communities of elements. The existence of these large groups suggests that phage–plasmids are ancient. Their gene repertoires are large, the average element is larger than an average phage or plasmid, and they include slightly more homologs to phages than to plasmids. We analyzed the pangenomes and the genetic organization of each group of phage–plasmids and found the key phage genes to be conserved and co-localized within distinct groups, whereas genes with homologs in plasmids are much more variable and include most accessory genes. Phage–plasmids are a sizeable fraction of the sequenced plasmids (∼7%) and phages (∼5%), and could have key roles in bridging the genetic divide between phages and other mobile genetic elements.

Induction and Genome Analysis of HY01, a Newly Reported Prophage from an Emerging Shrimp Pathogen Vibrio campbellii

Vibrio campbellii is an emerging aquaculture pathogen that causes luminous vibriosis in farmed shrimp. Although prophages in various aquaculture pathogens have been widely reported, there is still limited knowledge regarding prophages in the genome of pathogenic V. campbellii. Here, we describe the full-genome sequence of a prophage named HY01, induced from the emerging shrimp pathogen V. campbellii HY01. The phage HY01 was induced by mitomycin C and was morphologically characterized as long tailed phage. V. campbellii phage HY01 is composed of 41,772 bp of dsDNA with a G+C content of 47.45%. A total of 60 open reading frames (ORFs) were identified, of which 31 could be predicted for their biological functions. Twenty seven out of 31 predicted protein coding regions were matched with several encoded proteins of various Enterobacteriaceae, Pseudomonadaceae, Vibrionaceae, and other phages of Gram-negative bacteria. Interestingly, the comparative genome analysis revealed that the phage HY01 was only distantly related to Vibrio phage Va_PF430-3_p42 of fish pathogen V. anguillarum but differed in genomic size and gene organization. The phylogenetic tree placed the phage together with Siphoviridae family. Additionally, a survey of Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) spacers revealed two matching sequences between phage HY01 genome and viral spacer sequence of Vibrio spp. The spacer results combined with the synteny results suggest that the evolution of V. campbellii phage HY01 is driven by the horizontal genetic exchange between bacterial families belonging to the class of Gammaproteobacteria.

Prophage-Related Gene VpaChn25_0724 Contributes to Cell Membrane Integrity and Growth of Vibrio parahaemolyticus CHN25

Vibrio parahaemolyticus is a leading seafood-borne pathogen that can cause acute gastroenteritis and even death in humans. In aquatic ecosystems, phages constantly transform bacterial communities by horizontal gene transfer. Nevertheless, biological functions of prophage-related genes in V. parahaemolyticus remain to be fully unveiled. Herein, for the first time, we studied one such gene VpaChn25_0724 encoding an unknown hypothetical protein in V. parahaemolyticus CHN25. This gene deletion mutant ΔVpaChn25_0724 was constructed by homologous recombination, and its complementary mutant ΔVpaChn25_0724-com was also obtained. The ΔVpaChn25_0724 mutant exhibited a sever defect in growth and swimming motility particularly at lower temperatures. Biofilm formation and cytotoxicity capacity of V. parahaemolyticus CHN25 was significantly lowered in the absence of VpaChn25_0724. Comparative secretomic analysis revealed an increase in extracellular proteins of ΔVpaChn25_0724, which likely resulted from its damaged cell membrane. Comparison of transcriptome data showed twelve significantly altered metabolic pathways in ΔVpaChn25_0724, suggesting inactive transport and utilization of carbon sources, repressed energy production and membrane biogenesis in ΔVpaChn25_0724. Comparative transcriptomic analysis also revealed several remarkably down-regulated key regulators in bacterial gene regulatory networks linked to the observed phenotypic variations. Overall, the results here facilitate better understanding of biological significance of prophage-related genes remaining in V. parahaemolyticus.

Closely Related Vibrio alginolyticus Strains Encode an Identical Repertoire of Caudovirales-Like Regions and Filamentous Phages

Many filamentous vibriophages encode virulence genes that lead to the emergence of pathogenic bacteria. Most genomes of filamentous vibriophages characterized up until today were isolated from human pathogens. Despite genome-based predictions that environmental Vibrios also contain filamentous phages that contribute to bacterial virulence, empirical evidence is scarce. This study aimed to characterize the bacteriophages of a marine pathogen, Vibrio alginolyticus (Kiel-alginolyticus ecotype) and to determine their role in bacterial virulence. To do so, we sequenced the phage-containing supernatant of eight different V. alginolyticus strains, characterized the phages therein and performed infection experiments on juvenile pipefish to assess their contribution to bacterial virulence. We were able to identify two actively replicating filamentous phages. Unique to this study was that all eight bacteria of the Kiel-alginolyticus ecotype have identical bacteriophages, supporting our previously established theory of a clonal expansion of the Kiel-alginolyticus ecotype. We further found that in one of the two filamentous phages, two phage-morphogenesis proteins (Zot and Ace) share high sequence similarity with putative toxins encoded on the Vibrio cholerae phage CTXΦ. The coverage of this filamentous phage correlated positively with virulence (measured in controlled infection experiments on the eukaryotic host), suggesting that this phage contributes to bacterial virulence.

High diversity of Vibrio spp. associated with different ecological niches in a marine aquaria system and description of Vibrio aquimaris sp. nov

The aim of the study was to characterise the diversity and niche-specific colonization of Vibrio spp. in a marine aquaria system by a cultivation-dependent approach. A total of 53 Vibrio spp. isolates were cultured from different ecological niches in a marine aquarium including microplastic (MP) and sandy sediment particles (12 weeks after added sterile to the system), detritus, and the surrounding aquarium water. Based on the 16S rRNA gene sequence phylogeny and multilocus sequence analysis (MLSA) the isolates were assigned to seven different phylotypes. Six phylotypes were identified by high probability to the species level. The highest phylotype diversity was cultured from detritus and water (six out of seven phylotypes), while only two phylotypes were cultured from MP and sediment particles. Genomic fingerprinting indicated an even higher genetic diversity of Vibrio spp. at the strain (genotype) level. Again, the highest diversity of genotypes was recovered from detritus and water while only few partially particle-type specific genotypes were cultured from MP and sediment particles. Phylotype V-2 formed an independent branch in the MLSA tree and could not be assigned to a described Vibrio species. Isolates of this phylotype showed highest 16S rRNA gene sequence similarity to type strains of Vibrio japonicus (98.5%) and Vibrio caribbeanicus (98.4%). A representative isolate, strain THAF100^T, was characterised by a polyphasic taxonomic approach and Vibrio aquimaris sp. nov., with strain THAF100^T (=DSM 109633^T=LMG 31434^T=CIP 111709^T) as type strain, is proposed as novel species.

Looking for the hidden: Characterization of lysogenic phages in potential pathogenic Vibrio species from the North Sea

The incidence of potentially pathogenic Vibrio species in the marine environment around Europe, is correlated with the increase of surface seawater temperature. Despite their importance, little is known about the trigger factors of potential outbreak-causing strains in this region. As prophages may compose a major reservoir of virulence traits in marine ecosystems, this study aims to identify and characterize the genomes of lysogenic Vibrio phages exemplarily from the North Sea. Therefore, 31 isolates from potentially pathogenic Vibrio species from the North Sea were screened for inducible prophages with mitomycin C. From them, one V. cholerae isolate and 40% V. parahaemolyticus isolates carried inducible prophages. Three lysogenic phages were selected for genomic characterization. The phage vB_VpaM_VP-3212 (unclassified Myoviridae) has a genome with a length of 36.81 Kbp and 55 CDS were identified. This lysogenic phage of V. parahaemolyticus contains genes related to replicative transposition mechanism, such as transposase and mobile elements similar to Mu-like viruses. The phage vB_VpaP_VP-3220 (Podoviridae, unclassified Nona33virus) has a genome length of 58,14 Kbp and contains 63 CDS. This V. parahaemolyticus phage probably uses a headful (pac) packaging replication mechanism. The phage vB_VchM_VP-3213 (unclassified Myoviridae) has a genome with a length of 41 Kbp and 63 CDS were identified, including integrase and Xer system for lysogenic recombination. This lysogenic phage of V. cholerae has similar genomic features as lambdoid phages. Although no pathogenicity genes were identified, their similarity among other phage genomes indicates that these phages can affect the development of pathogenic Vibrio strains in marine environments.

Large Phenotypic and Genetic Diversity of Prophages Induced from the Fish Pathogen Vibrio anguillarum

Vibrio anguillarum is a marine pathogenic bacterium that causes vibriosis in fish and shellfish. Although prophage-like sequences have been predicted in V. anguillarum strains, many are not characterized, and it is not known if they retain the functional capacity to form infectious particles that can infect and lysogenize other bacterial hosts. In this study, the genome sequences of 28 V. anguillarum strains revealed 55 different prophage-related elements. Chemical and spontaneous induction allowed a collection of 42 phage isolates, which were classified in seven different groups according to a multiplex PCR assay. One shared prophage sequence, p41 (group III), was present in 17 V. anguillarum strains, suggesting that this specific element is very dynamically exchanged among V. anguillarum populations. Interestingly, the host range of genetically identical phages was highly dependent on the strains used for proliferation, indicating that phenotypic properties of phages were partly regulated by the host. Finally, experimental evidence displayed that the induced phage ɸVa_90-11-287_p41 was able to lysogenize V. anguillarum strain Ba35, and subsequently spontaneously become released from the lysogenized cells, demonstrating an efficient transfer of the phage among V. anguillarum strains. Altogether, the results showed large genetic and functional diversity and broad distribution of prophages in V. anguillarum, and demonstrated the potential of prophages as drivers of evolution in V. anguillarum strains.

A Host-Produced Quorum-Sensing Autoinducer Controls a Phage Lysis-Lysogeny Decision

Vibrio cholerae uses a quorum-sensing (QS) system composed of the autoinducer 3,5-dimethylpyrazin-2-ol (DPO) and receptor VqmA (VqmA_Vc), which together repress genes for virulence and biofilm formation. vqmA genes exist in Vibrio and in one vibriophage, VP882. Phage-encoded VqmA (VqmA_Phage) binds to host-produced DPO, launching the phage lysis program via an antirepressor that inactivates the phage repressor by sequestration. The antirepressor interferes with repressors from related phages. Like phage VP882, these phages encode DNA-binding proteins and partner antirepressors, suggesting that they, too, integrate host-derived information into their lysis-lysogeny decisions. VqmA_Phage activates the host VqmA_Vc regulon, whereas VqmA_Vc cannot induce phage-mediated lysis, suggesting an asymmetry whereby the phage influences host QS while enacting its own lytic-lysogeny program without interference. We reprogram phages to activate lysis in response to user-defined cues. Our work shows that a phage, causing bacterial infections, and V. cholerae, causing human infections, rely on the same signal molecule for pathogenesis.

Comparative genomics for non-O1/O139 Vibrio cholerae isolates recovered from the Yangtze River Estuary versus V. cholerae representative isolates from serogroup O1

Vibriocholerae, which is autochthonous to estuaries worldwide, can cause human cholera that is still pandemic in developing countries. A number of V. cholerae isolates of clinical and environmental origin worldwide have been subjected to genome sequencing to address their phylogenesis and bacterial pathogenesis, however, little genome information is available for V. cholerae isolates derived from estuaries, particularly in China. In this study, we determined the complete genome sequence of V. cholerae CHN108B (non-O1/O139 serogroup) isolated from the Yangtze River Estuary, China and performed comparative genome analysis between CHN108B and other eight representative V. cholerae isolates. The 4,168,545-bp V. cholerae CHN108B genome (47.2% G+C) consists of two circular chromosomes with 3,691 predicted protein-encoding genes. It has 110 strain-specific genes, the highest number among the eight representative V. cholerae whole genomes from serogroup O1: there are seven clinical isolates linked to cholera pandemics (1937-2010) and one environmental isolate from Brazil. Various mobile genetic elements (such as insertion sequences, prophages, integrative and conjugative elements, and super-integrons) were identified in the nine V. cholerae genomes of clinical and environmental origin, indicating that the bacterium undergoes extensive genetic recombination via lateral gene transfer. Comparative genomics also revealed different virulence and antimicrobial resistance gene patterns among the V. cholerae isolates, suggesting some potential virulence factors and the rising development of resistance among pathogenic V. cholerae. Additionally, draft genome sequences of multiple V. cholerae isolates recovered from the Yangtze River Estuary were also determined, and comparative genomics revealed many genes involved in specific metabolism pathways, which are likely shaped by the unique estuary environment. These results provide additional evidence of V. cholerae genome plasticity and will facilitate better understanding of the genome evolution and pathogenesis of this severe water-borne pathogen worldwide.

Characterization of a lytic vibriophage VP06 of Vibrio parahaemolyticus

Vibrio parahaemolyticus is a human enteropathogenic bacterium and is also pathogenic to shrimp and finfish. In a search for a biocontrol agent for V. parahaemolyticus and other pathogenic Vibrio species, a lytic phage VP06 was isolated from oyster using V. parahaemolyticus as the host. VP06 is a Siphoviridae phage with a polyhedral head and a long tail. The genome sequence of VP06 was 75,893 nucleotides in length and the G + C content was 49%; a total of 101 CDSs were identified in VP06, of which 39 exhibited functional domains/motifs. The genomic sequence of VP06 is similar to those of a lytic Vibrio vulnificus phage SSP002 and a temperate V. parahaemolyticus phage vB_VpaS_MAR10, although VP06 has distinct features in the CDS arrangement and 14 unique CDSs. Phylogenetic analysis revealed that VP06, SSP002 and vB_VpaS_MAR10 belong to a novel genus cluster of Siphoviridae phages. This phage lysed 28.1% of various Vibrio strains, and the efficiency of plating method revealed that VP06 was highly effective in lysing strains of Vibrio alginolyticus, Vibrio azureus, Vibrio harveyi and V. parahaemolyticus. The properties of VP06, including its broad range of hosts and resistance to environmental stresses, indicate that it may be a candidate biocontrol agent.

Tripartite species interaction: eukaryotic hosts suffer more from phage susceptible than from phage resistant bacteria

Background

Evolutionary shifts in bacterial virulence are often associated with a third biological player, for instance temperate phages, that can act as hyperparasites. By integrating as prophages into the bacterial genome they can contribute accessory genes, which can enhance the fitness of their prokaryotic carrier (lysogenic conversion). Hyperparasitic influence in tripartite biotic interactions has so far been largely neglected in empirical host-parasite studies due to their inherent complexity. Here we experimentally address whether bacterial resistance to phages and bacterial harm to eukaryotic hosts is linked using a natural tri-partite system with bacteria of the genus Vibrio, temperate vibriophages and the pipefish Syngnathus typhle. We induced prophages from all bacterial isolates and constructed a three-fold replicated, fully reciprocal 75 × 75 phage-bacteria infection matrix.

Results

According to their resistance to phages, bacteria could be grouped into three distinct categories: highly susceptible (HS-bacteria), intermediate susceptible (IS-bacteria), and resistant (R-bacteria). We experimentally challenged pipefish with three selected bacterial isolates from each of the three categories and determined the amount of viable Vibrio counts from infected pipefish and the expression of pipefish immune genes. While the amount of viable Vibrio counts did not differ between bacterial groups, we observed a significant difference in relative gene expression between pipefish infected with phage susceptible and phage resistant bacteria.

Conclusion

These findings suggest that bacteria with a phage-susceptible phenotype are more harmful against a eukaryotic host, and support the importance of hyperparasitism and the need for an integrative view across more than two levels when studying host-parasite evolution.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0930-2) contains supplementary material, which is available to authorized users.

Binding Specificities of the Telomere Phage ϕKO2 Prophage Repressor CB and Lytic Repressor Cro

Temperate bacteriophages possess a genetic switch which regulates the lytic and lysogenic cycle. The genomes of the temperate telomere phages N15, PY54, and ϕKO2 harbor a primary immunity region (immB) comprising genes for the prophage repressor (cI or cB), the lytic repressor (cro) and a putative antiterminator (q). The roles of these products are thought to be similar to those of the lambda proteins CI (CI prophage repressor), Cro (Cro repressor), and Q (antiterminator Q), respectively. Moreover, the gene order and the location of several operator sites in the prototype telomere phage N15 and in ϕKO2 are reminiscent of lambda-like phages. We determined binding sites of the ϕKO2 prophage repressor CB and lytic repressor Cro on the ϕKO2 genome in detail by electrophoretic mobility shift assay (EMSA) studies. Unexpectedly, ϕKO2 CB and Cro revealed different binding specificities. CB was bound to three O_R operators in the intergenic region between cB and cro, two O_L operators between cB and the replication gene repA and even to operators of N15. Cro bound exclusively to the 16 bp operator site O_R3 upstream of the ϕKO2 prophage repressor gene. The ϕKO2 genes cB and cro are regulated by several strong promoters overlapping with the O_R operators. The data suggest that Cro represses cB transcription but not its own synthesis, as already reported for PY54 Cro. Thus, not only PY54, but also phage ϕKO2 possesses a genetic switch that diverges significantly from the switch of lambda-like phages.

Two novel temperate bacteriophages co-existing in Aeromonas sp. ARM81 - characterization of their genomes, proteomes and DNA methyltransferases

Aeromonas species are causative agents of a wide spectrum of diseases in animals and humans. Although these bacteria are commonly found in various environments, little is known about their phages. Thus far, only one temperate Aeromonas phage has been characterized. Whole-genome sequencing of an Aeromonas sp. strain ARM81 revealed the presence of two prophage clusters. One of them is integrated into the chromosome and the other was maintained as an extrachromosomal, linear plasmid-like prophage encoding a protelomerase. Both prophages were artificially and spontaneously inducible. We separately isolated both phages and compared their genomes with other known viruses. The novel phages show no similarity to the previously characterized Aeromonas phages and might represent new evolutionary lineages of viruses infecting Aeromonadaceae. Apart from the comparative genomic analyses of these phages, complemented with their structural and molecular characterization, a functional analysis of four DNA methyltransferases encoded by these viruses was conducted. One of the investigated N6-adenine-modifying enzymes shares sequence specificity with a Dam-like methyltransferase of its bacterial host, while another one is non-specific, as it catalyzes adenine methylation in various sequence contexts. The presented results shed new light on the diversity of Aeromonas temperate phages.

Replication and Maintenance of Linear Phage-Plasmid N15

The lambdoid phage N15 of Escherichia coli is very unusual among temperate phages in that its prophage is not integrated into the chromosome but is a linear plasmid molecule with covalently closed ends (telomeres). Upon infection, the phage DNA circularizes via cohesive ends, and then a special phage enzyme of the tyrosine recombinase family, protelomerase, cuts at another site and joins the ends, forming hairpin telomeres of the linear plasmid prophage. Replication of the N15 prophage is initiated at an internally located ori site and proceeds bidirectionally, resulting in the formation of duplicated telomeres. The N15 protelomerase cuts them, generating two linear plasmid molecules with hairpin telomeres. Stable inheritance of the plasmid prophage is ensured by a partitioning operon similar to the F factor sop operon. Unlike the F centromere, the N15 centromere consists of four inverted repeats dispersed in the genome. The multiplicity and dispersion of centromeres are required for efficient partitioning of a linear plasmid. The centromeres are located in the N15 genome regions involved in phage replication and control of lytic development, and binding of partition proteins at these sites regulates these processes. The family of N15-like linear phage-plasmids includes lambdoid phages ɸKO2 and pY54, as well as Myoviridae phages ΦHAP-1, VHML, VP882, Vp58.5, and vB_VpaM_MAR of marine gamma-proteobacteria. The genomes of these phages contain similar protelomerase genes, lysogeny control modules, and replication genes, suggesting that these phages may belong to a group diverged from a common ancestor.

Genome characterization of a novel vibriophage VpKK5 (Siphoviridae) specific to fish pathogenic strain of Vibrio parahaemolyticus

Vibrio parahaemolyticus has long been known pathogenic to shrimp but only recently it is also reported pathogenic to tropical cultured marine finfish. Traditionally, bacterial diseases in aquaculture are often treated using synthetic antibiotics but concern due to side effects of these chemicals is elevating hence, new control strategies which are both environmental and consumer friendly, are urgently needed. One promising control strategy is the bacteriophage therapy. In this study, we report the isolation and characterization of a novel vibriophage (VpKK5), belonging to the family Siphoviridae that was specific and capable of complete lysing the fish pathogenic strain of V. parahaemolyticus. The VpKK5 exhibited short eclipse and latent periods of 24 and 36 min, respectively, but with a large burst size of 180 pfu/cell. The genome analysis revealed that the VpKK5 is a novel bacteriophage with the estimated genome size of 56,637 bp and has 53.1% G + C content. The vibriophage has about 80 predicted open reading frames consisted of 37 complete coding sequences which did not match to any protein databases. The analysis also found no lysogeny and virulence genes in the genome of VpKK5. With such genome features, we suspected the vibriophage is novel and could be explored for phage therapy against fish pathogenic strains of V. parahaemolyticus in the near future.

The diverse genetic switch of enterobacterial and marine telomere phages

Temperate bacteriophages possess a genetic switch which regulates the lytic and lysogenic cycle. The genomes of the enterobacterial telomere phages N15, PY54 and ϕKO2 harbor a primary immunity region (immB) comprising genes for the prophage repressor, the lytic repressor and a putative antiterminator, similar to CI, Cro and Q of lambda, respectively. Moreover, N15 and ϕKO2 contain 3 related operator (OR) sites between cI and cro, while only one site (OR3) has been detected in PY54. Marine telomere phages possess a putative cI gene but not a cro-like gene. Instead, a gene is located at the position of cro, whose product shows some similarity to the PY54 ORF42 product, the function of which is unknown. We have determined the transcription start sites of the predicted repressor genes of N15, PY54, ϕKO2 and of the marine telomere phage VP58.5. The influence of the genes on phage propagation was analyzed in E. coli, Y. enterocolitica and V.parahaemolyticus. We show that the repressors and antiterminators of N15, ϕKO2 and PY54 exerted their predicted activities. However, while the proteins of both N15 and ϕKO2 affected lysis and lysogeny by N15, they did not affect PY54 propagation. On the other hand, the respective PY54 proteins exclusively influenced the propagation of this phage. The immB region of VP58.5 contains 2 genes that revealed prophage repressor activity, while a lytic repressor gene could not be identified. The results indicate an unexpected diversity of the growth regulation mechanisms in these temperate phages.

From cholera to corals: Viruses as drivers of virulence in a major coral bacterial pathogen

Disease is an increasing threat to reef-building corals. One of the few identified pathogens of coral disease is the bacterium Vibrio coralliilyticus. In Vibrio cholerae, infection by a bacterial virus (bacteriophage) results in the conversion of non-pathogenic strains to pathogenic strains and this can lead to cholera pandemics. Pathogenicity islands encoded in the V. cholerae genome play an important role in pathogenesis. Here we analyse five whole genome sequences of V. coralliilyticus to examine whether virulence is similarly driven by horizontally acquired elements. We demonstrate that bacteriophage genomes encoding toxin genes with homology to those found in pathogenic V. cholerae are integrated in V. coralliilyticus genomes. Virulence factors located on chromosomal pathogenicity islands also exist in some strains of V. coralliilyticus. The presence of these genetic signatures indicates virulence in V. coralliilyticus is driven by prophages and other horizontally acquired elements. Screening for pathogens of coral disease should target conserved regions in these elements.

Different Ancestries of R Tailocins in Rhizospheric Pseudomonas Isolates

Bacterial genomes accommodate a variety of mobile genetic elements, including bacteriophage-related clusters that encode phage tail-like protein complexes playing a role in interactions with eukaryotic or prokaryotic cells. Such tailocins are unable to replicate inside target cells due to the lack of a phage head with associated DNA. A subset of tailocins mediate antagonistic activities with bacteriocin-like specificity. Functional characterization of bactericidal tailocins of two Pseudomonas putida rhizosphere isolates revealed not only extensive similarity with the tail assembly module of the Pseudomonas aeruginosa R-type pyocins but also differences in genomic integration site, regulatory genes, and lytic release modules. Conversely, these three features are quite similar between strains of the P. putida and Pseudomonas fluorescens clades, although phylogenetic analysis of tail genes suggests them to have evolved separately. Unlike P. aeruginosa R pyocin elements, the tailocin gene clusters of other pseudomonads frequently carry cargo genes, including bacteriocins. Compared with P. aeruginosa, the tailocin tail fiber sequences that act as specificity determinants have diverged much more extensively among the other pseudomonad species, mostly isolates from soil and plant environments. Activity of the P. putida antibacterial particles requires a functional lipopolysaccharide layer on target cells, but contrary to R pyocins from P. aeruginosa, strain susceptibilities surpass species boundaries.

Unique and conserved genome regions in Vibrio harveyi and related species in comparison with the shrimp pathogen Vibrio harveyi CAIM 1792

Vibrio harveyi CAIM 1792 is a marine bacterial strain that causes mortality in farmed shrimp in north-west Mexico, and the identification of virulence genes in this strain is important for understanding its pathogenicity. The aim of this work was to compare the V. harveyi CAIM 1792 genome with related genome sequences to determine their phylogenic relationship and explore unique regions in silico that differentiate this strain from other V. harveyi strains. Twenty-one newly sequenced genomes were compared in silico against the CAIM 1792 genome at nucleotidic and predicted proteome levels. The proteome of CAIM 1792 had higher similarity to those of other V. harveyi strains (78%) than to those of the other closely related species Vibrio owensii (67%), Vibrio rotiferianus (63%) and Vibrio campbellii (59%). Pan-genome ORFans trees showed the best fit with the accepted phylogeny based on DNA-DNA hybridization and multi-locus sequence analysis of 11 concatenated housekeeping genes. SNP analysis clustered 34/38 genomes within their accepted species. The pangenomic and SNP trees showed that V. harveyi is the most conserved of the four species studied and V. campbellii may be divided into at least three subspecies, supported by intergenomic distance analysis. blastp atlases were created to identify unique regions among the genomes most related to V. harveyi CAIM 1792; these regions included genes encoding glycosyltransferases, specific type restriction modification systems and a transcriptional regulator, LysR, reported to be involved in virulence, metabolism, quorum sensing and motility.

Strategies of Vibrio parahaemolyticus to acquire nutritional iron during host colonization

Iron is an essential element for the growth and development of virtually all living organisms. As iron acquisition is critical for the pathogenesis, a host defense strategy during infection is to sequester iron to restrict the growth of invading pathogens. To counteract this strategy, bacteria such as Vibrio parahaemolyticus have adapted to such an environment by developing mechanisms to obtain iron from human hosts. This review focuses on the multiple strategies employed by V. parahaemolyticus to obtain nutritional iron from host sources. In these strategies are included the use of siderophores and xenosiderophores, proteases and iron-protein receptor. The host sources used by V. parahaemolyticus are the iron-containing proteins transferrin, hemoglobin, and hemin. The implications of iron acquisition systems in the virulence of V. parahaemolyticus are also discussed.

The Molecular Switch of Telomere Phages: High Binding Specificity of the PY54 Cro Lytic Repressor to a Single Operator Site

Temperate bacteriophages possess a molecular switch, which regulates the lytic and lysogenic growth. The genomes of the temperate telomere phages N15, PY54 and ϕKO2 harbor a primary immunity region (immB) comprising genes for the prophage repressor, the lytic repressor and a putative antiterminator. The roles of these products are thought to be similar to those of the lambda proteins CI, Cro and Q, respectively. Moreover, the gene order and the location of several operator sites in the prototype telomere phage N15 and in ϕKO2 are also reminiscent of lambda-like phages. By contrast, in silico analyses revealed the presence of only one operator (O_R3) in PY54. The purified PY54 Cro protein was used for EMSA studies demonstrating that it exclusively binds to a 16-bp palindromic site (O_R3) upstream of the prophage repressor gene. The O_R3 operator sequences of PY54 and ϕKO2/N15 only differ by their peripheral base pairs, which are responsible for Cro specificity. PY54 cI and cro transcription is regulated by highly active promoters initiating the synthesis of a homogenious species of leaderless mRNA. The location of the PY54 Cro binding site and of the identified promoters suggests that the lytic repressor suppresses cI transcription but not its own synthesis. The results indicate an unexpected diversity of the growth regulation mechanisms in lambda-related phages.

Characterization of a large novel phage-like plasmid in Salmonella enterica serovar Typhimurium

Salmonella enterica serovar Typhimurium is a food-borne pathogen and a leading cause of gastroenteritis in humans. Recently, we sequenced a phage-type DT108 strain (L945) and found reads with high similarity to both Salmonella typhi strain CT18 plasmid pHCM2 and bacteriophage SSU5. In this study, we completely sequenced the novel phage-like plasmid which was designated as pSTM_Φ. The presence of this phage-like plasmid was examined in a collection of 284 Salmonella Typhimurium isolates using PCR of the parB gene and only one other isolate (L946) was found to carry the phage-like plasmid suggesting that it is infrequently present amongst Salmonella Typhimurium isolates. pSTM_Φ is a circular phage-like plasmid of 107.7 kb encoding 132 coding regions (ORFs) with the majority of the ORFs encoding hypothetical proteins. Comparative analysis with other closely related phage-like plasmids and the SSU5 phage revealed that there were four divergent lineages of phage-like plasmids found in the family of Enterobacteriaceae. In conclusion, pSTM_Φ is a new member of an emerging family of phage-like plasmids.

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

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

Bacteriophage remediation of bacterial pathogens in aquaculture: a review of the technology

Bacteriophages have been proposed as an alternative to antibiotic usage and several studies on their application in aquaculture have been reported. This review highlights progress to date on phage therapies for the following fish and shellfish diseases and associated pathogens: hemorrhagic septicemia (Aeromonas hydrophila) in loaches, furunculosis (Aeromonas salmonicida) in trout and salmon, edwardsiellosis (Edwardsiella tarda) in eel, columnaris disease (Flavobacterium columnare) in catfish, rainbow trout fry syndrome or cold water disease (Flavobacterium psychrophilum) in trout and salmon, lactococcosis (Lactococcus spp.) in yellowtail, ulcerative skin lesions (Pseudomonas aeruginosa) in freshwater catfish, bacterial hemorrhagic ascites disease (Pseudomonas plecoglossicida) in ayu fish, streptococcosis (Streptococcus iniae) in flounder, and luminescent vibriosis (Vibrio harveyi) in shrimp. Information is reviewed on phage specificity, host resistance, routes of administration, and dosing of fish and shellfish. Limitations in phage research are described and recommended guidelines are provided for conducting future phage studies involving fish and shellfish.

Ribosomally encoded antibacterial proteins and peptides from Pseudomonas

Members of the Pseudomonas genus produce diverse secondary metabolites affecting other bacteria, fungi or predating nematodes and protozoa but are also equipped with the capacity to secrete different types of ribosomally encoded toxic peptides and proteins, ranging from small microcins to large tailocins. Studies with the human pathogen Pseudomonas aeruginosa have revealed that effector proteins of type VI secretion systems are part of the antibacterial armamentarium deployed by pseudomonads. A novel class of antibacterial proteins with structural similarity to plant lectins was discovered by studying antagonism among plant-associated Pseudomonas strains. A genomic perspective on pseudomonad bacteriocinogeny shows that the modular architecture of S pyocins of P. aeruginosa is retained in a large diversified group of bacteriocins, most of which target DNA or RNA. Similar modularity is present in as yet poorly characterized Rhs (recombination hot spot) proteins and CDI (contact-dependent inhibition) proteins. Well-delimited domains for receptor recognition or cytotoxicity enable the design of chimeric toxins with novel functionalities, which has been applied successfully for S and R pyocins. Little is known regarding how these antibacterials are released and ultimately reach their targets. Other remaining issues concern the identification of environmental triggers activating these systems and assessment of their ecological impact in niches populated by pseudomonads.

Recognition of Protein-coding Genes Based on Z-curve Algorithms

Recognition of protein-coding genes, a classical bioinformatics issue, is an absolutely needed step for annotating newly sequenced genomes. The Z-curve algorithm, as one of the most effective methods on this issue, has been successfully applied in annotating or re-annotating many genomes, including those of bacteria, archaea and viruses. Two Z-curve based ab initio gene-finding programs have been developed: ZCURVE (for bacteria and archaea) and ZCURVE_V (for viruses and phages). ZCURVE_C (for 57 bacteria) and Zfisher (for any bacterium) are web servers for re-annotation of bacterial and archaeal genomes. The above four tools can be used for genome annotation or re-annotation, either independently or combined with the other gene-finding programs. In addition to recognizing protein-coding genes and exons, Z-curve algorithms are also effective in recognizing promoters and translation start sites. Here, we summarize the applications of Z-curve algorithms in gene finding and genome annotation.

A Brief Review: The Z-curve Theory and its Application in Genome Analysis

In theoretical physics, there exist two basic mathematical approaches, algebraic and geometrical methods, which, in most cases, are complementary. In the area of genome sequence analysis, however, algebraic approaches have been widely used, while geometrical approaches have been less explored for a long time. The Z-curve theory is a geometrical approach to genome analysis. The Z-curve is a three-dimensional curve that represents a given DNA sequence in the sense that each can be uniquely reconstructed given the other. The Z-curve, therefore, contains all the information that the corresponding DNA sequence carries. The analysis of a DNA sequence can then be performed through studying the corresponding Z-curve. The Z-curve method has found applications in a wide range of areas in the past two decades, including the identifications of protein-coding genes, replication origins, horizontally-transferred genomic islands, promoters, translational start sides and isochores, as well as studies on phylogenetics, genome visualization and comparative genomics. Here, we review the progress of Z-curve studies from aspects of both theory and applications in genome analysis.

What does the talking?: quorum sensing signalling genes discovered in a bacteriophage genome

The transfer of novel genetic material into the genomes of bacterial viruses (phages) has been widely documented in several host-phage systems. Bacterial genes are incorporated into the phage genome and, if retained, subsequently evolve within them. The expression of these phage genes can subvert or bolster bacterial processes, including altering bacterial pathogenicity. The phage phiCDHM1 infects Clostridium difficile, a pathogenic bacterium that causes nosocomial infections and is associated with antibiotic treatment. Genome sequencing and annotation of phiCDHM1 shows that despite being closely related to other C. difficile myoviruses, it has several genes that have not been previously reported in any phage genomes. Notably, these include three homologs of bacterial genes from the accessory gene regulator (agr) quorum sensing (QS) system. These are; a pre-peptide (AgrD) of an autoinducing peptide (AIP), an enzyme which processes the pre-peptide (AgrB) and a histidine kinase (AgrC) that detects the AIP to activate a response regulator. Phylogenetic analysis of the phage and C. difficile agr genes revealed that there are three types of agr loci in this species. We propose that the phage genes belonging to a third type, agr3, and have been horizontally transferred from the host. AgrB and AgrC are transcribed during the infection of two different strains. In addition, the phage agrC appears not to be confined to the phiCDHM1 genome as it was detected in genetically distinct C. difficile strains. The discovery of QS gene homologs in a phage genome presents a novel way in which phages could influence their bacterial hosts, or neighbouring bacterial populations. This is the first time that these QS genes have been reported in a phage genome and their distribution both in C. difficile and phage genomes suggests that the agr3 locus undergoes horizontal gene transfer within this species.

Complete genome sequence of Vibrio parahaemolyticus bacteriophage vB_VpaM_MAR

Vibrio parahaemolyticus is a major pathogen that is mainly associated with seafood and is a global food safety issue. Our objective was to isolate and completely sequence a specific phage against this bacterium. Phage vB_VpaM_MAR is able to lyse 76% of the V. parahaemolyticus strains tested. MAR belongs to the Myoviridae family and has a genome comprised of double-stranded DNA with a size of 41,351 bp, a G+C content of 51.3%, and 62 open reading frames (ORFs). Bioinformatic analysis showed that phage MAR is closely related to Vibrio phages VHML, VP58.5, and VP882 and Halomonas aquamarina phage ΦHAP-1.

Linear chromosome-generating system of Agrobacterium tumefaciens C58: protelomerase generates and protects hairpin ends

Agrobacterium tumefaciens C58, the pathogenic bacteria that causes crown gall disease in plants, harbors one circular and one linear chromosome and two circular plasmids. The telomeres of its unusual linear chromosome are covalently closed hairpins. The circular and linear chromosomes co-segregate and are stably maintained in the organism. We have determined the sequence of the two ends of the linear chromosome thus completing the previously published genome sequence of A. tumefaciens C58. We found that the telomeres carry nearly identical 25-bp sequences at the hairpin ends that are related by dyad symmetry. We further showed that its Atu2523 gene encodes a protelomerase (resolvase) and that the purified enzyme can generate the linear chromosomal closed hairpin ends in a sequence-specific manner. Agrobacterium protelomerase, whose presence is apparently limited to biovar 1 strains, acts via a cleavage-and-religation mechanism by making a pair of transient staggered nicks invariably at 6-bp spacing as the reaction intermediate. The enzyme can be significantly shortened at both the N and C termini and still maintain its enzymatic activity. Although the full-length enzyme can uniquely bind to its product telomeres, the N-terminal truncations cannot. The target site can also be shortened from the native 50-bp inverted repeat to 26 bp; thus, the Agrobacterium hairpin-generating system represents the most compact activity of all hairpin linear chromosome- and plasmid-generating systems to date. The biochemical analyses of the protelomerase reactions further revealed that the tip of the hairpin telomere may be unusually polymorphically capable of accommodating any nucleotide.

A rebuttal to the comments on the genome order index and the Z-curve

Background

Elhaik, Graur and Josic recently commented on the genome order index (S) and the Z-curve (Elhaik et al. Biol Direct 2010, 5: 10). S is a quantity defined as S = a² + c² + g² + t², where a, c, g and t denote corresponding base frequencies. The Z-curve is a three dimensional curve that represents a DNA sequence in the manner that each can be uniquely reconstructed given the other. Elhaik et al. made 4 major claims. 1) In the previous mapping system with the regular tetrahedron, calculation of the radius of the inscribed sphere is "a mathematical error". 2) S follows an exponential distribution and is narrowly distributed with a range of (0.25 - 0.33). 3) Based on the Chargaff's second parity rule (PR2), "S is equivalent to H [Shannon entropy]" and they are derivable from each other. 4) Z-curve "suffers from over dimensionality", because based on the analysis of 235 bacterial genomes, x and y components contributed only less than 1% of the variance and therefore "would be of little use".

Results

1) Elhaik et al. mistakenly neglected the parameter 4/3 when calculating the radius of the inscribed sphere. 2) The exponential distribution of S is a restatement of our previous conclusion, and the range of (0.25 - 0.33) only paraphrases the previously suggested S range (0.25 -1/3). 3) Elhaik et al. incorrectly disregard deviations from PR2 by treating the deviations as 0 altogether, reduce S and H, both having 4 variables, a, c, g and t, into functions of one single variable, a only, and apply this treatment to all DNA sequences as the basis of their "demonstration", which is therefore invalid. 4) Elhaik et al. confuse numeral smallness with biological insignificance, and disregard the distributions of purine/pyrimidine and amino/keto bases (x and y components), the variations of which, although can be less than that of GC content, contain rich information that is important and useful, such as in locating replication origins of bacterial and archaeal genomes, and in studies of gene recognition in various species.

Conclusion

Elhaik et al. confuse S (a single number) with Z-curve (a series of 3D coordinates), which are distinct. To use S as a case study of Z-curve, by itself, is invalid. S and H are neither equivalent nor derivable from each other. The criticisms of Elhaik, Graur and Josic are wrong.

Reviewers

This article was reviewed by Erik van Nimwegen.

N15: the linear phage-plasmid

The lambdoid phage N15 of Escherichia coli is very unusual among temperate phages in that its prophage is not integrated into chromosome but is a linear plasmid molecule with covalently closed ends. Upon infection the phage DNA circularises via cohesive ends, then phage-encoded enzyme, protelomerase, cuts at an inverted repeat site and forms hairpin ends (telomeres) of the linear plasmid prophage. Replication of the N15 prophage is initiated at an internally located ori site and proceeds bidirectionally resulting in formation of duplicated telomeres. Then the N15 protelomerase cuts duplicated telomeres generating two linear plasmid molecules with hairpin telomeres. Stable inheritance of the plasmid prophage is ensured by partitioning operon similar to the F factor sop operon. Unlike F sop, the N15 centromere consists of four inverted repeats dispersed in the genome. The multiplicity and dispersion of centromeres are required for efficient partitioning of a linear plasmid. The centromeres are located in N15 genome regions involved in phage replication and control of lysogeny, and binding of partition proteins at these sites regulates these processes. Two N15-related lambdoid Siphoviridae phages, φKO2 in Klebsiella oxytoca and pY54 in Yersinia enterocolitica, also lysogenize their hosts as linear plasmids, as well as Myoviridae marine phages VP882 and VP58.5 in Vibrio parahaemolyticus and ΦHAP-1 in Halomonas aquamarina. The genomes of all these phages contain similar protelomerase genes, lysogeny modules and replication genes, as well as plasmid-partitioning genes, suggesting that these phages may belong to a group diverged from a common ancestor.

A new group of cosmopolitan bacteriophages induce a carrier state in the pandemic strain of Vibrio parahaemolyticus

A clonal population of pathogenic Vibrio parahaemolyticus O3 : K6 serovar has spread in coastal waters, causing outbreaks worldwide since 1996. Bacteriophage infection is one of the main factors affecting bacterial strain concentration in the ocean. We studied the occurrence and properties of phages infecting this V. parahaemolyticus pandemic strain in coastal waters. Analysing 143 samples, phages were found in 13. All isolates clustered in a closely related group of podophages with at least 90% nucleotide sequence identity in three essential genes, despite distant geographical origins. These bacteriophages were able to multiply on the V. parahaemolyticus pandemic strain, but the impact on host concentration and subsequent growth was negligible. Infected bacteria continued producing the phage but were not lysogenized. The phage genome of prototype strain VP93 is 43 931 nucleotides and contains 337 bp direct terminal repeats at both ends. VP93 is the first non-Pseudomonas phage related to the PhiKMV-like subgroup of the T7 supergroup. The lack of a major effect on host growth suggests that these phages exert little control on the propagation of the pandemic strain in the environment. This form of phage growth can be modelled if phage-sensitive and -resistant cells that convert to each other with a high frequency are present in clonal cultures of pandemic V. parahaemolyticus.