Exploring Codon Adjustment Strategies towards Escherichia coli-Based Production of Viral Proteins Encoded by HTH1, a Novel Prophage of the Marine Bacterium Hypnocyclicus thermotrophus

Marine viral sequence space is immense and presents a promising resource for the discovery of new enzymes interesting for research and biotechnology. However, bottlenecks in the functional annotation of viral genes and soluble heterologous production of proteins hinder access to downstream characterization, subsequently impeding the discovery process. While commonly utilized for the heterologous expression of prokaryotic genes, codon adjustment approaches have not been fully explored for viral genes. Herein, the sequence-based identification of a putative prophage is reported from within the genome of Hypnocyclicus thermotrophus, a Gram-negative, moderately thermophilic bacterium isolated from the Seven Sisters hydrothermal vent field. A prophage-associated gene cluster, consisting of 46 protein coding genes, was identified and given the proposed name Hypnocyclicus thermotrophus phage H1 (HTH1). HTH1 was taxonomically assigned to the viral family Siphoviridae, by lowest common ancestor analysis of its genome and phylogeny analyses based on proteins predicted as holin and DNA polymerase. The gene neighbourhood around the HTH1 lytic cassette was found most similar to viruses infecting Gram-positive bacteria. In the HTH1 lytic cassette, an N-acetylmuramoyl-L-alanine amidase (Amidase_2) with a peptidoglycan binding motif (LysM) was identified. A total of nine genes coding for enzymes putatively related to lysis, nucleic acid modification and of unknown function were subjected to heterologous expression in Escherichia coli. Codon optimization and codon harmonization approaches were applied in parallel to compare their effects on produced proteins. Comparison of protein yields and thermostability demonstrated that codon optimization yielded higher levels of soluble protein, but codon harmonization led to proteins with higher thermostability, implying a higher folding quality. Altogether, our study suggests that both codon optimization and codon harmonization are valuable approaches for successful heterologous expression of viral genes in E. coli, but codon harmonization may be preferable in obtaining recombinant viral proteins of higher folding quality.

Comprehensive in silico survey of the Mycolicibacterium mobilome reveals an as yet underexplored diversity

The mobilome plays a crucial role in bacterial adaptation and is therefore a starting point to understand and establish the gene flow occurring in the process of bacterial evolution. This is even more so if we consider that the mobilome of environmental bacteria can be the reservoir of genes that may later appear in the clinic. Recently, new genera have been proposed in the family Mycobacteriaceae, including the genus Mycolicibacterium, which encompasses dozens of species of agricultural, biotechnological, clinical and ecological importance, being ubiquitous in several environments. The current scenario in the Mycobacteriaceae mobilome has some bias because most of the characterized mycobacteriophages were isolated using a single host strain, and the few plasmids reported mainly relate to the genus Mycobacterium. To fill in the gaps in these issues, we performed a systematic in silico study of these mobile elements based on 242 available genomes of the genus Mycolicibacterium. The analyses identified 156 putative plasmids (19 conjugative, 45 mobilizable and 92 non-mobilizable) and 566 prophages in 86 and 229 genomes, respectively. Moreover, a contig was characterized by resembling an actinomycete integrative and conjugative element (AICE). Within this diversity of mobile genetic elements, there is a pool of genes associated with several canonical functions, in addition to adaptive traits, such as virulence and resistance to antibiotics and metals (mercury and arsenic). The type-VII secretion system was a common feature in the predicted plasmids, being associated with genes encoding virulent proteins (EsxA, EsxB, PE and PPE). In addition to the characterization of plasmids and prophages of the family Mycobacteriaceae, this study showed an abundance of these genetic elements in a dozen species of the genus Mycolicibacterium.

Examination of Staphylococcus aureus Prophages Circulating in Egypt

Staphylococcus aureus infections are of growing concern given the increased incidence of antibiotic resistant strains. Egypt, like several other countries, has seen alarming increases in methicillin-resistant S. aureus (MRSA) infections. This species can rapidly acquire genes associated with resistance, as well as virulence factors, through mobile genetic elements, including phages. Recently, we sequenced 56 S. aureus genomes from Alexandria Main University Hospital in Alexandria, Egypt, complementing 17 S. aureus genomes publicly available from other sites in Egypt. In the current study, we found that the majority (73.6%) of these strains contain intact prophages, including Biseptimaviruses, Phietaviruses, and Triaviruses. Further investigation of these prophages revealed evidence of horizontal exchange of the integrase for two of the prophages. These Egyptian S. aureus prophages are predicted to encode numerous virulence factors, including genes associated with immune evasion and toxins, including the Panton-Valentine leukocidin (PVL)-associated genes lukF-PV/lukS-PV. Thus, prophages are likely to be a major contributor to the virulence of S. aureus strains in circulation in Egypt.

