Prophage Tracer: precisely tracing prophages in prokaryotic genomes using overlapping split-read alignment

The life cycle of temperate phages includes a lysogenic cycle stage when the phage integrates into the host genome and becomes a prophage. However, the identification of prophages that are highly divergent from known phages remains challenging. In this study, by taking advantage of the lysis-lysogeny switch of temperate phages, we designed Prophage Tracer, a tool for recognizing active prophages in prokaryotic genomes using short-read sequencing data, independent of phage gene similarity searching. Prophage Tracer uses the criterion of overlapping split-read alignment to recognize discriminative reads that contain bacterial (attB) and phage (attP) att sites representing prophage excision signals. Performance testing showed that Prophage Tracer could predict known prophages with precise boundaries, as well as novel prophages. Two novel prophages, dsDNA and ssDNA, encoding highly divergent major capsid proteins, were identified in coral-associated bacteria. Prophage Tracer is a reliable data mining tool for the identification of novel temperate phages and mobile genetic elements. The code for the Prophage Tracer is publicly available at https://github.com/WangLab-SCSIO/Prophage_Tracer.     

Identification of prophages in bacterial genomes by dinucleotide relative abundance difference

Background

Prophages are integrated viral forms in bacterial genomes that have been found to contribute to interstrain genetic variability. Many virulence-associated genes are reported to be prophage encoded. Present computational methods to detect prophages are either by identifying possible essential proteins such as integrases or by an extension of this technique, which involves identifying a region containing proteins similar to those occurring in prophages. These methods suffer due to the problem of low sequence similarity at the protein level, which suggests that a nucleotide based approach could be useful.

Methodology

Earlier dinucleotide relative abundance (DRA) have been used to identify regions, which deviate from the neighborhood areas, in genomes. We have used the difference in the dinucleotide relative abundance (DRAD) between the bacterial and prophage DNA to aid location of DNA stretches that could be of prophage origin in bacterial genomes. Prophage sequences which deviate from bacterial regions in their dinucleotide frequencies are detected by scanning bacterial genome sequences. The method was validated using a subset of genomes with prophage data from literature reports. A web interface for prophage scan based on this method is available at http://bicmku.in:8082/prophagedb/dra.html. Two hundred bacterial genomes which do not have annotated prophages have been scanned for prophage regions using this method.

Conclusions

The relative dinucleotide distribution difference helps detect prophage regions in genome sequences. The usefulness of this method is seen in the identification of 461 highly probable loci pertaining to prophages which have not been annotated so earlier. This work emphasizes the need to extend the efforts to detect and annotate prophage elements in genome sequences.     

High resolution assembly and characterization of genomes of Canadian isolates of Salmonella Enteritidis

Background

There is a need to characterize genomes of the foodborne pathogen, Salmonella enterica serovar Enteritidis (SE) and identify genetic information that could be ultimately deployed for differentiating strains of the organism, a need that is yet to be addressed mainly because of the high degree of clonality of the organism. In an effort to achieve the first characterization of the genomes of SE of Canadian origin, we carried out massively parallel sequencing of the nucleotide sequence of 11 SE isolates obtained from poultry production environments (n = 9), a clam and a chicken, assembled finished genomes and investigated diversity of the SE genome.

Results

The median genome size was 4,678,683 bp. A total of 4,833 chromosomal genes defined the pan genome of our field SE isolates consisting of 4,600 genes present in all the genomes, i.e., core genome, and 233 genes absent in at least one genome (accessory genome). Genome diversity was demonstrable by the presence of 1,360 loci showing single nucleotide polymorphism (SNP) in the core genome which was used to portray the genetic distances by means of a phylogenetic tree for the SE isolates. The accessory genome consisted mostly of previously identified SE prophage sequences as well as two, apparently full- sized, novel prophages namely a 28 kb sequence provisionally designated as SE-OLF-10058 (3) prophage and a 43 kb sequence provisionally designated as SE-OLF-10012 prophage.

Conclusions

The number of SNPs identified in the relatively large core genome of SE is a reflection of substantial diversity that could be exploited for strain differentiation as shown by the development of an informative phylogenetic tree. Prophage sequences can also be exploited for SE strain differentiation and lineage tracking. This work has laid the ground work for further studies to develop a readily adoptable laboratory test for the subtyping of SE.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-713) contains supplementary material, which is available to authorized users.     

MetaGeneAnnotator: Detecting Species-Specific Patterns of Ribosomal Binding Site for Precise Gene Prediction in Anonymous Prokaryotic and Phage Genomes

Recent advances in DNA sequencers are accelerating genome sequencing, especially in microbes, and complete and draft genomes from various species have been sequenced in rapid succession. Here, we present a comprehensive gene prediction tool, the MetaGeneAnnotator (MGA), which precisely predicts all kinds of prokaryotic genes from a single or a set of anonymous genomic sequences having a variety of lengths. The MGA integrates statistical models of prophage genes, in addition to those of bacterial and archaeal genes, and also uses a self-training model from input sequences for predictions. As a result, the MGA sensitively detects not only typical genes but also atypical genes, such as horizontally transferred and prophage genes in a prokaryotic genome. In this paper, we also propose a novel approach for analyzing the ribosomal binding site (RBS), which enables us to detect species-specific patterns of the RBSs. The MGA has the ingenious RBS model based on this approach, and precisely predicts translation starts of genes. The MGA also succeeds in improving prediction accuracies for short sequences by using the adapted RBS models (96% sensitivity and 93% specificity for 700 bp fragments). These features of the MGA expedite wide ranges of microbial genome studies, such as genome annotations and metagenome analyses.Key words: bioinformatics, gene-finding, prokaryote, phage, ribosomal binding site     

Regional Specificity of Illegitimate Recombination by the Translocatable Amplicillin-Resistance Element Tn1 in the Genome of Phage P22

Insertions of the translocatable ampicillin-resistance element Tn1 were selected in the genome of the temperate Salmonella phage P22 by growing the phage on hosts carrying the resistance plasmid RP4. Insertions of Tn1 into phage P22 are rare (10(-10) per phage) and nonrandomly distributed in the P22 genome. They are found mainly in the vicinity of the P22 ant gene. Insertions within the ant gene are found at many (at least 15) genetically separable sites, are found equally frequently in both orientations and cause irreversible loss of gene function. Some insertions in ant appear to be associated with an adjecent deletion. Prophage deletions were derived from P22::Tn1 phages by two methods. Low multiplicity transductants have nonrandomly distributed endpoints. One end is at or very near the site of the Tn1 insertion, and the other is in the vicinity of gene 12; however, there are many genetically distinguishable endpoints within gene 12. Prophage deletions selected as survivors of induction of a P22Ap mnt-ts lysogen have similarly nonrandom endpoints, with the Tn1-distal end frequently near the ant gene, as well as gene 12. Physical analysis of several prophage deletions suggests that the Tn1 is intact to the resolution of DNA electron microscopy and that the deletions begin at the end of the Tn1 insertion. These results suggest that illegitimate recombination associated with Tn1 shows regional specificity (i.e., preference for some large areas of the P22 genome over other areas), but that within these regions is quite nonspecific.     

Isolation, characterization and mapping of the N15 phasmid insertion mutations]

Prophage of N15 temperate bacteriophage is stably maintained in Escherichia coli lysogens as a 46.33 kb linear plasmid. Using different transposons we obtained 18 insertion mutants of the N15 plasmid prophage. They were analysed for plaque formation ability, stability of the plasmid state and lysogenic conversion. Restriction mapping of the insertions allowed us to localize on the map the regions necessary for lytic growth and to map the lysogenic conversion gene. A recombinant phage encoding two antibiotic resistance genes was obtained. The phage contains an additional 4.77 kb DNA fragment (over 10% of the N15 genome).     

DiffTool: building,visualizing and querying protein clusters

DiffTool is a resource to build and visualize protein clusters computed from a sequence database. The package provides a clustering tool to construct protein families according to sequence similarities and a web interface to query the corresponding clusters. A subtractive genome analysis tool selects protein families specific for a genome or a group of genomes. For each protein cluster, DiffTool includes access to sequences, coloured multiple alignments and phylogenetic trees. AVAILABILITY: A cluster database built from yeast and complete prokaryotic genomes is queryable at http://bioweb.pasteur.fr/seqanal/difftool. All the Perl sources are freely available to non-profit organizations upon request.     

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 excisive recombination are discussed.     

DNA Sequence Heterogeneity of Campylobacter jejuni CJIE4 Prophages and Expression of Prophage Genes