Characterization of Clinical and Carrier Streptococcus agalactiae and Prophage Contribution to the Strain Variability

Streptococcus agalactiae (group B Streptococcus, GBS) represents a leading cause of invasive bacterial infections in newborns and is also responsible for diseases in older and immunocompromised adults. Prophages represent an important factor contributing to the genome plasticity and evolution of new strains. In the present study, prophage content was analyzed in human GBS isolates. Thirty-seven prophages were identified in genomes of 20 representative sequenced strains. On the basis of the sequence comparison, we divided the prophages into eight groups named A–H. This division also corresponded to the clustering of phage integrase, even though several different integration sites were observed in some relative prophages. Next, PCR method was used for detection of the prophages in 123 GBS strains from adult hospitalized patients and from pregnancy screening. At least one prophage was present in 105 isolates (85%). The highest prevalence was observed for prophage group A (71%) and satellite prophage group B (62%). Other groups were detected infrequently (1–6%). Prophage distribution did not differ between clinical and screening strains, but it was unevenly distributed in MLST (multi locus sequence typing) sequence types. High content of full-length and satellite prophages detected in present study implies that prophages could be beneficial for the host bacterium and could contribute to evolution of more adapted strains.

Genomic analysis of 40 prophages located in the genomes of 16 carbapenemase-producing clinical strains of Klebsiella pneumoniae

Klebsiella pneumoniae is the clinically most important species within the genus Klebsiella and, as a result of the continuous emergence of multi-drug resistant (MDR) strains, the cause of severe nosocomial infections. The decline in the effectiveness of antibiotic treatments for infections caused by MDR bacteria has generated particular interest in the study of bacteriophages. In this study, we characterized a total of 40 temperate bacteriophages (prophages) with a genome range of 11.454-84.199 kb, predicted from 16 carbapenemase-producing clinical strains of K. pneumoniae belonging to different sequence types, previously identified by multilocus sequence typing. These prophages were grouped into the three families in the order Caudovirales (27 prophages belonging to the family Myoviridae, 10 prophages belonging to the family Siphoviridae and 3 prophages belonging to the family Podoviridae). Genomic comparison of the 40 prophage genomes led to the identification of four prophages isolated from different strains and of genome sizes of around 33.3, 36.1, 39.6 and 42.6 kb. These prophages showed sequence similarities (query cover >90 %, identity >99.9 %) with international Microbe Versus Phage (MVP) (http://mvp.medgenius.info/home) clusters 4762, 4901, 3499 and 4280, respectively. Phylogenetic analysis revealed the evolutionary proximity among the members of the four groups of the most frequently identified prophages in the bacterial genomes studied (33.3, 36.1, 39.6 and 42.6 kb), with bootstrap values of 100 %. This allowed the prophages to be classified into three clusters: A, B and C. Interestingly, these temperate bacteriophages did not infect the highest number of strains as indicated by a host-range assay, these results could be explained by the development of superinfection exclusion mechanisms. In addition, bioinformatic analysis of the 40 identified prophages revealed the presence of 2363 proteins. In total, 59.7 % of the proteins identified had a predicted function, mainly involving viral structure, transcription, replication and regulation (lysogenic/lysis). Interestingly, some proteins had putative functions associated with bacterial virulence (toxin expression and efflux pump regulators), phage defence profiles such as toxin-antitoxin modules, an anti-CRISPR/Cas9 protein, TerB protein (from terZABCDE operon) and methyltransferase proteins.

A novel uncultured marine cyanophage lineage with lysogenic potential linked to a putative marine Synechococcus 'relic' prophage

Marine cyanobacteria are important contributors to primary production in the ocean and their viruses (cyanophages) affect the ocean microbial communities. Despite reports of lysogeny in marine cyanobacteria, a genome sequence of such temperate cyanophages remains unknown although genomic analysis indicate potential for lysogeny in certain marine cyanophages. Using assemblies from Red Sea and Tara Oceans metagenomes, we recovered genomes of a novel uncultured marine cyanophage lineage, which contain, in addition to common cyanophage genes, a phycobilisome degradation protein NblA, an integrase and a split DNA polymerase. The DNA polymerase forms a monophyletic clade with a DNA polymerase from a genomic island in Synechococcus WH8016. The island contains a relic prophage that does not resemble any previously reported cyanophage but shares several genes with the newly identified cyanophages reported here. Metagenomic recruitment indicates that the novel cyanophages are widespread, albeit at low abundance. Here, we describe a novel potentially lysogenic cyanophage family, their abundance and distribution in the marine environment.

Comparative Genomics and Characterization of the Late Promoter pR' from Shiga Toxin Prophages in Escherichia coli

Shiga-toxin producing Escherichia coli (STEC) causes human illness ranging from mild diarrhea to death. The bacteriophage encoded stx genes are located in the late transcription region, downstream of the antiterminator Q. The transcription of the stx genes is directly under the control of the late promoter pR’, thus the sequence diversity of the region between Q and stx, here termed the pR’ region, may affect Stx toxin production. Here, we compared the gene structure of the pR’ region and the stx subtypes of nineteen STECs. The sequence alignment and phylogenetic analysis suggested that the pR’ region tends to be more heterogeneous than the promoter itself, even if the prophages harbor the same stx subtype. Furthermore, we established and validated transcriptional fusions of the pR’ region to the DsRed reporter gene using mitomycin C (MMC) induction. Finally, these constructs were transformed into native and non-native strains and examined with flow cytometry. The results showed that induction levels changed when pR’ regions were placed under different regulatory systems. Moreover, not every stx gene could be induced in its native host bacteria. In addition to the functional genes, the diversity of the pR’ region plays an important role in determining the level of toxin induction.

Relating Phage Genomes to Helicobacter pylori Population Structure: General Steps Using Whole-Genome Sequencing Data

The review uses the Helicobacter pylori, the gastric bacterium that colonizes the human stomach, to address how to obtain information from bacterial genomes about prophage biology. In a time of continuous growing number of genomes available, this review provides tools to explore genomes for prophage presence, or other mobile genetic elements and virulence factors. The review starts by covering the genetic diversity of H. pylori and then moves to the biologic basis and the bioinformatics approaches used for studding the H. pylori phage biology from their genomes and how this is related with the bacterial population structure. Aspects concerning H. pylori prophage biology, evolution and phylogeography are discussed.

The genome sequence of Escherichia coli tailed phage D6 and the diversity of Enterobacteriales circular plasmid prophages

The temperate Escherichia coli bacteriophage D6 can exist as a circular plasmid prophage, and we report here its 91,159bp complete genome sequence. It is a distant relative of the well-studied phage P1, but it is sufficiently different that it typifies a previously undescribed tailed phage type or cluster. Examination of the database of bacterial genome sequences revealed that phage P1 and D6 prophage plasmids are common in the Enterobacteriales, and in addition, previously described Salmonella phage SSU5 represents a different type of temperate tailed phage with a circular plasmid prophage that is also very common in this host order. This analysis also discovered additional divergent clusters of putative circular plasmid prophages within the two larger P1 and SSU5 groups (superclusters) that inhabit the Enterobacteriales as well as bacteria in several other orders in the Gamma-proteobacteria class. Very few of these sequences are annotated as putative prophages.

Genomics of Clostridium taeniosporum, an organism which forms endospores with ribbon-like appendages

Clostridium taeniosporum, a non-pathogenic anaerobe closely related to the C. botulinum Group II members, was isolated from Crimean lake silt about 60 years ago. Its endospores are surrounded by an encasement layer which forms a trunk at one spore pole to which about 12–14 large, ribbon-like appendages are attached. The genome consists of one 3,264,813 bp, circular chromosome (with 26.6% GC) and three plasmids. The chromosome contains 2,892 potential protein coding sequences: 2,124 have specific functions, 147 have general functions, 228 are conserved but without known function and 393 are hypothetical based on the fact that no statistically significant orthologs were found. The chromosome also contains 101 genes for stable RNAs, including 7 rRNA clusters. Over 84% of the protein coding sequences and 96% of the stable RNA coding regions are oriented in the same direction as replication. The three known appendage genes are located within a single cluster with five other genes, the protein products of which are closely related, in terms of sequence, to the known appendage proteins. The relatedness of the deduced protein products suggests that all or some of the closely related genes might code for minor appendage proteins or assembly factors. The appendage genes might be unique among the known clostridia; no statistically significant orthologs were found within other clostridial genomes for which sequence data are available. The C. taeniosporum chromosome contains two functional prophages, one Siphoviridae and one Myoviridae, and one defective prophage. Three plasmids of 5.9, 69.7 and 163.1 Kbp are present. These data are expected to contribute to future studies of developmental, structural and evolutionary biology and to potential industrial applications of this organism.

The gut virome of the protochordate model organism, Ciona intestinalis subtype A

The identification of host-specific bacterial and viral communities associated with diverse animals has led to the concept of the metaorganism, which defines the animal and all of its associated microbes as a single unit. Here we sequence the viruses found in the gut (i.e., the gut virome) of the marine invertebrate model system, Ciona intestinalis subtype A, in samples collected one year apart. We present evidence for a host-associated virome that is distinct from the surrounding seawater and contains some temporally-stable members. Comparison of gut tissues before and after starvation in virus-free water enabled the differentiation between the Ciona-specific virome and transient viral communities associated with dietary sources. The Ciona gut viromes were dominated by double-stranded DNA tailed phages (Order Caudovirales) and sequence assembly yielded a number of complete circular phage genomes, most of which were highly divergent from known genomes. Unique viral communities were found in distinct gut niches (stomach, midgut and hindgut), paralleling the compartmentalization of bacterial communities. Additionally, integrase and excisionase genes, including many that are similar to prophage sequences within the genomes of bacterial genera belonging to the Ciona core microbiome, were prevalent in the viromes, indicating the active induction of prophages within the gut ecosystem. Knowledge of the gut virome of this model organism lays the foundation for studying the interactions between viruses, bacteria, and host immunity.

Evolution of a zoonotic pathogen: investigating prophage diversity in enterohaemorrhagic Escherichia coli O157 by long-read sequencing

Enterohaemorrhagic Escherichia coli (EHEC) O157 is a zoonotic pathogen for which colonization of cattle and virulence in humans is associated with multiple horizontally acquired genes, the majority present in active or cryptic prophages. Our understanding of the evolution and phylogeny of EHEC O157 continues to develop primarily based on core genome analyses; however, such short-read sequences have limited value for the analysis of prophage content and its chromosomal location. In this study, we applied Single Molecule Real Time (SMRT) sequencing, using the Pacific Biosciences long-read sequencing platform, to isolates selected from the main sub-clusters of this clonal group. Prophage regions were extracted from these sequences and from published reference strains. Genome position and prophage diversity were analysed along with genetic content. Prophages could be assigned to clusters, with smaller prophages generally exhibiting less diversity and preferential loss of structural genes. Prophages encoding Shiga toxin (Stx) 2a and Stx1a were the most diverse, and more variable compared to prophages encoding Stx2c, further supporting the hypothesis that Stx2c-prophage integration was ancestral to acquisition of other Stx types. The concept that phage type (PT) 21/28 (Stx2a+, Stx2c+) strains evolved from PT32 (Stx2c+) was supported by analysis of strains with excised Stx-encoding prophages. Insertion sequence elements were over-represented in prophage sequences compared to the rest of the genome, showing integration in key genes such as stx and an excisionase, the latter potentially acting to capture the bacteriophage into the genome. Prophage profiling should allow more accurate prediction of the pathogenic potential of isolates.

Induction and characterization of a lysogenic bacteriophage of Lactococcus garvieae isolated from marine fish species

This study investigated the presence of prophages in Lactococcus garvieae isolated from several marine fish species in Japan. Representative strains of 16 bacterial genotypes (S1-S16) selected from more than 400 L. garvieae isolates were used to induce lysogenic bacteriophages. These strains were treated with 500 ng mL(-1) freshly prepared mitomycin C. A cross-spotting assay was performed to validate the lysogenic and indicator strains. The lysogenic strains were selected for isolation and concentration of the phages. Phage DNA was digested with EcoRI for biased sinusoidal field gel electrophoresis analysis. Polymerase chain reaction (PCR) was used to detect integrated prophage DNA. Of the 16 representative bacterial genotypes, 12 strains integrated prophages as indicated by the PCR assay, and 10 phages were detected and isolated using two indicator bacterial strains. Analysis of genomic DNA showed that these phages were homologous and named as PLgT-1. Transmission electron microscopy revealed that the morphology of PLgT-1 was consistent with the virus family Siphoviridae. PCR analysis of the prophage DNA revealed that all of the S1 genotype strains were lysogenic (30/30), but none of the S16 genotype strains were lysogenic (0/30). This is the first study to investigate lysogenic bacteriophages from L. garvieae.

Predominance of Single Prophage Carrying a CRISPR/cas System in "Candidatus Liberibacter asiaticus" Strains in Southern China

“Candidatus Liberibacter asiaticus” (CLas) is an uncultureable α-proteobacterium associated with citrus Huanglongbing (HLB, yellow shoot disease), a highly destructive disease affecting citrus production worldwide. HLB was observed in Guangdong Province of China over a hundred years ago and remains endemic there. Little is known about CLas biology due to its uncultureable nature. This study began with the genome sequence analysis of CLas Strain A4 from Guangdong in the prophage region. Within the two currently known prophage types, Type 1 (SC1-like) and Type 2 (SC2-like), A4 genome contained only a Type 2 prophage, CGdP2, namely. An analysis on CLas strains collected in Guangdong showed that Type 2 prophage dominated the bacterial population (82.6%, 71/86). An extended survey covering five provinces in southern China also revealed the predominance of single prophage (Type 1 or Type 2) in the CLas population (90.4%, 169/187). CLas strains with two and no prophage types accounted for 7.2% and 2.8%, respectively. In silico analyses on CGdP2 identified a CRISPR (clustered regularly interspaced short palindromic repeats)/cas (CRISPR-associated protein genes) system, consisting of four 22 bp repeats, three 23 bp spacers and 9 predicted cas. Similar CRISPR/cas systems were detected in all 10 published CLas prophages as well as 13 CLas field strains in southern China. Both Type 1 and Type 2 prophages shared almost identical sequences in spacer 1 and 3 but not spacer 2. Considering that the function of a CRISPR/cas system was to destroy invading DNA, it was hypothesized that a pre-established CLas prophage could use its CRISPR/cas system guided by spacer 1 and/or 3 to defeat the invasion of the other phage/prophage. This hypothesis explained the predominance of single prophage type in the CLas population in southern China. This is the first report of CRISPR/cas system in the “Ca. Liberibacter” genera.

Characterization of Functional Prophages in Clostridium difficile

Bacteriophages (phages) are present in almost, if not all ecosystems. Some of these bacterial viruses are present as latent "prophages," either integrated within the chromosome of their host, or as episomal DNAs. Since prophages are ubiquitous throughout the bacterial world, there has been a sustained interest in trying to understand their contribution to the biology of their host. Clostridium difficile is no exception to that rule and with the recent release of hundreds of bacterial genome sequences, there has been a growing interest in trying to identify and classify these prophages. Besides their identification in bacterial genomes, there is also growing interest in determining the functionality of C. difficile prophages, i.e., their capacity to escape their host and reinfect a different strain, thereby promoting genomic evolution and horizontal transfer of genes through transduction, for example of antibiotic resistance genes. There is also some interest in using therapeutic phages to fight C. difficile infections.The objective of this chapter is to share with the broader C. difficile research community the expertise we developed in the study of C. difficile temperate phages. In this chapter, we describe a general "pipeline" comprising a series of experiments that we use in our lab to identify, induce, isolate, propagate, and characterize prophages. Our aim is to provide readers with the necessary basic tools to start studying C. difficile phages.

Complete genome analysis of two new bacteriophages isolated from impetigo strains of Staphylococcus aureus

Exfoliative toxin A (ETA)-coding temperate bacteriophages are leading contributors to the toxic phenotype of impetigo strains of Staphylococcus aureus. Two distinct eta gene-positive bacteriophages isolated from S. aureus strains which recently caused massive outbreaks of pemphigus neonatorum in Czech maternity hospitals were characterized. The phages, designated ϕB166 and ϕB236, were able to transfer the eta gene into a prophageless S. aureus strain which afterwards converted into an ETA producer. Complete phage genome sequences were determined, and a comparative analysis of five designed genomic regions revealed major variances between them. They differed in the genome size, number of open reading frames, genome architecture, and virion protein patterns. Their high mutual sequence similarity was detected only in the terminal regions of the genome. When compared with the so far described eta phage genomes, noticeable differences were found. Thus, both phages represent two new lineages of as yet not characterized bacteriophages of the Siphoviridae family having impact on pathogenicity of impetigo strains of S. aureus.

Genomic analysis of a phage and prophage from a Bacillus thuringiensis strain

Bacteriophages have been found to be the most abundant and also potentially most diverse biological entities on Earth. In the present study, Bacillus phages were isolated rapidly and shown to have a high degree of diversity. The genomes of a newly isolated phage, phiCM3, and a prophage, proCM3, from the Bacillus thuringiensis strain YM-03 were sequenced and characterized. Comparative genome analysis showed that the phiCM3 genome is highly similar to the genomes of eight other Bacillus phages and seven of these phages were classified as the Wβ group of phages. Analysis of the differential evolution of the genes in the Wβ-group phages indicated that the genes encoding the antirepressor and tail fibre protein were more highly conserved than those encoding the major capsid protein, DNA replication protein, and RNA polymerase σ factor, which might have diverged to acquire mechanisms suitable for survival in different microbial hosts. Genome analysis of proCM3 revealed that proCM3 might be a defective phage because of mutations in the minor structural protein, and it was not inducible by mitomycin C treatment. The proCM3 genome was similar to those of two lytic Bacillus phages in sequence, but had a different genomic structure, composed of three regions in a different order. These data suggest that the three phages might have had a common ancestor and that genome rearrangement might have occurred during evolution. The findings of this study enrich our current knowledge of Bacillus phage diversity and evolution, especially for the Wβ-group and TP21-L-like phages, and may help the development of practical applications of Bacillus phages.

A modular view of the bacteriophage genomic space: identification of host and lifestyle marker modules

Bacteriophage genomes can be regarded as an ensemble of modules which are accessible to the whole phage population via recombination. The time spent by prophages in the bacterial host provides them with the opportunity to exchange modules with other prophages or infecting phages. Here we analyze the modular structure of a set of 457 phages and 760 prophages extracted from completely sequenced bacterial genomes using the ACLAME database and its associated tools. We identified 91 modules of proteins with similar phylogenetic profiles. Of these, 25 and 6 are associated with temperate and virulent phages, respectively; 57 are restricted to a host or small group of hosts; and 55 could be annotated with a phage function. We use the transposable phages as a study case and show how the inclusion of prophages allows us to unveil new types of genome organization (i.e. novel module combinations) and obtain insight into the host range for this particular group, highlighting the utility of prophage prediction to better characterize phage diversity.

'Ca. Liberibacter asiaticus' carries an excision plasmid prophage and a chromosomally integrated prophage that becomes lytic in plant infections

Huanglongbing (HLB), also known as citrus greening, is a lethal disease of citrus caused by several species of 'Candidatus Liberibacter', a psyllid-transmitted, phloem-limited, alpha proteobacteria. 'Ca. Liberibacter asiaticus' is widespread in Florida citrus. The recently published 'Ca. L. asiaticus' psy62 genome, derived from a psyllid, revealed a prophage-like region of DNA in the genome, but phage have not been associated with 'Ca. L. asiaticus' to date. In the present study, shotgun sequencing and a fosmid DNA library of curated 'Ca. L. asiaticus' UF506, originally derived from citrus symptomatic for HLB, revealed two largely homologous, circular phage genomes, SC1 and SC2. SC2 encoded putative adhesin and peroxidase genes that had not previously been identified in 'Ca. L. asiaticus' and which may be involved in lysogenic conversion. SC2 also appeared to lack lytic cycle genes and replicated as a prophage excision plasmid, in addition to being found integrated in tandem with SC1 in the UF506 chromosome. By contrast, SC1 carried suspected lytic cycle genes and was found in nonintegrated, lytic cycle forms only in planta. Phage particles associated with 'Ca. L. asiaticus' were found in the phloem of infected periwinkles by transmission electron microscopy. In psyllids, both SC1 and SC2 were found only as prophage.

Viruses in the faecal microbiota of monozygotic twins and their mothers

Viral diversity and life cycles are poorly understood in the human gut and other body habitats. Phages and their encoded functions may provide informative signatures of a human microbiota and of microbial community responses to various disturbances, and may indicate whether community health or dysfunction is manifest after apparent recovery from a disease or therapeutic intervention. Here we report sequencing of the viromes (metagenomes) of virus-like particles isolated from faecal samples collected from healthy adult female monozygotic twins and their mothers at three time points over a one-year period. We compared these data sets with data sets of sequenced bacterial 16S ribosomal RNA genes and total-faecal-community DNA. Co-twins and their mothers share a significantly greater degree of similarity in their faecal bacterial communities than do unrelated individuals. In contrast, viromes are unique to individuals regardless of their degree of genetic relatedness. Despite remarkable interpersonal variations in viromes and their encoded functions, intrapersonal diversity is very low, with >95% of virotypes retained over the period surveyed, and with viromes dominated by a few temperate phages that exhibit remarkable genetic stability. These results indicate that a predatory viral-microbial dynamic, manifest in a number of other characterized environmental ecosystems, is notably absent in the very distal intestine.

Subclassification and targeted characterization of prophage-encoded two-component cell lysis cassette

Bacteriophage induced lysis of host bacterial cell is mediated by a two component cell lysis cassette comprised of holin and lysozyme. Prophages are integrated forms of bacteriophages in bacterial genomes providing a repertoire for bacterial evolution. Analysis using the prophage database (http://bicmku.in:8082) constructed by us showed 47 prophages were associated with putative two component cell lysis genes. These proteins cluster into four different subgroups. In this process, a putative holin (essd) and endolysin (ybcS), encoded by the defective lambdoid prophage DLP12 was found to be similar to two component cell lysis genes in functional bacteriophages like p21 and P1. The holin essd was found to have a characteristic dual start motif with two transmembrane regions and C-terminal charged residues as in class II holins. Expression of a fusion construct of essd in Escherichia coli showed slow growth. However, under appropriate conditions, this protein could be over expressed and purified for structure function studies. The second component of the cell lysis cassette, ybcS, was found to have an N-terminal SAR (Signal Arrest Release) transmembrane domain. The construct of ybcS has been over expressed in E.coli and the purified protein was functional, exhibiting lytic activity against E.coli and Salmonella typhi cell wall substrate. Such targeted sequence- structure-function characterization of proteins encoded by cryptic prophages will help understand the contribution of prophage proteins to bacterial evolution.

Morphological and genetic diversity of temperate phages in Clostridium difficile

Eight temperate phages were characterized after mitomycin C induction of six Clostridium difficile isolates corresponding to six distinct PCR ribotypes. The hypervirulent C. difficile strain responsible for a multi-institutional outbreak (NAP1/027 or QCD-32g58) was among these prophage-containing strains. Observation of the crude lysates by transmission electron microscopy (TEM) revealed the presence of three phages with isometric capsids and long contractile tails (Myoviridae family), as well as five phages with long noncontractile tails (Siphoviridae family). TEM analyses also revealed the presence of a significant number of phage tail-like particles in all the lysates. Southern hybridization experiments with restricted prophage DNA showed that C. difficile phages belonging to the family Myoviridae are highly similar and most likely related to previously described prophages phiC2, phiC5, and phiCD119. On the other hand, members of the Siphoviridae phage family are more genetically divergent, suggesting that they originated from distantly related ancestors. Our data thus suggest that there are at least three genetically distinct groups of temperate phages in C. difficile; one group is composed of highly related myophages, and the other two groups are composed of more genetically heterogeneous siphophages. Finally, no gene homologous to genes encoding C. difficile toxins or toxin regulators could be identified in the genomes of these phages using DNA hybridization. Interestingly, each unique phage restriction profile correlated with a specific C. difficile PCR ribotype.

Phage_Finder: automated identification and classification of prophage regions in complete bacterial genome sequences

Phage_Finder, a heuristic computer program, was created to identify prophage regions in completed bacterial genomes. Using a test dataset of 42 bacterial genomes whose prophages have been manually identified, Phage_Finder found 91% of the regions, resulting in 7% false positive and 9% false negative prophages. A search of 302 complete bacterial genomes predicted 403 putative prophage regions, accounting for 2.7% of the total bacterial DNA. Analysis of the 285 putative attachment sites revealed tRNAs are targets for integration slightly more frequently (33%) than intergenic (31%) or intragenic (28%) regions, while tmRNAs were targeted in 8% of the regions. The most popular tRNA targets were Arg, Leu, Ser and Thr. Mapping of the insertion point on a consensus tRNA molecule revealed novel insertion points on the 5′ side of the D loop, the 3′ side of the anticodon loop and the anticodon. A novel method of constructing phylogenetic trees of phages and prophages was developed based on the mean of the BLAST score ratio (BSR) of the phage/prophage proteomes. This method verified many known bacteriophage groups, making this a useful tool for predicting the relationships of prophages from bacterial genomes.

Comparative genomics and transcriptional analysis of prophages identified in the genomes of Lactobacillus gasseri, Lactobacillus salivarius, and Lactobacillus casei

Lactobacillus gasseri ATCC 33323, Lactobacillus salivarius subsp. salivarius UCC 118, and Lactobacillus casei ATCC 334 contain one (LgaI), four (Sal1, Sal2, Sal3, Sal4), and one (Lca1) distinguishable prophage sequences, respectively. Sequence analysis revealed that LgaI, Lca1, Sal1, and Sal2 prophages belong to the group of Sfi11-like pac site and cos site Siphoviridae, respectively. Phylogenetic investigation of these newly described prophage sequences revealed that they have not followed an evolutionary development similar to that of their bacterial hosts and that they show a high degree of diversity, even within a species. The attachment sites were determined for all these prophage elements; LgaI as well as Sal1 integrates in tRNA genes, while prophage Sal2 integrates in a predicted arginino-succinate lyase-encoding gene. In contrast, Lca1 and the Sal3 and Sal4 prophage remnants are integrated in noncoding regions in the L. casei ATCC 334 and L. salivarius UCC 118 genomes. Northern analysis showed that large parts of the prophage genomes are transcriptionally silent and that transcription is limited to genome segments located near the attachment site. Finally, pulsed-field gel electrophoresis followed by Southern blot hybridization with specific prophage probes indicates that these prophage sequences are narrowly distributed within lactobacilli.

Prophage-like elements in bifidobacteria: insights from genomics, transcription, integration, distribution, and phylogenetic analysis

So far, there is only fragmentary and unconfirmed information on bacteriophages infecting the genus Bifidobacterium. In this report we analyzed three prophage-like elements that are present in the genomes of Bifidobacterium breve UCC 2003, Bifidobacterium longum NCC 2705, and Bifidobacterium longum DJO10A, designated Bbr-1, Bl-1, and Blj-1, respectively. These prophagelike elements exhibit homology with genes of double-stranded DNA bacteriophages spanning a broad phylogenetic range of host bacteria and are surprisingly closely related to bacteriophages infecting low-G+C bacteria. All three prophage-like elements are integrated in a tRNA(Met) gene, which appears to be reconstructed following phage integration. Analysis of the distribution of this integration site in many bifidobacterial species revealed that the attB sites are well conserved. The Blj-1 prophage is 36.9 kb long and was induced when a B. longum DJO10A culture was exposed to mitomycin C or hydrogen peroxide. The Bbr-1 prophage-like element appears to consist of a noninducible 28.5-kb chimeric DNA fragment composed of a composite mobile element inserted into prophage-like sequences, which do not appear to be widely distributed among B. breve strains. Northern blot analysis of the Bbr-1 prophage-like element showed that large parts of its genome are transcriptionally silent. Interestingly, a gene predicted to encode an extracellular beta-glucosidase carried within the Bbr-1 prophage-like element was shown to be transcribed.

GIL16, a new gram-positive tectiviral phage related to the Bacillus thuringiensis GIL01 and the Bacillus cereus pBClin15 elements

One of the most notable characteristics of Tectiviridae resides in their double-layer coats: the double-stranded DNA is located within a flexible lipoprotein vesicle covered by a rigid protein capsid. Despite their apparent rarity, tectiviruses have an extremely wide distribution compared to other phage groups. Members of this family have been found to infect gram-negative (PRD1 and relatives) as well as gram-positive (Bam35, GIL01, AP50, and phiNS11) hosts. Several reports have shown that tectiviruses infecting gram-negative bacteria are closely related, whereas no information is currently available on the genetic relationship among those infecting gram-positive bacteria. The present study reports the sequence of GIL16, a new isolate originating from Bacillus thuringiensis, and a genetic comparison of this isolate with the tectiviral bacteriophages Bam35 and GIL01, which originated from B. thuringiensis serovars Alesti and Israelensis, respectively. In contrast to PRD1 and its relatives, these are temperate bacteriophages existing as autonomous linear prophages within the host cell. Mutations in a particular motif in both the GIL01 and GIL16 phages are also shown to correlate with a switch to the lytic cycle. Interestingly, both bacterial viruses displayed narrow, yet slightly different, host spectrums. We also explore the hypothesis that pBClin15, a linear plasmid hosted by the Bacillus cereus reference strain ATCC 14579, is also a prophage. Sequencing of its inverted repeats at both extremities and a comparison with GIL01 and GIL16 emphasize its relationship to the Tectiviridae.

Identification of bacteriophage types and their carriage in Staphylococcus aureus

Conserved genomic sequences distinctive of Staphylococcus aureus phage types 3A, 11, 77, 187 and Twort, representative of phage serogroups A, B, F, L and D, were identified and characterized. PCR primers designed for the above sequences were used for development of a multiplex PCR assay which enabled us not only to classify all phages of the International Typing Set plus 16 additional phages, but also to detect prophages in S. aureus genomes. One to four different prophages were unambiguously detected in experimentally lysogenized S. aureus strains, and substantial variation in prophage content was found in 176 S. aureus clinical strains of different provenance. In addition, by using a comparative genomics approach, all the prophages in the S. aureus genomes sequenced to date could be revealed and classified.

The Bacillus thuringiensis linear double-stranded DNA phage Bam35, which is highly similar to the Bacillus cereus linear plasmid pBClin15, has a prophage state

Bam35, a 15-kbp double-stranded DNA phage, infects Bacillus thuringiensis. Recently, sequencing of the related Bacillus cereus revealed a 15.1-kbp linear plasmid, pBClin15. We show that pBClin15 closely resembles Bam35 and demonstrate conversion of Bam35 to a prophage. This state is common, as several B. thuringiensis strains release Bam35-related viruses.