Campylobacter jejuni carry temperate bacteriophages that can affect the biology or virulence of the host bacterium. Known effects include genomic rearrangements and resistance to DNA transformation. C. jejuni prophage CJIE1 shows sequence variability and variability in the content of morons. Homologs of the CJIE1 prophage enhance both adherence and invasion to cells in culture and increase the expression of a specific subset of bacterial genes. Other C. jejuni temperate phages have so far not been well characterized. In this study we describe investigations into the DNA sequence variability and protein expression in a second prophage, CJIE4. CJIE4 sequences were obtained de novo from DNA sequencing of five C. jejuni isolates, as well as from whole genome sequences submitted to GenBank by other research groups. These CJIE4 DNA sequences were heterogenous, with several different insertions/deletions (indels) in different parts of the prophage genome. Two variants of a 3–4 kb region inserted within CJIE4 had different gene content that distinguished two major conserved CJIE4 prophage families. Additional indels were detected throughout the prophage. Detection of proteins in the five isolates characterized in our laboratory in isobaric Tags for Relative and Absolute Quantitation (iTRAQ) experiments indicated that prophage proteins within each of the two large indel variants were expressed during growth of the bacteria on Mueller Hinton agar plates. These proteins included the extracellular DNase associated with resistance to DNA transformation and prophage repressor proteins. Other proteins associated with known or suspected roles in prophage biology were also expressed from CJIE4, including capsid protein, the phage integrase, and MazF, a type II toxin-antitoxin system protein. Together with the results previously obtained for the CJIE1 prophage these results demonstrate that sequence variability and expression of moron genes are both general properties of temperate bacteriophages in C. jejuni.     

Prophage genomics.

The majority of the bacterial genome sequences deposited in the National Center for Biotechnology Information database contain prophage sequences. Analysis of the prophages suggested that after being integrated into bacterial genomes, they undergo a complex decay process consisting of inactivating point mutations, genome rearrangements, modular exchanges, invasion by further mobile DNA elements, and massive DNA deletion. We review the technical difficulties in defining such altered prophage sequences in bacterial genomes and discuss theoretical frameworks for the phage-bacterium interaction at the genomic level. The published genome sequences from three groups of eubacteria (low- and high-G+C gram-positive bacteria and gamma-proteobacteria) were screened for prophage sequences. The prophages from Streptococcus pyogenes served as test case for theoretical predictions of the role of prophages in the evolution of pathogenic bacteria. The genomes from further human, animal, and plant pathogens, as well as commensal and free-living bacteria, were included in the analysis to see whether the same principles of prophage genomics apply for bacteria living in different ecological niches and coming from distinct phylogenetical affinities. The effect of selection pressure on the host bacterium is apparently an important force shaping the prophage genomes in low-G+C gram-positive bacteria and gamma-proteobacteria.     

Genetic control of prophage induction in haemophilus influenzae after exposure to psoralen plus near-UV light.

Prophage S2 could be induced by psoralen plus near-UV light (PNUV) from a wild-type strain of Haemophilus influenzae, from UV light-sensitive strains uvr-1 and uvr-2 and PNUV-sensitive strains PSO1 amd PSO7, but not from a recombination-deficient strain, rec-1. The levels of prophage induction were comparable in the wild type and an ATP-dependent DNase-deficient strain, KW31, even though the PNUV-induced degradation in the latter strain was considerably lower. Prophage induction could be observed even with chloramphenicol present before, during, and 30 min after PNUV treatment.     

Fast algorithms for large-scale genome alignment and comparison

We describe a suffix-tree algorithm that can align the entire genome sequences of eukaryotic and prokaryotic organisms with minimal use of computer time and memory. The new system, MUMmer 2, runs three times faster while using one-third as much memory as the original MUMmer system. It has been used successfully to align the entire human and mouse genomes to each other, and to align numerous smaller eukaryotic and prokaryotic genomes. A new module permits the alignment of multiple DNA sequence fragments, which has proven valuable in the comparison of incomplete genome sequences. We also describe a method to align more distantly related genomes by detecting protein sequence homology. This extension to MUMmer aligns two genomes after translating the sequence in all six reading frames, extracts all matching protein sequences and then clusters together matches. This method has been applied to both incomplete and complete genome sequences in order to detect regions of conserved synteny, in which multiple proteins from one organism are found in the same order and orientation in another. The system code is being made freely available by the authors.     

Genome Analysis of the Fruiting Body-Forming Myxobacterium Chondromyces crocatus Reveals High Potential for Natural Product Biosynthesis

Here, we report the complete genome sequence of the type strain of the myxobacterial genus Chondromyces, Chondromyces crocatus Cm c5. It presents one of the largest prokaryotic genomes featuring a single circular chromosome and no plasmids. Analysis revealed an enlarged set of tRNA genes, along with reduced pressure on preferred codon usage compared to that of other bacterial genomes. The large coding capacity and the plethora of encoded secondary metabolite biosynthetic gene clusters are in line with the capability of Cm c5 to produce an arsenal of antibacterial, antifungal, and cytotoxic compounds. Known pathways of the ajudazol, chondramide, chondrochloren, crocacin, crocapeptin, and thuggacin compound families are complemented by many more natural compound biosynthetic gene clusters in the chromosome. Whole-genome comparison of the fruiting-body-forming type strain (Cm c5, DSM 14714) to an accustomed laboratory strain which has lost this ability (nonfruiting phenotype, Cm c5 fr−) revealed genetic changes in three loci. In addition to the low synteny found with the closest sequenced representative of the same family, Sorangium cellulosum, extensive genetic information duplication and broad application of eukaryotic-type signal transduction systems are hallmarks of this 11.3-Mbp prokaryotic genome.     

Contrasting rates of evolution in Pm3 loci from three wheat species and rice

The Pm3 gene from wheat confers resistance against powdery mildew and recent studies have shown that it is a member of a multigene family in the wheat genome. We compared genomic sequences ranging from 178 to 332 kb containing six Pm3-like genes and five gene fragments from orthologous loci in the A genome of wheat at three different ploidy levels. We found that the wheat Pm3 loci display an extremely dynamic evolution where sequence conservation is minimal between species and basically limited to very short sequences containing the genetic markers that define the orthology. The Pm3-like genes and their up- and downstream regions were reshuffled by multiple rearrangements, resulting in a complex mosaic of conserved and unique sequences. Comparison with rice showed that the known wheat Pm3-like genes represent only one branch of a large superfamily with several clusters in rice and suggests the presence of additional similar genes in the wheat genome. Estimates of divergence times and transposable-element insertions indicate that the Pm3 locus in wheat has undergone more drastic changes in its recent evolution than its counterpart in rice. This indicates that loci containing homologous resistance gene analogs can evolve at highly variable speeds in different species.     

Comparative Genomics of a Helicobacter pylori Isolate from a Chinese Yunnan Naxi Ethnic Aborigine Suggests High Genetic Divergence and Phage Insertion

Helicobacter pylori is a common pathogen correlated with several severe digestive diseases. It has been reported that isolates associated with different geographic areas, different diseases and different individuals might have variable genomic features. Here, we describe draft genomic sequences of H. pylori strains YN4-84 and YN1-91 isolated from patients with gastritis from the Naxi and Han populations of Yunnan, China, respectively. The draft sequences were compared to 45 other publically available genomes, and a total of 1059 core genes were identified. Genes involved in restriction modification systems, type four secretion system three (TFS3) and type four secretion system four (TFS4), were identified as highly divergent. Both YN4-84 and YN1-91 harbor intact cag pathogenicity island (cagPAI) and have EPIYA-A/B/D type at the carboxyl terminal of cagA. The vacA gene type is s1m2i1. Another major finding was a 32.5-kb prophage integrated in the YN4-84 genome. The prophage shares most of its genes (30/33) with Helicobacter pylori prophage KHP30. Moreover, a 1,886 bp transposable sequence (IS605) was found in the prophage. Our results imply that the Naxi ethnic minority isolate YN4-84 and Han isolate YN1-91 belong to the hspEAsia subgroup and have diverse genome structure. The genome has been extensively modified in several regions involved in horizontal DNA transfer. The important roles played by phages in the ecology and microevolution of H. pylori were further emphasized. The current data will provide valuable information regarding the H. pylori genome based on historic human migrations and population structure.     

Mechanism of Transfection with Deoxyribonucleic Acid from the Temperate Bacillus Bacteriophage φ105

Bacteriophage phi105 is a temperate phage for the transformable Bacillus subtilis 168. The infectivity of deoxyribonucleic acid (DNA) extracted from mature phi105 phage particles, from bacteria lysogenic for phi105 (prophage DNA), and from induced lysogenic bacteria (vegetative DNA) was examined in the B. subtilis transformation system. About one infectious center was formed per 10(8) mature DNA molecules added to competent cells, but single markers could be rescued from mature DNA by a superinfecting phage at a 10(3)- to 10(4)-fold higher frequency. Single markers in mature DNA were inactivated at an exponential rate after uptake by a competent cell. Prophage and vegetative DNA gave about one infectious center per 10(3) molecules added to competent cells. Infectious prophage DNA entered competent cells as a single molecule; it gave a majority of lytic responses. Single markers in sheared prophage DNA were inactivated at the same rate as markers in mature DNA. Prophage DNA was dependent on the bacterial rec-1 function for its infectivity, whereas vegetative DNA was not. The mechanism of transfection of B. subtilis with viral DNA is discussed, and a model for transfection with phi105 DNA is proposed.